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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2022

Elucidating CO Oxidation Pathways on Rh Atoms and Clusters on the “29” Cu 2 O/Cu(111) Surface

We report single-atom catalysts have attracted a great deal of attention due to their distinct reactivity and potential for cost savings. However, despite the wealth of literature in recent years, identifying the exact nature of the active sites and associated reaction mechanisms remains challenging in many cases. Herein, we take a surface science approach to understand how Rh single atoms and small clusters behave on the thin film “29” Cu 2 O grown on Cu(111). We find that in contrast to Pt, which is present solely as single atoms on the “29” Cu 2 O surface, Rh atoms and clusters coexist and each enable low-temperature CO oxidation, but via different pathways. Specifically, the single Rh atoms produce CO 2 at 444 K via a Mars van Krevelen mechanism whereas the Rh clusters can also dissociate CO, as demonstrated via isotope labeling, and liberate CO 2 at 313 K. Density functional theory (DFT) calculations quantify the energetics of these different pathways and demonstrate that only extended Rh is capable of CO dissociation. Low-temperature scanning tunneling microscopy (STM) reveals that unlike Pt atoms on the same surface, which stay atomically dispersed, the distribution of Rh structures is dependent on pretreatment conditions. DFT calculations reveal the greater tendency of Rh atoms to cluster than Pt, and STM image simulations confirm the active sites. Ambient pressure X-ray photoelectron spectroscopy studies on the same single crystal model systems demonstrate that 1% of a monolayer of Rh on the “29” Cu 2 O thin film significantly accelerates its reduction by CO at 400 K, thus confirming the ultrahigh vacuum surface science findings. Together, these results illustrate how well-defined single crystal experiments are useful in building structure–function relationships that elucidate the reactivity of different ensemble sizes with a level of detail beyond what is possible with high surface area catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemically Mediated Au–C(sp 2 ) Anchors for Molecular Electronics

Terminal anchor groups play a key role in the stability and electronic properties of molecular junctions. Single molecule junctions typically consist of two preinstalled terminal anchors linking organic molecules to metal electrodes. Here, in this work, we show that p -terphenyl derivatives containing only a single terminal anchor show conductance features similar to junctions with two preinstalled terminal anchors. A set of p -terphenyl derivatives with one terminal anchor was prepared using automated chemical synthesis and characterized using single molecule electronics experiments, molecular dynamics (MD) simulations, bulk electrochemistry and spectroscopy, and nonequilibrium Green’s function-density functional theory (NEGF-DFT) calculations. Our results show that 4-amino- p -terphenyl (PPP) and related analogs exhibit a well-defined high conductance state that is diminished or absent in other p-terphenyl derivatives lacking a preinstalled amine terminal anchor or fluorine or methyl substitutions at the terminal para position. However, a low conductance state is observed in all amino- p -terphenyl derivatives with one preinstalled anchor due to molecular junctions formed by noncovalent dimeric π–π stacking interactions. The observed high conductance state diminishes upon the addition of reducing agents and is restored upon the addition of an oxidizing agent. Our results suggest that the high conductance state arises due to Au–C(sp 2 ) bond formation facilitated by a single electron oxidation event at the electrode surface. A series of control experiments with different anchor groups shows that primary amines play a key role in forming Au–C bonds for molecular junctions. Overall, these results suggest that Au–C bond formation gives rise to high conductance pathways in organic molecules containing only one preinstalled terminal anchor. Insights from this work can be leveraged in the design of molecular electronic devices, particularly in understanding the mechanisms of molecular binding and junction formation.

charge transport↗
DOE OSTI2021

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

Modeling of Branched (L, T and Y) Carbon Nanotubes

Models for connecting two or three carbon nanotubes (CNT) using topological defects (i.e., pentagons and heptagons) are presented for the characterization of experimentally observed L, T and Y CNT junctions. The effects of the separation and orientation of the topological defects on the structures and energetics of these junctions are investigated using the nonlocal density function theory (DFT) and semi-empirical molecular orbital (AM1) calculations, and the Brenner empirical potential molecular mechanics simulations. The potential applications of L, Y and T CNT junctions in nanoelectronic devices are also discussed.

Han, Jie↗
DOE OSTI2026

Ab Initio Bulk Free Energy Surface of Proper Ferroelectrics

We report a systematic and accurate approach for deriving the bulk free energy surface (FES), a function of temperature, polarization, and strain, from the first-principles density functional theory (DFT) of proper ferroelectrics. The core of our approach is the metadynamics algorithm that extracts the polarization dependence of the FES from all-atom molecular dynamics simulations without an a priori ansatz. The rest of the FES is derived from the metadynamics trajectories that span the relevant phase space. We demonstrate our approach in the case of lead titanate. The errors across the phase transition, due to DFT numerics, all-atom molecular dynamics, and free energy evaluation by enhanced sampling, can be systematically controlled and are of the order of 1 meV/atom. The accuracy of the resulting ab initio FES is only limited by the adopted functional approximation of DFT.

Xie, Pinchen [Lawrence Berkeley National Laborator↗
DOE OSTI2023

Active Learning Framework

Machine learning (ML) of interatomic potentials show great promise to accelerate scientific simulation, e.g., by emulating expensive computations at a high accuracy but much reduced computational cost. Training datasets are calculated from computationally expensive ab initio quantum mechanics methods, density functional theory (DFT). Trained on this data, an ML model can be very successful in predicting energy and forces for new atomic configurations. A critical factor is the quality and diversity of the training dataset. Thus, a highly automated approach to dataset construction based on active learning framework is designed suitable for material physics. The active learning scheme begins with fully randomized atomic configurations. Then, many Molecular Dynamics (MD) trajectories are simulated using current ML potentials, where each MD trajectory is initialized to a random disordered configuration. The temperature is varied in order to diversify the sampled configuration during these simulations. The variance of predictions for eight neural networks within an ensemble is analyzed to determine whether the model is operating as expected. This helps in determining whether collecting more data would be helpful to the model by checking the ensemble variance is greater than the threshold. In this case, the MD trajectory is terminated and the final atomic configuration is placed on a queue (SQL database) for DFT calculations and added to training dataset. Periodically, ML model is retrained to the updated training model. This Active Learning loop is iterated until the cost of MD simulations becomes prohibitively expensive. The MD simulations will hopefully be sufficiently robust to support nucleation after many active learning iterations. In this sense, active learning scheme must automatically discover the important low energy and nonequilibrium physics.

Nebgen, Benjamin↗
DOE OSTI2025

Accuracy of Kohn–Sham density functional theory for warm- and hot-dense matter equation of state

We study the accuracy of Kohn–Sham density functional theory (DFT) for warm- and hot-dense matter (WDM and HDM). Specifically, considering a wide range of systems, we perform accurate ab initio molecular dynamics simulations with temperature-independent local/semilocal density functionals to determine the equations of state at compression ratios of 3x–7x and temperatures near 1 MK. We find very good agreement with path integral Monte Carlo benchmarks, while having significantly smaller error bars and smoother data, demonstrating the accuracy of DFT for the study of WDM and HDM at such conditions. In addition, using a Δ-machine learned force field scheme, we confirm that the DFT results are insensitive to the choice of exchange-correlation functional, whether local, semilocal, or nonlocal.

Suryanarayana, Phanish (ORCID:0000000151720049)↗
DOE OSTI2021

Modeling dark- and light-induced crystal structures and single-crystal optical absorption spectra of ruthenium-based complexes that undergo SO 2 -linkage photoisomerization

A family of coordination complexes of the type [Ru(SO 2 )(NH 3 ) 4 X] m+ Y n – (m, n = 1 or 2) exhibit optical switching capabilities in their single-crystal states. This striking effect is caused by the light-induced formation of SO 2 -linkage photoisomers, which are metastable if kept at suitably cool temperatures. We modeled the dark- and light-induced states of these large crystalline complexes via plane-wave (PW)- and molecular-orbital (MO)-based density functional theory (DFT) and time-dependent DFT in order to calculate their structural and optical properties; the calculated results are compared with experimental data. We show that the PW-DFT-based periodic models replicate the structural properties of these complexes more effectively than the MO-DFT-based molecular-fragment models, observing only small deviations in key bond lengths relative to the experimentally derived crystal structures. The periodic models were also found to more effectively simulate trends seen in experimental optical absorption spectra, with optical absorbance and coverage of the visible region increasing with the formation of the photoinduced geometries. The contribution of the metastable photoisomeric species more heavily focuses on the lower-energy end of the spectra. Spectra generated from the molecular-fragment models are limited by the geometry of the fragment used and the number of excited-state roots considered in those calculations. In general, periodic models outperform the molecular-fragment models owing to their ability to better appreciate the periodic phenomena that are present in these crystalline materials as opposed to MO approaches, which are finite methods. We thus demonstrate that PW-DFT-based periodic models should be considered as a more than viable method for simulating the optical and electronic properties of these single-crystal optical switches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗
DOE OSTI2023

Development of a Reactive Force Field for Simulating Photoinitiated Acrylate Polymerization

Light-driven and photo-curable polymer based additive manufacturing (AM) has enormous potential due to its excellent resolution and precision. Acrylated radical chain-growth polymerized resins are widely used in photopolymer AM due to their fast kinetics, and often serve as a departure point for developing other resin materials for photopolymer-based AM technologies. For successful control of the photopolymer resins, the molecular basis of the acrylate free-radical polymerization has to be understood in detail. We present an optimized reactive force field (ReaxFF) for molecular dynamics (MD) simulations of acrylate polymer resins that captures radical polymerization thermodynamics and kinetics. The force field is trained against an extensive training set including density functional theory (DFT) calculations of reaction pathways along the radical polymerization from methyl acrylate to methyl butyrate, bond dissociation energies, and structures and partial charges of several molecules and radicals. We also found that it was critical to train the force field against an incorrect, nonphysical reaction pathway observed in simulations that used parameters not optimized for acrylate polymerization. As a result, the parameterization process utilizes a parallelized search algorithm, and the resulting model can describe polymer resin formation, crosslinking density, conversion rate, and residual monomers of the complex acrylate mixtures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structurally Driven Environmental Degradation of Friction in MoS2 Films

Abstract We report an investigation of the friction mechanisms of MoS 2 thin films under changing environments and contact conditions using a variety of computational and experimental techniques. Molecular dynamics simulations were used to study the effects of water and molecular oxygen on friction and bonding of MoS 2 lamellae during initial sliding. Characterization via photoelectron emission microscopy (PEEM) and Kelvin probe force microscopy (KPFM) were used to determine work function changes in shear modified material within the top few nanometers of MoS 2 wear scars. The work function was shown to change with contact conditions and environment, and shown by density functional theory (DFT) calculations and literature reports to be correlated with lamellae size and thickness of the basally oriented surface layer. Results from nanoscale simulations and macroscale experiments suggest that the evolution of the friction behavior of MoS 2 is linked primarily to the formation or inhibition of a basally oriented, molecularly thin surface film with long-range order.

Curry, John F.↗
DOE OSTI2020

Microwave-induced selective decomposition of cellulose: Computational and experimental mechanistic study

Understanding microwave-material interactions will help facilitate the utilization of microwave technology in gasification and renewable energy production. In this study, cellulose was used as a model compound to simulate the organic matter in biomass, and its catalytic decomposition under a microwave (MW) field was studied to identify structural changes from the reaction. The study was conducted using a MW source coupled to a fixed-bed gas-flow reactor, mass spectroscopic, and Fourier transform infrared spectroscopy post-reaction analysis. Zeolite 13X was chosen as a microwave absorber to study the catalytic enhancement of the decomposition of cellulose. Density functional theory (DFT) was used to gain insights into the molecular transformations occurring in the presence of a static electric field, which was used to simulate the electric field component of the microwave electromagnetic radiation. Theoretical calculations demonstrated that both the positive and negative portion of the electric field interact with the permanent dipoles of the cellulose effecting the decomposition through the glycosidic bond breaking mechanism. The theoretical result was verified using infrared spectroscopic analysis of the pure cellulose during microwave heating. The theoretical calculations suggest only the positive electric field component may be active for the cellulose decomposition in the presence of Zeolite 13X. Physically mixing Zeolite 13X with the cellulose led to a significant enhancement in the decomposition rate of the glycosidic bond. Zeolite 13X enhanced the glycosidic O-C decomposition at lower MW power (lower temperatures), whereas the O-H functional group required higher MW power (higher temperature) for its decomposition. The DFT study coupled with the reaction studies revealed that the electric field polarizability is dependent upon both the direction and the orientation of the cellulose. Gas products revealed that applying 250 W of MW power led to the production of CO, H 2 , along with some CO 2 , CH 4 , and benzene at 305 °C. Finally, reaction under 500, 750, and 1000 W of power at constant temperature (305 °C) revealed that higher power led to the complete decomposition of cellulose to mostly CO, H 2 , and CO 2 .

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

Identification of Carbonyl Species on Palladium Supported on Ceria in Complex Microenvironments

Herein, we present a systematic comparison between Pd carbonyl (Pd-CO) species, specifically over Pd/CeO 2 based catalysts, observed during isothermal adsorption and in several prototypical catalytic reactions to identify and understand CO adsorption on palladium-ceria based catalysts. Pd-CO is observed via DRIFTS to probe the gas-solid conditions, while ATR-IR is used to probe the affinity of Pd-CO under more complex solvated gas-solid-liquid conditions to discern the influence of the microenvironments for carbonyl adsorption. Here, we explore the presence of Pd-CO under several reactive environments, including CO adsorption, CO 2 + H 2 , CO + H 2 , CH 4 + CO 2 and CO under gas-solid-liquid media, highlighting reactions with notable Pd-CO formation. The differences between palladium carbonyls and carbonate species show that carbonyl species are much more affected via a shifting of the peak position than carbonates, which remain static irrespective of the immediate chemical environment. By following the rate of CO accumulation via K-M mode DRIFTS, we observe migration from linear, 2095 cm -1 , to bridge site, 1978 cm -1 , as a function of time under a static CO atmosphere. With the use of DFT, we discerned changes in Pd-carbonyl stretches due to both coverage effects of CO under simulated reaction conditions and temperature effects. Regardless of whether CO is formed as an intermediate or a reactant, the competitive adsorption of *H and *CO affects the binding strength of *CO at all temperatures, with low temperature favoring atop binding and high temperature favoring the more stable FCC Pd-CO site.

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

Final Technical Report for project entitled Highly Active, Durable, and Ultra-Low PGM NSTF Thin Film ORR Catalysts and Support

In this project, the objective was to develop new oxygen reduction reaction (ORR) electrocatalysts for proton exchange membrane fuel cells (PEMFCs) which could exceed all of the Department of Energy (DOE) 2020 targets listed in DE-FOA-0001224, Subtopic 1b, Table I. The expected outcome was development of one or more electrocatalysts which are substantially improved over the current state-of-the-art in terms of overall activity, durability, and cost, and are suitable for automotive traction and stationary fuel cell applications. This report summarizes this project’s progress towards meeting the stated objectives and includes a summary of key findings and conclusions. The development activities towards new highly active and durable thin film electrocatalysts has led to the discovery of several electrocatalysts which approached or exceeded several DOE 2020 targets for activity, durability, and performance in proton exchange membrane fuel cell (membrane electrode assemblies (MEAs). In this project, electrocatalyst development focused on systematic physical and electrochemical characterization of electrocatalyst activity, durability, and performance in MEAs as a function of electrocatalyst fabrication, compositional and structural variables. The project approach was to develop relationships between the catalyst functional responses (activity, durability and performance) with catalyst physical properties and catalyst fabrication methods. Additionally, the thin film catalysts were integrated onto the unique and durable 3M Nanostructured Thin Film (NSTF) support, consisting of arrays of self-assembled organic crystalline whiskers. The catalysts were evaluated against the project targets via extensive electrochemical characterization in MEA format at 3M, advanced structural and compositional microscopy at Oak Ridge National Laboratory (ORNL), and atomic structure analysis via XAFS at Argonne National Laboratory (ANL). The extensive electrochemical and physical characterization resulted in development of several trends which correlate the catalysts’ electrochemical properties to their physical properties, including composition, structure, and method of fabrication. The development was guided by density functional theory (DFT) modeling at Purdue University and kinetic Monte Carlo (kMC) modeling at Johns Hopkins University. The catalyst simulations provided key insights into the observed experimental activity and durability trends, and additionally were utilized to assess new electrocatalyst concepts prior to or coincident with physical catalyst development. This development has led to several electrocatalyst candidates with activity and durability which approach or exceed DOE targets. One class of catalysts, based on nm-scale thin layers of Pt on Ir, met or exceeded 6 of the 6 DOE targets the project addressed.

08 HYDROGEN↗
DOE OSTI2025

An atomic cluster expansion (ACE) potential for water under extreme conditions

We present a machine learning interatomic potential for water designed to capture its complex multiphase behavior, including both molecular and superionic ice phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to enable high-fidelity molecular dynamics simulations of water under extreme conditions, for pressures up to 100 GPa and for temperatures between 500 and 6000 K. A diverse range of configurations was generated through ab initio molecular dynamics (AI-MD) simulations, covering insulating and superionic ice phases, liquid water, and dissociated plasma phase. We demonstrate that the H 2 O ACE potential accurately reproduces experimental and DFT predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of water, including the transition from molecular fluid to the appropriate solid ice phases, and the superionic ice phases. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of water under extreme thermodynamic conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Molecular simulation using transfer-learned potentials for the disordered nanoscale structure of nitrogen-doped nanoporous carbons

Machine learning (ML)-based molecular dynamics (MD) simulations of the formation of a class of N-doped nanoporous carbons are performed to assess their disordered partially graphitized nanoscale structure. The study is motivated by the effectiveness of so-called nitrogen assembly carbons (NACs) for catalysis applications. Benchmark simulations for pure-C disordered graphitic systems reveal the importance of reliably capturing the vdW component of the potentials in order to accurately describe the tendency for layering of disordered graphene-like sheets. In our modeling, this is achieved by a transfer learning strategy incorporating features of the energetics from the optB88-vdW DFT functional into potentials initially trained with a less expensive functional, thereby providing a superior description of the pure-C systems. Generation from MD simulations of realistic partially graphitized structures is significantly more challenging for N-doped versus for pure C systems. However, such structures are achieved by a tailored MD simulation protocol mimicking the experimental synthesis process and in particular incorporating an annealing and subsequent quenching stages. Simulated PXRD patterns effectively reproduce the features of experimental observations for NACs, including the appearance of a prominent but broad (002) peak at around 25, and the development of another weaker feature associated with in-layer ordering of mixed C-N graphene-like sheets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗