Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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 343 records · Page 19

Interpretation of experimental differential elastic scattering cross section for H/+/ + Ne.

The interatomic V(r) found by potential-model calculations is compared to an ab initio calculation of the intermolecular potential for NeH+ due to Peyerimhoff (1965). The reason for a rather significant difference in the two results is discussed, and a new method due to Remler (1971) involving Regge poles is applied. The method can be used to calculate the differential cross section, where the starting point for such calculations is based on intuition in light of the classical deflection function.

Bobbio, S. M.↗

Molten Salt Electrolytes: Computational Analysis of Transport, Electrochemical Properties and Designing New Mixtures

Electric aircraft propulsion has gained a traction over the last decade due to possible high-energy density electrochemistries with reasonable cycle life and identification of non-flammable chemistries that can guarantee much better safety. Commercial Li-ion batteries cannot reach such high capacities because of weight limitations that arise from the widely used intercalation electrodes. Also, use of organic electrolytes make them highly susceptible to fire on exposure to air and humidity. Currently, use of Lithium metal anode along with conversion electrochemistries such as Li-O2 and Li-S are being pursued to achieve such high energy densities. Safer inorganic solid-state electrolytes, recently discovered Water-in-Salt electrolytes, molten salt electrolytes are some of the alternatives being pursued for a safer/non-flammable battery. Of these safer alternatives, molten-salt electrolytes offer some attractive properties: liquid at operating temperatures resulting in better electrode-wetting, stable interface with Li-metal anodes and good ionic conductivity; their only drawback being high operating temperatures. In this talk, we will discuss some of the computational studies, driven via atomistic simulations, of the properties of molten-salt electrolytes including transport mechanism and interface stability. We will also discuss predicting electrochemical properties of these high-temperature electrolytes. Further, we will discuss a thermodynamic approach to predicting new molten-salt eutectics with lower melting points.

Molten Salt Electrolytes↗

Actinide opacities for modelling the spectra and light curves of kilonovae

ABSTRACT We extend previous ab initio calculations of lanthanide opacities to include a complete set of actinide opacities for use in the modelling of kilonova (KN) light curves and spectra. Detailed, fine-structure line features are generated using the configuration-interaction approach. These actinide opacities display similar trends to those observed for lanthanide opacities, such as the lighter actinides producing higher opacity than the heavier ones for relevant conditions in the dynamical ejecta. A line-binned treatment is employed to pre-compute opacity tables for 14 actinide elements (89 ≤ Z ≤ 102) over a grid of relevant temperatures and densities. These tabular opacities will be made publicly available for general usage in KN modelling. We demonstrate the usefulness of these opacities in KN simulations by exploring the sensitivity of light curves and spectra to different actinide abundance distributions that are predicted by different nuclear theories, as well as to different choices of ejecta mass and velocity. We find very little sensitivity to the two considered distributions, indicating that opacities for actinides with Z ≥ 99 do not contribute strongly. On the other hand, a single actinide element, protactinium, is found to produce faint spectral features in the far-infrared at late times (5–7 d post merger). More generally, we find that the choice of ejecta mass and velocity have the most significant effect on KN emission for this study.

79 ASTRONOMY AND ASTROPHYSICS↗

An ab initio molecular dynamics exploration of associates in Ba-Bi liquid with strong ordering trends

Fictive associates are widely used to describe and model liquid phases with strong ordering trends. However, little evidence is known about the assumed associates in most cases. In the present work, an ab initio molecular dynamics (AIMD) study is employed to investigate the characters of the Ba-Bi liquid, in which associates have been assumed in existing thermodynamic modeling. It is found that in the Ba rich melt, the Bi atoms are almost completely surrounded by Ba atoms. The Bi-centered coordination polyhedrons are strongly associated to crystalline structures of Ba 5 Bi 3 and Ba 4 Bi 3 with a longer lifetime than other polyhedrons during the AIMD simulations. In addition, these Bi-centered polyhedrons in Ba rich melt connect with each other through vertex, edge, face, and/or bipyramid sharing to form medium range orders (MRO). In the Bi rich melt, the Ba-centered polyhedrons also form MROs, but they are both structurally and compositionally diverse with a shorter lifetime. These findings from AIMD study provide evidences that there exist a strongly ordering Ba 4 Bi 3 associate and a weakly ordering BaBi 3 associate in the Ba-Bi liquid. Here, the predicted enthalpy of mixing in the liquid agrees well with the results by the CALPHAD modeling in the literature.

36 MATERIALS SCIENCE↗

Endo-Fullerene and Doped Diamond Nanocrystallite Based Models of Qubits for Solid-State Quantum Computers

Models of encapsulated 1/2 nuclear spin H-1 and P-31 atoms in fullerene and diamond nanocrystallite, respectively, are proposed and examined with ab-initio local density functional method for possible applications as single quantum bits (qubits) in solid-state quantum computers. A H-1 atom encapsulated in a fully deuterated fullerene, C(sub 20)D(sub 20), forms the first model system and ab-initio calculation shows that H-1 atom is stable in atomic state at the center of the fullerene with a barrier of about 1 eV to escape. A P-31 atom positioned at the center of a diamond nanocrystallite is the second model system, and 3 1P atom is found to be stable at the substitutional site relative to interstitial sites by 15 eV, Vacancy formation energy is 6 eV in diamond so that substitutional P-31 atom will be stable against diffusion during the formation mechanisms within the nanocrystallite. The coupling between the nuclear spin and weakly bound (valance) donor electron coupling in both systems is found to be suitable for single qubit applications, where as the spatial distributions of (valance) donor electron wave functions are found to be preferentially spread along certain lattice directions facilitating two or more qubit applications. The feasibility of the fabrication pathways for both model solid-state qubit systems within practical quantum computers is discussed with in the context of our proposed solid-state qubits.

Park, Seongjun↗

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗

Electronic Structure and Mechanical Properties of Solvated Montmorillonite Clay Using Large-Scale DFT Method

Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.

36 MATERIALS SCIENCE↗

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗

A model of the Schumann-Runge continuum of O2

Ab initio calculations of the B3Sigma sub u(-) potential energy curve of O2 are reported together with the transition moment connecting the B3Sigma sub u(-) and the ground electronic states. Photodissociation cross sections are presented over the wavelength region 127-152 nm. Small adjustments to the potential energy curve and the transition moment lead to a model of the absorption which agrees well with the experimental cross sections and which reproduces the structural features observed on the high-energy edge. The shoulder at 135.6 nm is attributed to the lowest 3Pi sub u state and the auxiliary maxima to transitions into the B3Sigma sub u(-) state.

Allison, A. C.↗

Ab initio multishell valence-space Hamiltonians and the island of inversion

In the shell-model framework, valence-space Hamiltonians connecting multiple major-oscillator shells are of key interest for investigating the physics of neutron-rich nuclei, which have been the subject of intense experimental activity for decades. Here we present an extension of the ab initio valence-space in-medium similarity renormalization group, which allows the derivation of such Hamiltonians nonperturbatively. Starting from initial two- and three-nucleon forces from chiral effective field theory, we then calculate properties of nuclei in the important island-of-inversion region above oxygen, so far unexplored with ab initio methods. In this work, our results in the neon and magnesium isotopes indicate the importance of neutron excitation from the $\textit{sd}$ to $\textit{pf}$ shells and ground states dominated by intruder configurations around $\textit{N}$ = 20, consistent with the conclusions from phenomenological studies. We also benchmark the excitation spectrum of 16 O with coupled-cluster theory, finding generally good agreement, and discuss implications for ground-state energies and charge radii in oxygen and calcium isotopes. Finally we outline the proper procedure for treating the longstanding issue of center-of-mass contamination, and show that with a particular choice of valence space, these spurious states can be removed successfully.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ab initio studies on structural and thermodynamic properties of magnetic Fe

The present work systematically investigates the total energy, phonon spectra, and thermodynamic properties of different polymorphs of pure Fe, i.e., FCC, BCC, and HCP, with the ab initio approach, considering various magnetic configurations. In general, the calculated energy vs. volume curves and phonon spectra agree well with previous calculations and the experimental data. In addition, their thermodynamic properties are estimated by the quasiharmonic approximation (QHA). Specifically, a superposition approach based on the latest Zentropy theory was utilized to predict magnetic transition temperatures and thermodynamic properties of pure Fe. With the ensemble of the partition function considering the multiplicity of each magnetic microstate, the current work successfully reproduced the Curie/Néel temperature and the Schottky anomaly of heat capacity in FCC, BCC, and HCP Fe purely based on the ab initio input, which exhibits good agreement with the experimental data and CALPHAD modeling.

36 MATERIALS SCIENCE↗

Uncertainty quantification of collective nuclear observables from the chiral potential parametrization

We perform an uncertainty estimate of quadrupole moments and B(E2) transition rates that inform nuclear collectivity. In particular, we study the low-lying states of 6 Li and 12 C using the ab initio symmetry-adapted no-core–shell model. For a narrow standard deviation of approximately 1% on the low-energy constants which parametrize high-precision chiral potentials, we find output standard deviations in the collective observables ranging from approximately 3%–6%. The results mark the first step towards a rigorous uncertainty quantification of collectivity in nuclei that aims to account for all sources of uncertainty in ab initio descriptions of challenging collective and clustering observables.

ab initio↗

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗

Ab-Initio Study of Stacking Fault Segregation Behavior in Ni-based Superalloys at High Temperatures

Ni-based superalloys demonstrate extraordinary creep properties compared to traditional metallic systems due to the formation of L12 (γ’) precipitates. Despite experimental evidence of alloying element segregation to the stacking faults, prior ab-initio efforts indicated repulsion of Nb from superlattice intrinsic stacking faults (SISFs) . Due to the importance of this element in key strengthening mechanisms, we developed a new technique to account for multi-component segregation behavior at temperature combining ab-initio molecular dynamics (AIMD) and Monte-Carlo (MC) simulations. Using this technique, we replicated experimentally observed segregation behavior of Nb, Co, and Ti to the SISF in L12 (Figure 1). Furthermore, we investigated the co-segregation behaviors of individual pairs of alloying elements to determine the interaction effects which lead to this segregation behavior in the full alloy system. By simulating the reduced-chemistry models using our approach, we reconciled previous ab-initio calculations demonstrating repulsion of Nb from the SISF and the experimentally observed segregation behavior.

Molecular Dynamics↗

Understanding Adsorption and Reactions at Aqueous Oxide Interfaces with Neural Network Potential Molecular Dynamics

Chemical processes at metal oxide−water interfaces are of central importance in geochemistry, biology, and energy technologies. A better understanding of these processes would allow us to make a significant step toward optimizing and controlling them, which could in turn lead to broader impacts. Computational modeling is indispensable to accomplishing this task because complexity and disorder often make it difficult to extract atomistic information from experiments. Balancing computational cost and accuracy, simulation schemes based on efficient machine learning representations of the potential energy surface (PES) predicted by ab initio calculations have become increasingly popular over the past decade. In particular, several studies have demonstrated the ability of machine learning models to accurately reproduce the complex ab initio PESs of aqueous oxide interfaces, allowing simulations of systems and processes that are not accessible with ab initio methods. In this Account, we review our recent efforts to understand adsorption processes and reactions at aqueous oxide interfaces using deep potential molecular dynamics (DPMD), a simulation scheme employing deep neural networks (DNNs), which has proven to be quite successful in accurately describing many different systems in the condensed phase. After summarizing the DPMD methodology, we first review our work on the acid−base chemistry of oxide surfaces in contact with water, a fundamental characteristic that controls proton transfer and surface charge at the interface. We focus on the aqueous interface of rutile IrO 2 , an oxide material thus far considered the best catalyst for the oxygen evolution reaction (OER). We show that this interface is characterized by a large fraction of dissociated water and a strong Brønsted acidity of the surface sites, in good agreement with the experimentally measured value of the point of zero proton charge. In our second example, we investigate how the adsorption of organic species from ambient air or water affects the structure and wettability of the aqueous interfaces of TiO 2 , a prototypical photocatalytic material. This is a question that is relevant to understanding the UV-induced hydrophilicity of TiO 2 surfaces, a property at the basis of self-cleaning windows and related applications. Specifically focusing on formic and acetic acids, the two most common atmospheric organic acids, our simulations reveal that these acids control the wettability of TiO 2 largely through acid−base chemistry at the interface rather than chemisorption on the oxide surface, a finding that could help improve the design of self-cleaning surfaces and photocatalytic devices. Finally, we review our recent study of methanol at TiO 2 −water interfaces, a system whose interest is largely motivated by the role of methanol in enhancing photocatalytic hydrogen evolution on TiO 2 . Our simulations provide mechanistic insights into the coupled roles of the organic adsorbate and water at the TiO 2 interface, with implications for how methanol enhances the activity of H 2 evolution.

adsorption↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Planar Systems for Quantum Information

This project aims to develop two‐dimensional (2D) moiré materials as a quantum simulator to implement model Hamiltonians and their phase diagrams. Progress in quantum information science (QIS) requires the development of advanced quantum materials systems. The rich family of layered van der Waals materials and their heterostructures present opportunities to create previously unrealized types of applications for QIS. Specifically, when two layers of van der Waals materials are overlaid with a small twist angle or/and lattice mismatch, a moiré superlattice with a period of about ten nanometers is formed. This provides a periodic trapping potential for electrons. Electrons can tunnel between the traps and repel each other by their mutual Coulomb interactions. The platform of 2D moiré materials provides many attractive features, including tunability of length and energy scales, charge density, and even lattice symmetry. It presents new possibilities for realizing quantum simulation of the many-body physics in a solid-state platform. This integrated team of six investigators seeks to develop relevant theoretical treatments to link ab-initio studies of 2D moiré materials to model Hamiltonians and to evaluate correlated phases predicted by these model Hamiltonians in the relevant regimes. On the experimental side, the team aims to develop methods to realize a homogeneous and highly controlled potential landscape for the electrons and to initiate, protect, and measure their quantum many-body states.

36 MATERIALS SCIENCE↗