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

Band gap renormalization, carrier mobilities, and the electron-phonon self-energy in crystalline naphthalene

Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. Here, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and the resulting scattering lifetimes of electronic states near the valence- and conduction-band edges employing density functional theory. Further, our calculated phonon renormalization of the GW-corrected quasiparticle band structure predicts a fundamental band gap of 5 eV for naphthalene at room temperature, in good agreement with experiments. From our calculated phonon-induced electron lifetimes, we obtain the temperature-dependent mobilities of electrons and holes in good agreement with experimental measurements at room temperature. Finally, we show that an approximate energy self-consistent computational scheme for the electron-phonon self-energy leads to the prediction of strong satellite bands in the electronic band structure. We find that a single calculation of the self-energy can reproduce the self-consistent results of the band gap renormalization and electrical mobilities for naphthalene, provided that the on-the-mass-shell approximation is used, i.e., if the self-energy is evaluated at the bare eigenvalues.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

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

Predicted thermophysical properties of UN, PuN and (U,Pu)N

Molecular dynamics and density functional theory simulations are used to predict the lattice and electronic contributions of thermophysical properties for UN, PuN, and mixed (U,Pu)N systems. The properties predicted include the lattice parameter, linear thermal expansion, enthalpy, and specific heat capacity, as a function of temperature. The simulation predictions for high temperature specific heat capacity are compared against experimental measurements to understand the behavior, and why differences in the experimental measurements are observed. The influence of adding U vacancies, N interstitials, and Pu to UN is also examined. For this, a new PuN potential parameter set is developed and used with the Kocevski UN potential, enabling the dynamics of mixed (U,Pu)N systems to be studied. How defects impact the thermophysical properties is important for understanding fuel behavior under different reactor conditions, and these mechanistic predictions can be used to support fuel performance codes where data is scarce.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Photoemission and dynamical mean field theory study of electronic correlations in a t 2 g 5 metal SrRhO 3 thin film

Perovskite rhodates are characterized by intermediate strengths of both electronic correlation as well as spin-orbit coupling (SOC) and usually behave as moderately correlated metals. A recent publication on epitaxial SrRhO 3 thin films reported a bad-metallic behavior and suggested the occurrence of antiferromagnetism below 100 K. In this work, we have further studied this SrRhO 3 thin film by hard x-ray photoemission spectroscopy and found a very small density of states (DOS) at the Fermi level $(E_F)$, which is consistent with the previously reported bad-metallic behavior. However, the absence of DOS persists up to room temperature, which contradicts the explanation of antiferromagnetic transition at ~100 K. We also employed electronic structure calculations within the framework of density functional theory and dynamical mean-field theory. In contrast to the photoemission results, our calculations indicate metallic behaviors of both bulk SrRhO 3 and the SrRhO 3 thin film, and a stronger correlation effect was observed in the thin film than that in the bulk. The calculated uniform magnetic susceptibility is substantially larger in the thin film than that in the bulk. We also investigated the role of SOC and found only a moderate modulation of the band structure. Hence SOC is not expected to significantly affect the electronic correlation in SrRhO 3 . Extrinsic effects of finite-thickness effects beyond our calculations and localization effects may play important roles to induce the negligible spectral weight at $(E_F)$.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding Solid Electrolyte Interphase Nucleation and Growth on Lithium Metal Surfaces

Experiments and theory are needed to decode the exact structure and distribution of components of a passivation layer formed at the anode surface of Li metal batteries, known as the Solid Electrolyte Interphase (SEI). Due to the inherent dynamic behavior as well as the lithium reactivity, the SEI structure and its growth mechanisms are still unclear. This study uses molecular simulation and computational chemistry tools to investigate the initial nucleation and growth dynamics of LiOH and Li2O that provide us with thermodynamics and structural information about the nucleating clusters of each species. Following the most favorable pathways for the addition of each of the components to a given nascent SEI cluster reveals their preferential nucleation mechanisms and illustrates different degrees of crystallinity and electron density distribution that are useful to understand ionic transport through SEI blocks.

25 ENERGY STORAGE↗
DOE OSTI2022

Low-energy interband transition in the infrared response of the correlated metal SrVO 3 in the ultraclean limit

We studied the low-energy electronic response of the prototypical correlated metal SrVO 3 in the ultraclean and disordered limit using infrared spectroscopy and density functional theory plus dynamical mean field theory calculations (DFT+DMFT). A strong optical excitation at 70 meV is observed in the optical response of the ultraclean samples but is hidden by the low-energy Drude-like response from intraband excitations in the more disordered samples. DFT+DMFT calculations reveal that this optical excitation originates from interband transitions between the bands split by orbital off-diagonal hopping, which has often been ignored in cubic systems, such as SrVO 3 . A memory function analysis of the optical data shows that this interband transition can lead to deviations of optical self-energy from the expected Fermi-liquid behavior. Our findings demonstrate that analysis schemes employed to extract many-body effects from optical spectra may be oversimplified to study the true electronic ground state and that improvements in material quality can guide efforts to refine theoretical approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Charge correlations and magnetoelastic coupling in intercalated transition metal dichalcogenides

The large van der Waals gap in transition metal dichalcogenides (TMDs) offers an avenue to tune the ground state of 2D materials through the intercalation of magnetic atoms. Here, we investigate the charge correlations in Fe 1/3 ⁢TaS 2 , Co 1/3 ⁢TaS 2 , and Fe 0.35 ⁢NbS 2 by combining angle-resolved photoemission spectroscopy (ARPES), x-ray scattering, magnetometry, and density functional theory (DFT). We find that, while short-range charge fluctuations develop in Ta-based compounds, Fe 0.35 ⁢NbS 2 exhibits long-range charge order which is strongly coupled with magnetic order and tunable by external magnetic field. Our electronic structure analysis reveals that intercalation reconstructs the Fermi surface via charge transfer and band renormalization, yet does not generate the nesting conditions compatible with the observed ordering vectors. Complementary phonon calculations further exclude a conventional electron-phonon origin of charge order. Together, these results establish magnetoelastic coupling as the dominant mechanism behind charge ordering in Fe 0.35 ⁢NbS 2 and highlight the contrasting role of Nb and Ta hosts in stabilizing correlated ground states in intercalated TMDs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formation of pyrophosphates across grain boundaries induces the formation of mismatched but oriented interfaces in silver phosphate polypods

The interfaces and their misfit defects determine the materials properties for a wide range of applications, such as electronic devices, photocatalysis, and mechanical engineering, etc. However, current understanding of atomic interfacial structures is limited. Here we discover a special interfacial structure, mismatched but oriented interface via two distinct facets. Transmission and scanning electron microscopy results suggest that, in Ag 3 PO 4 polypods structures, interfaces of {1 0 0} and {1 1 0}, ({1 0 0}/(1 1 0}), {1 0 0}/{1 1 1}, {1 1 0}/{1 1 1}, and {1 0 0}/{1 0 0}, etc., have a certain orientation relationship, corresponding to the energy minima and coincident site lattice of interfacial atoms as demonstrated by molecular dynamics simulations. Density functional theory demonstrates that the formation of pyrophosphate and/or phosphates rotation to bond across the interface compensate the lattice mismatch at the interfaces, as well as deformations of Ag-O bonds. Furthermore, our work opens up a new avenue for a much wider range of interfacial structures, allow for a higher diversity of structures, and shine light on tailoring crystal structures, morphologies, and the resulting properties.

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

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Controlled Formation of Stacked Si Quantum Dots in Vertical SiGe Nanowires

In this work, we demonstrate the ability to fabricate vertically stacked Si quantum dots (QDs) within SiGe nanowires with QD diameters down to 2 nm. These QDs are formed during high-temperature dry oxidation of Si/SiGe heterostructure pillars, during which Ge diffuses along the pillars’ sidewalls and encapsulates the Si layers. Continued oxidation results in QDs with sizes dependent on oxidation time. The formation of a Ge-rich shell that encapsulates the Si QDs is observed, a configuration which is confirmed to be thermodynamically favorable with molecular dynamics and density functional theory. The type-II band alignment of the Si dot/SiGe pillar suggests that charge trapping on the Si QDs is possible, and electron energy loss spectra show that a conduction band offset of at least 200 meV is maintained for even the smallest Si QDs. Our approach is compatible with current Si-based manufacturing processes, offering a new avenue for realizing Si QD devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Insights into Spontaneous Solid Electrolyte Interphase Formation at Magnesium Metal Anode Surface from Ab Initio Molecular Dynamics Simulations

Spontaneous chemical reactivity at multivalent (Mg, Ca, Zn, Al) electrode surfaces is critical to solid electrolyte interphase (SEI) formation, and hence, directly affects the longevity of batteries. Here, we report an investigation of the reactivity of 0.5 M Mg(TFSI)2 in 1,2-dimethoxyethane (DME) solvent at a Mg(0001) surface using ab initio molecular dynamics (AIMD) simulations and detailed Bader charge analysis. Based on the simulations, the initial degradation reactions of the electrolyte strongly depend on the structure of the Mg(TFSI)2 species near the anode surface. At the surface, the dissociation of Mg(TFSI)2 species occurs via cleavage of the N-S bond for the solvent separated ion pair (SSIP) and via cleavage of the C-S bond for the contact ion pair (CIP) configuration. In the case of the CIP, both TFSI anions undergo spontaneous bond dissociation reactions to form atomic O, C, S, F, and N species adsorbed on the surface of the Mg anode. These products indicate that the initial SEI layer formed on the surface of the pristine Mg anode consists of a complex mixture of multiple components such as oxides, carbides, sulfides, fluorides, and nitrides. We believe that the atomic level insights gained from these simulations will lay the groundwork for the rational design of tailored and functional interphases that are critical for the success of multivalent battery technology.

Agarwal, Garvit↗
DOE OSTI2022

Understanding Hydrogenation Chemistry at MgB 2 Reactive Edges from Ab Initio Molecular Dynamics

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB 2 without a priori assumption of reaction pathways. Focusing on highly reactive ($10\overline{1}0$) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH 4 – units and key chemical intermediates, involving H 2 dissociation, generation of functionalities and molecular complexes containing BH 2 and BH 3 motifs, and B–B bond breaking. Further, the genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B–H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg 2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB 2 -based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

08 HYDROGEN↗
DOE OSTI2025

Coordinated Cation Transport in Ti 3 C 2 T x MXene Membranes

Membrane nanofiltration is an attractive strategy for the selective recovery of high-demand metals from wastewater and brine. Effective sieving of ions in aqueous environments will require precise control over membranes’ nanochannel size and chemistry. Ti 3 C 2 T x MXene is an environmentally stable 2D material that can be processed into laminar membranes containing nanoscale interlayer spaces. The MXene interlayer environment depends on the ion species and amount of water intercalated between MXene sheets, and it is the major factor governing permeation and selectivity through MXene membranes. Coordinated ion–ion and ion-interlayer dynamics in the presence of complex mixtures can impact ion permeability and selectivity. Herein, we observe strong competitive effects between different cations (Li + , Na + , and Ca 2+ ) in binary mixtures, resulting in reduced selectivity when compared with single-salt permeability ratios. X-ray diffraction, molecular dynamics, and density functional theory simulations support the conclusion that cations with stronger attraction to MXene flakes can preferentially occupy the MXene nanochannels and hinder other ions via charge or size exclusion. In conclusion, elucidation of ion transport behavior in MXene under complex conditions will allow for more rational design of efficient ion-sieving membranes.

2D Materials↗
DOE OSTI2022

Mapping 1D Confined Electromagnetic Edge States in 2D Monolayer Semiconducting MoS 2 Using 4D-STEM

Here, four-dimensional (4D) scanning transmission electron microscopy is used to study the electric fields at the edges of 2D semiconducting monolayer MoS2. Sub-nanometer 1D features in the 2D electric field maps are observed at the outermost region along zigzag edges and also along nanowire MoS-terminated MoS 2 edges. Atomic-scale oscillations are detected in the magnitude of the 1D electromagnetic edge state, with spatial variations that depend on the specific periodic edge reconstructions. Electric field reconstructions, along with integrated differential phase contrast reconstructions, reveal the presence of low Z number atoms terminating many of the uniform edges, which are difficult to detect by annular dark field scanning transmission electron microscopy due to its limited dynamic range. Density functional theory calculations support the formation of periodic 1D edge states and also show that enhancement of the electric field magnitude can occur for some edge terminations. The experimentally observed electric fields at the edges are attributed to the absence of an opposing electric field from a nearest neighbor atom when the electron beam propagates through the 2D monolayer and interacts. These results show the potential of 4D-STEM to map the atomic scale structure and fluctuations of electric fields around edge atoms with different bonding states than bulk atoms in 2D materials, beyond conventional imaging.

2D materials↗
DOE OSTI2020

Probing Heterogeneous Charge Distributions at the α-Al 2 O 3 (0001)/H 2 O Interface

Unlike metal or semiconductor electrodes, the surface charge resulting from the protonation or deprotonation of insulating mineral oxides is highly localized and heterogeneous in nature. In this work the Stark active C≡N stretch of potassium thiocyanate is used as a molecular probe of the heterogeneity of the interfacial electrostatic potential at the α-Al 2 O 3 (0001)/H 2 O interface. Vibrational sum frequency generation (vSFG) measurements performed in the OH stretching region suggest that thiocyanate species organize interfacial water similarly to halide ions. Changes in the electrostatic potential are then tracked via Stark shifts of the vibrational frequency of the thiocyanate stretch. Our vSFG measurements show that we can simultaneously measure the vSFG response of SCN – ions experiencing charged and neutral surface sites. We assign local potentials of +308 and -154 mV to positively and negatively charged aluminol groups that are present at pH = 4 and pH = 10, respectively. Thiocyanate anions at positively charged surface sites and negatively charged surface sites and those participating in contact ion pairing adopt similar orientations and are oppositely oriented relative to thiocyanate ions near neutral surface sites. All four species followed Langmuir adsorption isotherms. Density functional theory–molecular dynamics (DFT-MD) simulations of SCN – near the neutral α-Al 2 O 3 (0001)/H 2 O interface show that the vSFG response in the C≡N stretch region originates from a SCN–H–O–Al complex, suggesting the surface site specificity of these experiments. To our knowledge this is the first spectroscopic measurement of local potentials associated with a heterogeneously charged surface. The ability to probe the evolution of local charges in situ could provide vital insight into many industrial, electrochemical, and geochemically relevant interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗