Calculs ab initio de structures électroniques et de leur dépendance en température avec la méthode GW

Cette thèse porte sur le calcul de structures électroniques dans les solides. À l'aide de la théorie de la fonctionnelle de densité, puis de la théorie des perturbations à N-corps, on cherche à calculer la structure de bandes des matériaux de façon aussi précise et efficace que possible. Dans un premier temps, les développements théoriques ayant mené à la théorie de la fonctionnelle de densité (DFT), puis aux équations de Hedin sont présentés. On montre que l'approximation GW constitue une méthode pratique pour calculer la self-énergie, dont les résultats améliorent l'accord de la structure de bandes avec l'expérience par rapport aux calculs DFT. On analyse ensuite la performance des calculs GW dans différents oxydes transparents, soit le ZnO, le SnO2 et le SiO2. Une attention particulière est portée aux modèles de pôle de plasmon, qui permettent d'accélérer grandement les calculs GW en modélisant la matrice diélectrique inverse. Parmi les différents modèles de pôle de plasmon existants, celui de Godby et Needs s'avère être celui qui reproduit le plus fidèlement le calcul complet de la matrice diélectrique inverse dans les matériaux étudiés. La seconde partie de la thèse se concentre sur l'interaction entre les vibrations des atomes du réseau cristallin et les états électroniques. Il est d'abord montré comment le couplage électron-phonon affecte la structure de bandes à température finie et à température nulle, ce qu'on nomme la renormalisation du point zéro (ZPR). On applique ensuite la méthode GW au calcul du couplage électron-phonon dans le diamant. Le ZPR s'avère être fortement amplifié par rapport aux calculs DFT lorsque les corrections GW sont appliquées, améliorant l'accord avec les observations expérimentales.

Vertically coupled double quantum rings at zero magnetic field

Within local-spin-density functional theory, we have investigated the ¿dissociation¿ of few-electron circular vertical semiconductor double quantum ring artificial molecules at zero magnetic field as a function of interring distance. In a first step, the molecules are constituted by two identical quantum rings. When the rings are quantum mechanically strongly coupled, the electronic states are substantially delocalized, and the addition energy spectra of the artificial molecule resemble those of a single quantum ring in the few-electron limit. When the rings are quantum mechanically weakly coupled, the electronic states in the molecule are substantially localized in one ring or the other, although the rings can be electrostatically coupled. The effect of a slight mismatch introduced in the molecules from nominally identical quantum wells, or from changes in the inner radius of the constituent rings, induces localization by offsetting the energy levels in the quantum rings. This plays a crucial role in the appearance of the addition spectra as a function of coupling strength particularly in the weak coupling limit.

Theory for the optical properties of small Hg_n clusters

The static and dynamical polarizabilities of the Hg-dimer are calculated by using a Hubbard Hamiltonian to describe the electronic structure. The Hamiltonian is diagonalized exactly within a subspace of second-quantized electronic states from which only multiply ionized atomic configurations have been excluded. With this approximation we can describe the most important electronic transitions including the effect of charge fluctuations. We analyze the polarizability as a function of the intraatomic Coulomb interaction which represents the repulsion between electrons. We obtain that this interaction results in strong electronic correlations in the excited states and increases the first excitation energy of the dimer by 0.8 eV in comparison to a calculation which neglects correlations, resulting in a better agreement with the experiment.

Femtosecond time-resolved molecular multiphoton ionization and fragmentation of Na_2

We present a comparison between experimental and theoretical results for pump/probe multiphoton ionizing transitions of the sodium dimer, initiated by femtosecond laser pulses. It is shown that the motion of vibrational wavepackets in two electronic states is probed simultaneously and their dynamics is reflected in the total Na^+_2 ion signal which is recorded as a function of the time delay between pump and probe pulse. The time dependent quantum calculations demonstrate that two ionization pathways leading to the same final states of the molecularion exist: one gives an oscillating contribution to the ion signal, the other yields a constant background. From additional measurements of the Na^+ -transient photofragmentation spectrum it is deduced that another ionization process leading to different final ionic states exists. The process includes the excitation of a doubly excitedbound Rydberg state. This conclusion is supported by the theoretical simulation.

A variational method for the calculation of spin-rovibronic energy levels of any triatomic molecule in an electronic triplet state

The triatomic spin-rovibronic variational code RVIB3 has been extended to include the effect of two uncoupled electrons, for both (3)Sigma(-) and (3)Pi (Renner-Teller) electronic states. The spin-orbital-rotational kinetic energy is included in the usual way, via terms (J+L+S). The phenomenological terms AL.S and lambda 2/3(3S(z)(2)) are introduced to reproduce the 3 spin-orbit and spin-spin splittings, respectively. Calculations are performed to evaluate the spin-rovibronic energy levels of CCO (X) over tilde (3) Sigma(-) and CCO (A) over tilde (3) Pi for which the Born-Oppenheimer potentials are derived from high-accuracy ab initio calculations.

Theoretical calculation of absorption intensities of C2H and C2D

The theory of dipole-allowed absorption intensities in triatomic molecules is presented for systems with three close-lying electronic states of doublet multiplicity. Its derivation is within the framework of a recently developed variational method [CARTER, S., HANDY, N. C., PUZZARINI, C., TARRONI, R., and PALMIERI, P., 2000, Molec. Phys., 98,1967]. The method has been applied to the calculation of the infrared absorption spectrum of the C2H radical and its deuterated isotopomer for energies up to 10000 cm(-1) above the ground state, using highly accurate ab initio diabatic potential energy and dipole moment surfaces. The calculated spectra agree very well with those recorded experimentally in a neon matrix [FORNEY, D., JACOX, M. E., and THOMPSON, W. E., 1995, J. molee. Spectrosc., 170, 178] and assignments in the high energy region of the IR spectra are proposed for the first time.

Electronic structure and chemical bonding in W(2) molecule

The electronic structure of the lowest-lying electronic states of W(2) were investigated at the CASPT2 level. The ground state is a X(1)Sigma(+)(g) state, followed by the a(3)Delta(u), b(3)Sigma(+)(u) and A(1)Delta(u) electronic states. Seven low-lying Omega-states were computed: (1)0(g)(+), (2)3(u), (3)2(u), (4)1(u), (5)0(u)(-), (6)1(u), and (7)2(u), with the ground state corresponding to the (1)0(g)(+)(X(1)Sigma(+)(g)) state. Comparison with the other VIB transition metal group dimers indicates a common pattern of electronic structure and spectroscopic properties. (C) 2010 Elsevier B.V. All rights reserved.

Structural and spectroscopic properties of the diazocarbene radical (CNN) and its ions CNN(+) and CNN(-): a high-level theoretical investigation

The diazocarbene radical, CNN, and the ions CNN(+) and CNN(-) were investigated at a high level of theory. Very accurate structural parameters for the states X(3)Sigma(-) and A(3)Pi of CNN, and X(2)Pi of both CNN(+) and CNN(-) were obtained with the UCCSD(T) method using correlated-consistent basis functions with extrapolations to the complete basis set limit, with valence only and also with all electrons correlated. Harmonic and anharmonic frequencies were obtained for all species and the Renner parameter and average frequencies evaluated for the Pi states. At the UCCSD(T)/CBS(T-5) level of theory, Delta(f)H(0 K) = 138.89 kcal/mol and Delta(f)H(298 K) = 139.65 kcal/mol were obtained for diazocarbene; for the ionization potential and the electron affinity of CNN, 10.969 eV (252.95 kcal/mol), and 1.743 eV (40.19 kcal/mol), respectively, are predicted. Geometry optimization was also carried out with the CASSCF/MRCI/CBS(T-5) approach for the states X(3)Sigma(-) A(3)Pi, and a(1)Delta of CNN, and with the CASSCF/MRSDCI/aug-cc-pVTZ approach for the states b(1)Sigma(+), c(1)Pi, d(1)Sigma(-), and B(3)Sigma(-), and excitation energies (T(e)) evaluated. Vertical energies were calculated for 15 electronic states, thus improving on the accuracy of the five transitions already described, and allowing for a reliable overview of a manifold of other states, which is expected to guide future spectroscopic experiments. This study corroborates the experimental assignment for the vertical transition X (3)Sigma(-) <- E (3)Pi.

Excited states of the CaAl molecule: An MRCI study

Accurate potential energy curves, dissociation energies and spectroscopic constants for several low-lying doublet and quartet electronic states of CaAl were investigated using the CASSCF/MRCI methodology, and the cc-pVQZ basis set. Our results represent an improvement over a previous theoretical description, and also characterizes new higher excited states not previously investigated, thus confirming the assignment of four excited states investigated experimentally. With the theoretical transition moment functions, transition probabilities and radiative lifetimes were estimated via Einstein spontaneous emission coefficients. (c) 2008 Elsevier B. V. All rights reserved.

A theoretical study of the binding and electronic spectrum of the Mo-2 molecule

Multiconfigurational SCF and second-order perturbation theory have been employed to study seven low-lying singlet and triplet electronic states of the Mo-2 molecule. The bond order of the ground state has been analyzed based on the effective bond order (EBO), indicating that a fully developed sextuple bond is formed between the two Mo atoms. The experimentally observed excited states a(3)Sigma(+)(u) and A(1)Sigma(+)(u) have been determined and the so-called (3)Lambda excited state identified as the b(3)Sigma(+)(u) state, in agreement with experimental expectations. (C) 2007 Elsevier B.V. All rights reserved.

Theoretical study on the molecular mechanism for the reaction of VO2+ with C2H4

The complex reaction between VO2+ ((1)A(1)/(3)A) and C2H4 (Ag-1(g)/(3)A(1)) to yield VO+ ((1)Delta/(3)Sigma) and CH3CHO ('A'/(3)A) has been studied by means of B3LYP/6-31G* and B3LYP/6-311G(2d,p) calculations. The structures of all reactants, products, intermediates, and transition structures of this reaction have been optimized and characterized at the fundamental singlet and first excited triplet electronic states. Crossing points are localized, and possible spin inversion processes are discussed by means of the intrinsic reaction coordinate approach. Relevant stationary points along the most favorable reaction pathways have been studied at the CCSD/6-311G(2d,p)//B3LYP/6-311G(2d,p) calculation level. The theoretical results allow the development of thermodynamic and kinetic arguments about the reaction pathways of the title process. In the singlet state, the first step is the barrierless obtention of a reactant complex associated with the formation of a V-C bond, while in the triplet state a three-membered ring addition complex with the V bonded to the two C atoms is obtained. Similar behavior is found in the exit channels: the product complexes can be formed from isolated products without barriers. The reactant and product complexes are the most stable stationary points in the singlet and triplet electronic states. From the singlet state reactant complex, two reaction pathways are posssible to reach the triplet state product complex. (i) A mechanism in which a hydrogen transfer process is the first and rate limiting step and the second step is an oxygen transfer between vanadium and carbon atoms with a concomitant change in the spin state. The crossing point between singlet and triplet spin states is not kinetically relevant because it takes place at a later stage occurring in the exit channel. (ii) A mechanism in which the first stage renders a four-membered ring between vanadyl cation and the ethylene fragment and an oxygencarbon bond is formed; on going from this minimum to the second transition structure, associated with a carbon-vanadium bond breaking process, the crossing point between singlet and triplet spin states is reached. The final step is the hydrogen transfer between both carbon atoms to yield the product complex. In this case the spin change opens a lower barrier pathway. The transition structures with larger values of relative energies for both reactive channels of VO2+ ((1)A(1)) + C2H4 (Ag-1) --> VO+ ((3)Sigma) + CH3CHO ((1)A') present similar energies, and the two reaction pathways can be considered as competitive.

A theoretical study on the gas phase reactions of the anions NbO3-, NbO5-, and NbO2(0H)(2)(-) with H2O and O-2

The potential energy surfaces at the singlet (s) and the triplet (t) electronic states associated with the gas-phase ion/molecule reactions of NbO3-, NbO5-, and NbO2(OH)(2)(-) with H2O and O-2 have been investigated by means of DFT calculations at the B3LYP level. An analysis of the results points out that the most favorable reactive channel comprises s-NbO3- reacting with H2O to give an ion-molecule complex s-NbO3(H2O)without a barrier. From this minima, an intramolecular hydrogen transfer takes place between the incoming water molecule and an oxygen atom of the NbO3- fragment to render the most stable minimum, s-NbO2(OH)(2)(-). This oxyhydroxide system reacts with O-2 along a barrierless process to obtain the triplet t-NbO4(OH)(2)(-)-A intermediate, and the crossing point, CP1, between s and t electronic states has been characterized. The next step is the hydrogen-transfer process between the oxygen atom of a hydroxyl group and the one adjacent oxygen atom to render a minimum with the two OH groups near each other, t-NbO4(OH)(2)(-)-B. From this point, the last hydrogen migration takes place, to obtain the product complex, t-NbO5(H2O)(-), that can be connected with the singlet separated products, s-NbO5- and H2O. Therefore, a second crossing point, CP2, has been localized. The nature of the chemical bonding of the key minima (NbO3-, NbO2(OH)(2)(-), NbO4(OH)(2)(-)-B, and NbO5-) in both electronic states of the reaction and an interaction with O-2 has been studied by topological analysis of Becke-Edgecombe electron-localization function (ELF) and atoms-in-molecules (AIM) methodology. The niobium-oxygen interactions are characterized as unshared-electron (ionic) interactions and some oxygen-oxygen interactions as protocovalent bonds.

Periodic study on the structural and electronic properties of bulk, oxidized and reduced SnO2(110) surfaces and the interaction with O-2

The structural and electronic properties of bulk and both oxidized and reduced SnO2(110) surfaces as well as the adsorption process of O-2 on the reduced surface have been investigated by periodic DFT calculations at B3LYP level. The lattice parameters, charge distribution, density of states and band structure are reported for the bulk and surfaces. Surface relaxation effects have been explicitly taken into account by optimizing slab models of nine and seven atomic layers representing the oxidized and reduced surfaces, respectively. The conductivity behavior of the reduced SnO2(110) surface is explained by a distribution of the electrons in the electronic states in the band gap induced by oxygen vacancies. Three types of adsorption approaches of O-2 on the four-fold tin at the reduced SuO(2)(110) surface have been considered. The most exothermic channel corresponds to the adsorption of O-2 parallel to the surface and to the four-fold tin row, and it is believed to be associated with the formation of a peroxo O-2(2-) species. The chemisorption of O-2 on reduced SnO2(110) surface causes a significant depopulation of states along the band gap and it is shown to trap the electrons in the chemisorbed complex producing an electron-depleted space-charge layer in the inner surface region of the material in agreement with some experimental evidences. (C) 2002 Elsevier B.V. B.V. All rights reserved.

BOUND HYDROGENIC STATES IN A MAGNETIC FIELD: A FINITE-DIFFERENCE APPROACH

A finite-difference scheme is used to calculate bound electronic states of an electron in a hydrogen atom subject to a magnetic field. The numerical results are in good agreement with exact results, in the absence of the magnetic field, and with a two-parameters variational calculation, when the magnetic field is applied.

Study of the decay B-s(0)->(DsDs(*))-D-(*)

We report a study of the decay B-s(0)->(DsDs(*))-D-(*) using a data sample corresponding to 1.3 fb(-1) of integrated luminosity collected by the D0 experiment in 2002-2006 during run II of the Fermilab Tevatron collider. One D-s((*)) meson was partially reconstructed in the decay D-s ->phi mu nu, and the other D-s((*)) meson was identified using the decay D-s ->phi pi where no attempt was made to distinguish D-s and D-s(*) states. For the branching fraction Br(B-s(0)->(DsDs(*))-D-(*)) we obtain a 90% C.L. range [0.002,0.080] and central value 0.039(-0.017)(+0.019)(stat)(-0.015)(+0.016)(syst). This was subsequently used to make the most precise estimate of the width difference Delta Gamma(CP)(s) in the B-s(0)-(B)over bar(s)(0) system: Delta Gamma(CP)(s)/Gamma(s)=0.079(-0.035)(+0.038)(stat)(-0.030)(+0.031)(syst).