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.

Facilitating Realistic and Efficient Applications of Ab Initio Electronic Dynamics for Polyatomic Systems

Thesis (Ph.D.)--University of Washington, 2016-06

Towards Accurate and Efficient Description of Excited States

Thesis (Ph.D.)--University of Washington, 2016-06

Investigating routes toward atomic layer deposition of silicon carbide: Ab initio screening of potential silicon and carbon precursors

Silicon carbide (SiC) is a promising material for electronics due to its hardness, and ability to carry high currents and high operating temperature. SiC films are currently deposited using chemical vapor deposition (CVD) at high temperatures 1500–1600 °C. However, there is a need to deposit SiC-based films on the surface of high aspect ratio features at low temperatures. One of the most precise thin film deposition techniques on high-aspect-ratio surfaces that operates at low temperatures is atomic layer deposition (ALD). However, there are currently no known methods for ALD of SiC. Herein, the authors present a first-principles thermodynamic analysis so as to screen different precursor combinations for SiC thin films. The authors do this by calculating the Gibbs energy ΔGΔG of the reaction using density functional theory and including the effects of pressure and temperature. This theoretical model was validated for existing chemical reactions in CVD of SiC at 1000 °C. The precursors disilane (Si2H6), silane (SiH4), or monochlorosilane (SiH3Cl) with ethyne (C2H2), carbontetrachloride (CCl4), or trichloromethane (CHCl3) were predicted to be the most promising for ALD of SiC at 400 °C.

The catalytic role of an isolated-Ti atom in the hydrogenation of Ti-doped Al(001) surface : An ab initio density functional theory calculation

Recent work [S. Chaudhuri, J.T. Muckerman, J. Phys. Chem. B 109 (2005) 6952] reported that two Ti-substituted atoms on an Al(0 0 1) surface can form a catalytically active site for the dissociation of H2, but the diffusion barrier of atomic H away from Ti site is as high as 1.57 eV. By using ab initio density functional calculations, we found that two hydrogen molecules can dissociate on isolated-Ti atom doped Al(0 0 1) surface with small activation barriers (0.21 and 0.235 eV for first and second H2, respectively). Additionally, the diffusion barrier of atomic H away from Ti site is also moderate (0.47 eV). These results contribute further towards understanding the improved kinetics observed in recycling of hydrogen with Ti-doped NaAlH4.

The role of V2O5 on the dehydrogenation and hydrogenation in magnesium hydride : An ab initio study

Ab initio density functional theory calculations are performed to study the experimentally observed catalytic role of V2O5 in the recycling of hydrogen in magnesium hydride. We find that the Mg–H bond length becomes elongated when MgH2 clusters are positioned on single, two, and three coordinated oxygen sites (O1, O2, and O3) on the V2O5(001) surface. Molecular hydrogen is predicted to spontaneously form at the hole site on the V2O5(001) surface. Additionally, the activation barrier for the dissociation of hydrogen on V-doped Mg(0001) surface is 0.20 eV, which is only 1/5 of that on pure Mg(0001) surface. Our results indicate that oxygen sites on the V2O5(001)surface and the V dopant in Mg may be important facilitators for dehydrogenation and rehydrogenation, respectively. The understanding gained here will aid in the rational design and development of Mg-based hydrogen storage materials.

Electron capture of tetracyanoethylene oxide in the gas phase. Rearrangement of the parent radical anion to form [(NC)(3)C](-). A joint experimental and ab initio study

Capture of an electron by tetracyanoethylene oxide can initiate a number of decomposition pathways. One of these decompositions yields [(NC)3C]− as the ionic product. Ab initio calculations (at the B3LYP/6-31+G∗ level of theory) indicate that the formation of [(NC)3C]− is initiated by capture of an electron into the LUMO of tetracyanoethylene oxide to yield the anion radical [(NC)2C–O–C(CN)2]−· that undergoes internal nucleophilic substitution to form intermediate [(NC)3C–OCCN]−·. This intermediate dissociates to form [(NC)3C]− (m/z 90) as the ionic product. The radical (NC)3C· has an electron affinity of 4.0 eV (385 kJ mol−1). Ab initio calculations show that [(NC)3C]− is trigonal planar with the negative charge mainly on the nitrogens. A pictorial representation of this structure is the resonance structure formed from three degenerate contributing structures (NC)2–CCN−. The other product of the reaction is nominally (NCCO)·, but there is no definitive experimental evidence to indicate whether this radical survives intact, or decomposes to NC· and CO. The overall process [(NC)2C–O–C(CN)2]−· → [(NC)3C]− + (NCCO)· is calculated to be endothermic by 21 kJ mol−1 with an overall barrier of 268 kJ mol−1.

Ab initio molecular dynamics simulations of excited hydrogen halides and methyl halides

Using excited-state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field approach, we study the mechanism of photodissociation in terms of time evolution of structure, kinetic energy, charges and potential energy for the first excited state of hydrogen halides and methyl halides. Although the hydrogen halides and methyl halides are similar in the photodissociation mechanism, their dynamics are slightly different. The presence of the methyl group causes delay in photodissociation as compared to hydrogen halides.

Ab initio electronic structures of rhombohedral and cubic HgXO3 (X = Ti, Pb)

First-principles calculations were performed for orthorhombic HgO, rhombohedral and cubic phases of HgTiO3 (HTO) and HgPbO3 (HPO). The calculations show that in the rhombohedral phase HTO is a direct gap insulator with a gap of ~1.6 eV. The rhombohedral phase of HPO, on the other hand, shows a weak metallic character. The results provide an explanation for the electrical properties of these compounds. The cubic phases of HTO and HPO are invariably metallic in nature, thereby suggesting that for HTO the rhombohedral–cubic transition must also be accompanied by a change in the electrical state. Examination of the electronic density of states of these systems revealed no significant on-site mixing of Hg 5d and Hg 6s states in any of these materials.

Quantum Chemical Studies on Tropospheric Nucleation Mechanisms Involving Sulfuric Acid

Nucleation is the first step of the process by which gas molecules in the atmosphere condense to form liquid or solid particles. Despite the importance of atmospheric new-particle formation for both climate and health-related issues, little information exists on its precise molecular-level mechanisms. In this thesis, potential nucleation mechanisms involving sulfuric acid together with either water and ammonia or reactive biogenic molecules are studied using quantum chemical methods. Quantum chemistry calculations are based on the numerical solution of Schrödinger's equation for a system of atoms and electrons subject to various sets of approximations, the precise details of which give rise to a large number of model chemistries. A comparison of several different model chemistries indicates that the computational method must be chosen with care if accurate results for sulfuric acid - water - ammonia clusters are desired. Specifically, binding energies are incorrectly predicted by some popular density functionals, and vibrational anharmonicity must be accounted for if quantitatively reliable formation free energies are desired. The calculations reported in this thesis show that a combination of different high-level energy corrections and advanced thermochemical analysis can quantitatively replicate experimental results concerning the hydration of sulfuric acid. The role of ammonia in sulfuric acid - water nucleation was revealed by a series of calculations on molecular clusters of increasing size with respect to all three co-ordinates; sulfuric acid, water and ammonia. As indicated by experimental measurements, ammonia significantly assists the growth of clusters in the sulfuric acid - co-ordinate. The calculations presented in this thesis predict that in atmospheric conditions, this effect becomes important as the number of acid molecules increases from two to three. On the other hand, small molecular clusters are unlikely to contain more than one ammonia molecule per sulfuric acid. This implies that the average NH3:H2SO4 mole ratio of small molecular clusters in atmospheric conditions is likely to be between 1:3 and 1:1. Calculations on charged clusters confirm the experimental result that the HSO4- ion is much more strongly hydrated than neutral sulfuric acid. Preliminary calculations on HSO4- NH3 clusters indicate that ammonia is likely to play at most a minor role in ion-induced nucleation in the sulfuric acid - water system. Calculations of thermodynamic and kinetic parameters for the reaction of stabilized Criegee Intermediates with sulfuric acid demonstrate that quantum chemistry is a powerful tool for investigating chemically complicated nucleation mechanisms. The calculations indicate that if the biogenic Criegee Intermediates have sufficiently long lifetimes in atmospheric conditions, the studied reaction may be an important source of nucleation precursors.

Ab-initio molecular orbital studies on C---H…X hydrogen bonded systemsstar

C---H…X hydrogen bonded systems are studied by the STO-3G method. The proton donor ability of carbon is analysed in terms of its hybridization states and the substituents.

AB-initio study of alpha-silyl and beta-silyl substituent effects on vinyl radicals

The geometries of alpha- and beta-silyl substituted vinyl radicals and of alpha,beta-disilylvinyl radical have been optimised with the STO-3G and the STO-3G* basis sets. The relative stabilities of various conformers have been determined at the UMP2/6-31G* level. The stabilisation of vinyl radicals through alpha-silyl substitution is larger than that due to corresponding alkyl groups. The presence of an alpha-silyl group also leads to a tendency towards linearisation of the vinyl radical centre and a corresponding reduction in the inversion barrier. In marked contrast, the beta-silyl effect is negligible. The geometric, conformational and energetic consequences are insignificant. Overall, the silyl substituent effect at vinyl radicals is very different from that computed earlier for the vinyl cations, but qualitatively similar to that found in carbanions.

Use of L matrix to obtain reliable ab initio force constants of polyatomic molecules: ethylene as a test system

Direct use of experimental eigenvalues of the vibrational secular equation on to the ab initio predicted eigenvector space is suggested as a means of obtaining a reliable set of intramolecular force constants. This method which we have termed RECOVES (recovery in the eigenvector space) is computationally simple and free from arbitrariness. The RECOVES force constants, by definition, reproduce the experimental vibrational frequencies of the parent molecule exactly. The ab initio calculations were carried out for ethylene as a test molecule and the force constants obtained by the present procedure also correctly predict the vibrational frequencies of the deuterated species. The RECOVES force constants for ethylene are compared with those obtained by using the SQM procedure.

Electron donor-acceptor complex of ICl with diethyl ether - He I photoelectron spectroscopy and ab initio molecular orbital study

The He I photoelectron spectrum of the diethyl ether-ICl complex has been obtained. The oxygen orbitals are shifted to higher binding energies and that of ICl to lower binding energies owing to complex formation. Ab initio molecular orbital (MO) calculations of the complex molecule showed that the bonding is between the sigma-type lone pair of oxygen and the I atom and that the complex has C-2v symmetry. The binding energy of the complex is computed to be 8.06 kcal mol(-1) at the MP2/3-21G* level. The orbital energies obtained from the photoelectron spectra of the complex are compared and assigned with orbital energies obtained by MO calculations. Natural bond orbital analysis (NBO) shows that charge transfer is from the sigma-type oxygen lone pair to the iodine atom and the magnitude of charge transfer is 0.0744 e.