Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Microscopic Description of the Mechanism of Transition between the 2H and 4H Polytypes of Silicon Carbide

Abstract—: The mechanism of displacement of one close-packed SiC layer from one minimum position to another on the example of SiC polytype transition 2H → 4H has been studied by ab initio methods. It has been shown that the intermediate state with monoclinic symmetry Cm greatly facilitates this displacement breaking it into two stages. Initially, the Si atom chiefly moves, only then—mainly the C atom. In this case, the Si–C bond is significantly tilted in comparison with the initial position, which allows the reducing of the compression of the SiC bonds in the (11 2-bar 0) plane. Two transition states of this process, which also possess the Cm symmetry, have been computed. It has been found that the height of the activation barrier of the process of moving the close-packed layer of SiC from one position to another is equal to 1.8 eV. The energy profile of this movement has been calculated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Techniques for Polytypic Transformations in Silicon Carbide

Abstract—: Two main polytype transformations in silicon carbide, namely, 2H → 6H and 3C → 6H, have been studied by ab initio methods. It has been shown that the intermediate phases with trigonal symmetry P3m1 and monoclinic symmetry Cm make it much easier to move the close-packed layers in such transitions by breaking them up into separate stages. It has been found that these two polytype transformations proceed in completely different ways. The links being relocated noticeably tilted compared to their initial position at the transition 2H → 6H, which allows the compression of the SiC links in the plane (11 2-bar 0). The transition 3C → 6H is carried out through the formation of Si–Si and C–C auxiliary links, living for a short time and helping densely packed layers to swap places. As a result, the activation barrier of the transformation 2H → 6H (1.7 eV/atom) is significantly less than the activation barrier of the transformation 3C → 6H (3.6 eV/atom), which means that the second transition should occur at the temperatures by 750–800°C higher than the first one. The energy profiles of this polytypic transformations, as well as the geometry of all intermediate and transition phases have been calculated. It has been shown that all the transition states have monoclinic symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

The Effect of Hydrogen on Fluctuation Embrittlement of Aluminum

In this paper, we describe the main processes occurring during the generation of vacancies in aluminum in the presence of hydrogen by ab initio methods using the SCAN functional. Hydrogen is shown to reduce the generation energy of vacancies from 2.8 to 0.8 eV. In this case, the eight hydrogen atoms located in the tetrahedral lattice voids around single aluminum atom greatly facilitate its movement to the interstitial site. The dependence of the activation energy of hydrogen embrittlement of aluminum on the hydrogen concentration in aluminum and temperature is calculated based on the kinetic strength concept. Hydrogen is shown to reduce the time of aluminum destruction only if its mole fraction in aluminum is greater than the critical level (~ 3 × 10{sup –4} at T = 293 K).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

A Theoretical Investigation of the Structure and Reactivity of the Molecular Constituents of Oil Sand and Oil Shale

We used a variety of small organic models of asphaltenes to investigate the molecular mechanism for the high temperature decomposition that would take place as part of the oil refinery process. We determined that the decomposition is initiated via four different types of hydrogen migration reactions. According to the energetics of the reactions, the dominant 1,2-H shift mechanism involves two competitive product channels, namely, C 2 H 2 + CH 2 CS and CS + CH 3 CCH. The minor channels include the formation of CS + CH 2 CCH 2 , H 2 S + C 4 H 2 , HCS + CH 2 CCH, CS + CH 2 CHCH, H + C 4 H 3 S, and HS + C 4 H 3 . We also investigated the alkyl substitution effect by exploring the decomposition pathways of models with alkyl arms. The energetics of such systems were very similar to that for unsubstituted model compounds, which suggests that asphaltene alkylation may not play a significant role in the decomposition of asphaltene compounds. This work was published in the Journal of Physical Chemistry A 2011, 115, 2882-2891. A MECHANISTIC STUDY OF THE 2-THIENYLMETHYL + HO2 RADICAL RECOMBINATION REACTION Radicals are molecules which contain single electrons. They are very reactive. Radical recombination reactions are important in the combustion of fuel oils. Shale oil contains radicals. We used quantum mechanics to explore the reactivity of shale oil model radical compounds. Seventeen product channels corresponding to either addition/elimination or direct hydrogen abstraction were characterized. Direct hydrogen abstraction proceeds via a weakly bonded complex, which leads to 2-methylthiophene, 2-methylene-2,3-dihydrothiophene or 2-methylene-2,5-dihydrothiophene depending upon the 2-thienylmethyl radical reaction site. The addition pathway for the two radical reactants is barrierless with the formation of three adducts, as distinguished by HO 2 reaction at three different sites on the 2-thienylmethyl radical. The addition is exothermic by 37 ~ 55 kcal mol-1 relative to the entrance channel. These excess energies are available to promote further decomposition or rearrangement of the adducts that lead to nascent products such as H, OH, H 2 O and CH 2 O. The reaction surfaces are characterized by relatively low barriers (most are lower than 10 kcal mol-1). Based upon a careful analysis of the overall barrier heights and reaction exothermicities, the formation of O2, OH and H2O is likely to be an important pathway in the radical recombination reactions of 2-thienylmethyl + HO 2 . This work was published in the Journal of Physical Chemistry A, 2011, 115, 14546-14557. REACTION OF THIOPHENE AND METHYLTHIOPHENE WITH SINGLET AND TRIPLET MOLECULAR OXYGEN Mechanisms for the reaction of thiophene and 2-methylthiophene with molecular oxygen on both the triplet and singlet potential energy surfaces (PESs) were investigated using ab initio methods. Thiophene and 2-methylthiophene where shown to react with O 2 via two types of mechanisms; namely, direct hydrogen abstraction and addition/elimination. The barriers for reaction with triplet oxygen are all significantly large (i.e., > 30 kcal mol-1), which indicates that the direct oxidation of thiophene by ground state oxygen might be important only in high temperature processes. Reaction of thiophene with singlet oxygen via a 2+4 cycloaddition leading to endoperoxides is the most favorable channel. Moreover, it was found that alkylation of the thiophene ring (i.e., methyl-substituted thiophene) is capable of lowering the barrier height for the addition pathway. The implication of the current theoretical results may shed new light on the initiation mechanisms for combustion of asphaltenes. This work was published in the Journal of Physical Chemistry A, 2012 116, 4934-4946. JAHN-TELLER STABILIZATION IN POSS CATIONS We have a long standing interest in polyoligomeric silsesquioxane (POSS) molecules. 1-2 These molecules have recently been used as advanced surface coatings for photovoltaic devices and have potential as molecular-based energy storage devices as well as magnetically controllable liquid marbles. 3-5 We have been investigating the small molecule encapsulation properties of POSS and discovered some interesting symmetry breaking processes that need to be better understood in order to use POSS in advanced materials. We have investigated this symmetry breaking mechanism in POSS monocations Si8O12(C(CH3)3)8+ and Si8O12Cl8+, using density functional theory (DFT) and group theory. Under Oh symmetry, these ions possess 2T2g and 2Eg electronic states, respectively, and undergo different symmetry breaking mechanisms. The ground states of Si 8 O 12 (C(CH 3 ) 3 ) 8 + and Si 8 O 12 Cl 8 + belong to the C 3v and D 4h point groups and are characterized by Jahn-Teller stabilization energies of 3959 and 1328 cm-1, respectively, at the B3LYP/def2-SVP level of theory. The symmetry distortion mechanism in Si 8 O 12 Cl 8 + is Jahn-Teller type, whereas in Si 8 O 12 (C(CH 3 ) 3 ) 8 + the distortion is a combination of both Jahn-Teller and pseudo-Jahn-Teller effects. The distortion force acting in Si 8 O 12 (C(CH 3 ) 3 ) 8 + is mainly localized on one Si-(tert-butyl) group while in Si 8 O 12 Cl 8 + it is distributed over the oxygen atoms. The main distortion forces acting on the Si8O12 core arise from the coupling between the electronic state and the vibrational modes; identified as 9t 2g +1e g +3a 2u for the Si 8 O 12 (C(CH 3 ) 3 ) 8 + and 1e g +2e g for Si 8 O 12 Cl 8 + . This work was published in the Journal of Physical Chemistry A, 2015, 119, 4237-4243.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

A Theoretical Analysis on the Dynamic Behavior and Electromechanical Response of PZT

Piezoelectric materials are used as a power source in stockpile components for safety and reliability assurances. The objective of this project is to gain insights into the fundamental pressure-induced behavior and electromechanical response of lead zirconate titanate (PZT). Specifically, to establish the basis for an accurate, physics-based material model for robust model-based design to rapidly optimize PZT-based materials and components for performance studies. Established ab-initio methods are used to interrogate and understand the dynamic behavior of PZT as a function of composition (50/50, 65/35, 80/20, 95/5) and dopants (La, Nb) to overcome the costly and time-consuming experimental methodologies. New cold curves for pure and doped single-crystal PZT are obtained as the reference equation-of-state (EOS). A negligible change in the pressure responses was observed for the systems and strains studied. Dielectric and piezoelectric responses of pure and doped single-crystal PZT were also calculated as a function of pressure. For undoped PZT, there is a clear and orderly pressure and composition dependence. For doped PZT, there is a significant increase in the responses, but the behavior is very disordered and inconclusive.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Studying electroweak few-body observables in chiral effective field theory

The use of nuclei to study electroweak probes is becoming increasingly relevant experimentally. The success of dark matter and neutrino experiments strongly depends on the ability to control nuclear effects in order to extract the fundamental parameters associated with external probes. Therefore, reliable theoretical calculations of nuclear structure and reactions, with well-controlled errors, are crucial for the success of experimental efforts. Currently, chiral effective field theory ($\chi$EFT) coupled with {\it ab-initio} methods represents one of the best approaches that fulfills these requirements. To use this approach as a tool for studying fundamental physics, it is essential to validate it against experimental data for which the calculations are well under control, such as the elastic scattering of electrons on nuclei. In this proceeding, I will present recent developments in the fitting of electromagnetic currents derived using $\chi$EFT and the calculation of electromagnetic form factors of light nuclei. The results of these calculations demonstrate the strength of the theory in describing the interaction of nuclei with electromagnetic probes over a broad range of momentum transfers and highlight the robustness of $\chi$EFT for analyzing future experimental data aimed at extracting fundamental parameters.

Gnech, Alex [Old Dominion Univ., Norfolk, VA (Unit↗
DOE OSTI2021

Quantum Calculations of VX Ammonolysis and Hydrolysis Pathways via Hydrated Lithium Nitride

Recently, lithium nitride (Li 3 N) has been proposed as a chemical warfare agent (CWA) neutralization reagent for its ability to produce nucleophilic ammonia molecules and hydroxide ions in aqueous solution. Quantum chemical calculations can provide insight into the Li 3 N neutralization process that has been studied experimentally. Here, we calculate reaction-free energies associated with the Li3N-based neutralization of the CWA VX using quantum chemical density functional theory and ab initio methods. We find that alkaline hydrolysis is more favorable to either ammonolysis or neutral hydrolysis for initial P-S and P-O bond cleavages. Reaction-free energies of subsequent reactions are calculated to determine the full reaction pathway. Notably, products predicted from favorable reactions have been identified in previous experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

How Accurate Can Crystal Structure Predictions Be for High-Energy Molecular Crystals?

Molecular crystals have shallow potential energy landscapes, with multiple local minima separated by very small differences in total energy. Predicting molecular packing and molecular conformation in the crystal generally requires ab initio methods of high accuracy, especially when polymorphs are involved. We used dispersion-corrected density functional theory (DFT-D) to assess the capabilities of an evolutionary algorithm (EA) for the crystal structure prediction (CSP) of well-known but challenging high-energy molecular crystals (HMX, RDX, CL-20, and FOX-7). While providing the EA with the experimental conformation of the molecule quickly re-discovers the experimental packing, it is more realistic to start instead from a naïve, flat, or neutral initial conformation, which reflects the limited experimental knowledge we generally have in the computational design of molecular crystals. By doing so, and using fully flexible molecules in fully variable unit cells, we show that the experimental structures can be predicted in fewer than 20 generations. Nonetheless, one must be aware that some molecular crystals have naturally hindered evolutions, requiring as many attempts as there are space groups of interest to predict their structures, and some may require the accuracy of all-electron calculations to discriminate between closely ranked structures. To save resources in this computationally demanding process, we showed that a hybrid xTB/DFT-D approach could be considered in a subsequent study to push the limits of CSP beyond 200+ atoms and for cocrystals.

42 ENGINEERING↗
DOE OSTI2016

SPECTROSCOPIC CONSTANTS OF THE X{sup 1}Σ{sup +} AND 1{sup 3}Π STATES OF AlO{sup +}

Using both standard and explicitly correlated ab initio methods in conjunction with several atomic basis sets, the ground state of AlO(X{sup 2}Σ{sup +}) and the two lowest electronic states of AlO{sup +} ({sup 1}Σ{sup +} and {sup 3}Π) are investigated. Potential energy curves for these species are mapped, which are incorporated later to solve the nuclear motion problem. Benchmark computations on AlO(X{sup 2}Σ{sup +}) are used to determine the reliability of the theoretical methods and basis sets used for an accurate description of aluminum oxide compounds. The electronic ground state of AlO{sup +} is X{sup 1}Σ{sup +}, followed by the low-lying 1{sup 3}Π state. For both cationic electronic states, a set of spectroscopic parameters are recommended that may help in the identification of this ion in laboratory and astrophysical media. An accurate estimation of the adiabatic ionization energy of AlO, AIE = 9.70 eV, is also reported.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Pressure Distortion of the H{sub 2}–He Collision-induced Absorption at the Photosphere of Cool White Dwarf Stars

Collision-induced absorption (CIA) from molecular hydrogen is a dominant opacity source in the atmosphere of cool white dwarfs. It results in a significant flux depletion in the near-IR and IR parts of their spectra. Because of the extreme conditions of helium-rich atmospheres (where the density can be as high as a few g cm{sup −3}), this opacity source is expected to undergo strong pressure distortion and the currently used opacities have not been validated at such extreme conditions. To check the distortion of the CIA opacity, we applied state-of-the-art ab initio methods of computational quantum chemistry to simulate the CIA opacity at high densities. The results show that the CIA profiles are significantly distorted above densities of 0.1 g cm{sup −3} in a way that is not captured by the existing models. The roto-translational band is enhanced and shifted to higher frequencies as an effect of the decrease of the interatomic separation of the H{sub 2} molecule. The vibrational band is blueward shifted and split into Q {sub R} and Q {sub P} branches, separated by a pronounced interference dip. Its intensity is also substantially reduced. The distortions result in a shift of the maximum of the absorption from 2.3 μm to 3–7 μm, which could potentially explain the spectra of some very cool, helium-rich white dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Radiative lifetimes of the second negative system of O2(+)

The dipole moment of the A2Pi(u)-X2Pi(g) transition of O2(+) is calculated as a function of internuclear distance using ab initio methods. The band absorption oscillator strengths and band transition probabilities of the second negative system are derived and the resulting lifetimes are compared with experimental data. The high-lying v double prime levels of the ground state may decay into low-lying v prime levels of the excited state. The corresponding radiative lifetimes are calculated.

Wetmore, R. W.↗
NASA NTRS1989

On the binding energy of He(n)+, for n = 2-7

The atomization energy of small positively charged He clusters has been studied using ab initio methods that include correlation. The atomization energy increases monotonically with cluster size. The inclusion of correlation changes the nature of the bonding and hence changes the structure of the most stable cluster.

Rosi, Marzio↗
NASA NTRS1992

Ab initio study of the molecular structure and vibrational spectrum of nitric acid and its protonated forms

The equilibrium structures, harmonic vibrational frequencies, IR intensities, and relative energetics of HNO3 and its protonated form H2NO3+ were investigated using double-zeta plus polarization and triple-zeta plus polarization basis sets in conjunction with high-level ab initio methods. The latter include second-order Moller-Plesset perturbation theory, the single and double excitation coupled cluster (CCSD) methods, a perturbational estimate of the effects of connected triple excitations (CCSD(T)), and the self-consistent field. To determine accurate energy differences CCSD(T) energies were computed using large atomic natural orbital basis sets. Four different isomers of H2NO3+ were considered. The lowest energy form of protonated nitric acid was found to correspond to a complex between H2O and NO2+, which is consistent with earlier theoretical and experimental studies.

Lee, Timothy J.↗
NASA NTRS1993

Calculations on the competition between association and reaction for C3H(+) + H2

The ab initio results presently obtained for the potential energy surface of competing associative and reactive ion-molecule processes of the reactants H2 and C3H(+) show that the linear ion is able to directly access the deep potential well of the propargyl ion H2CCCH(+). Using the minimum energy potential pathway and properties of the stationary point structures determined via the ab initio methods, calculation results are obtained for both the association-rate coefficient for C3H3(+) production as a function of density, and the branching ratio between the propargyl and cyclic structures of the ion.

Maluendes, Sergio A.↗
NASA NTRS1996

Optimized Structures and Proton Affinities of Fluorinated Dimethyl Ethers: An Ab Initio Study

Ab initio methods have been used to investigate the proton affinity and the geometry changes upon protonation for the molecules (CH3)2O, (CH2F)2O, (CHF2)2O, and (CF3)2O. Geometry optimizations were performed at the MP2/3-2 I G level, and the resulting geometries were used for single-point energy MP2/6-31G calculations. The proton affinity calculated for (CH3)2O was 7 Kjoule/mole from the experimental value, within the desired variance of +/- 8Kjoule/mole for G2 theory, suggesting that the methodology used in this study is adequate for energy difference considerations. For (CF3)20, the calculated proton affinity of 602 Kjoule/mole suggests that perfluorinated ether molecules do not act as Lewis bases under normal circumstances; e.g. degradation of commercial lubricants in tribological applications.

Orgel, Victoria B.↗
NASA NTRS1994

The Fraternal Twins of Quartet O4(+)

Eleven stationary geometries of quartet O4(+) have been studied by ab initio methods. The geometries were optimized at the Complete Active Space Self-Consistent Field (CASSCF) level of theory and the energies were calculated by the multiconfigurational second order pertubation method (CASPT2), using Double-Zeta Plus polarization (DZP), Triple-Zeta Plus Double Polarization (TZ2P), average Atomic Natural Orbital (ANO) (5s4p2d) and average ANO (6s5p3d2f) basis sets. The rectangular and trans-planar structures are found to be the most stable, with an energy barrier to conversion between the two at the threshold of dissociation. Both have a delocalized hole and are stable relative to separated 02 and 02(+)by 11.0 and 11.5 kcal/mol for the rectangular and the 2 trans-planar structure, respectively, compared with the experimentally deduced energy in the range of 9.2 to 10.8 kcal/mol. The adiabatic ionization potentials of 04 and 02 are computed to be 11.67 and 12.21 eV, while experimental values are 11.66 and 12.07 eV, respectively. The vibrational frequencies have been computed for all degrees of freedom at the CASSCF level of theory. Symmetry breaking is found to be a particular problem in the computation of the antisymmetric stretch frequency for the delocalized structures at the CASSCF level of theory. Attempts to rectify these problems using the restricted active space self-consistent field (RASSCF) method leads to additional difficulties, but further analysis yields insight into the symmetry breaking and problems with earlier calculations. Finally, a nonorthogonal configuration interaction (CI) calculation based on the interaction of localized CASSCF wave functions using the Complete Active Space State Interation (CASSI) method leads to a balanced treatment of the antisymmetric stretch which is free from symmetry breaking. The study explains the four most prominent absorption frequencies observed in the partially unassigned IR spectrum of O4(+) isolated in solid neon as the antisymmetric OO stretch, and the combination band of the symmetric and antisymmetric 00 stretch of both the rectangular and trans-planar structures.

Lindh, Roland↗
NASA NTRS1996

AB Initio Characterization of MgCCH, MgCCH(+), and MgC2, and Pathways to their Formation in the Interstellar Medium

A study of Mg-bearing compounds has been performed in order to determine molecular properties which are critical for planning new astronomical searches and laboratory studies. The primary focus of the work is on MgCCH, MgCCH(+), and the isomers of MgC2. Only MgCCH has been identified in laboratory studies. Additional calculations have been carried out on MgH, MgNC, MgCN, and their cations in an effort to evaluate pathways to the formation of MgCCH and MgCCH(+) in the InterStellar Medium (ISM) or in circumstellar envelopes. Correlated ab initio methods and correlation-consistent basis sets have been employed. Properties including structures, rotational constants, dipole moments, and harmonic frequencies are reported. A transition state between linear MgCC and cyclic MgC2 has been characterized and was found to yield a minimal barrier (approx. 0.5 kcal/mole), indicating easy interconversion to the cyclic form. Direct reactions in the ISM between Mg or Mg(+) and HCCH are precluded by energetic considerations, but a number of ion- molecule or neutral-neutral exchange reactions between CCH and various Mg-containing species offer plausible pathways to MgCCH or MgCCH(+). Weakly bound MgH may react with CCH to form MgCCH, but MgH has not been detected. Both MgNC and MgCN have been observed, but reactions with CCH are slightly endothermic by 1-3 kcal/mole. Although MgH(+), MgNC(+), and MgCN(+) have not been detected, their reactions with CCH to form MgCCH(+) are all exothermic. With only a small barrier separating linear MgCC and cyclic MgC2, the dissociative recombination of MgCCH(+) with an electron is expected to yield cyclic MgC2, and regenerate Mg and CCH. New astronomical searches for MgCCH, MgCCH(+), cyclic MgC2, MgNC(+), and MgCN(+) will provide further insight into organo-magnesium astrochemistry.

Woon, David E.↗
NASA NTRS1997

Can Chlorine Anion Catalyze the Reaction fo HOCl with HCl?

The reaction of HOCl + HCl -> Cl2 + H20 in the presence of Cl has been studied using ab initio methods. This reaction has been shown to have a high activation barrier of 46.5 kcal/mol. The chlorine anion, Cl- is found to catalyze the reaction, viz. two mechanisms. The first involves Cl- interacting through the concerted four-center transition state of the neutral reaction. The other mechanism involves the formation of a HCl-HOCl-Cl- intermediate which dissociates into Cl2 + Cl- + H20. The steps are found to have no barriers. The overall exothermicity is 15.5 kcal/mol.

Richardson, S. L.↗