Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic density functional theory”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Solvent oligomerization pathways facilitated by electrolyte additives during solid-electrolyte interphase formation

The solid–electrolyte interphase (SEI) layer formation is known to play an important role in determining the lifetime of lithium-ion batteries. A thin, stable SEI layer is linked to overall improved battery performance and longevity, however, the factors and mechanisms that lead to optimal SEI morphology and composition are not well understood. Inclusion of electrolyte additives (fluoroethylene carbonate, FEC; and vinylene carbonate, VC) is often necessary for improving SEI characteristics. To understand how these electrolyte additives impact SEI formation, we employed molecular dynamics (MD) and density functional theory (DFT) simulations to study the reaction networks and oligomerization pathways, respectively, for three systems containing ethylene carbonate (EC), a lithium ion, and FEC or VC. MD simulations suggest radical oligomerization pathways analogous to traditional oligomerization with nucleophilic alkoxide species via SN1 reaction mechanisms. Both S N 1 and S N 2 mechanisms were studied for all three systems using DFT. Oligomerization reactions were studied with both a standard alkoxide species and a ring-opened EC radical as the nucleophiles and EC, FEC, and VC as the electrophiles. For all cases, FEC and VC exhibited lower free energy barriers and more stable adducts when compared with EC. We conclude that one of the role of additives is to modify the oligomerization process of EC by introducing branching points (FEC) or termination points (VC).

25 ENERGY STORAGE↗
DOE OSTI2026

Revealing EDL-driven reduction mechanisms in binary, ternary, and quaternary fluorinated electrolytes via an integrated MD–DFT–ML framework

Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.

DFT-MD-ML workflow↗
DOE OSTI2023

A High-Speed Rotational Diamond Anvil Cell for In Situ Analysis of Hierarchical Microstructural Evolution of Metallic Alloys during Extreme Shear Deformation

High speed shear deformation is ubiquitous in engineering applications, ranging from material processing methods such as friction stir processing/extrusion and in tribological contacts. However, analyzing the microstructural evolution of materials while they are undergoing high speed shear deformation have been a long-standing challenge. This led to predominant reliance on ex situ microscopy before and after shear deformation. But ex situ microscopy lacks the ability to analyze dynamic and transient hierarchical microstructural evolution mechanisms that could occur during shear deformation of materials. Therefore, to better understand the dynamic mechanisms of mass and energy transfer in materials under shear deformation, we developed a first of its kind high-speed rotational diamond anvil cell (HS-RDAC) for synchrotron-based in situ high-energy x-ray diffraction (XRD). We studied the time resolved lattice strain evolution, XRD peak broadening and changes in spatial variation of shear deformation induced alloying in pure metal and metal alloy sheets and powder mixture using the HS-RDAC. These in situ results were combined with detailed ex situ microstructural characterization before and after the shear deformation using transmission electron microscopy and atom probe tomography, which revealed the different stages of evolution of a shear deformation induced hierarchical nanostructure. Multiscale computational simulations including computational fluid dynamics, crystal plasticity, molecular dynamic simulation and density functional theory uncovered the mechanisms behind morphological changes, evolution of defect structures and changes in driving force for shear deformation induced intermixing. In conclusion, this in situ HS-RDAC capability, in combination with ex situ microstructural characterization and computational simulations, can provide new insights into the hierarchical microstructural evolution pathway during shear deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Carbon-carbon interactions in warm dense titanium carbide

Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.

Ab initio molecular dynamics↗
DOE OSTI2020

Correlation strength and orbital differentiation across the phase diagram of plutonium metal

We compare the trends on the strength of electronic correlations across the different phases of elemental Pu focusing on its site and orbital dependence, using a combination of density functional theory (DFT) and dynamical mean-field theory (DMFT) calculations within the vertex corrected one crossing approximation. We find that Pu-5f states are more correlated in δ -Pu, followed by some crystallographic sites in α and β phases. In addition, we observe that Pu-5f 5/2 and Pu-5f 7/2 orbital differentiation is a general feature of this material, as is site differentiation in the low-symmetry phases. The Pu-5f 5/2 states show Fermi liquid like behavior, whereas the Pu-5f 7/2 states remaining incoherent down to very low temperatures. We correlate the correlation strength in the different phases to their structure and the Pu-5 f occupancy of their crystallographic sites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

DFT + DMFT study of dopant effects in the heavy-fermion compound CeCoIn 5

Here we study the dopant-induced inhomogeneity effect on the electronic properties of heavy fermion CeCoIn 5 using a combined approach of density functional theory (DFT) and dynamical mean-field theory (DMFT). The inhomogeneity of the hybridization between Ce-4f and conduction electrons is introduced to impose the inequivalent Ce atoms with respect to the dopant. From the DFT to the DFT + DMFT results, we demonstrate a variation of the hybridization strength depending on the hole or electron doping. A drastic asymmetric mass renormalization could be reproduced in the DFT + DMFT calculation. Finally, the calculated Kondo temperature reflects the different development of the heavy quasiparticle states, depending on the dopant.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Role of orbital selectivity on crystal structures and electronic states in BiMnO 3 and LaMnO 3 perovskites

Correlated oxides, such as BiMnO 3 and LaMnO 3 , show a complex interplay of electronic correlations and crystal structure exhibiting multiple first-order phase transitions, some without a clear order parameter. The quantitative theoretical description of this temperature-dependent electronic-structural interplay in the vicinity of a Mott transition is still a challenge. Here we address this issue by simultaneously considering both structural and electronic degrees of freedom, within a self-consistent density functional theory with embedded dynamical mean field theory. Our results show the existence of novel electronic states characterized by the coexistence of insulating, semimetallic, and metallic orbitals. This state is in agreement with resonant x-ray scattering. Here we also show that electronic entropy plays a decisive role in both electronic and structural phase transitions. By self-consistent determination of both the electronic state and the corresponding crystal structure, we show that the temperature evolution of these phases can be quantitatively explained from first principles, thus demonstrating the predictive power of the theoretical method for both the structural and the electronic properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Understanding the lattice thermal conductivity of SrTiO 3 from an ab initio perspective

Here, we present a detailed analysis of the structure dependence of the lattice thermal conductivity of SrTiO 3 . We have used both ab initio molecular dynamic simulations and density functional theory calculations to decouple the effect of different structural distortions on the thermal conductivity. We have identified two main mechanisms for tuning the thermal conductivity when a distortion is applied. First, the modification of the acoustic-modes energy dispersion when a change in the lattice parameters is imposed and second, the low energy polar modes. In particular and counterintuitively, we have found that an increase in the angle of the oxygen octahedral rotations increases the thermal conductivity due to its coupling to these polar modes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reversible phase transformations between Pb nanocrystals and a viscous liquid-like phase

Phase transformations have been a prominent topic of study for both fundamental and applied science. Solid-liquid reaction–induced phase transformations can be hard to characterize, and the transformation mechanisms are often not fully understood. Here, we report reversible phase transformations between a metal (Pb) nanocrystal and a viscous liquid-like phase unveiled by in situ liquid cell transmission electron microscopy. The reversible phase transformations are obtained by modulating the electron current density (between 1000 and 3000 electrons Å –2 s –1 ). The metal-organic viscous liquid-like phase exhibits short-range ordering with a preferred Pb-Pb distance of 0.5 nm. Assisted by density functional theory and molecular dynamics calculations, we show that the viscous liquid-like phase results from the reactions of Pb with the CH 3 O fragments from the triethylene glycol solution under electron beam irradiation. Such reversible phase transformations may find broad implementations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal Neutron Scattering Research and Development at Oak Ridge National Laboratory [Slides]

This presentation provides an overview of the project, specifically describing how neutron cross sections can be broken up into several regions. It also discusses how thermal scattering law (TSL) $\textit{S(α, β)}$ is related to dynamic structure factor $\textit{S(q, E)}$ and how $\textit{S(q, E)}$ can be broken up into 2 primary components, incoherent and coherent, and how it can be calculated in several ways, including density functional theory (DFT) and molecular dynamics (MT). Additionally, it states that ENDF currently stores coherent elastic (Bragg edges), incoherent elastic and incoherent inelastic. Also discussed are the experimental capabilities, evaluations, and uncertainty quantification of the project.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Bidentate Lewis bases are preferred for passivation of MAPbI 3 surfaces: A time-domain ab initio analysis

We report chemical passivation is a predominant approach to inhibit intrinsic defects responsible for electron–hole recombination in perovskite solar cells. Using time-domain density functional theory combined with nonadiabatic molecular dynamics, we demonstrate defect passivation by separation of electrons and holes in metal halide perovskites (MHPs), and show that bidentate ligands exhibit the best performance. Defects in traditional semiconductors create deep midgap states, and passivation eliminates these states. In contrast, common defects produce no deep midgap states in MHPs. Instead, defects localize electrons and holes around defect sites, enhancing electron–hole interaction. Defect passivation in MHPs acts to separate charges, decreasing electron–hole and charge–phonon interactions, and increasing charge lifetimes. Bidentate ligands work best, because they can passivate both unsaturated chemical bonds created due to vacancy defects. Bidentate ligands with spatially separated binding sites are preferred, since they provide better match to the sparse inorganic lattice of MHPs. Similar to the Lewis base ligands, water also acts as a ligand and extends charge lifetimes, with the oxygen atom donating its lone electron pair to the defect site. However, water accelerates chemical degradation of MHPs. The binding energy of most ligands is larger than that of water, and therefore, the ligands displace water, increase MHP stability and prolong carrier lifetimes. The established physical mechanism of defect passivation, and the specific principles guiding the choice of passivating molecules advance our understanding of the exceptional properties of MHPs and suggest routes for further improvement of MHP performance in solar energy and optoelectronic applications.

36 MATERIALS SCIENCE↗
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 tungsten doping on nonradiative electron–hole recombination in monolayer MoSe 2 with Se vacancies

Two-dimensional transition metal dichalcogenides (TMDs) are receiving significant attention due to their excellent electronic and optoelectronic properties. The material quality is greatly affected by defects that are inevitably generated during material synthesis. Focusing on chalcogenide vacancies, which constitute the most common defect, we use the state-of-the-art simulation methodology developed in our group to demonstrate that W doping of MoSe 2 with Se vacancies reduces charge carrier losses by two mechanisms. First, W doping makes the formation of double Se vacancies unfavorable, while it is favorable in undoped MoSe 2 . Second, if a Se vacancy is present, the charge carrier lifetimes are extended in the W-doped MoSe 2 . Combining ab initio real-time time-dependent density functional theory with nonadiabatic molecular dynamics, the simulations show that the nonradiative carrier losses in the presence of Se vacancies proceed by sub-10 ps electron trapping and relaxation down the manifold of trap states, followed by a 100 ps recombination of trapped electrons with free holes. The electron–vibrational energy exchange is driven by both in-plane and out-of-plane vibrational motions of the MoSe 2 layer. Additionally, the atomistic studies advance our understanding of the influence of defects on charge carrier properties in TMDs and guide improvements of material quality and development of TMD applications.

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

Band Structure-Driven Design of a α-CsPbI 3 Ammonia Sensor for Industrial Applications

We investigate the defect-dependent electronic structure and gas-sensing potential of cubic α-CsPbI 3 using first-principles density functional theory and nonadiabatic molecular dynamics. Among the intrinsic defects, interstitials, vacancies, antisites, and switches studied, the I Pb and Pb I antisite defects exhibit transition energy levels near the middle of the band gap, thus functioning as deep traps. Short-term adsorption of ammonia selectively modifies the electronic structure, coordinating with Pb at Pb I sites and Cs at I Pb sites, significantly altering recombination pathways. Detailed analysis reveals that NH 3 reduces anharmonicity at I Pb defects, enabling enhanced recombination at elevated temperatures, while trap-assisted recombination dominates at room temperature. Other analytes, including CH 3 NH 2 and NO 2 , show negligible impact on the band gap or recombination dynamics, highlighting the potential selectivity of NH 3 interactions. Ab initio nonadiabatic molecular dynamics simulations at 300 K and 600 K further demonstrate temperature-dependent modulation of carrier lifetimes, with NH 3 accelerating recombination at ambient conditions and suppressing certain pathways at higher temperatures. These findings suggest that α-CsPbI 3 can serve as a selective and sensitive ammonia sensor over a broad temperature range and offer insights for ammonia detection under industrially relevant conditions.

ammonia sensor↗
DOE OSTI2023

Atomistic simulations for investigation of substrate effects on lipid in-source fragmentation in secondary ion mass spectrometry

In beam-based ionization methods, the substrate plays an important role on the desorption mechanism of molecules from surfaces. Both the specific orientation that a molecule adopts at a surface and the strength of the molecule-surface interaction can greatly influence desorption processes, which in turn will affect the ion yield and the degree of in-source fragmentation of a molecule. In the beam-based method of secondary ion mass spectrometry (SIMS), in-source fragmentation can be significant and molecule specific due to the hard ionization method of using a primary ion beam for molecule desorption. To investigate the role of the substrate on orientation and in-source fragmentation, we have used atomistic simulations—molecular dynamics in combination with density functional theory calculations—to explore the desorption of a sphingolipid (palmitoylsphingomyelin) from a model surface (gold). We then compare SIMS data from this model system to our modeling findings. Using this approach, we found that the combined adsorption and binding energy of certain bonds associated with the headgroup fragments (C3H8N+, C5H12N+, C5H14NO+, and C5H15PNO4+) was a good predictor for fragment intensities (as indicated by relative ion yields). This is the first example where atomistic simulations have been applied in beam-based ionization of lipids, and it presents a new approach to study biointerfacial lipid ordering effects on SIMS imaging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗