Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗
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 OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The lion, the witch, and the wormhole: ensemble averaging the symmetric product orbifold

Abstract We consider the ensemble average of two dimensional symmetric product orbifold CFTs Sym N (𝕋 D ) over the Narain moduli space. We argue for a bulk dual given byNcopies of an abelian Chern-Simons theory coupled to topological gravity, endowed with a discrete gauge symmetry exchanging theNcopies. As a check of this proposal, we calculate the ensemble average of various partition and correlation functions of the symmetric product orbifold theory and compare the resulting expressions to gauge theory quantities in the bulk. We comment on the ensemble average of the tensionless string partition function on AdS 3 × S 3 × 𝕋 4 by considering the specific case ofD= 4 with the addition of supersymmetry.

Physics↗
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↗