Search NASA⌕ Search

SEARCH · Search NASA

Results for “Exchange correlation functionals”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Free-energy orbital-free density functional theory: recent developments, perspective, and outlook

By summarizing the constraint-based development of orbital-free free-energy density functional approximations, we provide a perspective on progress over the last 15 years, the limitations of existing functionals, and the challenges awaiting resolution. We outline the chronology of the development of noninteracting and exchange-correlation free-energy orbital-free functionals and summarize the theoretical basis of existing local density approximation, second-order approximation, generalized gradient approximation (GGA), and meta-GGAs. We discuss limitations and challenges such as problems with thermodynamic derivatives, free-energy nonadditivity and the closely related issue of all-electron versus valence-only local pseudo-potential performance.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Scalable learning of potentials to predict time-dependent Hartree–Fock dynamics

We propose a framework to learn the time-dependent Hartree–Fock (TDHF) inter-electronic potential of a molecule from its electron density dynamics. Although the entire TDHF Hamiltonian, including the inter-electronic potential, can be computed from first principles, we use this problem as a testbed to develop strategies that can be applied to learn a priori unknown terms that arise in other methods/approaches to quantum dynamics, e.g., emerging problems such as learning exchange–correlation potentials for time-dependent density functional theory. We develop, train, and test three models of the TDHF inter-electronic potential, each parameterized by a four-index tensor of size up to 60 × 60 × 60 × 60. Two of the models preserve Hermitian symmetry, while one model preserves an eight-fold permutation symmetry that implies Hermitian symmetry. Across seven different molecular systems, we find that accounting for the deeper eight-fold symmetry leads to the best-performing model across three metrics: training efficiency, test set predictive power, and direct comparison of true and learned inter-electronic potentials. All three models, when trained on ensembles of field-free trajectories, generate accurate electron dynamics predictions even in a field-on regime that lies outside the training set. To enable our models to scale to large molecular systems, we derive expressions for Jacobian-vector products that enable iterative, matrix-free training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗
DOE OSTI2024

Efficient simulations of charge density waves in the transition metal Dichalcogenide TiSe 2

Charge density waves (CDWs) in transition metal dichalcogenides are the subject of growing scientific interest due to their rich interplay with exotic phases of matter and their potential technological applications. Here, using density functional theory with advanced meta-generalized gradient approximations (meta-GGAs) and linear response time-dependent density functional theory (TDDFT) with state-of-the-art exchange-correlation kernels, we investigate the electronic, vibrational, and optical properties in 1T-TiSe 2 with and without CDW. In both bulk and monolayer TiSe 2 , the electronic bands and phonon dispersions in either normal or CDW (semiconducting) phase are described well via meta-GGAs, which separate the valence and conduction bands just as HSE06 does but with significantly more computational feasibility. The experimentally observed humps of electron energy loss spectroscopy are successfully reproduced in TDDFT. Our work opens the door to simulating these complexities in CDW compounds from first principles by revealing meta-GGAs as an accurate low-cost alternative to HSE06.

Electronic properties and materials↗
DOE OSTI2025

Atomistic characterization of hydration-dependent fuel cell ionomer nanostructure: validation by vibrational spectroscopy

The development of Nafion alternatives for fuel cells and electrolyzers requires a fundamental understanding of hydration-dependent ion-exchange site acid/base chemistry. We present here reactive force field (ReaxFF) molecular dynamics (MD) simulations of Nafion at varying molar water/ion-exchange-site ratios (λ), which we correlate to our experimental and density functional theory-based vibrational spectra. ReaxFF describes the formation and breaking of covalent bonds, enabling simulations of proton exchange between sulfonic acid/sulfonate groups and water/hydronium. Our MD simulations determine the λ-dependent equilibrium proportions of protonated and deprotonated sites. We find that protonated sites persist across all λ from 0 to 20, challenging the widely accepted notion that all sites are ionized above a threshold λ value (e.g., 3 or 4). Our simulations generate hundreds of realistic exchange site environments, the characterization of which are based on 6 Å radii sulfur-centered ‘inner-spheres’. These inner- and outer-sphere regions elucidate an interplay of stereoelectronic factors that influence protonation states, including the number of inner-sphere waters (Λ). Our simulations produce broad λ-dependent distributions of Λ values, representing non-uniform exchange site hydration. Moreover, we demonstrate that these distributions are specific to both protonated and deprotonated sites, with significant overlap between the two distributions for all non-zero λ. These distributions underpin the IR spectra of hydrated membranes, with each exchange site contributing an IR spectrum characteristic of its protonation state and Λ. We expect that these nanostructural characterizations of Nafion exchange sites will contribute to the development of new ionomers.

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

Out-of-time-order correlators bridge classical transport and quantum dynamics

The out-of-time-order correlator (OTOC) has emerged as a central tool for quantifying decoherence across wide-ranging physical platforms. Here, we demonstrate its direct measurement in a classical ensemble using nuclear magnetic resonance with a modulated gradient spin echo sequence and extend the method into a multidimensional correlation to track exchange phenomena. Position is encoded through magnetic field gradients and momentum through the velocity autocorrelation function, enabling experimental access to OTOCs for proton motion confined within the self-similar lattice of the metal–organic framework MOF-808. Here, water confined to specified geometries within the MOF pores gives rise to spatially distinct diffusive eigenmodes with characteristic relative entropies. We demonstrate that periodic radio frequency driving combined with gradient modulation yields entropy evolution through the selection of distinct diffusion modes. Frequency-resolved diffusion spectra connect these entropy dynamics to classical heat exchange laws, revealing how operational features of quantum systems are mirrored in confined, macroscopic spin ensembles.

Fricke, Sophia N. [University of California, Berke↗
DOE OSTI2025

Continuum of magnetic excitations in the Kitaev honeycomb iridate D 3 LiIr 2 O 6

Inelastic neutron scattering (INS) measurements of powder D 3 ( 7 Li)( 193 Ir) 2 O 6 reveal low energy magnetic excitations with a scattering cross-section that is broad in ∣Q∣ and energy transfer. The magnetic nature of the excitation spectrum is demonstrated by longitudinally polarized neutron scattering. The total magnetic moment of 1.8(4)μ B /Ir inferred from the observed magnetic scattering cross-section is consistent with the effective moment inferred from magnetic susceptibility data and expectations for the J eff = 1/2 single ion state. The rise in the dynamic correlation function $\mathcal{S}(Q,w)$ for ℏω < 5 meV can be described by a simple model assuming nearest-neighbor anisotropic spin exchange, such as that found in the Kitaev model. Exchange disorder associated with the D site likely plays an important role in stabilizing the low T quantum fluctuating state.

Halloran, Thomas [Johns Hopkins Univ., Baltimore, ↗
DOE OSTI2026

Collapsing of spin-split Fermi surfaces in the magnetic Weyl semimetal candidate EuAuSb

EuAuSb is a magnetic Weyl semimetal candidate that exhibits in-plane ferromagnetic ordering and a helical magnetic order across the plane at temperatures below 3.9 K. These magnetic characteristics can drive an interaction-induced reconstruction of the Fermi surface. Here, we report a rare instance of Fermi surface splitting and collapsing driven by exchange coupling, as revealed by strongly temperature-dependent quantum oscillation frequencies in EuAuSb. Shubnikov–de Haas measurements reveal two distinct quantum oscillation frequencies (𝐹 𝛼 = 80 T and 𝐹 𝛽 = 165 T) at 2 K. As the temperature increases toward 40 K, the oscillation peaks at these two frequencies gradually merge into a single peak centered around 128 T. This behavior signifies the collapse of the spin-split Fermi surfaces, driven by the weakening of short-range magnetic correlations and the exchange coupling between conduction electrons and local magnetic moments as temperature increases. Furthermore, our results are further supported by density functional theory calculations, angle-resolved photoemission spectroscopy, and magnetization measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Gluon Sivers function from forward exclusive χ c 1 photoproduction on unpolarized protons

Exclusive production of a χ c 1 axial vector quarkonia in photon-proton scattering at high energies requires a C -odd t -channel exchange. In the limit of vanishing momentum transfer this occurs either via the exchange of a photon, the Primakoff process, where the spin of the proton does not change. For axial-vector meson production, as a consequence of the Landau-Yang theorem, the Primakoff cross section is finite as t → 0 . Alternatively, a C -odd spin dependent Odderon can be exchanged, which involves a spin flip of the proton. The resulting cross section is related to the square of the collinear trigluon correlator or the k ⊥ -moment of the gluon Sivers function. Using two models for the gluon Sivers function from the literature we compute the ratio of Sivers to Primakoff cross sections and the angular coefficient λ θ governing the angular distribution of the χ c 1 → J / ψ + γ decay as functions of x . We point out that these observables constrain the magnitude of the gluon Sivers function at small x which could be accessed in electron-proton scattering and ultraperipheral proton-proton and nucleus-proton collisions. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Local spin structure in the layered van der Waals materials MnPS 𝑥 ⁢Se 3−𝑥

Two-dimensional (2D) layered materials, whether in bulk form or reduced to just a single layer, have potential applications in spintronics and capacity for advanced quantum phenomena. A prerequisite for harnessing these opportunities lies in gaining a comprehensive understanding of the spin behavior in 2D materials. The low dimensionality motivates an understanding of the spin correlations over a wide length scale, from local to long-range order. In this context, we focus on the magnetism in bulk MnPSe 3 and MnPS 3 , 2D layered van der Waals antiferromagnetic semiconductors. These materials have similar honeycomb Mn layers and magnetic ordering temperatures, but distinct spin orientations and exchange interactions. Here, we utilize neutron scattering to gain deeper insights into the local magnetic structures and spin correlations in the paramagnetic and ordered phases by systematically investigating a MnPS 𝑥 ⁢Se 3−𝑥 (𝑥 = 0, 1, 1.5, 2, 3) series of powder samples using total neutron scattering measurements. By employing magnetic pair distribution function analysis, we unraveled the short-range magnetic correlations in these materials and explored how the nonmagnetic anion S/Se mixing impacts the magnetic correlations. The results reveal that the magnetism can be gradually tuned through alteration of the nonmagnetic S/Se content, which tunes the atomic structure. The change in magnetic structure is also accompanied by a control of the magnetic correlation length within the 2D honeycomb layers. Complimentary inelastic neutron scattering measurements allowed a quantification of the change in the magnetic exchange interactions for the series and further highlighted the gradual evolution of spin interactions in the series MnPS 𝑥 ⁢Se 3−𝑥 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Chemical Design of Spin Frustration to Realize Topological Spin Glasses

Patterning spins to generate collective behavior is at the core of condensed matter physics. Physicists develop techniques, including the fabrication of magnetic nanostructures and precision layering of materials specifically to engender frustrated lattices. As chemists, we can access such exotic materials through targeted chemical synthesis and create new lattice types by chemical design. Here, in this study, we introduce a new approach to induce magnetic frustration on a modified honeycomb lattice through a competition of alternating antiferromagnetic (AFM) and ferromagnetic (FM) nearest-neighbor interactions. By subtly modulating these two types of interactions through facile synthetic modifications, we created two systems: (1) a topological spin glass and (2) a frustrated spin-canted magnet with low-temperature exchange bias. To design this unconventional magnetic lattice, we used a metal–organic framework (MOF) platform, Ni 3 (pymca) 3 X 3 (NipymcaX where pymca = pyrimidine-2-carboxylato and X = Cl, Br). We isolated two MOFs, NipymcaCl and NipymcaBr, featuring canted Ni 2+ -based moments. Despite this similarity, differences in the single-ion anisotropies of the Ni 2+ spins result in distinct magnetic properties for each material. NipymcaCl is a topological spin glass, while NipymcaBr is a rare frustrated magnet with low-temperature exchange bias. Density functional theory calculations and Monte Carlo simulations on the NipymcaX lattice support the presence of magnetic frustration as a result of alternating AFM and FM interactions. Our calculations enabled us to determine the ground-state spin configuration and the distribution of spin–spin correlations relative to paradigmatic kagomé and triangular lattices. This modified honeycomb lattice is similar to the electronic Kekulé-O phase in graphene and provides a highly tunable platform to realize unconventional spin physics.

Exchange bias↗
DOE OSTI2025

Facet-Dependent Doping and Dopant-Dependent Faceting in Si-Doped GaAsSb Nanowires

Here, we correlate the spatial distributions of Si, Sb, and rotational twins in Si-doped GaAs 1-x Sb x nanowires. GaAs 1-x Sb x nanowires were grown epitaxially on Si(111) substrates by tuning process conditions to achieve repeated nucleation of rotational twins and growth along the [111]B direction; dilute Sb and Si fluxes were chosen to create a sufficient twin density to achieve high yield while avoiding growth of the wurtzite phase. While the impact of Si and Sb on twin density and nanowire growth rate has been previously reported, the facet-dependent incorporation of these species has not been established. Scanning transmission electron microscopy was used to confirm that Sb incorporates preferentially on the (111)B facets relative to {1 ̅1 ̅0} facets prior to nucleation of a rotational twin. With periodic twinning, this facet dependence leads to alternating regions of enriched and depleted Sb concentration attributed to a growth rate-dependent Sb-As-exchange mechanism. Atom probe tomography measurements establish that while Si doping is not perturbed by twinning on (111)B facets, Si and Sb concentrations are anti-correlated for growth on non-(111)B facets. Density functional theory calculations underpin a thermodynamic model that explains the observed anisotropies in dopant incorporation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring interlayer coupling in the twisted bilayer PtTe 2

We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

Performance of Diffusion Monte Carlo Calculations for Predicting the Relative Energies of Quinoidal and Nonquinoidal Species

Coupled cluster singles and doubles with perturbative triples [CCSD(T)] and single determinant fixed-node diffusion Monte Carlo (SD-DMC) have emerged as two of the most useful methods for providing benchmark reaction and interaction energies of chemical systems without strong static correlation. The errors in DMC energies are dominated by an inexact description of the nodal surfaces for electron exchange. One of the main approaches to addressing the fixed-node error is to use multideterminant (MD) trial wave functions. We consider here the energy differences between pairs of related molecules with aromatic and quinoidal structures as well as between quinoidal isomers. Quinoidal systems tend to have some diradical character, leading one to anticipate that SD-DMC calculations may face challenges in accurately describing their energetics. The MD trial wave functions were generated from the complete active space calculations. A comparison is made with the predictions of well-converged CCSD(T) calculations.

basis sets↗
DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (ACS Appl. Mater. Interfaces 2026, 18, 15654-15664). The calculations model the interfacial energetics of two MnTe polymorphs — NiAs-MnTe (hexagonal, alpha phase) and ZnS-MnTe (cubic, gamma phase) — on InP(111) substrates with two surface terminations: In-terminated InP(111)A and P-terminated InP(111)B. This gives four interface configurations: NiAs on In-terminated (experimentally observed), NiAs on P-terminated (computed for comparison), ZnS on In-terminated (computed for comparison), and ZnS on P-terminated (experimentally observed). The dataset is organized into four calculation types, each covering all four polymorph/termination combinations: (i) Slabs: Pristine MnTe/InP heterostructure slabs used to compute total energies and interface energy densities (Eint) for all four configurations, as reported in Fig. 6 of the main text. (ii) Disorder: Same slab geometries with a P_Te + Te_P antisite defect pair introduced near the interface, used to assess chemical intermixing effects on interface stability (Fig. S8, SI). (iii) Strain: Pristine slab calculations with in-plane lattice parameters strained by -1% and +1% relative to the InP lattice constant, used to evaluate strain-dependent interface energetics (Fig. S9, SI). (iv) Charge_Density: Single-point calculations on the full heterostructure, the isolated InP slab, and the isolated MnTe slab at fixed geometry, used to compute differential charge density plots showing interfacial charge accumulation and depletion as a function of surface termination (Fig. S10, SI). Each calculation folder contains INCAR, KPOINTS, POSCAR, CONTCAR, OUTCAR, and POTCAR_info.txt (PAW potential information, excluding the full POTCAR due to VASP licensing restrictions). The calculations were performed using VASP 6.4.3 with PBE exchange-correlation, PAW potentials, a Hubbard correction of Ueff = 5 eV on Mn d-states, and A-type AFM spin initialization.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Pressure‐Induced Structural and Magnetic Evolution in Layered Antiferromagnet YbMn 2 Sb 2

Electronic states under pressure exhibit unconventional spin and charge dynamics that provide a powerful route to uncover exotic phases in quantum materials. Here, we present the structural, magnetic, and electronic evolution of YbMn 2 Sb 2 under pressure. Single-crystal X-ray diffraction reveals a pressure-induced structural transition from the space group trigonal P$\overline{3}$m1 to the monoclinic P2 1 /m phase near 3.5 GPa, which remains stable up to 10 GPa. Magnetization measurements display an anomalously weak net magnetic moment and the absence of Curie–Weiss behavior up to 400 K, suggesting the formation of short-range Mn moment pairs that cancel macroscopically and subsequently evolve into long-range order upon cooling. Temperature-dependent resistivity shows semiconducting behavior with a transition at ∼119 K at ambient pressure, while pressure induces a dramatic suppression of resistance and the emergence of metallic-like temperature dependence, stabilized beyond 5 GPa. This pressure-driven semiconductor-metal transition is consistent with our density functional theory calculations, confirming the closing of the band gap under compression. Neutron diffraction under pressure identifies an incommensurate magnetic structure with antiparallel correlations between paired spins. Together, these results demonstrate how pressure-driven structural tuning and competing exchange interactions stabilize unconventional magnetic states in this low-dimensional magnetic semiconductor.

antiferromagnetic semiconductor↗
DOE OSTI2024

Excitons in 2D Magnets

Theoretical abilities that we built: In our QSGW formalism, the electronic eigenfunctions are calculated self-consistently in presence of the ladder e-h vertex corrections to the screened Coulomb exchange. Charge and self-energies are recalculated iterated until the desired tolerance is achieved in the one-particle Green's function. Further, to include the multi-determinantal (spin-flip structure of multiplets) nature of many body correlations we combine QSGW with DMFT. Below we explore various classes of strongly correlated magnetic systems where excitons are sufficiently described within either the QSGW or the QSGW +DMFT approach. In the process, we establish the digrammatic requirements for a minimum-sufficient theory for describing excitons in large classes of materials.

2D magnets↗
DOE OSTI2025

Effect of Ionic Core on the Properties of Non-Isothermal Plasmas

This study investigates the effects of ionic cores on non-isothermal plasmas using a novel ion-ion interaction potential that incorporates screening effects from both ion cores and exchange-correlation interactions. Our findings indicate that with increasing distance, the effective potential approaches a Yukawa-like screening potential, while at shorter distances, strong electron binding weakens screening. The different values of the cutoff radius ' cut and core edge steepness α significantly influence the potential behavior and radial distribution functions (RDFs). Higher coupling parameters (Γ i ) strengthen the electron-ion interactions, leading to deeper potential wells and more pronounced non-ideality corrections. Increasing Θ decreases the absolute values of non-ideality corrections, indicating fewer interactions in the system. A larger cutoff radius r cut due to weaker screening effects. As Γ i at a fixed parameter also reduces corrections increases, non-ideality corrections grow, reflecting stronger coupling. The results show the importance of taking into account the ion core effects in dense non-isothermal plasma research.

dense plasmas↗