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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effective Storage of Electrons in Water by the Formation of Highly Reduced Polyoxometalate Clusters

Aqueous solutions of polyoxometalates (POMs) have been shown to have potential as high-capacity energy storage materials due to their potential for multi-electron redox processes, yet the mechanism of reduction and practical limits are currently unknown. Herein, we explore the mechanism of multi-electron redox processes that allow the highly reduced POM clusters of the form {MO 3 } y to absorb y electrons in aqueous solution, focusing mechanistically on the Wells-Dawson structure X 6 [P 2 W 18 O 62 ], which comprises 18 metal centers and can uptake up to 18 electrons reversibly (y = 18) per cluster in aqueous solution when the countercations are lithium. This unconventional redox activity is rationalized by density functional theory, molecular dynamics simulations, UV-vis, electron paramagnetic resonance spectroscopy, and small-angle X-ray scattering spectra. These data point to a new phenomenon showing that cluster protonation and aggregation allow the formation of highly electron-rich meta-stable systems in aqueous solution, which produce H 2 when the solution is diluted. Finally, we show that this understanding is transferrable to other salts of [P 5 W 30 O 110 ] 15- and [P 8 W 48 O 184 ] 40- anions, which can be charged to 23 and 27 electrons per cluster, respectively.

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

Discovery of Stacking Heterogeneity, Layer Buckling, and Residual Water in COF-999-NH 2 and Implications on CO 2 Capture

Covalent organic frameworks (COFs), with their modular architectures and tunable functionalities, provide a versatile platform to design sorbents for the direct capture of CO 2 from air. Here, for this work, we combined density functional theory, molecular dynamics, and grand canonical Monte Carlo simulations with experiment to understand structural factors for furthering COF-999-NH 2 ’s performance as the precursor to COF-999 for direct air CO 2 capture. Small energy differences among laterally shifted stackings suggest intrinsic stacking heterogeneity. The simulations show pronounced layer buckling coupled to extensive amine–nitrile hydrogen bonding and persistent pore water, which initiates undesired polymerization and undermines uptake. The predicted presence of water is confirmed by subsequent experiments. These insights point to a single, actionable design rule: exclude retained water by introducing hydrophobic pore environments to maximize the CO 2 capture efficiency.

adsorption↗
DOE OSTI2024

A topological Hund nodal line antiferromagnet

The interplay of topology, magnetism, and correlations gives rise to intriguing phases of matter. In this study, through state-of-the-art angle-resolved photoemission spectroscopy, density functional theory, and dynamical mean-field theory calculations, we visualize a fourfold degenerate Dirac nodal line at the boundary of the bulk Brillouin zone in the antiferromagnet YMn 2 Ge 2 . We further demonstrate that this gapless, antiferromagnetic Dirac nodal line is enforced by the combination of magnetism, space-time inversion symmetry, and nonsymmorphic lattice symmetry. The corresponding drumhead surface states traverse the whole surface Brillouin zone. YMn 2 Ge 2 thus serves as a platform to exhibit the interplay of multiple degenerate nodal physics and antiferromagnetism. Interestingly, the magnetic nodal line displays a d-orbital dependent renormalization along its trajectory in momentum space, thereby manifesting Hund's coupling. Our findings offer insights into the effect of electronic correlations on magnetic Dirac nodal lines, leading to an antiferromagnetic Hund nodal line.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electrical conductivities of (Mg,Fe)O at extreme pressures and implications for planetary magma oceans

During planet formation, planets undergo many impacts that can generate magma oceans. When these crystallize, part of the magma densifies via iron enrichment and migrates to the core–mantle boundary, forming an iron-rich basal magma ocean (BMO). The BMO could generate a dynamo in early Earth and super-Earths if the electrical conductivity of the BMO, which is thought to be sensitive to its Fe content, is sufficiently high. To test this hypothesis, here we conduct laser-driven shock experiments on ferropericlase (Mg x ,Fe 1−x )O (0.95 ≤ x ≤ 1) as an Fe-rich BMO analogue, perform density functional theory molecular dynamics simulations on MgO and calculate the long-term evolution of super-Earths. We find that the d.c. conductivities of MgO and (Mg,Fe)O are indistinguishable between 467 GPa and 1,400 GPa, despite previous predictions. Here, we predict that super-Earths larger than 3–6 Earth masses can produce BMO-driven dynamos that are almost one order of magnitude stronger than core-driven dynamos for several billion years.

Exoplanets↗
DOE OSTI2020

Surface softening in palladium nanoparticles: effects of a capping agent on vibrational properties

The presence of a capping agent (CTAB) on Pd nanoparticles produces a strong static disorder in the surface region. This results in a surface softening, which contributes to an overall increase in the Debye–Waller coefficient measured by X-ray powder diffraction. Here, molecular dynamics and density functional theory simulations show that the adsorption-induced surface disorder is strong enough to overcome the effects of nanoparticle size and shape.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Unexpected electrochemical behavior of an anolyte redoxmer in flow battery electrolytes: solvating cations help to fight against the thermodynamic–kinetic dilemma

Redoxmers are redox-active molecules that can store energy in electrolytes for redox flow batteries (RFBs), and their electrochemical properties are significantly affected by the choice of supporting electrolytes. Herein, we use 2,1,3-benzothiadiazole (BzNSN) as a model system to scrutinize the supporting electrolyte impact. By systemically varying the components of supporting salts, BzNSN not only shows substantial redox potential shifts but also exhibits varying electrochemical stabilities. Specifically, changing the size of cations can effectively alter the coordination between the supporting salt and BzNSN species. From Li + , Na + , K + , to NEt 4 + , the redox potential of BzNSN shifts negatively, from -1.63 V to -1.82 V vs. Ag/Ag + . Molecular dynamics and density functional theory simulations revealed that smaller cations, like Li + , are closer to the charged BzNSN when coordinated, implying stronger coordination, while larger cations, like K + and NEt 4 + , are farther away. Interestingly, the large cation electrolytes also lead to much improved electrochemical stability, evidenced by the extraordinarily enhanced kinetic lifetime from electron paramagnetic resonance measurement. This study demonstrates the first example of tuning an anolyte redoxmer toward a concurrent improvement of lowered redox potentials AND enhanced calendar lives via solvation means, which is usually constrained by the thermodynamic–kinetic relation.

2,1,3-benzothiadiazole↗
DOE OSTI2022

The solution structures and relative stability constants of lanthanide–EDTA complexes predicted from computation

Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ~7 and ~11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ~11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ~7 and ~11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. As a result, increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.

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

Explaining the structure sensitivity of Pt and Rh for aqueous-phase hydrogenation of phenol

Phenol is an important model compound to understand the thermocatalytic (TCH) and electrocatalytic hydrogenation (ECH) of biomass to biofuels. Although Pt and Rh are among the most studied catalysts for aqueous-phase phenol hydrogenation, the reason why certain facets are active for ECH and TCH is not fully understood. Herein, we identify the active facet of Pt and Rh catalysts for aqueous-phase hydrogenation of phenol and explain the origin of the size-dependent activity trends of Pt and Rh nanoparticles. Phenol adsorption energies extracted on the active sites of Pt and Rh nanoparticles on carbon by fitting kinetic data show that the active sites adsorb phenol weakly. We predict that the turnover frequencies (TOFs) for the hydrogenation of phenol to cyclohexanone on Pt(111) and Rh(111) terraces are higher than those on (221) stepped facets based on density functional theory modeling and mean-field microkinetic simulations. The higher activities of the (111) terraces are due to lower activation energies and weaker phenol adsorption, preventing high coverages of phenol from inhibiting hydrogen adsorption. We measure that the TOF for ECH of phenol increases as the Rh nanoparticle diameter increases from 2 to 10 nm at 298 K and -0.1 V vs the reversible hydrogen electrode, qualitatively matching prior reports for Pt nanoparticles. The increase in experimental TOFs as Pt and Rh nanoparticle diameters increase is due to a larger fraction of terraces on larger particles. These findings clarify the structure sensitivity and active site of Pt and Rh for the hydrogenation of phenol and will inform the catalyst design for the hydrogenation of bio-oils.

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

Effects of carbon concentration on the local atomic structure of amorphous GST

Ge-Sb-Te (GST) alloys are leading phase-change materials for data storage due to the fast phase transition between amorphous and crystalline states. Ongoing research aims at improving the stability of the amorphous phase to improve retention. This can be accomplished by the introduction of carbon as a dopant to Ge 2 Sb 2 Te 5 , which is known to alter the short- and mid-range structure of the amorphous phase and form covalently bonded C clusters, both of which hinder crystallization. The relative importance of these processes as a function of C concentration is not known. Here we used molecular dynamics simulation based on density functional theory to study how carbon doping affects the atomic structure of GST-C. Carbon doping results in an increase in tetrahedral coordination, especially of Ge atoms, and this is known to stabilize the amorphous phase. We observe an unexpected, non-monotonous trend in the number of tetrahedral bonded Ge with the amount of carbon doping. Our simulations show an increase in the number of tetrahedral bonded Ge up to 5 at.% C, after which the number saturates and begins to decrease above 14 at.% C. The carbon atoms aggregate into clusters, mostly in the form of chains and graphene flakes, leaving less carbon to disrupt the GST matrix at higher carbon concentrations. Different degrees of carbon clustering can explain divergent experimental results for recrystallization temperature for carbon doped GST.

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

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

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

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

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

Thermal exchange-correlation functionals: Capturing quantum electron behavior in warm, dense plasmas

We summarize and give perspective upon recent progress in developing non-empirical constraint-based thermal (i.e., free energy) exchange-correlation (XC) density functionals essential for accurate description of the quantum behavior of electrons in warm, dense plasmas. After delineating the critical role of ground-state functionals for zero-temperature, time-dependent DFT, we outline the underpinnings of local density approximation, generalized gradient approximation (GGA), and meta-GGA XC free-energy functionals. Two basic thermalization principles for upgrading ground-state XC functionals to successful thermal ones are emphasized. Then, we turn to a long-standing challenge, assessment of the accuracy of well-founded functionals. Unlike the ground state, there are a few exact results for large T and P. An exception is path integral Monte Carlo (PIMC) data for dense H/D and He plasmas. For those, we did ab initio molecular dynamics simulations under selected thermodynamic conditions employing five thermal XC functionals: two approximate thermal GGAs, fully thermal GGA, an approximate meta-GGA, and fully thermal meta-GGA. Comparisons with the PIMC data show that functionals thermalized by augmenting a non-thermal functional with a lower-level thermal contribution are inferior to functionals with thermal XC and spatial inhomogeneity effects taken into account at the same level of refinement. We believe this and similar evidence should be convincing to the high-energy density physics community of the necessity of use of proper thermal XC functionals in simulation studies of finite-temperature quantum effects in warm, dense plasmas.

Ab-initio molecular dynamics↗
DOE OSTI2025

Impact of structural distortions on the correlated electronic structure of orbital-selective Mott insulating Na 3 ⁢Co 2 ⁢SbO 6 under strain

Na 3 Co 2 SbO 6 is a promising candidate to realize the Kitaev spin liquid phase since the large Kitaev spin exchange interaction is tunable via the change in electronic structure, such as the trigonal crystal field splitting (Δ TCF ). Here, in this study, we show that the uncorrelated electronic structure of Na 3 Co 2 SbO 6 is rather insensitive to the strain effect due to the low crystal symmetry accompanied by oxygen displacements and the presence of Sb s orbitals. Using density functional theory plus dynamical mean field theory, we find that the correlated electronic structure of Na 3 Co 2 SbO 6 is an orbital-selective Mott insulating state where the trigonal a 1g orbital is insulating due to the nearly full occupation, while other d orbitals behave as typical Mott insulators, resulting in the effective tunability of Δ TCF under the strain effect. The sign change of Δ TCF can occur as the in-plane tensile strain is applied, and the Kitaev spin liquid phase could possibly be realized due to the strongly suppressed Δ TCF under tensile strain. Our results show that the local Co-site symmetry and dynamical correlation effects will play an important role in engineering the novel magnetic phase in this and related materials.

Nguyen, Nam [Univ. of Illinois, Chicago, IL (Unite↗
DOE OSTI2020

DFT Plus DMFT Study of Spin-Charge-Lattice Coupling in Covalnt LaCoO{siub 3}.

We study energetics and the nature of both homogeneous and mixed spin (MS) states in LaCoO3 incorporating structural changes of the crystal volume expansion and the Co-O bond disproportionation (BD) during the spin-state transition using the density functional theory plus dynamical mean field theory (DFT + DMFT) method. DFT + DMFT predicts that energetics of both excited spin states are almost the same while DFT + U calculations of the same structures energetically favor the MS states and produce various metastable solutions whose energetics depend sensitively on final spin states. Within DFT + DMFT, the homogeneous spin state in the expanded crystal volume shows the multiconfigurational nature with non-negligible occupancy probabilities of both high spin (HS) and low spin (LS) states along with d(6) and d(7) charge configurations indicating the dynamically fluctuating nature of spin and charge states due to the Co-O covalency. The nature of the MS state under the BD structure reveals that Co sites with the long Co-O bonds develop a Mott insulating state and favor HS with a d(6) configuration, while more covalent Co sites with the short Co-O bonds occupy more LS states with a d(7) configuration and behave as a band insulator, as a result, charge ordering is induced in the BD structure from the spin-state ordering. We also find that both energetics and electronic structure sensitively depend on the Co-O covalency effect, which can be tuned by changing the double counting potential and the resulting d occupancy (N-d), and N-d close to 6.7 is consistent with the nature of the spin-state transition. Our results show that structural changes during the spin-state transition can play an important role in understanding energetics and electronic structure of LaCoO3.

Park, Hyowon↗
DOE OSTI2020

Compositional phase stability of correlated electron materials within DFT+DMFT

Predicting the compositional phase stability of strongly correlated electron materials is an outstanding challenge in condensed matter physics, requiring precise computations of total energies. In this work, we employ the density functional theory plus dynamical mean-field theory (DFT+DMFT) formalism to address local correlations due to transition metal $\textit{d}$ electrons on compositional phase stability in the prototype rechargeable battery cathode material Li x CoO 2 , and detailed comparisons are made with the simpler DFT+$\textit{U}$ approach (i.e., the Hartree-Fock solution of the DMFT impurity problem). Local interactions are found to strongly impact the energetics of the band insulator LiCoO 2 , most significantly via the $E_g$ orbitals, which are partially occupied via hybridization with $\textit{O p}$ states. We find CoO 2 and Li 1/2 CoO 2 to be moderately correlated Fermi liquids with quasiparticle weights of 0.6–0.8 for the $T_{2g}$ states, which are most impacted by the interactions. Further, as compared to DFT+$\textit{U}$, DFT+DMFT considerably dampens the increase in total energy as $\textit{U}$ is increased, which indicates that dynamical correlations are important to describe this class of materials despite the relatively modest quasiparticle weights. Unlike DFT+$\textit{U}$, which can incorrectly drive Li x CoO 2 toward spurious phase separating or charge-ordered states, DFT+DMFT correctly captures the system's phase stability and does not exhibit a strong charge-ordering tendency. Most importantly, the error within DFT+$\textit{U}$ varies strongly as the composition changes, challenging the common practice of artificially tuning $\textit{U}$ within DFT+$\textit{U}$ to compensate the errors of Hartree-Fock. DFT+DMFT predicts the average intercalation voltage decreases relative to DFT, opposite to the result of DFT+$\textit{U}$, which would yield favorable agreement with experiment in conjunction with the overprediction of the voltage by the strongly constrained and appropriately normed DFT functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Orbital selective spin waves in detwinned NaFeAs

The existence of orbital-dependent electronic correlations has been recognized as an essential ingredient to describe the physics of iron-based superconductors. NaFeAs, a parent compound of iron-based superconductors, exhibits a tetragonal-to-orthorhombic lattice distortion below T s ≈ 60 K, forming an electronic nematic phase with two 90° rotated (twinned) domains, and orders antiferromagnetically below T N ≈ 42 K. We use inelastic neutron scattering to study spin waves in uniaxial pressure-detwinned NaFeAs. By comparing the data with combined density functional theory and dynamical mean-field theory calculations, we conclude that spin waves up to an energy scale of E crossover ≈ 100 meV are dominated by d yz -d yz intraorbital scattering processes, which have the twofold (C 2 ) rotational symmetry of the underlying lattice. On the other hand, the spin wave excitations above E crossover , which have approximately fourfold (C 4 ) rotational symmetry, arise from the d xy -d xy intraorbital scattering that controls the overall magnetic bandwidth in this material. In addition, we find that the low-energy (E ≈ 6 meV) spin excitations change from approximate C 4 to C 2 rotational symmetry below a temperature T* (>T s ), while spin excitations at energies above E crossover have approximate C 4 rotational symmetry and are weakly temperature dependent. Here, these results are consistent with angle-resolved photoemission spectroscopy measurements, where the presence of a uniaxial strain necessary to detwin NaFeAs also raises the onset temperature T* of observable orbital-dependent band splitting to above T s , thus supporting the notion of orbital selective spin waves in the nematic phase of iron-based superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hund's metal physics: From SrNiO 2 to LaNiO 2

We study the normal state electronic structure of the recently discovered infinite-layer nickelate superconductor, Nd 1-x Sr x NiO 2 , using density functional theory plus dynamical mean-field theory calculations. Starting with the multiorbital compound SrNiO 2 , our calculations show that despite large charge-carrier doping from SrNiO 2 to LaNiO 2 , the Ni-3d total occupancy is barely changed due to the decreased hybridization with the occupied oxygen-2p states and increased hybridization with the unoccupied La-5d states. Thus, using SrNiO 2 as a reference, La 1-x Sr x NiO 2 is naturally and conclusively found to be a multiorbital electronic system with characteristic Hund's metal behaviors, such as metallicity, the importance of high-spin configurations, tendency towards orbital differentiation, and the absence of magnetism in regimes which are ordered according to static mean-field theories. Finally, our results are in good agreement with the existing spectroscopic studies and make an essential step towards an understanding of the electronic structures of Nd 1-x Sr x NiO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗