Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Temperature dependence of emission measure in solar X-ray plasmas. 1: Non-flaring active regions

X-ray and ultraviolet line emission from hot, optically thin material forming coronal active regions on the sun may be described in terms of an emission measure distribution function, Phi (T). A relationship is developed between line flux and Phi (T), a theory which assumes that the electron density is a single-valued function of temperature. The sources of error involved in deriving Phi (T) from a set of line fluxes are examined in some detail. These include errors in atomic data (collisional excitation rates, assessment of other mechanisms for populating excited states of transitions, element abundances, ion concentrations, oscillator strengths) and errors in observed line fluxes arising from poorly - known instrumental responses. Two previous analyses are discussed in which Phi (T) for a non-flaring active region is derived. A least squares method of Batstone uses X-ray data of low statistical significance, a fact which appears to influence the results considerably. Two methods for finding Phi (T) ab initio are developed. The coefficients are evaluated by least squares. These two methods should have application not only to active-region plasmas, but also to hot, flare-produced plasmas.

Phillips, K. J. H.↗
NASA NTRS2000

Improvement of the Database on the 1.13-microns Band of Water Vapor

Corrections have recently been reported (Giver et al.) on the short-wave (visible and near-infrared) line intensities of water vapor that were catalogued in the spectroscopic database known as HITRAN. These updates have been posted on www.hitran.com, and are being used to reanalyze the polar stratospheric absorption in the 0.94 microns band as observed in POAM. We are currently investigating additional improvement in the 1.13 microns band using data obtained by us with an absorption path length of 1.107 km and 4 torr of water vapor and the ab initio line list of Partridge and Schwenke (needs ref). We are proposing the following four types of improvement of the HITRAN database in this region: 1) HITRAN has nearly 200 lines in this region without proper assignments of rotational quantum levels. Nearly all of them can now be assigned. 2) We have measured positions of the observable H2O-17 and H2O-18 lines. These lines in HITRAN currently have approximate positions based upon rather aged computations. 3) Some additional lines are observed and assigned which should be included in the database. 4) Corrections are necessary for the lower state energies E" for the HITRAN lines of the 121-010 "hot" band.

Giver, Lawrence P.↗
NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2017

First-Principles-Based Computational Study on Nucleation and Growth Mechanisms of U on Mo(110) Surface Solvated in an Eutectic LiCl-KCl Molten Salt

In nuclear engineering and industry electrochemical methods have been intensively employed in pyro-processing technology to recycle used nuclear fuels for substantial alleviation of the concern with the safe treatment of high-level radioactive nuclear wastes. Pyro-processing is based on sequential electrochemical reactions at high temperature (typically above 700 K at ambient pressure) to enhance the efficiency for decoupling nuclear fuels from unfavorable fission fragments. Using a molten salt in electrorefining process, at which actinides can be selectively deposited in solid cathode from ionic states, gives several advantages over an aqueous solution: higher thermal conductivity (it is liquid phase at high temperature, molten salt plays a role as medium at high temperature over 773 K in electrochemical reaction, lower kinetic barrier in ionic diffusion mobility and better stability against undesired side reactions with radioactive materials. Therefore, the volume of high-level nuclear wastes to be, otherwise, disposed into secured sites is dramatically reduced. Most of all it conforms to nuclear non-proliferation treaty. Electrochemical reactions in a molten salt medium are, however, much more complicated than in aqueous solution rendering the accurate understanding of thermodynamic and kinetic properties formidably challenging. It is represented by extensive researches in experiments and computations to identify underlying mechanisms on how U and other fission fragments behave in the electrorefining system. Regardless of the researches there is a general agreement that the initial stage of the nucleation and growth process of the U deposit is critical for determining the properties of final product and the electrorefining efficiency. In spite of the focused studies it still remains unclear. Theoretical and computational researches have been often carried out with various empirical parameter for interatomic potential energy, which are unknown variables until numerical fits to experimental measurements or just assumption is made. Thus, empirical simulations sometimes predicted outcomes are inconsistent depending on materials of interest with experimental measurements. This may be largely because the accuracy of an empirical interatomic potential energy models critically depends on the input parameters and thus, can be limited in transferability of the results. Furthermore, it may be devoid of fundamental physics or chemistry. This situation requires special caution when an electrorefining system is studied to recycle spent nuclear fuels because of the extreme difficulty in reliable acquisition of thermodynamic and kinetic data from direct experimental measurements. In this presentation, we extensively utilized first principles density functional theory (DFT) calculations and ab-initio molecular dynamic (AIMD) simulations to figure out underlying mechanisms of the initial stage of nucleation and growth behaviors of U on Mo(110) surface in a eutectic LiCl-KCl molten salt at T = 773 K. To rigorously investigate the effect of the molten salt electrolyte we compared the results with those calculated under vacuum condition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Five ab initio potential energy and dipole moment surfaces for hydrated NaCl and NaF. I. Two-body interactions

We report full-dimensional, ab initio-based potentials and dipole moment surfaces for NaCl, NaF, Na{sup +}H{sub 2}O, F{sup −}H{sub 2}O, and Cl{sup −}H{sub 2}O. The NaCl and NaF potentials are diabatic ones that dissociate to ions. These are obtained using spline fits to CCSD(T)/aug-cc-pV5Z energies. In addition, non-linear least square fits using the Born-Mayer-Huggins potential are presented, providing accurate parameters based strictly on the current ab initio energies. The long-range behavior of the NaCl and NaF potentials is shown to go, as expected, accurately to the point-charge Coulomb interaction. The three ion-H{sub 2}O potentials are permutationally invariant fits to roughly 20 000 coupled cluster CCSD(T) energies (awCVTZ basis for Na{sup +} and aVTZ basis for Cl{sup −} and F{sup −}), over a large range of distances and H{sub 2}O intramolecular configurations. These potentials are switched accurately in the long range to the analytical ion-dipole interactions, to improve computational efficiency. Dipole moment surfaces are fits to MP2 data; for the ion-ion cases, these are well described in the intermediate- and long-range by the simple point-charge expression. The performance of these new fits is examined by direct comparison to additional ab initio energies and dipole moments along various cuts. Equilibrium structures, harmonic frequencies, and electronic dissociation energies are also reported and compared to direct ab initio results. These indicate the high fidelity of the new PESs.

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

Reliable electrochemical phase diagrams of magnetic transition metals and related compounds from high-throughput ab initio calculations

Abstract Magnetic transition metals (mTM = Cr, Mn, Fe, Co, and Ni) and their complex compounds (oxides, hydroxides, and oxyhydroxides) are highly important material platforms for diverse technologies, where electrochemical phase diagrams with respect to electrode potential and solution pH can be used to effectively understand their corrosion and oxidation behaviors in relevant aqueous environments. Many previous decades-old mTM–Pourbaix diagrams are inconsistent with various direct electrochemical observations, because experimental complexities associated with extracting reliable free energies of formation (Δ f G ) lead to inaccuracies in the data used for modeling. Here, we develop a high-throughput simulation approach based on density-functional theory (DFT), which quickly screens structures and compounds using efficient DFT methods and calculates accurate Δ f G values, using high-level exchange-correlation functions to obtain ab initio Pourbaix diagrams in comprehensive and close agreement with various important electrochemical, geological, and biomagnetic observations reported over the last few decades. We also analyze the microscopic mechanisms governing the chemical trends among the Δ f G values and Pourbaix diagrams to further understand the electrochemical behaviors of mTM-based materials. Last, we provide probability profiles at variable electrode potential and solution pH to show quantitatively the likely coexistence of multiple-phase areas and diffuse phase boundaries.

Huang, Liang-Feng (ORCID:0000000329375183)↗
DOE OSTI2016

Structural and vibrational study of 2-MethoxyEthylAmmonium Nitrate (2-OMeEAN): Interpretation of experimental results with ab initio molecular dynamics

In this work we report an analysis of the bulk phase of 2-methoxyethylammonium nitrate based on ab initio molecular dynamics. The structural and dynamical features of the ionic liquid have been characterized and the computational findings have been compared with the experimental X-ray diffraction patterns, with infrared spectroscopy data, and with the results obtained from molecular dynamics simulations. The experimental infrared spectrum was interpreted with the support of calculated vibrational density of states as well as harmonic frequency calculations of selected gas phase clusters. Particular attention was addressed to the high frequency region of the cation (ω > 2000 cm{sup −1}), where the vibrational motions involve the NH{sub 3}+ group responsible for hydrogen bond formation, and to the frequency range 1200-1400 cm{sup −1} where the antisymmetric stretching mode (ν{sub 3}) of nitrate is found. Its multiple absorption lines in the liquid arise from the removal of the degeneracy present in the D{sub 3h} symmetry of the isolated ion. Our ab initio molecular dynamics leads to a rationalization of the frequency shifts and splittings, which are inextricably related to the structural modifications induced by a hydrogen bonding environment. The DFT calculations lead to an inhomogeneous environment.

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

Communication: Mode specific quantum dynamics of the F + CHD{sub 3} → HF + CD{sub 3} reaction

The mode specific reactivity of the F + CHD{sub 3} → HF + CD{sub 3} reaction is investigated using an eight-dimensional quantum dynamical model on a recently developed ab initio based full-dimensional potential energy surface. Our results indicate prominent resonance structures at low collision energies and absence of an energy threshold in reaction probabilities. It was also found that excitation of the C–D stretching or CD{sub 3} umbrella mode has a relatively small impact on reactivity. On the other hand, the excitation of the C–H vibration (v{sub 1}) in CHD{sub 3} is shown to significantly increase the reactivity, which, like several recent quasi-classical trajectory studies, is at odds with the available experimental data. Possible sources of the disagreement are discussed.

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

Elucidating the Solution-Phase Structure and Behavior of 8-Hydroxyquinoline Zinc in DMSO

Here, the solution-phase structure and electronic relaxation dynamics of zinc bis-8-hydroxyquinoline [Zn(8HQ) 2 ] in dimethyl sulfoxide (DMSO) were examined using a broad array of spectroscopic techniques, complimented by ab initio calculations of molecular structure. The ground-state structure was determined using extended X-ray absorption fine structure (EXAFS) data collected on the Zn Kedge and diffusion ordered spectroscopy (DOSY) NMR The complex was found to be monomeric and octahedral, with two bidentate 8-hydroxyquinolate ligands and two DMSO molecules coordinated to the zinc through oxygen atoms. Electronic relaxation dynamics were examined with ultrafast transient absorption spectroscopy and complementary density functional calculations. Electronic relaxation was observed to proceed through both singlet and triplet pathways. This solution-phase data provides a deeper physical understanding of the behavior of this molecule, which has a variety of uses such as sensing, OLEDs, and biological applications.

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

Quartic force field-derived vibrational frequencies and spectroscopic constants for the isomeric pair SNO and OSN and isotopologues

The SNO and OSN radical isomers are likely to be of significance in atmospheric and astrochemistry, but very little is known about their gas phase spectroscopic properties. State-of-the-art ab initio composite quartic force fields are employed to analyze the rovibrational features for both systems. Comparison to condensed-phase experimental data for SNO has shown that the 1566.4 cm{sup −1} ν{sub 1} N–O stretch is indeed exceptionally bright and likely located in this vicinity for subsequent gas phase experimental analysis. The OSN ν{sub 1} at 1209.4 cm{sup −1} is better described as the antisymmetric stretch in this molecule and is also quite bright. The full vibrational, rotational, and rovibrational data are provided for SNO and OSN and their single {sup 15}N, {sup 18}O, and {sup 34}S isotopic substitutions in order to give a more complete picture as to the chemical physics of these molecules.

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

Lepton scattering from Ar 40 and Ti 48 in the quasielastic peak region

Neutron and proton spectral functions of 40 Ar, 40 Ca, and 48 Ti isotopes are computed using the ab initio self-consistent Green's function approach. The resulting radii and charge distributions are in good agreement with available experimental data. The spectral functions of Ar and Ti are then utilized to calculate inclusive (e,e') cross sections within a factorization scheme and are found to correctly reproduce the recent Jefferson Lab measurements. Based on these successful agreements, the weak charged and neutral current double-differential cross sections for neutrino- 40 Ar scattering are predicted in the quasielastic region. Finally, results obtained by replacing the (experimentally inaccessible) neutron spectral distribution of 40 Ar with the (experimentally accessible) proton distribution of 48 Ti are compared and the accuracy of this approximation is assessed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Comparative critical mass calculations for NNL and ENDF/B-VIII.0 - Zirconium hydride thermal neutron scattering laws

Zirconium hydride (ZrH{sub x}) is a moderator material for TRIGA reactors and historical space reactor systems, such as SNAP-10A. Thermal neutron scattering laws (TSL) for two phases of this material, δ and ε, have been previously evaluated by Naval Nuclear Laboratory (NNL) and submitted to the National Nuclear Data Center (NNDC) for inclusion in the US national ENDF/B-VIII.1 nuclear data library. In contrast to the current ENDF/B-VIII.0 TSL evaluations, which consider only a single phase, the new evaluations are derived from separate ab initio calculations for both phases and include coherent elastic effects of the zirconium sublattice. To estimate the impact of these changes to the TSL evaluation of this material, comparative critical mass calculations were performed with MC21 for homogenous mixtures of high-enriched uranium (HEU) and ZrH{sub x} in bare and water reflected sphere configurations. These calculations yield an impact on the estimated critical mass as a function of {sup 235}U loading density with maximum differences as large as 1% - 5% for over-moderated thermal spectrum systems. Consequently, the NNL TSL evaluations are anticipated to have a small impact on criticality calculations of thermal reactor systems regardless of the material phase. Nevertheless, characteristic differences exist in the predicted thermal spectra as function of energy for the two sets of TSL evaluations, which are attributed to difference in the underlying phonon density of states of hydrogen bound in ZrH{sub x}. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Modeling non-harmonic behavior of materials from experimental inelastic neutron scattering and thermal expansion measurements

Based on thermodynamic principles, we derive expressions quantifying the non-harmonic vibrational behavior of materials, which are rigorous yet easily evaluated from experimentally available data for the thermal expansion coefficient and the phonon density of states. These experimentally-derived quantities are valuable to benchmark first-principles theoretical predictions of harmonic and non-harmonic thermal behaviors using perturbation theory, ab initio molecular-dynamics, or Monte-Carlo simulations. In this study, we illustrate this analysis by computing the harmonic, dilational, and anharmonic contributions to the entropy, internal energy, and free energy of elemental aluminum and the ordered compound FeSi over a wide range of temperature. Our results agree well with previous data in the literature and provide an efficient approach to estimate anharmonic effects in materials.

thermal expansion↗
DOE OSTI2015

Three-dimensional potential energy surface of Ar–CO

A three-dimensional intermolecular potential energy surface of the Ar–CO complex has been determined by fitting most of the previously reported spectroscopic data, where observed transition frequencies by microwave, millimeter-wave, submillimeter-wave, and infrared spectroscopy were reproduced simultaneously within their experimental accuracies. A free rotor model Hamiltonian considering all the freedom of motions for an atom-diatom system was applied to calculate vibration-rotation energies. A three-dimensional potential energy surface obtained by ab initio calculations at the CCSD(T)-F12b/aug-cc-pV5Z level of theory was parameterized by a model function consisting of 46 parameters. They were used as initial values for the least-squares analysis of the experimental data. A total of 20 parameters were optimized to reproduce all the spectroscopic data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Structural evolution of low-temperature liquid GaIn eutectic alloy

Structural changes of liquid Galn eutectic alloy upon cooling from about 373 K to 242 K are studied by using in situ synchrotron radiation X-ray diffraction (XRD) combined with the ab initio molecular dynamics simulations (AIMD). The experimental and theoretical results show that no abnormal structural change is observed in the supercooled liquid region prior to crystallization. The XRD data exhibit that the liquid can crystallize into a mixture of orthorhombic Ga-I and body-centered tetragonal (bct)-In phases. Also, the AIMD calculations simulate the solidification process, forming a bct-like structure locally similar to bct-Ga-III and bct-In phases. The aggregation of In atoms forms In-rich zones, acting as the processor of crystallization. Our results provide an in-depth understanding of liquid Galn alloys and will trigger more interests to other metallic liquids. (C) 2019 Elsevier B.V. All rights reserved.

36 MATERIALS SCIENCE↗
DOE OSTI2023

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗