Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density response”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electronic density response of warm dense matter

Matter at extreme temperatures and pressures—commonly known as warm dense matter (WDM)—is ubiquitous throughout our Universe and occurs in astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are probed in x-ray Thomson scattering (XRTS) experiments and are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear response theory and nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques, including path-integral Monte Carlo simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrary complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Green’s Function Perspective on the Nonlinear Density Response of Quantum Many-Body Systems

We derive equations of motion for higher order density response functions using the theory of thermodynamic Green’s functions. We also derive expressions for the higher order generalized dielectric functions and polarization functions. Moreover, we relate higher order response functions and higher order collision integrals within the Martin–Schwinger hierarchy. We expect our results to be highly relevant to the study of a variety of quantum many-body systems such as matter under extreme temperatures, densities, and pressures.

Density Functional Theory↗
DOE OSTI2024

Consistency between reflection momentum-resolved electron energy loss spectroscopy and optical spectroscopy measurements of the long-wavelength density response of Bi 2 Sr 2 CaCu 2 O 8 + x

The density fluctuation spectrum captures many fundamental properties of strange metals. Using momentum-resolved electron energy-loss spectroscopy (M-EELS), we recently showed that the density response of the strange metal Bi 2 Sr 2 CaCu 2 O 8+x (Bi-2212) at large momentum, q, exhibits a constant-in-frequency continuum reminiscent of the marginal Fermi liquid (MFL) hypothesis of the late 1980s. However, reconciling this observation with infrared (IR) optics experiments, which show a well-defined plasmon excitation at q ~ 0, has been challenging. Here we report M-EELS measurements of Bi-2212 using 4× improved momentum resolution, allowing us to reach the optical limit. For momenta q < 0.04 r.l.u., the M-EELS data show a plasmon feature that is quantitatively consistent with IR optics. For q > 0.04 r.l.u., the spectra become incoherent with an MFL-like, constant-in-frequency form. Here we speculate that, at finite frequency, ω, and nonzero q, some attribute of this Planckian metal randomizes the probe electron, causing it to lose information about its own momentum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Outer Radiation Belt Flux and Phase Space Density Response to Sheath Regions: Van Allen Probes and GPS Observations

Turbulent and compressed sheath regions preceding interplanetary coronal mass ejections strongly impact electron dynamics in the outer radiation belt. Changes in electron flux can occur on timescales of tens of minutes, which are unlikely to be captured by a two-satellite mission. The recently released Global Positioning System (GPS) data set generally has shorter revisit times (at L ~ 4–8) owing to the large number of satellites in the constellation equipped with energetic particle detectors. Investigating electron fluxes at energies from 140 keV to 4 MeV and sheaths observed in 2012–2018, we show that the flux response to sheaths on a timescale of 6 hr, previously reported from Van Allen Probes (RBSP) data, is reproduced by GPS measurements. Furthermore, GPS data enables derivation of the response on a timescale of 30 min, which further confirms that the energy and L-shell dependent changes in electron flux are associated with the impact of the sheath. Sheath-driven loss is underestimated over longer timescales as the electrons recover during the ejecta. We additionally show the response of electron phase space density (PSD), which is a key quantity in identifying non-adiabatic loss from the system and electron energization through wave-particle interactions. The PSD response is calculated from both RBSP and GPS data for the 6 hr timescale, as well as from GPS data for the 30 min timescale. The response is divided based on the geoeffectiveness of the sheaths revealing that electrons are effectively accelerated only during geoeffective sheaths, while loss commonly occurs during all sheaths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum Monte Carlo calculations of electron scattering from 12 C in the short-time approximation

The short-time approximation is a method introduced to evaluate electroweak nuclear response for systems with A ≥ 12, extending the reach of first-principle many-body quantum Monte Carlo calculations. Using realistic two- and three-body nuclear interactions and consistent one- and two-body electromagnetic currents, we calculate longitudinal and transverse response densities and response functions of 12 C. Here, we compare the resulting cross sections with experimental data for electron-nucleus scattering, finding good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory

While time-dependent density functional theory has emerged as a method of choice for computing electronic spectra and response of molecules and materials, its reliability is hindered by the adiabatic approximation that is commonly made. In this work, we address one problematic aspect that arises from this approximation: an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. Here, we demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

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

Charge-Density-Specific Response of Grafted Polyelectrolytes to Electric Fields: Bending or Tilting?

The response of polyelectrolytes (PEs) to applied electric fields drive applications in energetics, diagnostics, materials development, and many more. Here we employ all-atom molecular dynamics (MD) simulations to probe the response of grafted PE brushes to axial electric fields. For PEs with large charge densities, the electric field triggers a left-right asymmetric distribution of counterions around the PE backbone: consequently, depending on the location (left or right), there is an unequal screening and an unequal force on the PE functional groups causing a bending-driven brush height reduction. However, for the weakly charged PEs, the electric field causes a uniform distribution of the counterions across the brush enforcing a uniform (without left-right asymmetry) partially unscreened PE brush: therefore, all the brush segments experience similar force causing a brush tilting-driven brush height reduction. As a result, such bending versus tilting responses is commensurate with the electric-field-driven increase (decrease) in the flexibility of strongly (weakly) charged PEs.

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

Excited-State Densities from Time-Dependent Density Functional Response Theory

While the variational principle for excited-state energies leads to a route to obtaining excited-state densities from time-dependent density functional theory, relatively little attention has been paid to the quality of the resulting densities in real space obtained with different exchange-correlation functional approximations or how nonadiabatic approximations developed for energies of states of double-excitation character perform for their densities. Here we derive an expression directly in real space for the excited-state density, which includes the case of nonadiabatic kernels and consequently is able, for the first time, to yield densities of states of double-excitation character. Under some well-defined simplifications, we compare the performance of the local-density approximation and exact-exchange approximation, which are in a sense at the opposite extremes of the fundamental functional approximations, on local and charge-transfer excitations in one-dimensional model systems and show that the dressed Time-Dependent Density Functional Theory (TDDFT) approach gives good densities of double excitations.

approximation↗
DOE OSTI2021

Colossal Density-Driven Resistance Response in the Negative Charge Transfer Insulator MnS 2

A reversible density driven insulator to metal to insulator transition in high-spin MnS 2 is experimentally observed, leading with a colossal electrical resistance drop of $10^8$ Omega by 12 GPa. Density functional theory simulations reveal the metallization to be unexpectedly driven by previously unoccupied $S_{2}^{-2}$ $\sigma^{*}_{3p}$ antibonding states crossing the Fermi level. This is a unique variant of the charge transfer insulator to metal transition for negative charge transfer insulators having anions with an unsaturated valence. By 36 GPa the emergence of the low-spin insulating arsenopyrite ($P2_1/c$) is confirmed, and the bulk metallicity is broken with the system returning to an insulative electronic state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mitigating Green&#x27;s function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion-sign problem is crucial, as the signal-to-noise ratio decreases exponentially with Euclidean time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parameterized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. We only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross-sections of larger nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Mitigating Green's function Monte Carlo signal-to-noise problems using contour deformations

The Green's function Monte Carlo (GFMC) method provides accurate solutions to the nuclear many-body problem and predicts properties of light nuclei starting from realistic two- and three-body interactions. Controlling the GFMC fermion sign problem is crucial, as the signal-to-noise ratio decreases exponentially with imaginary time, requiring significant computing resources. Inspired by similar scenarios in lattice quantum field theory and spin systems, in this work, we employ integration contour deformations to improve the GFMC signal-to-noise ratio. Machine learning techniques are used to select optimal contours with minimal variance from parametrized families of deformations. As a proof of principle, we consider the deuteron binding energies and Euclidean density response functions. Here, we only observe mild signal-to-noise improvement for the binding energy case. On the other hand, we achieve an order of magnitude reduction of the variance for Euclidean density response functions, paving the way for computing electron- and neutrino-nucleus cross sections of larger nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

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

Analysis and visualization of energy densities. II. Insights from linear-response time-dependent density functional theory calculations

Inspired by the analysis of Kohn–Sham energy densities by Nakai and coworkers, we extended the energy density analysis to linear-response time-dependent density functional theory (LR-TDDFT) calculations. Using ethylene–tetrafluoroethylene and oxyluciferin–water complexes as examples, distinctive distribution patterns were demonstrated for the excitation energy densities of local excitations (within a molecular fragment) and charge-transfer excitations (between molecular fragments). It also provided a simple way to compute the effective energy of both hot carriers (particle and hole) from charge-transfer excitations via an integration of the excitation energy density over the donor and acceptor grid points.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Spin response and neutrino mean free path in neutron matter

The longitudinal and transverse density- and spin-density response functions in pure neutron matter (PNM) are derived over a wide range of densities within the time dependent local spin density approximation. The underlying density functional was derived from an Auxiliary Field Diffusion Monte Carlo computation of the equation of state of unpolarized and fully spin polarized PNM. In order to assess the dependence of the results on the specific underlying nucleon–nucleon Hamiltonian, we used both the phenomenological Argonne AV8' + UIX force, and local chiral forces up to next-to-next-to-leading-order (N 2 LO). Finally, the resulting response function has then been applied to the study of the neutrino mean free path in PNM.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

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

Collective excitations and quantum incompressibility in electron-hole bilayers

In this study, we apply quantum continuum mechanics to the calculation of the excitation spectrum of a coupled electron-hole bilayer. The theory expresses excitation energies in terms of ground-state intra- and interlayer pair correlation functions, which are available from quantum Monte Carlo calculations. The final formulas for the collective modes deduced from this approach coincide with the formulas obtained in the “quasilocalized particle approximation” by Kalman et al. [G. Kalman, V. Valtchinov, and K. I. Golden, Phys. Rev. Lett. 82 , 3124 (1999) ], and likewise the theory predicts the existence of gapped excitations in the charged channel, with the gap arising from electron-hole correlation. An immediate consequence of the gap is that the static density-density response function of the charged channel vanishes as q 2 for wave vector q → 0 , rather than linearly in q , as commonly expected. In this sense, the system is incompressible . This feature, which has no analog in the classical electron-hole plasma, is consistent with the existence of an excitonic ground state and implies the existence of a discontinuity in the chemical potential of electrons and holes when the numbers of electrons and holes are equal. It should be experimentally observable by monitoring the densities of electrons and holes in response to potentials that attempt to change these densities in opposite directions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗