Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal–Organic Frameworks

Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO 2 . Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal–Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. In conclusion, a strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.

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

Building Toward the Future in Chemical and Materials Simulation with Accessible and Intelligently Designed Web Applications

Over the last few decades, significant progress has been made in the development and use of electronic structure and other molecular simulation methods. As these methods become more mature and are able to simulate larger and more complex chemical simulations, the need for improvement in scientific visualization, molecular builders, simplified input to simulation methods, and the development of new approaches and languages to describe simulations, along with workflows to carry them out, becomes more apparent. In this chapter, we describe our recent efforts in developing a prototype open-source computational tool called Arrows that combines NWChem, SQL and NoSQL databases, email, web APIs, and web applications in a way that make molecular and materials modeling accessible to all scientists and engineers. At the same time, because of its simplified input, it provides a framework for expert users to carry out large numbers of calculations and run complex workflows.

2↗
DOE OSTI2022

Characterizing particle-like charge-migration dynamics with high-order harmonic sideband spectroscopy

Here, we introduce high-order harmonic sideband spectroscopy and show that it can be a robust probe of charge migration (CM) in a halogenated carbon-chain molecule. We simulate both the CM and high-order harmonic generation (HHG) dynamics using ab initio time-dependent density-functional theory. We find that CM dynamics initiated along the molecular backbone induces sidebands in the HHG spectrum driven by a delayed laser pulse that is polarized perpendicular to the molecular axis. Monitoring the spectrum as either the HHG laser frequency or the relative delay is scanned allows for the extraction of detailed information about the time-domain characteristics of the CM process.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Influence of Ag metal dispersion on the thermal conversion of ethanol to butadiene over Ag-ZrO 2 /SiO 2 catalysts

Atomistic scale models were developed and coupled with experimental investigation to deliver a functional understanding of catalytic activity and selectivity in the conversion of ethanol to 1,3-butadiene over Ag/ZrO 2 /SiO 2 . A detailed evaluation of the structural and electronic properties of the resultant catalyst models led to the identification of critical active sites of the catalyst. More importantly, the extent of Ag dispersion on the SiO 2 support and relative proximity to ZrO 2 were found to vary with the oxidation state of Ag and local coordination environment (Ag-O SiO2 ) allowing for critical control of ethanol conversion towards butadiene or ethylene. Simulations revealed that less dispersed or clustered Ag contain predominantly Ag 0 charge state and promote conversion of ethanol to ethylene. The well-dispersed Ag/ZrO 2 /SiO 2 catalyst instead contain a larger fraction of cationic Ag δ+ , and predominantly promote ethanol dehydrogenation and subsequent production of butadiene. The theoretical insights drawn were validated and confirmed experimentally using TEM, XRD and reactivity measurements demonstrating the effect of Ag dispersion on the selectivity of ethanol conversion.

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

Three Jahn-Teller States of Matter in Spin-Crossover System Mn(taa)

Three high-spin phases recently discovered in the spin-crossover system Mn(taa) are identified through analysis by a combination of first-principles calculations and Monte Carlo simulation as a low-temperature Jahn-Teller ordered (solid) phase, an intermediate-temperature dynamically correlated (liquid) phase, and an uncorrelated (gas) phase. In particular, the Jahn-Teller liquid phase arises from competition between mixing with low-spin impurities, which drive the disorder, and intermolecular strain interactions. Additionally, the latter are a key factor in both the spin-crossover phase transition and the magnetoelectric coupling. Jahn-Teller liquids may exist in other spin-crossover materials and materials that have multiple equivalent Jahn-Teller axes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Selective Hydration of Rutile TiO 2 as a Strategy for Site-Selective Atomic Layer Deposition

In this report the feasibility of a site-selective hydration strategy that enables site-selective atomic layer deposition (ALD) is investigated among four rutile TiO 2 facets [(110), (100), (101) and (001)] and their most prevalent step edges. First-principles simulations of asymmetric slab models were utilized to create accurate representations of pristine terrace and step edge sites. The adsorption free energies for molecular and dissociative adsorption of H 2 O were calculated to evaluate this strategy as a viable route to step edge selectivity. We predict that selective hydroxylation is possible on the 110 and 001 step edges and further computationally evaluate three metalorganic ALD precursors for their compatibility with the selective hydration strategy. Experimental evidence for delayed nucleation of ALD on rutile (001), (110), and (100) TiO2 single crystals corroborates predictions of the dehydration of the surface and suggests the possibility of site-selective ALD.

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

CH 4 steam reforming on Pt + Pd/Al 2 O 3 monolith: impact of Mn 0.5 Fe 2.5 O 4 spinel addition

The performance of Pt + Pd/Al 2 O 3 monolith catalysts was analyzed under a feed (H 2 + CO + CH 4 + NO + O 2 + H 2 O + CO 2 ) simulating the exhaust gas of stoichiometric natural gas vehicles (NGVs). The flow reactor results show enhanced low temperature (<400 °C) CH 4 oxidation activity under cyclic lean/rich feed modulation and with Mn 0.5 Fe 2.5 O 4 (MFO) spinel addition to the platinum group metal (PGM) catalyst. In the absence of the MFO spinel, the PGM catalyst can achieve full CH 4 conversion, but the presence of spinel in the catalyst limits the high temperature (>500 °C) CH 4 conversion. The depletion of gas phase O 2 , as indicated by the detection of products CO and H 2 in the reactor effluent at higher temperature and CH 4 conversion, portends the emergence of steam reforming of CH 4 (SRM). To investigate the cause of the detrimental impact of spinel under the high temperature, net rich regime, we show through a combination of SRM experiments and post-reaction catalyst characterization that SRM is negatively impacted by the accumulation of Mn and Fe spinel species on the PGM catalyst. This is due to the migration of the base metal species and potential encapsulation and blockage of PGM active sites. The latter is supported by density functional theory (DFT) calculations that indicate favorable metal oxide decoration of active Pt step sites. Lastly, a zoned catalyst design is developed that mitigates the detrimental metals migration while exploiting the methane conversion enhancement afforded by spinel oxides under oxygen rich conditions.

03 NATURAL GAS↗
DOE OSTI2025

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

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

Evaluation of Sb-Nd and Te-Nd phases within the U-Zr fuel matrix and their interactions with HT9 alloy

Antimony (Sb) and tellurium (Te) were investigated as potential additives for U-10Zr (wt.%) metallic fuel to limit the fuel-cladding chemical interaction (FCCI) with HT-9 alloy. Neodymium (Nd) was utilized to simulate the formation of lanthanide-based solid fission products which are known to play a detrimental role in FCCI. Fuel alloys of U-Zr-Sb-Nd and U-Zr-Te-Nd were evaluated in their annealed condition and compared against their as-cast conditions. Isothermal diffusion couple experiments were performed between U-Zr-Nd, U-Zr-Sb-Nd, and U-Zr-Te-Nd against the cladding alloy HT9 to evaluate the effectiveness of the additives to stabilize Nd within the fuel alloys, as well as investigate the interaction regions that form between the different fuel alloys and HT9. Further, SbNd and Sb 3 Nd 4 , and TeNd are found to be the primary neodymium-based phases formed in the U-Zr-Sb-Nd and U-Zr-Te-Nd alloys, respectively. The zirconium-based phase, Zr 2 Sb, is also found to form within the former alloy. All phases were found to remain stable through the diffusion experiments and exhibited no interaction with HT9 constituent elements. Preferential interaction between Nd with additivities Te and Sb compared to constituting elements in HT9 was further verified based on density functional theory (DFT) calculated enthalpy of mixing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Efficient 3D Molecular Design with an E(3) Invariant Transformer VAE

This work introduces a three-dimensional (3D) invariant graph-to-string transformer variational autoencoders (VAE) (Vagrant) for generating molecules with accurate density functional theory (DFT)-level properties. Vagrant learns to model the joint probability distribution of a 3D molecular structure and its properties by encoding molecular structures into a 3D-aware latent space. Directed navigation through this latent space implicitly optimizes the 3D structure of a molecule, and the latent embedding can be used to condition a generative transformer to predict the candidate structure as a one-dimensional (1D) sequence. Additionally, we introduce two novel sampling methods that exploit the latent characteristics of a VAE to improve performance. Here, we show that our method outperforms comparable 3D autoregressive and diffusion methods for predicting quantum chemical property values of novel molecules in terms of both sample quality and computational efficiency.

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

Effects of charged interfaces on electrolyte decomposition at the lithium metal anode

Lithium–Sulfur batteries are promising candidates to substitute conventional Li-ion batteries due to their higher energy density and reduced cost. However, several challenges related to the reactivity of the lithium metal anode have prevented this technology from becoming broadly commercialized. Lithium's high reactive nature leads to the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. Thus, a comprehensive understanding of how the SEI film is formed is crucial to help elucidate improvements for this battery technology. In this work, we use density functional theory (DFT) based computational methods to investigate the effect of charged interfaces on the electrolyte degradation at the lithium anode. Here, several electrolyte mixtures were considered including 1,2-dimethoxyethane (DME) and 1,3-dioxolane (DOL) as solvents, and lithium bis(fluorosulfonyl)imide (LiFSI) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as salts. It is found that the extent of salt decomposition is higher when the interface is charged. In addition, two types of solvent reduction mechanisms are identified: C–O bond cleavage and radical attacks. Charge transfer and evolution analysis are characterized in detail. Finally, solvent decomposition reactions are found to become more thermodynamically favorable and their activation barriers to be diminished under the effect of constant potentials and electric fields.

25 ENERGY STORAGE↗
DOE OSTI2020

(Photo)Electrocatalytic CO 2 Reduction at the Defective Anatase TiO 2 (101) Surface

Excessive carbon dioxide (CO 2 ) emissions by combustion of fossil fuels are linked to global warming and rapid climate change. One promising route to lowering the concentration of CO 2 in the atmosphere is to reduce it to useful small molecules via photoelectrocatalytic hydrogenation, which would enable solar energy storage with a zero-carbon emission cycle and perform a more efficient separation of the photogenerated electron and hole pair than pure photocatalysis. Indeed, photoelectrocatalytic CO 2 reduction has been an intense focus of research. Using the density functional theory (DFT), we studied the CO 2 reduction reaction on the defective anatase TiO 2 (101) surface, at both the solvent/catalyst and the electrolyte/catalyst interfaces. The analysis of the electronic structure of the surface shows a contrast between the solvent/catalyst and the electrolyte/catalyst interfaces, which results in the two corresponding catalytic cycles being distinct. Furthermore, our study explains at the electronic and mechanistic levels why methanol is the main product in the presence of the electrolyte and why the overpotential is not only controlled by the reaction process but also by the diffusion process.

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

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2024

VOTCA: multiscale frameworks for quantum and classical simulations in soft matter

Many physical phenomena in liquids and soft matter are multiscale by nature and can involve processes with quantum and classical degrees of freedom occurring over a vast range of length- and timescales. Examples range from structure formation processes of complex polymers or even polymer blends (Svaneborg & Everaers, 2023) on the classical side to charge and energy transport and conversion processes (Lee et al., 2019) involving explicit electronic and, therefore, quantum information. The Versatile Object-oriented Toolkit for Coarse-graining Applications (VOTCA) provides multiscale frameworks built on a comprehensive set of methods for the development of classical coarse-grained potentials (VOTCA-CSG) as well as state-of-the art excited state electronic structure methods based on density-functional and many-body Green’s function theories, coupled in mixed quantum-classical models and used in kinetic network models (VOTCA-XTP).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Insights from Computational Studies on the Anisotropic Volume Change of Li x NiO 2 at High States of Charge (x < 0.25)

The need for high-capacity Li-ion battery cathodes has favored the increase of Ni content in commercial battery cells. However, at high states of charge (SOCs), Ni-rich materials undergo a phase transition and volume collapse with deleterious effects on battery performance. It is uncertain whether this drastic volume change is caused by the phase transition or not. To provide more insight into the volume-phase transition relationship in the high Ni cathode LixNiO 2 , we performed density functional theory calculations, along with molecular dynamics simulations using machine learning potentials to calculate the temperature- and composition-dependent free energy differences between the suspected phases at high SOCs (x < 0.25). Herein, we find that the calculated free energy difference between the suspected phases containing different oxygen stacking sequences is small at room temperature. Furthermore, we find that the collapse of the layered LiNiO 2 c-lattice parameter at high SOCs is mainly due to the electronic depletion of the oxygen sublattice and the lack of screening from positive Li ions. The interactions between adjacent oxygen ions across an empty Li layer (NiO 2 ) are largely controlled by van der Waals interactions and are in fact similar regardless of the oxygen stacking, which explains the negligible free energy differences between O1 and O3 stacking in NiO 2 .

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

Design and Synthesis of PtPdNiCoMn High‐Entropy Alloy Electrocatalyst for Enhanced Alkaline Hydrogen Evolution Reaction: A Theoretically Supported Predictive Design Approach

Electrocatalytic hydrogen generation requires a multifunctional electrocatalyst with abundant active sites to drive multielectron transfer reactions. High entropy alloys (HEA) are five or more-elements with high configurational entropy are considered unique materials for next-generation electrocatalysts. Here, in this work, based on new screening guidelines for catalyst selections that combine density-functional theory calculated Gibbs formation-enthalpy with bond length and electronegativity variance, a novel HEA electrocatalyst consisting of five elements, namely, Pt, Pd, Ni, Co, and Mn has been designed. By simple room temperature electrodeposition, the designed catalyst is prepared and its hydrogen evolution reaction (HER) is explored and validated through experimental and theoretical approaches. The HEA demonstrated a superior HER activity with an overpotential of 22.6 mV at -10 mA cm -2 which outperforms Pt/C commercial catalyst. No evident degradation of the material is detected even after 100 hours of continuous operation under high current density. Moreover, the HEA has shown exceptional performance in harsh electrolyte conditions such as in simulated seawater and actual seawater. Remarkably, the density-functional theory calculated Gibbs formation-enthalpy is small (≈0 eV) compared to Pt/C placing the new HEA near the apex of Trasatti's model of Volcano plot, which is also suggestive of superior HER activity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

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

Hydride Accessibility and Reactivity in the Configurational and Stoichiometric Space of β-Ga 2 O 3 for CO 2 Hydrogenation

Understanding how surface species evolve under reaction conditions is essential for improving catalyst design for efficient CO 2 hydrogenation. This work combines systematic DFT calculations with grand canonical sampling to investigate the stability and reactivity of Ga–H species on β-Ga 2 O 3 across a range of reaction conditions. Initial DFT studies reveal that when Ga–H species are present, they facilitate formate formation via a low-barrier pathway, largely independent of the surface termination or hydrogen site. However, grand canonical sampling shows that under a broad range of reaction conditions─especially at high oxygen chemical potentials associated with high water content─Ga–H species are thermodynamically inaccessible. Furthermore, adsorbed water molecules can block reactive sites, inhibiting CO 2 activation even when hydrides are present. These findings suggest that the lack of accessible hydride species, rather than their intrinsic reactivity, could contribute to reduced catalytic performance of β-Ga 2 O 3 under more oxidizing, high-conversion conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗