Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio phasing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗

Theoretical Calculation of Jet Fuel Thermochemistry: Tetrahydrodicylopentadiene (JP10) Thermochemistry Using the CBS-QB3 and G3(MP2)//B3LYP Methods - 1

High-level ab initio calculations have been performed on the exo and endo isomers of gas-phase tetrahydrodicyclopentadiene (THDCPD), a principal component of the jet fuel JP10, using the Gaussian Gx and Gx(MPx) composite methods, as well as the CBS-QB3 method, and using a variety of isodesmic and homodesmotic reaction schemes. The impetus for this work is to help resolve large discrepancies existing between literature measurements of the formation enthalpy Delta (sub f)H deg (298) for exo-THDCPD. We find that use of the isodesmic bond separation reaction C10H16 + 14CH4 yields 12C2H6 yields results for the exo isomer (JP10) in between the two experimentally accepted values, for the composite methods G3(MP2), G3(MP2)//B3LYP, and CBS-QB3. Application of this same isodesmic bond separation scheme to gas-phase adamantane yields a value for Delta (sub f)H deg (298) within 5 kJ/mol of experiment. Isodesmic bond separation calculations for the endo isomer give a heat of formation in excellent agreement with the experimental measurement. Combining our calculated values for the gas-phase heat of formation with recent measurements of the heat of vaporization yields recommended values for Delta (sub f)H deg (298)liq of -126.4 and -114.7 kJ/mol for the exo and endo isomers, respectively.

Zehe, Michael J.↗

Ab-Initio Description and Prediction of Properties of Carbon-Based and Other Non-Metallic Materials

We have resolved the long-standing problem consisting of 30%-50% theoretical underestimates of the band gaps of non-metallic materials. We describe the Bagayoko, Zhao, and Williams (BZW) method that rigorously circumvents the basis-set and variational effect presumed to be a cause of these underestimates. We present ab-initio, computational results that are in agreement with experiment for diamond (C), silicon (Si), silicon carbides (3C-SiC and 4H-SiC), and other semiconductors (GaN, BaTiO3, AlN, ZnSe, ZnO). We illustrate the predictive capability of the BZW method in the case of the newly discovered cubic phase of silicon nitride (c-Si3N4) and of selected carbon nanotabes [(10,0), and (8,4)]. Our conclusion underscores the inescapable need for the BZW method in ab-initio calculations that employ a basis set in a variational approach. Current nanoscale trends amplify this need. We estimate that the potential impact of applications of the BZW method in advancing our understanding of nonmetallic materials, in informing experiment, and particularly in guiding device design and fabrication is simply priceless.

Bagayoko, D.↗

Structure of 2,4-Diaminopyrimidine - Theobromine Alternate Base Pairs

We report the structure of clusters of 2,4-diaminopyrimidine with 3,7-dimethylxanthine (theobromine) in the gas phase determined by IR-UV double resonance spectroscopy in both the near-IR and mid-IR regions in combination with ab initio computations. These clusters represent potential alternate nucleobase pairs, geometrically equivalent to guanine-cytosine. We have found the four lowest energy structures, which include the Watson-Crick base pairing motif. This Watson-Crick structure has not been observed by resonant two-photon ionization (R2PI) in the gas phase for the canonical DNA base pairs.

Gengeliczki, Zsolt↗

Calculations concerning the reaction C + H3(+) - CH(+) + H2

The ion-molecule reaction of C + H3(+) - CH(+) + H2 has not been studied in the laboratory but is thought to be important in the gas phase synthesis of organic molecules in dense interstellar clouds. This reaction was studied theoretically by performing quantum chemical ab initio calculations on the potential surface. It is found that there is no activation barrier to the reaction and that it proceeds smoothly to the first excited electronic state of CH(+). The rate coefficient as a function of temperature can then be estimated using the proper long-range potentials. The rate coefficient at 10 K is calculated to be 2.9 x 10 to the -9th cu cm/sec.

Talbi, D.↗

Ab initio Quantum Chemical Studies of Reactions in Astrophysical Ices. Reactions Involving CH3OH, CO2, CO, HNCO in H2CO/NH3/H2O Ices

While reactions between closed shell molecules generally involve prohibitive barriers in the gas phase, prior experimental and theoretical studies have demonstrated that some of these reactions are significantly enhanced when confined within an icy grain mantle and can occur efficiently at temperatures below 100 K with no additional energy processing. The archetypal case is the reaction of formaldehyde (H2CO) and ammonia (NH3) to yield hydroxymethylamine (NH2CH2OH). In the present work we have characterized reactions involving methanol (CH3OH), carbon dioxide (CO2), carbon monoxide (CO), and isocyanic acid (HNCO) in search of other favorable cases. Most of the emphasis is on CH3OH, which was investigated in the two-body reaction with one H2CO and the three-body reaction with two H2CO molecules. The addition of a second H2CO to the product of the reaction between CH3OH and H2CO was also considered as an alternative route to longer polyoxymethylene polymers of the -CH2O- form. The reaction between HNCO and NH3 was studied to determine if it can compete against the barrierless charge transfer process that yields OCN(-) and NH4(+). Finally, the H2CO + NH3 reaction was revisited with additional benchmark calculations that confirm that little or no barrier is present when it occurs in ice.

Woon, David E.↗

AB Initio Study of the Structure and Spectroscopic Properties of Halogenated Thioperoxy Radicals

Thioperoxy (XSO or XOS) radicals exist in a variety of chemical environments, and they have as a consequence drawn some interest. HSO, an important species in the chemistry of the troposphere, has been examined both experimentally. The halogenated (X = F, Cl or Br) peroxy species and isovalent thioperoxy species have been studied less, but they too are potentially interesting because oxidized sulfur species and halogen sources are present in the atmosphere. Learning the fate of XSO and XOS radicals is important to understanding the atmospheric oxidation chemistry of sulfur compounds. Of these, FSO and ClSO are particularly interesting because they have been directly detected spectroscopically. Recent studies in our laboratory on the photochemistry of thionyl halides (X2SO; where X = F or Cl) have suggested new ways to generate XSO species. The laser-induced photodissociation of thionyl fluoride, F2SO, at 193 nm and thionyl chloride, ClSO, at 248 nm is characterized by a radical mechanism, X2SO -> XSO + X. The structure of FSO has been characterized experimentally by Endo et cd. employing microwave spectroscopy. Using the unrestricted Hartree-Fock (UHF) self-consistent field (SCF) method, Sakai and Morokuma computed the electronic structure of the ground (sup 2)A" and the first excited (sup 2)A' states of FSO. Electron correlation was not taken into account in their study. In a laser photodissociation experiment, Huber et al. identified ClSO mass spectromctrically. ClSO has also been detected in low temperature matrices by EPR and in the gas phase by far IR laser magnetic resonance. Although the structure of FSO is known in detail, the only study, experimental or theoretical, of CISO has been an ab initio HFSCF study by Hinchliffe. Electron correlation corrections were also excluded from this study. In order to better understand the isomerization and dissociation dynamics of the radical species, we have performed ab initio correlated studies of the potential energy surfaces (PES) of ClSO and its isomer ClOS at the QCISD(T)/6-31 G* level of theory. For FSO and FOS, more extensive QCISD/6-31 1G(2df) calculations have been possible, and the results are summarized here.

Munoz, Luis A.↗

Calculations on the rate of the ion-molecule reaction between NH3(+) and H2

The rate coefficient for the ion-molecule reaction NH3(+) + H2 yields NH4(+) + H has been calculated as a function of temperature with the use of the statistical phase space approach. The potential surface and reaction complex and transition state parameters used in the calculation have been taken from ab initio quantum chemical calculations. The calculated rate coefficient has been found to mimic the unusual temperature dependence measured in the laboratory, in which the rate coefficient decreases with decreasing temperature until 50-100 K and then increases at still lower temperatures. Quantitative agreement between experimental and theoretical rate coefficients is satisfactory given the uncertainties in the ab initio results and in the dynamics calculations. The rate coefficient for the unusual three-body process NH3(+) + H2 + He yields NH4(+) + H + He has also been calculated as a function of temperature and the result found to agree well with a previous laboratory determination.

Herbst, Eric↗

Theoretical Studies of the Extra-terrestrial Chemistry of Biogenic Elements and Compounds

Results are presented on the following:(A) Ab initio quantum chemical studies of reactions in astrophysical ices.Theoretical electronic structure calculations were used to investigate reactions between formaldehyde (H2CO) and both hydrogen cyanide (HCN) and isocyanide (HNC) in search of other favorable reactions such as ammonia-formaldehyde addition, which was found in a recent theoretical study to be strongly enhanced when it occurs within cold ices.The present study examines further reactions between this product and H2CO in ices.(B) Heterogeneous hydrogenation of CO and H2CO on icy grain mantles.Formaldehyde (H2CO) and methanol (CH30H) are thought to be produced in the interstellar medium by the successive hydrogenation of carbon monoxide (CO) on grain surfaces. In the gas phase, the steps in which H adds to CO and H2CO possess modest barriers and are too inefficient to account for the observed abundances. Recent laboratory work has confirmed that formaldehyde and methanol are formed when H atoms are deposited on CO ice at 12 K. The present study employed ab initio quantum chemical calculations to investigate the impact of water ice on the sequential hydrogenation of CO.(C) Ice-bound condensed-phase reactions involving formic acid (HCOOH), methylenimine (CH2NH), hydrogen cyanide (HCN), hydrogen isocyanide (HNC), and ammonia ( 3) were investigated in order to characterize possible pathways to larger organic species that are efficient at the cold temperatures prevalent in cometary nuclei and the interstellar medium. (D) Pathways to glycine and other amino acids in ultraviolet-irradiated ices determined via quantum chemical modeling.(E) Photoionization in ultraviolet processing of astrophysical ice analogs at cryogenic temperatures.

Woon, David E.↗

The molecular composition of dense interstellar clouds

Presented in this paper is an ab initio chemical model for dense interstellar clouds that incorporates 598 grain surface reactions, with small grains providing the reaction area. Gas-phase molecules are depleted through collisions with grains. The abundances of 372 chemical species are calculated as a function of time and are found to be of sufficient magnitude to explain most observations. Peak abundances are achieved on time scales of the order of 100,000 to 1 million years, depending on cloud density and kinetic temperature. The reaction rates for ion-molecule chemistry are approximately the same, indicating that surface and gas-phase chemistry may be coupled in certain regions. The composition of grain mantles is shown to be a function of grain radius. In certain grain-size ranges, large molecules containing two or more heavy atoms are more predominant than lighter 'ices' - H2O, NH3, and CH4. It is possible that absorption due to these large molecules in the mantle may contribute to the observed 3-micron band in astronomical spectra.

Allen, M.↗

A reinvestigation of the rate of the C/+/ + H2 radiative association reaction

It is noted that new experimental results and statistical theories have prompted a reinvestigation of the rate coefficient of the interstellar reaction in which C(+) + H2 yields CH2(+) + h(nu) in the 10-100 K temperature range. The results presented here indicate a rate coefficient between 10 to the -16 and 10 to the -15 cu cm/s at all temperatures studied. In applying the modified thermal and phase space theories, it is expected that they will be as accurate as for radiative association as for three-body association, provided the ab initio value for the radiative decay rate is correct. It is expected that the calculated values of the rate coefficient will be accurate to within an order of magnitude and will have the correct temperature dependence.

Herbst, E.↗

The rate of the reaction between CN and C2H2 at interstellar temperatures

The rate coefficient for the important interstellar reaction between CN and C2H2 has been calculated as a function of temperature between 10 and 300 K. The potential surface for this reaction has been determined through ab initio quantum chemical techniques; the potential exhibits no barrier in the entrance channel but does show a small exit channel barrier, which lies below the energy of reactants. Phase-space calculations for the reaction dynamics, which take the exit channel barrier into account, show the same unusual temperature dependence as determined by experiment, in which the rate coefficient at first increases as the temperature is reduced below room temperature and then starts to decrease as the temperature drops below 50-100 K. The agreement between theory and experiment provides strong confirmation that the reaction occurs appreciably at cool interstellar temperatures.

NASA Discipline Exobiology↗

Toward an In-Depth Material Model for Cermet Nuclear Thermal Rocket Fuel Elements

The development and qualification of nuclear thermal propulsion (NTP) fuel element technologies would be aided by an in-depth model of material response and failure modes at operating conditions. Integrated computational materials engineering techniques have the potential to provide such a model, as demonstrated here through three case studies focused on a tungsten-uranium mononitride cermet fuel. The first case focuses on the erosion of tungsten (W, also named wolfram), a nominal coating/cladding material, in hot hydrogen. Ab initio techniques are used to calculate erosion rates and thermal expansion at NTP operating conditions. The second focuses on the stability of uranium mononitride (UN) fuels at high temperature and in the presence of hydrogen. Phase diagram techniques reveal potential instabilities and decomposition pathways at high hydrogen concentrations. The third focuses on using microstructure information to predict high temperature mechanical response and failure of tungsten, used in refractory cermet materials. Combined finite element and discrete dislocation dynamics techniques provide mechanical properties in agreement with experimental methods. The integration of these techniques for an all-encompassing material model is discussed.

cermet↗

Predicting Melt Properties Using Atomistic Simulations With A Highly Accurate Physically Informed Neural Network Interatomic Potential

The use of a recently developed machine learning (ML) interatomic potential for molecular dynamics simulations of aluminum melt properties will be presented. Such properties are critical for process modeling in additive manufacturing, including the melt pool size, solidification, and formation of solidification microstructures. Direct first-principles modeling of these processes is computationally prohibitive whereas simulations employing ML potentials combine the high accuracy of quantum-mechanical methods with high computational speeds. The physically-informed neural network (PINN) method used herein, integrates a high-dimensional regression implemented by an artificial neural network with a physics-based bond-order interatomic potential. PINN potentials can accurately reproduce many properties of aluminum in both crystalline-solid and liquid phases. We examine the accuracy of a PINN Al potential in predicting the density, self-diffusivity, viscosity, and the tension of the liquid surface and liquid-solid interfaces. Comparison with experimental data and ab initio molecular dynamics calculations shows very good agreement for all properties tested.

molecular dynamics↗

Ab initio Screening of Refractory Nitrides and Carbides for High Temperature Hydrogen Permeation Barriers

Density functional theory was used to screen eleven refractory materials – two pure metals, six nitrides, and three carbides–as high-temperature hydrogen permeation barriers to prevent hydrogen embrittlement. Activation energies were calculated for atomic hydrogen (H) diffusion into the first subsurface layer from the lowest energy surface of the high-temperature phase of each candidate material. The candidate barrier materials with the highest activation energies are h-BN, c-BN, HfN, and ZrN with predicted barriers of 3.25 eV, 3.23 eV, 3.14 eV, and 2.76 eV, respectively. Strain energies, Bader charges, and density of states were calculated for the diffusing H at the relaxed initial state and the transition state to provide insight into contributing factors to high energy barriers. The diffusing H atom in materials with the highest predicted barriers are protic. In addition, interstitial H atoms induce mid-gap states in the density of states of both BN polymorphs. The nitrogen retention of each nitride material at high temperatures was predicted using nitrogen vacancy formation energies with respect to gaseous nitrogen. Experimental evaluation of nitrogen retention in h-BN, ZrN, and TiN confirmed their resistance to nitrogen loss at 1773 K. However, of these nitrides, TiN is predicted to be the least stable. This work identifies multiple promising materials that are predicted to be effective hydrogen barriers at high temperatures and that are stable at temperatures above 2700 K, with BN predicted to perform best.

Density Functional Theory (DFT)↗

The Vaporization of B2O3(l) to B2O3(g) and B2O2(g)

The vaporization of B2O3 in a reducing environment leads to formation of both B2O3(g) and B2O2(g). While formation of B2O3(g) is well understood, many questions about the formation of B2O2(g) remain. Previous studies using B(s) + B2O3(l) have led to inconsistent thermodynamic data. In this study, it was found that after heating, B(s) and B2O3(l) appear to separate and variations in contact area likely led to the inconsistent vapor pressures of B2O2(g). To circumvent this problem, an activity of boron is fixed with a two-phase mixture of FeB and Fe2B. Both second and third law enthalpies of formation were measured for B2O2(g) and B2O3(g). From these the enthalpies of formation at 298.15 K are calculated to be -479.9 +/- 41.5 kJ/mol for B2O2(g) and -833.4 +/- 13.1 kJ/mol for B2O3(g). Ab initio calculations to determine the enthalpies of formation of B2O2(g) and B2O3(g) were conducted using the W1BD composite method and show good agreement with the experimental values.

Jacobson, Nathan S.↗

Molecular Simulations of Ultrafast Radiation Induced Melting at Metal-Semiconductor Interfaces

Understanding radiation induced ultrafast melting at material interfaces is essential in designing robust electronic devices for aviation/space applications and in laser machining. While it is difficult to achieve the spatial and temporal resolution required to quantify the phenomenon experimentally, simulations can provide the detailed mechanisms of the structural changes that happen during phase transition. In this work, we use molecular simulations to study the effect of radiation damage on silicon carbide (SiC) - tungsten (W) interfaces which is of interest in high power electronics. A multi-scale approach is involved wherein the reactions at the interfaces are quantified using ab-initio molecular dynamics (MD) simulations and classical MD simulations are employed to understand the structural and diffusional changes across the material interface. Finally, coarse-grained Lennard-Jones type models are used to study the larger scale mechanisms and structures obtained due to the induced damages. We show that the response of the material to radiation damage depends on factors such as energy of the incident radiation, thermal properties, and molecular structure of the material.

Ravichandran, Ashwin↗

Li-Doped Ionic Liquid Electrolytes: From Bulk Phase to Interfacial Behavior

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.

molecular dynamics↗