New Relativistic Atomic Natural Orbital Basis Sets for Lanthanide Atoms with Applications to the Ce Diatom and LuF3

New basis sets of the atomic natural orbital (ANO) type have been developed for the lanthanide atoms La-Lu. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies and some excitation energies. Computed ionization energies have an accuracy better than 0.1 eV in most cases. Two molecular applications are inluded as illustration: the cerium diatom and the LuF3 molecule. In both cases it is shown that 4f orbitals are not involved in the chemical bond in contrast to an earlier claim for the latter molecule.

Axial ligand influence on geometries, charge distributions and electronic structures of iron tetraazamacrocycle [Fe((II))TIM(X)(Y)](2+) complexes assessed by Density Functional Theory

We employed the Density Functional Theory along with small basis sets, B3LYP/LANL2DZ, for the study of FeTIM complexes with different pairs of axial ligands (CO, H(2)O, NH(3), imidazole and CH(3)CN). These calculations did not result in relevant changes of molecular quantities as bond lengths, vibrational frequencies and electronic populations supporting any significant back-donation to the carbonyl or acetonitrile axial ligands. Moreover, a back-donation mechanism to the macrocycle cannot be used to explain the observed changes in molecular properties along these complexes with CO or CH(3)CN. This work also indicates that complexes with CO show smaller binding energies and are less stable than complexes with CH(3)CN. Further, the electronic band with the largest intensity in the visible region (or close to this region) is associated to the transition from an occupied 3d orbital on iron to an empty pi* orbital located at the macrocycle. The energy of this Metal-to-Ligand Charge Transfer (MLCT) transition shows a linear relation to the total charge of the macrocycle in these complexes as given by Mulliken or Natural Population Analysis (NPA) formalisms. Finally, the macrocycle total charge seems to be influenced by the field induced by the axial ligands. (C) 2011 Elsevier Ltd. All rights reserved.

Efficiency of extrinsic and intrinsic charge-carrier photogeneration processes obtained from the steady-state photocurrent action spectra of poly(p-phenylene vinylene) derivatives

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Giant first-forbidden resonances

Recent experimental data on first-forbidden charge-exchange resonances are discussed in the framework of a schematic model. We also evaluate the screening of the weak coupling constants induced by both the giant resonances and the Δ-isobar. It is shown that the last effect does not depend on the multipolarity of the one-particle moment. Due to the same reason, the fraction of the reaction strength pushed up into the Δ-resonance region is always the same regardless of the quantum numbers carried by the excitation. Simple expressions are derived for the dependence of the excitation energies of the first-forbidden giant resonances on the mass number and isospin of the target. The model reproduces consistently both the Gamow-Teller and the first-forbidden resonances. © 1983.

Computer simulation study of collective dynamics in the glass former Ca(NO3)(2)center dot 4H(2)O

Time correlation functions of current fluctuations were calculated by molecular dynamics (MD) simulations in order to investigate sound waves of high wavevectors in the glass-forming liquid Ca(NO3)(2)center dot 4H(2)O. Dispersion curves, omega(k), were obtained for longitudinal (LA) and transverse acoustic (TA) modes, and also for longitudinal optic (LO) modes. Spectra of LA modes calculated by MD simulations were modeled by a viscoelastic model within the memory function framework. The viscoelastic model is used to rationalize the change of slope taking place at k similar to 0.3 angstrom(-1) in the omega(k) curve of acoustic modes. For still larger wavevectors, mixing of acoustic and optic modes is observed. Partial time correlation functions of longitudinal mass currents were calculated separately for the ions and the water molecules. The wavevector dependence of excitation energies of the corresponding partial LA modes indicates the coexistence of a relatively stiff subsystem made of cations and anions, and a softer subsystem made of water molecules. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751548]

Two-photon absorption in oxazole derivatives: An experimental and quantum chemical study

Experimental and theoretical studies on the two-photon absorption properties of two oxazole derivatives: 2,5-diphenyloxazole (PPO) and 2-(4-biphenylyI)-5-phenyl-1,3,4-oxadiazole (PBD) are presented. The two-photon absorption cross-section spectra were determined by means of the Z-scan technique, from 460 up to 650 nm, and reached peak values of 84 GM for PBD and 27 GM for PPO. Density Functional Theory and response function formalism are used to determine the molecular structures and the one- and two-photon absorption properties and to assist in the interpretation of the experimental results. The Polarizable Continuum Model in one-photon absorption calculations is used to estimate solvent effects. (C) 2011 Elsevier B.V. All rights reserved.

Efficiency of extrinsic and intrinsic charge-carrier photogeneration processes obtained from the steady-state photocurrent action spectra of poly(p-phenylene vinylene) derivatives

The efficiency of the charge-carrier photogeneration processes in poly(2,5-bis(3',7'-dimethyl-octyloxy)-1,4-phenylene vinylene) (OC(1)OC10-PPV) has been analyzed by the spectral response of the photocurrent of devices in ITO/polymer/Al structures. The symbatic response of the photocurrent action spectra of the OC1OC10-PPV devices, obtained for light-excitation through the ITO electrode and for forward bias, has been fitted using a phenomenological model which considers that the predominant transport mechanism under external applied electric field is the drift of photogenerated charge-carriers, neglecting charge-carrier diffusion. The proposed model takes into account that charge-carrier photogeneration occurs via intermediate stages of bounded pairs (excitonic states), followed by dissociation processes. Such processes result in two different contributions to the photoconductivity: The first one, associated to direct creation of unbound polaron pairs due to intrinsic photoionization; and the second one is associated to secondary processes like extrinsic photoinjection at the metallic electrodes. The results obtained from the model have shown that the intrinsic component of the photoconductivity at higher excitation energies has a considerably higher efficiency than the extrinsic one, suggesting a dependence on the photon energy for the efficiency of the photogeneration process.

Determination of degree of ionization of poly(allylamine hydrochloride) (PAH) and poly[1-[4-(3-carboxy‑4 hydroxyphenylazo)benzene sulfonamido]-1,2-ethanediyl, sodium salt] (PAZO) in layer-by-layer films using vacuum photoabsorption spectroscopy

Electrostatic and hydrophobic interactions govern most of the properties of supramolecular systems, which is the reason determining the degree of ionization of macromolecules has become crucial for many applications. In this paper, we show that highresolution ultraviolet spectroscopy (VUV) can be used to determine the degree of ionization and its effect on the electronic excitation energies of layer-by-layer (LbL) films of poly(allylamine hydrochloride) (PAH) and poly[1-[4-(3-carboxy-4 hydroxyphenylazo)- benzene sulfonamido]-1,2-ethanediyl, sodium salt] (PAZO). A full assignment of the VUV peaks of these polyelectrolytes in solution and in cast or LbL films could be made, with their pH dependence allowing us to determine the p'K IND. a' using the Henderson-Hasselbach equation. The p'K IND. a' for PAZO increased from ca. 6 in solution to ca. 7.3 in LbL films owing to the charge transfer from PAH. Significantly, even using solutions at a fixed pH for PAH, the amount adsorbed on the LbL films still varied with the pH of the PAZO solutions due to these molecular-level interactions. Therefore, the procedure based on a comparison of VUV spectra from solutions and films obtained under distinct conditions is useful to determine the degree of dissociation of macromolecules, in addition to permitting interrogation of interface effects in multilayer films.

Optische Spektroskopie an Fermium (Z=100)

Im Rahmen der vorliegenden Arbeit wurde erstmals Laser-Atomspektroskopie an einem Element durchgeführt, für das bisher keine atomaren Niveaus bekannt waren. Die Experimente wurden am Element Fermium mit der Ordnungszahl Z=100 mit der Resonanzionisationsspektroskopie (RIS) in einer Puffergaszelle durchgeführt. Verwendet wurde das Isotop 255Fm mit einer Halbwertszeit von 20.1 h, das im Hochflusskernreaktor des ORNL, Oak Ridge, USA, hergestellt wurde. Die von einem elektrochemischen Filament in das Argon-Puffergas bei einer Temperatur von 960(20)°C abgedampften Fm-Atome wurden mit Lasern in einem Zweistufenprozess resonant ionisiert. Dazu wurde das Licht eines Excimerlaser gepumpten Farbstofflasers für den ersten Anregungsschritt um die Wellenlänge 400 nm durchgestimmt. Ein Teil des Excimer (XeF) Laser Pumplichtes mit den Wellenlänge 351/353 nm wurde für die nicht-resonante Ionisation verwendet. Die Ionen wurden mit Hilfe elektrischer Felder aus der optischen Zelle extrahiert und nach einem Quadrupol Massenfilter mit einem Channeltron-Detektor massenselektiv nachgewiesen. Trotz der geringen Probenmenge von 2.7 x 10^10 eingesetzten Atomen wurden zwei atomare Resonanzen bei Energien von 25099.8(2) cm-1 und 25111.8(2) cm-1 gefunden und das Sättigungsverhalten dieser Linien gemessen. Es wurde ein theoretisches Modell entwickelt, dass sowohl das spektrale Profil der sättigungsverbreiterten Linien als auch die Sättigungskurven beschreibt. Durch Anpassung an die Messdaten konnten die partiellen Übergangsraten in den 3H6 Grundzustand Aki=3.6(7) x 10^6/s und Aki=3.6(6) x 10^6/s bestimmt werden. Der Vergleich der Niveauenergien und Übergangsraten mit Multikonfigurations Dirac-Fock Rechnungen legt die spektroskopische Klassifizierung der beobachteten Niveaus als 5f12 7s7p 5I6 und 5G6 Terme nahe. Weiterhin wurde ein Übergang bei 25740 cm-1 gefunden, der aufgrund der beobachteten Linienbreite von 1000 GHz als Rydbergzustand Zustand mit der Niveauenergie 51480 cm-1 interpretiert wurde und über einen Zweiphotonen Prozess angeregt werden kann. Basierend auf dieser Annahme wurde die Obergrenze für die Ionisationsenergie IP = 52140 cm-1 = 6.5 eV abgeschätzt. In den Messungen wurden Verschiebungen in den Zeitverteilungsspektren zwischen den mono-atomaren Ionen Fm+ und Cf+ und dem Molekül-Ion UO+ festgestellt und auf Driftzeitunterschiede im elektrischen Feld der gasgefüllten optischen Zelle zurückgeführt. Unter einfachen Modellannahmen wurde daraus auf die relativen Unterschiede Delta_r(Fm+,Cf+)/r(Cf+) ˜ -0.2 % und Delta_r(UO+,Cf+)/r(Cf+) ˜ 20 % in den Ionenradien geschlossen. Über die Bestimmung der Abnahme der Fm-a Aktivität des Filamentes auf der einen Seite und die Messung der Resonanzzählrate auf der anderen Seite, wurde die Nachweiseffizienz der Apparatur zu 4.5(3) x 10^-4 bestimmt. Die Nachweisapparatur wurde mit dem Ziel weiterentwickelt, Laserspektroskopie am Isotop 251Fm durchzuführen, das über die Reaktion 249Cf(a,2n)251Fm direkt in der optischen Zelle erzeugt werden soll. Das Verfahren wurde am chemischen Homolog Erbium getestet. Dabei wurde das Isotop 163Er über die Reaktion 161Dy(a,2n)163Er erzeugt und nach Resonanzionisation nachgewiesen. Die Nachweiseffizienz der Methode wurde zu 1 x 10^-4 bestimmt.

Towards a Reconstruction of Thermal Properties of Light Nuclei from Fusion - Evaporation reactions

This thesis work has been developed in the framework of a new experimental campaign, proposed by the NUCL-EX Collaboration (INFN III Group), in order to progress in the understanding of the statistical properties of light nuclei, at excitation energies above particle emission threshold, by measuring exclusive data from fusion-evaporation reactions. The determination of the nuclear level density in the A~20 region, the understanding of the statistical behavior of light nuclei with excitation energies ~3 A.MeV, and the measurement of observables linked to the presence of cluster structures of nuclear excited levels are the main physics goals of this work. On the theory side, the contribution to this project given by this work lies in the development of a dedicated Monte-Carlo Hauser-Feshbach code for the evaporation of the compound nucleus. The experimental part of this thesis has consisted in the participation to the measurement 12C+12C at 95 MeV beam energy, at Laboratori Nazionali di Legnaro - INFN, using the GARFIELD+Ring Counter(RCo) set-up, from the beam-time request to the data taking, data reduction, detector calibrations and data analysis. Different results of the data analysis are presented in this thesis, together with a theoretical study of the system, performed with the new statistical decay code. As a result of this work, constraints on the nuclear level density at high excitation energy for light systems ranging from C up to Mg are given. Moreover, pre-equilibrium effects, tentatively interpreted as alpha-clustering effects, are put in evidence, both in the entrance channel of the reaction and in the dissipative dynamics on the path towards thermalisation.

Search for K isomers in 252,254 No and 260 Sg and investigation of their nuclear structure

Study of K isomerism in the transfermium region around the deformed shells at N=152, Z=102, and N=162, Z=108 provides important information on the structure of heavy nuclei. Recent calculations suggest that the K-isomerism can enhance the stability of such nuclei against alpha emission and spontaneous fission. Nuclei showing K isomerism have neutron and proton orbitals with large spin projections on the symmetry axis which is due to multi quasiparticle states with aligned spins K. Quasi-particle states are formed by breaking pairs of nucleons and raising one or two nucleons in orbitals near the Fermi surface above the gap, forming high K (multi)quasi-particle states mainly at low excitation energies. Experimental examples are the recently studied two quasi-particle K isomers in 250,256-Fm, 254-No, and 270-Ds. Nuclei in this region, are produced with cross sections ranging from several nb up to µb, which are high enough for a detailed decay study. In this work, K isomerism in Sg and No isotopes was studied at the velocity filter SHIP of GSI, Darmstadt. The data were obtained by using a new data acquisition system which was developed and installed during this work. 252,254-No and 260-Sg were produced in fusion evaporation reactions of 48-Ca and 54-Cr projectiles with 206,208-Pb targets at beam energies close to the Coulomb barrier. A new K isomer was discovered in 252-No at excitation energy of 1.25 MeV, which decays to the ground state rotational band via gamma emission. It has a half-life of about 100 ms. The population of the isomeric state was about 20% of the ground state population. Detailed investigations were performed on 254-No in which two isomeric states (275 ms and 198 µs) were already discovered by R.-D. Herzberg, but due to the higher number of observed gamma decays more detailed information about the decay path of the isomers was obtained in the present work. In 260-Sg, we observed no statistically significant component with a half life different from that of the ground state. A comparison between experimental results and theoretical calculations of the single particle energies shows a fair agreement. The structure of the here studied nuclei is in particular important as single particle levels are involved which are relevant for the next shell closure expected to form the region of the shell stabilized superheavy elements at proton numbers 114, 120, or 126 and neutron number 184. K isomers, in particular, could be an ideal tool for the synthesis and study of these isotopes due to enhanced spontaneous fission life times which could result in higher alpha to spontaneous fission branching ratios and longer half lifes.

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.

Experimental and Calculated Spectra of π-Stacked Mild Charge-Transfer Complexes: Jet-Cooled Perylene·(Tetrachloroethene) n , n = 1,2

The S0 ↔ S1 spectra of the mild charge-transfer (CT) complexes perylene·tetrachloroethene (P·4ClE) and perylene·(tetrachloroethene)2 (P·(4ClE)2) are investigated by two-color resonant two-photon ionization (2C-R2PI) and dispersed fluorescence spectroscopy in supersonic jets. The S0 → S1 vibrationless transitions of P·4ClE and P·(4ClE)2 are shifted by δν = −451 and −858 cm–1 relative to perylene, translating to excited-state dissociation energy increases of 5.4 and 10.3 kJ/mol, respectively. The red shift is ∼30% larger than that of perylene·trans-1,2-dichloroethene; therefore, the increase in chlorination increases the excited-state stabilization and CT character of the interaction, but the electronic excitation remains largely confined to the perylene moiety. The 2C-R2PI and fluorescence spectra of P·4ClE exhibit strong progressions in the perylene intramolecular twist (1au) vibration (42 cm–1 in S0 and 55 cm–1 in S1), signaling that perylene deforms along its twist coordinate upon electronic excitation. The intermolecular stretching (Tz) and internal rotation (Rc) vibrations are weak; therefore, the P·4ClE intermolecular potential energy surface (IPES) changes little during the S0 ↔ S1 transition. The minimum-energy structures and inter- and intramolecular vibrational frequencies of P·4ClE and P·(4ClE)2 are calculated with the dispersion-corrected density functional theory (DFT) methods B97-D3, ωB97X-D, M06, and M06-2X and the spin-consistent-scaled (SCS) variant of the approximate second-order coupled-cluster method, SCS-CC2. All methods predict the global minima to be π-stacked centered coplanar structures with the long axis of tetrachloroethene rotated by τ ≈ 60° relative to the perylene long axis. The calculated binding energies are in the range of −D0 = 28–35 kJ/mol. A second minimum is predicted with τ ≈ 25°, with ∼1 kJ/mol smaller binding energy. Although both monomers are achiral, both the P·4ClE and P·(4ClE)2 complexes are chiral. The best agreement for adiabatic excitation energies and vibrational frequencies is observed for the ωB97X-D and M06-2X DFT methods.

The elusive S2 state, the S1/S2 splitting, and the excimer states of the benzene dimer

We observe the weak S 0 → S 2 transitions of the T-shaped benzene dimers (Bz)2 and (Bz-d 6)2 about 250 cm−1 and 220 cm−1 above their respective S 0 → S 1 electronic origins using two-color resonant two-photon ionization spectroscopy. Spin-component scaled (SCS) second-order approximate coupled-cluster (CC2) calculations predict that for the tipped T-shaped geometry, the S 0 → S 2 electronic oscillator strength f el (S 2) is ∼10 times smaller than f el (S 1) and the S 2 state lies ∼240 cm−1 above S 1, in excellent agreement with experiment. The S 0 → S 1 (ππ ∗) transition is mainly localized on the “stem” benzene, with a minor stem → cap charge-transfer contribution; the S 0 → S 2 transition is mainly localized on the “cap” benzene. The orbitals, electronic oscillator strengths f el (S 1) and f el (S 2), and transition frequencies depend strongly on the tipping angle ω between the two Bz moieties. The SCS-CC2 calculated S 1 and S 2 excitation energies at different T-shaped, stacked-parallel and parallel-displaced stationary points of the (Bz)2 ground-state surface allow to construct approximate S 1 and S 2 potential energy surfaces and reveal their relation to the “excimer” states at the stacked-parallel geometry. The f el (S 1) and f el (S 2) transition dipole moments at the C 2v -symmetric T-shape, parallel-displaced and stacked-parallel geometries are either zero or ∼10 times smaller than at the tipped T-shaped geometry. This unusual property of the S 0 → S 1 and S 0 → S 2 transition-dipole moment surfaces of (Bz)2 restricts its observation by electronic spectroscopy to the tipped and tilted T-shaped geometries; the other ground-state geometries are impossible or extremely difficult to observe. The S 0 → S 1/S 2 spectra of (Bz)2 are compared to those of imidazole ⋅ (Bz)2, which has a rigid triangular structure with a tilted (Bz)2 subunit. The S 0 → S 1/ S 2 transitions of imidazole-(benzene)2 lie at similar energies as those of (Bz)2, confirming our assignment of the (Bz)2 S 0 → S 2 transition.

A new set of relativistic screening constants for the screened hydrogenic model

AnewRelativisticScreenedHydrogenicModel has been developed to calculate atomic data needed to compute the optical and thermodynamic properties of high energy density plasmas. The model is based on anewset of universal screeningconstants, including nlj-splitting that has been obtained by fitting to a large database of ionization potentials and excitation energies. This database was built with energies compiled from the National Institute of Standards and Technology (NIST) database of experimental atomic energy levels, and energies calculated with the Flexible Atomic Code (FAC). The screeningconstants have been computed up to the 5p3/2 subshell using a Genetic Algorithm technique with an objective function designed to minimize both the relative error and the maximum error. To select the best set of screeningconstants some additional physical criteria has been applied, which are based on the reproduction of the filling order of the shells and on obtaining the best ground state configuration. A statistical error analysis has been performed to test the model, which indicated that approximately 88% of the data lie within a ±10% error interval. We validate the model by comparing the results with ionization energies, transition energies, and wave functions computed using sophisticated self-consistent codes and experimental data.