Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗
DOE OSTI2020

Aqueous Phase Aldol Condensation of Formaldehyde and Acetone on Anatase TiO2(101) Surface: A Theoretical Investigation

A mechanistic understanding of catalytic reactions at solid-liquid interface is limited both experimentally and theoretically but attracts much interest. Using density functional theory calculations (DFT) and ab initio molecular dynamics (AIMD) simulations, we investigated the effect of liquid water on a-H abstraction, C-C coupling, and dehydration steps of aldol condensation of formaldehyde and acetone on an anatase TiO2(101) surface. The existence of the aqueous phase lowered the Gibbs energy of activation of dehydration step pronouncedly from 187 to 74 kJ/mol through proton transfer mechanism, making the hydrogenation pathway more favorable in the aqueous phase. In contrast, a mixed route prevails in the vapor phase. This work provides insights into the effect of the bulk water through a proton transfer mechanism on the dehydrogenation, C-C coupling, and dehydration

anatase TiO2, formaldehyde and acetone, aldol cond↗
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

The role of excess vacancies in stabilizing solute clusters in low-alloy steels

The formation of solute clusters in irradiated low-alloy steels, such as reactor pressure vessel steels, is a critical cause of radiation hardening and embrittlement. However, the chemical interactions among solute elements and excess vacancies remain to be fully understood. This study employs density functional theory, cluster expansion, and lattice-based Monte Carlo simulations to examine the stability and morphology of nano-size coherent solute clusters with excess vacancies. The findings reveal that excess vacancies are crucial in stabilizing and promoting the growth of Mn-Ni-Si-vacancy clusters. A minimum of seven vacancies is required for stable nucleation and growth of these clusters. These Mn-Ni-Si clusters act as defect sinks, effectively trapping and absorbing mobile vacancies generated under irradiation. Furthermore, phosphorus (P) preferentially dissolves in Mn-Ni-Si clusters due to chemical coupling with vacancies. In conclusion, this study offers new insights into solute-vacancy interactions, enhancing our understanding of solute-defect cluster formation and embrittlement in low-alloy steels.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Switching of Hybrid Improper Ferroelectricity in Oxide Double Perovskites

In ABO 3 -type perovskite oxides with Pnma symmetry, rotation (Q R+ , a 0 a 0 c + ) and tilt (Q T , a – a – c 0 ) of BO 6 octahedra are the two primary order parameters. These order parameters establish an inherent trilinear coupling with anti-ferroelectric A-site displacement (Q AFE ) to form the low-symmetry phase. The symmetry is further lowered in double perovskite oxides (DPOs) due to A/A' cation ordering. It in turn makes these systems polar via hybrid improper ferroelectric mechanism, primarily driven by Q R+ and Q T . Naturally, it has been believed that functionalities such as polarization can also be switched by tuning these primary order parameters. However, mystery around finding switching mechanism still remains. Our study based on density functional theory calculations combined with finite-temperature molecular dynamics simulations shows that the polarization switching is a two-step process, driven by out-of-phase rotation (Q R– , a 0 a 0 c – when Q T = 0 or, a – a – b – when Q T ≠ 0). A series of polar DPOs such as KLnFeOsO 6 [Ln = Sm, Gd, Dy, Tm (lanthanides) and Y (rare earth)], all belonging to P2 1 symmetry, are considered in this investigation. The polarization switching P ($\overrightarrow{P}$) occurs at a very high temperature of ~1150 K through a phase transition, from a polar (P2 1 ) phase with $\overrightarrow{P}$(+) to $\overrightarrow{P}$(-) via a non-polar P4/n phase. The switching itself is metastable in nature. The switching (both polarization and spin state) is only observed for a very short period of time (~23 ps) that poses limitation on using such a mechanism in memory device realization. We demonstrate a concurrent heating–cooling procedure to overcome such shortcoming. In conclusion, simulations conducted at 600 K further imply that long lasting switching can be achieved, at least for 1.2 ns for 600 K, and ideally for an infinite time, if the material is heated just above the T c followed by rapid cooling to a temperature below T c .

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

Electron-Beam-Induced Molecular Plasmon Excitation and Energy Transfer in Silver Molecular Nanowires

In this work, we investigate (1) electron-beam-induced plasmon absorption spectra of Ag molecular nanowire dimers and (2) electron-beam-induced energy transfer between two nanowires. We employ linear-response time-dependent density functional theory (TDDFT) and real-time TDDFT methods to simulate the electron-beam-induced plasmonic excitations, dynamics, and corresponding electron energy loss spectrum for small models of a single molecular nanowire with four Ag atoms and for two Ag nanowires. An array of different relative orientations of nanowires and of different initial excitation conditions resulting from applying an electron beam at different positions with respect to the Ag nanowires is investigated. The results demonstrate (1) an electron beam can induce plasmonic excitations from the molecular Ag nanowire ground state to the excited states that are both optically allowed and forbidden, (2) a tunability for selective excitations that can be controlled by the position of a focused electron beam, and (3) kinetic and dynamic behaviors of time-dependent electron-beam-induced energy transfer between two Ag molecular nanowires depend on the position of the beam source and nanowire separation distance, providing insights into the spatial dependences of plasmonic couplings in nanowire arrays.

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

Structure and Stability of Pristine and Carboxylate-Covered Anatase TiO 2 (001) in Aqueous Environment

The interactions of carboxylic acids and water with TiO 2 surfaces are important in applications ranging from solar cells to biomedical devices. Here we focus on the aqueous interface with the minority (001) surface of anatase TiO 2 (A-001) and the behavior of adsorbed formic and acetic acid monolayers at this interface. We investigated the structure and stability of the pristine reconstructed and formic/acetic acid covered A-001 surfaces in contact with water using density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations. The (1 × 4) reconstruction of the pristine surface is found to be stable in aqueous environment, within the time scale of our simulation. Here, carboxylic acids adsorb in deprotonated bidentate (BD) form on A-001, with the dissociated proton transferred to a surface oxygen to form a bridging hydroxyl. Of the two possible configurations, BD bridging and BD straddling, of the adsorbed species, the latter is found to rapidly transform to a monodentate structure during our simulations. Further investigation of mixed acetate–formate monolayers on A-001 in water indicates that also BD bridging species can become unstable at the boundaries between formate and acetate-covered regions, transforming to a monodentate form that does not prevent water adsorption on the TiO 2 surface.

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

Methanol Synthesis from CO 2 Hydrogenation over a Potassium-Promoted Cu x O/Cu(111) ( x ≤ 2) Model Surface: Rationalizing the Potential of Potassium in Catalysis

Alkalis have been reported as a promotor in the heterogeneous catalysis, being able to enhance the activity and selectivity of catalysts. The effective utilization of alkalis in catalyst optimization requires the fundamental understanding of the underlying mechanism. In this work, we take a potassium (K)-modified Cu x O/Cu(111) (x ≤ 2) model surface as a case study to rationalize the nature of K during the carbon dioxide hydrogenation using combined density functional theory (DFT) calculation and the kinetic Monte Carlo (KMC) simulation. Our result demonstrates the significant tuning of selectivity from carbon monoxide to methanol on going from Cu(111) to K-modified Cu x O/Cu(111). The deposited K + stabilizes the Cu x O thin film under the reducing condition of carbon dioxide hydrogenation. More importantly, our study reveals that K + acts as an active center for selective tuning in the binding, an accelerator for charge transfer, and a mediator for the electron tunneling. As a result, the K-modified Cu x O/Cu(111) opens a methanediol [H 2 C(OH) 2 ]-mediated formate pathway to facilitate the selective conversion of carbon dioxide to methanol. Our study develops the intrinsic rules of design to tune the catalytic performance using alkali metals.

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

Insights into lithium ion deposition on lithium metal surfaces

Lithium metal is among the most promising anodes for the next generation of batteries due to its high theoretical energy density and high capacity. Challenges such as extreme reactivity and lithium dendrite formation have kept lithium metal anodes away from practical applications. However, the underlying mechanisms of Li ion deposition from the electrolyte solution onto the anode surface are still poorly understood due to their inherent complexity. In this work, density functional theory calculations and thermodynamic integration via constrained molecular dynamics simulations are conducted to study the electron and ion transfer between lithium metal slab and the electrolyte in absence of an external field. Here, we explore the effect of the solvent chemistry and structure, distance of the solvated complex from the surface, anion–cation separation, and concentration of Li-salts on the deposition of lithium ions from the electrolyte phase onto the surface. Ethylene carbonate (EC), 1,2-dimethoxyethane (DME), 1,3-dioxolane (DOL), and mixtures of them are used as solvents. These species compete with the salt anion and the Li cation for electron transfer from the surface. It is found that the structure and properties of the solvation shell around the lithium cation has a great influence on the ability of the cation to diffuse as well as on its surrounding electron environment. DME molecules allow easier motion of the lithium ion compared with EC and DOL molecules. The slow growth approach allows the study of energy barriers for the ion diffusion and desolvation during the deposition pathway. This method helps elucidating the underlying mechanisms on lithium-ion deposition and provides a better understanding of the early stages of Li nucleation.

25 ENERGY STORAGE↗
DOE OSTI2020

A Bayesian analysis of nuclear deformation properties with Skyrme energy functionals

In spite of numerous scientific and practical applications, there is still no comprehensive theoretical description of the nuclear fission process based solely on protons, neutrons and their interactions. Furthermore, the most advanced simulations of fission are currently carried out within nuclear density functional theory (DFT). In spite of being fully quantum-mechanical and rooted in the theory of nuclear forces, DFT still depends on a dozen or so parameters characterizing the energy functional. Calibrating these parameters on experimental data results in uncertainties that must be quantified for applications. This task is very challenging because of the high computational cost of DFT calculations for fission. In this paper, we use Gaussian processes to build emulators of DFT models in order to quantify and propagate statistical uncertainties of theoretical predictions for a range of nuclear deformations relevant to describing the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Surface Wettability Analysis from Adsorption Energy and Surface Electrical Charge

Surface wettability is determined by the attraction of a liquid phase to a solid surface. It is typically quantified by using contact angle measurements at mineral surfaces in the case of the flotation of mineral particles. Contact angle research to describe wettability has been investigated at different scales by sessile drop measurements, molecular dynamic simulation, and atomic force microscopy. In this study, the density functional theory (DFT) was employed for predicting the surface free energy and contact angles of a well-known hydrophobic phyllosilicate mineral talc and a well-known hydrophilic phyllosilicate mineral muscovite based on the calculated interfacial energy and surface charge. The results revealed that the predicted contact angle at the atomic scale was larger than the experimental value, and identified two interactions: electrostatic interaction and hydrogen bonding, between the hydrophilic muscovite surface and the water layer, while a water-exclusion zone of 3.346 Å was found between the hydrophobic talc surface and the first water layer. This investigation gives a new perspective for wettability determination at the atomic scale.

Geochemistry & Geophysics↗
DOE OSTI2022

A combined first principles study of the structural, magnetic, and phonon properties of monolayer CrI3

The first magnetic 2D material discovered, monolayer (ML) CrI3, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI3’s structural, electronic, and magnetic properties. Here, we leverage Density Functional Theory (DFT) and high-accuracy Diffusion Monte Carlo (DMC) simulations to predict lattice parameters, magnetic moments, and spin–phonon and spin–lattice coupling of ML CrI3. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI3’s ML geometry with DMC accuracy, yielding lattice parameters in good agreement with recently published STM measurements—an accomplishment given the ~10% variability in previous DFT-derived estimates depending upon the functional. Strikingly, we find that previous DFT predictions of ML CrI3’s magnetic spin moments are correct on average across a unit cell but miss critical local spatial fluctuations in the spin density revealed by more accurate DMC. DMC predicts that magnetic moments in ML CrI3 are 3.62 μB per chromium and -0.145 μB per iodine, both larger than previous DFT predictions. The large disparate moments together with the large spin–orbit coupling of CrI3’s I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI3, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI3 exhibits a substantial spin–phonon coupling of ~3.32 cm-1. Our work, thus, establishes many of ML CrI3’s key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials.

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

The total neutron cross section of liquid and solid ammonia

Ammonia is a material of interest for future neutron moderators at high-power sources due to its high hydrogen density, low melting point, and resistance to polymerization in an intense radiation field. Its performance in such applications cannot currently be calculated due to the absence of suitable computer models for the interaction of neutrons with ammonia under relevant conditions. In an effort to develop suitable scattering kernels for computer simulations of moderator performance, we have conducted a series of Density Functional Theory and Molecular Dynamics calculations of the molecular-level thermal properties of ammonia at various temperatures within both the solid and liquid phases. In this paper, we compare computer calculations for the energy-dependent total neutron cross section of ammonia, based on these models, to experimental measurements of those cross sections at temperatures of 221 K, 180 K, and 35 K. The experimental data were collected over an energy range from 0.1 meV to 10 eV using time-of-flight techniques at the Low Energy Neutron Source (LENS) facility at Indiana University. This comparison provides a first validation in the development of thermal scattering libraries for Monte Carlo source design simulations based on liquid and solid ammonia. In conclusion, we also provide some insights into where additional development of tools for creating such models may be needed.

Ammonia↗
DOE OSTI2021

Effect of Fluoroethylene Carbonate Additives on the Initial Formation of the Solid Electrolyte Interphase on an Oxygen-Functionalized Graphitic Anode in Lithium-Ion Batteries

The formation of a solid electrolyte interphase (SEI) at the electrode/electrolyte interface substantially affects the stability and lifetime of lithium-ion batteries (LIBs). One of the methods to improve the lifetime of LIBs is by the inclusion of additive molecules to stabilize the SEI. To understand the effect of additive molecules on the initial stage of SEI formation, in this work we compare the decomposition and oligomerization reactions of a fluoroethylene carbonate (FEC) additive on a range of oxygen-functionalized graphitic anodes to those of an ethylene carbonate (EC) organic electrolyte. A series of density functional theory (DFT) calculations augmented by ab initio molecular dynamics (AIMD) simulations reveal that EC decomposition on an oxygen-functionalized graphitic ($11\bar{2}0$) edge facet through a nucleophilic attack on an ethylene carbon site (C E ) of an EC molecule (S2 mechanism) is spontaneous during the initial charging process of LIBs. However, decomposition of EC through a nucleophilic attack on a carbonyl carbon (C C ) site (S1 mechanism) results in alkoxide species regeneration that is responsible for continual oligomerization along the graphitic surface. In contrast, FEC prefers to decompose through an S1 pathway, which does not promote alkoxide regeneration. Including FEC as an additive is thus able to suppress alkoxide regeneration and results in a smaller and thinner SEI layer that is more flexible toward lithium intercalation during the charging/discharging process. In addition, we find that the presence of different oxygen functional groups at the surface of graphite dictates the oligomerization products and the LiF formation mechanism in the SEI.

25 ENERGY STORAGE↗
DOE OSTI2023

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

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

Coverage, repulsion, and reactivity of hydrogen on High-Entropy alloys

Modeling hydrogen evolution reaction (HER) activity probability on IrPdPtRhRu(1 1 1) high-entropy alloys. Determining hydrogen coverages based on ligand effects and generalized hydrogen–hydrogen repulsion. The rate of H 2 formation is highly impacted by the level of hydrogen coverage on the catalyst surface. In search of optimal catalytic properties high-entropy alloys (HEA) are promising candidates that utilize the compositional space of multiple elements. Based on simulations of HEA model (1 1 1) surfaces with a range of hydrogen coverages, distributions of binding energies are used to construct a framework that approximates the probability that adsorbed hydrogen may lead to the formation of H 2 as a function of applied potential. By optimizing the alloy compositions for the highest activity probability at given potentials the best and most efficient catalyst candidates for HER can be identified. Treating hydrogen–hydrogen repulsion effects and binding energy separately, we find that the repulsion is larger for HEAs than for pure metals. Differing isotherm slopes in the mean adsorption and desorption energies demonstrate a possible hysteresis for hydrogen adsorption on HEAs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗