Quantum molecular dynamics and ab initio studies of the crystal structure of hydrogen and deuterium

Hydrogen isotopes play a critical role both in inertial and magnetic confinemen Nuclear Fusion. Since the preferent fuel needed for this technology is a mixture of deuterium and tritium. The study of these isotopes particularly at very low temperatures carries a technological interest in other applications. The present line promotes a deep study on the structural configuration that hydrogen and deuterium adopt at cryogenic temperatures and at high pressures. Typical conditions occurring in present Inertial Fusion target designs. Our approach is aims to determine the crystal structure characteristics, phase transitions and other parameters strongly correlated to variations of temperature and pressure.

Electronic structure and transport properties of atomic NiO spinvalves

Ab initio quantum transport calculations show that short NiO chains suspended in Ni nanocontacts present a very strong spin-polarization of the conductance.The generalized gradient approximation we use here predicts a similar polarization of the conductance as the one previously computed with non-local exchange, confirming the robustness of the result. Their use as nanoscopic spinvalves is proposed.

Mécanismes de transfert de proton d’une réaction acido-basique en phase aqueuse : une étude ab-initio

Les réactions de transfert de proton se retrouvent abondamment dans la nature et sont des processus cruciaux dans plusieurs réactions chimiques et biologiques, qui se produisent souvent en milieu aqueux. Les mécanismes régissant ces échanges de protons sont complexes et encore mal compris, suscitant un intérêt des chercheurs en vue d’une meilleure compréhension fondamentale du processus de transfert. Le présent manuscrit présente une étude mécanistique portant sur une réaction de transfert de proton entre un acide (phénol fonctionnalisé) et une base (ion carboxylate) en phase aqueuse. Les résultats obtenus sont basés sur un grand nombre de simulations de dynamique moléculaire ab-initio réalisées pour des systèmes de type « donneur-pont-accepteur », où le pont se trouve à être une unique molécule d’eau, permettant ainsi l’élaboration d’un modèle cinétique détaillé pour le système étudié. La voie de transfert principalement observée est un processus ultra-rapide (moins d’une picoseconde) passant par la formation d’une structure de type « Eigen » (H9O4+) pour la molécule d’eau pontante, menant directement à la formation des produits. Une seconde structure de la molécule d’eau pontante est également observée, soit une configuration de type « Zündel » (H5O2+) impliquant l’accepteur de proton (l’ion carboxylate) qui semble agir comme un cul-de-sac pour la réaction de transfert de proton.

Mécanismes de transfert de proton d’une réaction acido-basique en phase aqueuse : une étude ab-initio

Les réactions de transfert de proton se retrouvent abondamment dans la nature et sont des processus cruciaux dans plusieurs réactions chimiques et biologiques, qui se produisent souvent en milieu aqueux. Les mécanismes régissant ces échanges de protons sont complexes et encore mal compris, suscitant un intérêt des chercheurs en vue d’une meilleure compréhension fondamentale du processus de transfert. Le présent manuscrit présente une étude mécanistique portant sur une réaction de transfert de proton entre un acide (phénol fonctionnalisé) et une base (ion carboxylate) en phase aqueuse. Les résultats obtenus sont basés sur un grand nombre de simulations de dynamique moléculaire ab-initio réalisées pour des systèmes de type « donneur-pont-accepteur », où le pont se trouve à être une unique molécule d’eau, permettant ainsi l’élaboration d’un modèle cinétique détaillé pour le système étudié. La voie de transfert principalement observée est un processus ultra-rapide (moins d’une picoseconde) passant par la formation d’une structure de type « Eigen » (H9O4+) pour la molécule d’eau pontante, menant directement à la formation des produits. Une seconde structure de la molécule d’eau pontante est également observée, soit une configuration de type « Zündel » (H5O2+) impliquant l’accepteur de proton (l’ion carboxylate) qui semble agir comme un cul-de-sac pour la réaction de transfert de proton.

Theory of strongly correlated electron systems. I. Intersite coulomb interaction and the approximation of renormalized fermions in total energy calculations

The diagrammatic strong-coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many - electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one-to-one correspondence is found between the subset of such graphs for the many - electron GFs and the complete set of Feynman graphs of weak-coupling perturbation theory (WCPT) for single-electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many-electron quasi-particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self-energy and the exchange- correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LIDA (RF LDA), is obtained by introducing the spectral weights of the many electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self-interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LIDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part 11 of this study (this issue). (c) 2005 Wiley Periodicals, Inc.

Theory of strongly correlated electron systems. II. Including correlation effects into electronic structure calculations

We have previously shown that a division of the f-shell into two subsystems gives a better understanding of the cohesive properties as well the general behavior of lanthanide systems. In this article, we present numerical computations, using the suggested method. We show that the picture is consistent with most experimental data, e.g., the equilibrium volume and electronic structure in general. Compared with standard energy band calculations and calculations based on the self-interaction correction and LIDA + U, the f-(non-f)-mixing interaction is decreased by spectral weights of the many-body states of the f-ion. (c) 2005 Wiley Periodicals, Inc.

Ab initio quantum chemical studies on neutral-radical reactions of ethynyl (C2H) and cyano (CN) with unsaturated hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA's Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth's primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan's atmosphere, the ISM, and cold celestial bodies.<.

Ab Initio Quantum Chemical Studies on Neutral-Radical Reactions of Ethynyl (C2H) and Cyano (CN) with Unsaturated Hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA’s Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth’s primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan’s atmosphere, the ISM, and cold celestial bodies.

Unraveling effects of disorder on the electronic structure of SiO(2) from first principles

We present a first-principles systematic study of the electronic structure of SiO(2) including the crystalline polymorphs alpha quartz and beta cristobalite, and different types of disorder leading to the amorphous phase. We start from calculations within density functional theory and proceed to more sophisticated quasiparticle calculations according to the GW scheme. Our results show that different origins of disorder have also different impact on atomic and electronic-density fluctuations, which affect the electronic structure and, in particular, the size of the mobility gap in each case.

Accurate prediction of the Si/SiO(2) interface band offset using the self-consistent ab initio DFT/LDA-1/2 method

We use the density functional theory/local-density approximation (DFT/LDA)-1/2 method [L. G. Ferreira , Phys. Rev. B 78, 125116 (2008)], which attempts to fix the electron self-energy deficiency of DFT/LDA by half-ionizing the whole Bloch band of the crystal, to calculate the band offsets of two Si/SiO(2) interface models. Our results are similar to those obtained with a ""state-of-the-art"" GW approach [R. Shaltaf , Phys. Rev. Lett. 100, 186401 (2008)], with the advantage of being as computationally inexpensive as the usual DFT/LDA. Our band gap and band offset predictions are in excellent agreement with experiments.

Amorphous HfO(2) and Hf(1-x)Si(x)O via a melt-and-quench scheme using ab initio molecular dynamics

We have performed ab initio molecular dynamics simulations to generate an atomic structure model of amorphous hafnium oxide (a-HfO(2)) via a melt-and-quench scheme. This structure is analyzed via bond-angle and partial pair distribution functions. These results give a Hf-O average nearest-neighbor distance of 2.2 angstrom, which should be compared to the bulk value, which ranges from 1.96 to 2.54 angstrom. We have also investigated the neutral O vacancy and a substitutional Si impurity for various sites, as well as the amorphous phase of Hf(1-x)Si(x)O(2) for x=0.25, 0375, and 0.5.

Electronic structure of In(2)O(3) and Sn-doped In(2)O(3) by hard x-ray photoemission spectroscopy

The valence and core levels of In(2)O(3) and Sn-doped In(2)O(3) have been studied by hard x-ray photoemission spectroscopy (hv = 6000 eV) and by conventional Al K alpha (hv = 1486.6 eV) x-ray photoemission spectroscopy. The experimental spectra are compared with density-functional theory calculations. It is shown that structure deriving from electronic levels with significant In or Sn 5s character is selectively enhanced under 6000 eV excitation. This allows us to infer that conduction band states in Sn-doped samples and states at the bottom of the valence band both contain a pronounced In 5s contribution. The In 3d core line measured at hv = 1486.6 eV for both undoped and Sn-doped In(2)O(3) display an asymmetric lineshape, and may be fitted with two components associated with screened and unscreened final states. The In 3d core line spectra excited at hv = 6000 eV for the Sn-doped samples display pronounced shoulders and demand a fit with two components. The In 3d core line spectrum for the undoped sample can also be fitted with two components, although the relative intensity of the component associated with the screened final state is low, compared to excitation at 1486.6 eV. These results are consistent with a high concentration of carriers confined close to the surface of nominally undoped In(2)O(3). This conclusion is in accord with the fact that a conduction band feature observed for undoped In(2)O(3) in Al K alpha x-ray photoemission is much weaker than expected in hard x-ray photoemission.

Full-dimensional analytical ab initio potential energy surface of the ground state of HOI

Extensive ab initio calculations using a complete active space second-order perturbation theory wavefunction, including scalar and spin-orbit relativistic effects with a quadruple-zeta quality basis set were used to construct an analytical potential energy surface (PES) of the ground state of the [H, O, I] system. A total of 5344 points were fit to a three-dimensional function of the internuclear distances, with a global root-mean-square error of 1.26 kcal mol(-1). The resulting PES describes accurately the main features of this system: the HOI and HIO isomers, the transition state between them, and all dissociation asymptotes. After a small adjustment, using a scaling factor on the internal coordinates of HOI, the frequencies calculated in this work agree with the experimental data available within 10 cm(-1). (C) 2011 American Institute of Physics. [doi: 10.1063/1.3615545]