Wide-ranged multiphase equation of state for iron and model variations addressing uncertainties in high-pressure melting
Not Available
36 MATERIALS SCIENCE↗
SEARCH · Engineering Papers
Results for “initio molecular dynamics”
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.
Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications
Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Spectral-partitioned Kohn-Sham density functional theory
Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Transferable interatomic potential for aluminum from ambient conditions to warm dense matter
Not Available
36 MATERIALS SCIENCE↗
Low-frequency lattice phonons in halide perovskites explain high defect tolerance toward electron-hole recombination
Low-cost solution-based synthesis of metal halide perovskites (MHPs) invariably introduces defects in the system, which could form Shockley-Read-Hall (SRH) electron-hole recombination centers detrimental to solar conversion efficiency. Here, we investigate the nonradiative recombination processes due to native point defects in methylammonium lead halide (MAPbI 3 ) perovskites using ab initio nonadiabatic molecular dynamics within surface-hopping framework. Regardless of whether the defects introduce a shallow or deep band state, we find that charge recombination in MAPbI 3 is not enhanced, contrary to predictions from SRH theory. We demonstrate that this strong tolerance against defects, and hence the breakdown of SRH, arises because the photogenerated carriers are only coupled with low-frequency phonons and electron and hole states overlap weakly. Both factors appreciably decrease the nonadiabatic coupling. We argue that the soft nature of the inorganic lattice with small bulk modulus is key for defect tolerance, and hence, the findings are general to other MHPs.
Science & Technology - Other Topics↗
Reactivity of the Hydrated Electron
Hydrated electrons are created in virtually every radiation environment and in many photochemical or electrochemical environments where liquid water is present, so their reaction products and reaction rate constants are naturally important in applications. Thanks to the strong optical absorbance of (e - ) aq , these rate constants are easy to measure, and a large database has been accumulated. It is not generally appreciated that no working theory of hydrated electron reaction rates presently exists. We discuss key experimental observations of hydrated electron reactions in the context of recent progress in theoretical and simulation developments toward understanding them, made possible by ever increasing computational power.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation
An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.
36 MATERIALS SCIENCE↗
ORNL_AISD_NiNb
This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.
36 MATERIALS SCIENCE↗
Roles of Carburization and Oxidation on Transition-Metal Catalyzed Surface Reactions
Many industrial chemical transformation processes are based on reactions catalyzed by transition metals. In some cases the nano-size particles used as catalysts are synthesized simultaneously or right before the initiation of the main catalytic reaction. Synthesis treatments involve complex processes such as reductions followed by calcinations in oxygenated environments coupled to the presence of gases used for the main reaction. Synthesis details are crucial because they have important consequences on the catalyst shape, size, and chemical composition, and may also affect the structure and composition of the catalyst support. All these effects may alter the main catalyzed reaction in unexpected ways. This work focuses on the synthesis of the overall catalytic system that includes catalyst and support taking place in parallel with the initial stages of the main catalyzed reaction. Two important processes are characterized: carburization and oxidation of the transition metal catalyst during synthesis or initial stages of the main reaction. The objective is to understand how the coking mechanism coupled to surface oxidation and changes in the substrate influence the dry reforming of methane taken as a model reaction. The approach uses first principles computational methods to investigate the catalyst synthesis and characterization and catalytic reaction mechanisms. These studies are complemented with experimental characterization done by the group of Dr. Renu Sharma at the National Institute of Standards and Technology (unfunded collaborator) using high resolution transmission electron microscopy. Specific goals are to identify the role of the synthesis process on the initial nanocatalyst structure and composition, and to determine how that structure and composition may evolve after being exposed to the reactants atmosphere. These studies facilitate a systematic analysis of conditions where the nanoparticle composition can be optimized regarding activity and stability. Moreover, the fundamental study of catalyst/support interactions yields valuable insights for other catalytic processes and is a first step toward controlled nanocatalyst synthesis.
36 MATERIALS SCIENCE↗
Palladium at high pressure and high temperature: A combined experimental and theoretical study
Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid– solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.
36 MATERIALS SCIENCE↗
Final Scientific Report for SPIRAL (9/1/21-8/31/25)
The final scientific report for our DOE-funded project covering the period 9/1/21-8/31/25
36 MATERIALS SCIENCE↗
Ab Initio Phase Diagram of Chromium to 2.5 TPa
Chromium possesses remarkable physical properties such as hardness and corrosion resistance. Chromium is also a very important geophysical material as it is assumed that lighter Cr isotopes were dissolved in the Earth’s molten core during the planet’s formation, which makes Cr one of the main constituents of the Earth’s core. Unfortunately, Cr has remained one of the least studied 3d transition metals. In a very recent combined experimental and theoretical study (Anzellini et al., Scientific Reports, 2022), the equation of state and melting curve of chromium were studied to 150 GPa, and it was determined that the ambient body-centered cubic (bcc) phase of crystalline Cr remains stable in the whole pressure range considered. However, the importance of the knowledge of the physical properties of Cr, specifically its phase diagram, necessitates further study of Cr to higher pressure. In this work, using a suite of ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries, we obtain the theoretical phase diagram of Cr to 2.5 TPa. We calculate the melting curves of the two solid phases that are present on its phase diagram, namely, the lower-pressure bcc and the higher-pressure hexagonal close-packed (hcp) ones, and obtain the equation for the bcc-hcp solid–solid phase transition boundary. We also obtain the thermal equations of state of both bcc-Cr and hcp-Cr, which are in excellent agreement with both experimental data and QMD simulations. We argue that 2180 K as the value of the ambient melting point of Cr which is offered by several public web resources (“Wikipedia,” “WebElements,” “It’s Elemental,” etc.) is most likely incorrect and should be replaced with 2135 K, found in most experimental studies as well as in the present theoretical work.
equation of state↗
Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds
The organic compounds found in carbonaceous chondrite meteorites provide insight into primordial solar system chemistry. Evaluating the formation and decomposition mechanisms of meteoritic amino acids may aid our understanding of the origins of life and homochirality on Earth. The amino acid glycine is widespread in meteorites and other extraterrestrial environments; other amino acids, such as isovaline, are found with enantiomeric excesses in some meteorites. The relationship between meteoritic amino acids and other compounds with similar molecular structures, such as aliphatic monoamines and monocarboxylic acids is unclear; experimental results evaluating the decomposition of amino acids have produced inconclusive results about the preferred pathways, reaction intermediates, and if the conditions applied may be compatible with those occurring inside meteoritic parent bodies. In this work, we performed extensive tandem metadynamics, umbrella sampling, and committor analysis to simulate the neutral mild hydrothermal decomposition mechanisms of glycine and isovaline and put them into context for the origins of meteoritic organic compounds. Our ab initio simulations aimed to determine free energy profiles and decomposition pathways for glycine and isovaline. We found that under our modeled conditions, methylammonium, glycolic acid, and sec-butylamine are the most likely decomposition products. These results suggest that meteoritic aliphatic monocarboxylic acids are not produced from decomposition of meteoritic amino acids. Our results also indicate that the decomposition of L-isovaline prefers an enantioselective pathway resulting in the production of (S)-sec-butylamine.
Enhanced sampling↗
An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem
The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Ultrafast charge transfer coupled to quantum proton motion at molecule/metal oxide interface
Understanding how the nuclear quantum effects (NQEs) in the hydrogen bond (H-bond) network influence the photoexcited charge transfer at semiconductor/molecule interface is a challenging problem. By combining two kinds of emerging molecular dynamics methods at the ab initio level, the path integral–based molecular dynamics and time-dependent nonadiabatic molecular dynamics, and choosing CH 3 OH/TiO 2 as a prototypical system to study, we find that the quantum proton motion in the H-bond network is strongly coupled with the ultrafast photoexcited charge dynamics at the interface. The hole trapping ability of the adsorbed methanol molecule is notably enhanced by the NQEs, and thus, it behaves as a hole scavenger on titanium dioxide. The critical role of the H-bond network is confirmed by in situ scanning tunneling microscope measurements with ultraviolet light illumination. It is concluded the quantum proton motion in the H-bond network plays a critical role in influencing the energy conversion efficiency based on photoexcitation.
74 ATOMIC AND MOLECULAR PHYSICS↗
SSAIMS—Stochastic-Selection Ab Initio Multiple Spawning for Efficient Nonadiabatic Molecular Dynamics
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
CTRAMER: An open-source software package for correlating interfacial charge transfer rate constants with donor/acceptor geometries in organic photovoltaic materials
In this paper, we present CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory)—an opensource software package for calculating interfacial charge-transfer (CT) rate constants in organic photovoltaic (OPV) materials based on ab initio calculations and molecular dynamics simulations. The software is based on identifying representative donor/acceptor geometries within interfacial structures obtained from molecular dynamics simulation of donor/acceptor blends and calculating the corresponding Fermi’s golden rule CT rate constants within the framework of the linearized-semiclassical approximation. While the methods used are well established, the integration of these state-of-the-art tools originating from different disciplines to study photoinduced CT processes with explicit treatment of the environment, in our opinion, makes this package unique and innovative. The software also provides tools for investigating other observables of interest. After outlining the features and implementation details, the usage and performance of the software are demonstrated with results from an example OPV system.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials
Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.
Park, Sanghyun J. [Princeton University , , , ,]↗