Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First-principles modeling of plasmons in aluminum under ambient and extreme conditions

The theoretical understanding of plasmon behavior is crucial for an accurate interpretation of inelastic scattering diagnostics in many experiments. In this work, we highlight the utility of linear response time-dependent density functional theory (LR-TDDFT) as a first-principles framework for consistently modeling plasmon properties. We provide a comprehensive analysis of plasmons in aluminum from ambient to warm dense matter conditions and assess typical properties such as the dynamical structure factor, the plasmon dispersion, and the plasmon lifetime. We compare our results with scattering measurements and with other TDDFT results as well as models such as the random phase approximation, the Mermin approach, and the dielectric function obtained using static local field corrections of the uniform electron gas parametrized from path-integral Monte Carlo simulations. We conclude that results for the plasmon dispersion and lifetime are inconsistent between experiment and theories and that the common practice of extracting and studying plasmon dispersion relations is an insufficient procedure to capture the complicated physics contained in the dynamic structure factor in its full breadth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Origins of anisotropic transport in the electrically switchable antiferromagnet Fe 1/3 NbS 2

Recent experiments on the antiferromagnetic intercalated transition metal dichalcogenide Fe 1/3 NbS 2 have demonstrated reversible resistivity switching by application of orthogonal current pulses below its magnetic ordering temperature, making Fe 1/3 NbS 2 promising for spintronics applications. Here, we perform density functional theory calculations with Hubbard U corrections of the magnetic order, electronic structure, and transport properties of crystalline Fe 1/3 NbS 2 , clarifying the origin of the different resistance states. The two experimentally proposed antiferromagnetic ground states, corresponding to in-plane stripe and zigzag ordering, are computed to be nearly degenerate. In-plane cross sections of the calculated Fermi surfaces are anisotropic for both magnetic orderings, with the degree of anisotropy sensitive to the Hubbard U value. The in-plane resistance, computed within the Kubo linear response formalism using a constant relaxation time approximation, is also anisotropic, supporting a hypothesis that the current-induced resistance changes are due to a repopulating of antiferromagnetic domains. Our calculations indicate that the transport anisotropy of Fe 1/3 NbS 2 in the zigzag phase is reduced relative to stripe, consistent with the relative magnitudes of resistivity changes in experiment. Finally, our calculations reveal the likely directionality of the current-domain response, specifically, which domains are energetically stabilized for a given current direction.

36 MATERIALS SCIENCE↗
DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Changes in polarization dictate necessary approximations for modeling electronic deexcitation intensity: Application to x-ray emission

Accurate simulation of electronic excitations and deexcitations are critical for complementing complex spectroscopic experiments and can provide validation to theoretical approaches. Here, using a generalized framework, we contrast the accuracy and validity of orbital-constrained and linear-response approaches that build upon Kohn-Sham density functional theory (DFT) to simulate emission spectra of electronic origin and propose an efficient approximation, named many-body x-ray emission spectroscopy or MBXES, for simulating such processes. We show analytically as well as with computed examples that for electronic (de)excitation leading to an appreciable change in polarization (i.e., density rearrangement), the adiabatic approximation in a response-based formalism will be inadequate for the calculation of oscillator strength. Thus, such a change (e.g., in the net electrostatic dipole moment of a finite system) can be used as a metric for evaluating the applicability of the adiabatic response-based approach and can be particularly valuable in x-ray emission spectroscopy. On the other hand, MBXES, the flexible method introduced in this paper, can compute oscillator strengths accurately at a much lower computational expense on the basis of two DFT-based self-consistent field calculations. Using illustrative examples of emission spectra, the efficacy of the MBXES method is demonstrated by comparison with its parent theory, orbital-optimized DFT, and with experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Real-space Green's function approach for intrinsic losses in x-ray spectra

Intrinsic inelastic losses in x-ray spectra originate from excitations in an interacting electron system due to a suddenly created core-hole. These losses characterize the features observed in x-ray photoemission spectra (XPS), as well as many-body effects such as satellites and edge-singularities in x-ray absorption spectra (XAS). However, they are usually neglected in practical calculations. As shown by Langreth these losses can be treated within linear response in terms of a cumulant Green's function in momentum space. Here we present a complementary ab initio real-space Green's function generalization of the Langreth cumulant in terms of the dynamically screened core-hole interaction W c (ω) and the independent particle response function. Here we find that the cumulant kernel β⁡(ω) is analogous to XAS, but with the transition operator replaced by the core-hole potential with monopole selection rules. The behavior reflects the analytic structure of the loss function, with peaks near the zeros of the dielectric function, consistent with delocalized quasiboson excitations. The approach simplifies when W c (ω) is localized and spherically symmetric. In conclusion, illustrative results and comparisons are presented for the electron gas, sodium, and some early transition metal compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Extracting spinon self-energies from two-dimensional coherent spectroscopy

Two-dimensional coherent spectroscopy (2DCS) is a nonlinear spectroscopy technique capable of identifying whether apparent continua in linear response are made out of multiplets of sharp deconfined quasiparticles. This makes it a potent tool for experimental identification of fractionalized phases. Previous discussions have focused on limits where the quasiparticles in question are infinitely long lived. In this paper we discuss 2DCS in the regime where the fractionalized quasiparticles can themselves decay. We introduce a powerful path-integral-based approach, whereby the computation of nonlinear susceptibilities reduces to an efficient exercise in diagrammatic perturbation theory. Here, we apply this method to compute the 2DCS response of the one-dimensional transverse field Ising model, in the presence of integrability-breaking perturbations. We discuss aspects of the self-energy of the fractionalized quasiparticles that may be extracted via 2DCS, such as the momentum-dependent decay rate.

1-dimensional spin chains↗
DOE OSTI2023

Hydrodynamic relaxation of spin helices

Motivated by recent cold-atom experiments, here we study the relaxation of spin helices in quantum XXZ spin chains. The experimentally observed relaxation of spin helices follows scaling laws that are qualitatively different from linear-response transport. We construct a theory of the relaxation of such spin helices, combining generalized hydrodynamics with diffusive corrections and a generalized form of the local density approximation. Although helices are far from local equilibrium, our hydrodynamic approach reproduces the experimentally observed relaxational dynamics and also predicts the late-time relaxation, which is outside the experimentally accessible timescales. In particular, our theory explains the existence of temporal regimes with apparent anomalous diffusion, as well as the asymmetry between positive and negative anisotropy regimes at short and intermediate times.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Aluminum vacancy/sulfur complex in wurtzite AlN as an optically controllable spin qubit

Using our rational methodology, we reveal a defect in wurtzite AlN that can serve as an optically controllable spin qubit. It combines an Al vacancy and a S atom substituting the neighboring N atom (V Al⁢ S N ). Linear response GW and Bethe-Salpeter equation calculations guide us to find suitable ground and excited triplet and singlet states of V Al⁢ S N . The obtained optical spin-polarization cycle is similar to that observed in the negative nitrogen-vacancy (NV – ) center in diamond. Furthermore, the calculated optical oscillator strengths for V Al⁢ S N suggest that, in contrast to the NV – center, the optical emission in the singlet and triplet states have comparable rates, which is favorable for the qubit functionality.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Spin dynamics in the itinerant antiferromagnet SrCr 2 ⁢As 2

SrCr 2 ⁢As 2 is an itinerant antiferromagnet in the same structural family as the SrFe 2 ⁢As 2 high-temperature superconductors. Here, we report our calculations of exchange-coupling parameters 𝐽 𝑖⁢𝑗 for SrCr 2 ⁢As 2 using a static linear-response method based on first-principles electronic-structure calculations. We find that the dominant nearest-neighbor exchange coupling 𝐽 1 >0 is antiferromagnetic whereas the next-nearest-neighbor exchange coupling 𝐽 2 <0 is ferromagnetic with 𝐽 2 /𝐽 1 = −0.68, reinforcing the checkerboard in-plane magnetic structure. Thus, unlike other transition-metal arsenides based on Mn, Fe, or Co, we find no competing magnetic interactions in SrCr 2 ⁢As 2 , which aligns with experimental findings. Moreover, the orbital resolution of exchange interactions shows that 𝐽 1 and 𝐽 2 are dominated by direct exchange mediated by the Cr 𝑑 orbitals. To validate the calculations we conduct inelastic neutron-scattering measurements on powder samples that show steeply dispersive magnetic excitations arising from the magnetic Γ points and persisting up to energies of at least 175 meV. The spin-wave spectra are then modeled using the Heisenberg Hamiltonian with the theoretically calculated exchange couplings. In conclusion, the calculated neutron-scattering spectra are in good agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multi-messenger heavy-ion collision physics

We report that this work studies the production of direct photons in relativistic nuclear collisions, along with the production of hadrons. Radiation from the very first instants to the final moments of the evolution is included. The hybrid model used here describes all stages of relativistic heavy-ion collisions. Chronologically, those are an initial state reflecting the collision of nuclei described within the Color Glass Condensate effective theory; a pre-equilibrium phase based on non-equilibrium linear response; relativistic viscous hydrodynamics, and a hadronic afterburner. The effect of the preequilibrium phase on both photonic and hadronic observables is highlighted for the first time. The potential of photon observables - spectrum, differential elliptic and triangular flow - to reveal the chemical equilibration time is studied. Finally, we consider "small collision systems", including proton+nucleus collisions and collisions of light nuclei, as probed by hadronic and electromagnetic observables. We demonstrate how photon production can signal the formation of quark-gluon plasma in such small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Two-body weak currents in heavy nuclei

In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. Here we examine the systematic effects of two-body axial currents on Gamow-Teller strength and β-decay rates in heavy nuclei within energy-density functional theory. Using a Skyrme functional and the charge-changing finite amplitude method, we add the contributions of two-body currents to the usual one-body linear response in the Gamow-Teller channel, both exactly and though a densitymatrix expansion. The two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows. They can enhance individual low-lying transitions, however, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. We show that both these effects are related to changes in the total nucleon density as the number of neutrons increases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Critical and near-critical relaxation of holographic superfluids

We investigate the relaxation of holographic superfluids after quenches, when the end state is either tuned to be exactly at the critical point, or very close to it. By solving the bulk equations of motion numerically, we demonstrate that in the former case the system exhibits a power law falloff, as well as an emergent discrete scale invariance. The latter case is in the regime dominated by critical slowing down, and we show that there is an intermediate time range before the onset of late-time exponential falloff, where the system behaves similarly to the critical point with its power law falloff. We further postulate a phenomenological Gross-Pitaevskii-like equation (corresponding to model F of Hohenberg and Halperin) that is able to make quantitative predictions for the behavior of the holographic superfluid after near-critical quenches into the superfluid and normal phase. Intriguingly, all parameters of our phenomenological equation, which describes the nonlinear time evolution, may be fixed with information from the static equilibrium solutions and linear response theory. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Search for very-short-baseline oscillations of reactor antineutrinos with the SoLid detector

In this paper we report the first scientific result based on antineutrinos emitted from the BR2 reactor at SCK CEN. The SoLid experiment uses a novel type of highly granular detector whose basic detection unit combines two scintillators, polyvinyl toluene (PVT) and Li 6 F : ZnS ( Ag ) , to measure antineutrinos via their inverse- β -decay products. An advantage of PVT is its highly linear response as a function of deposited particle energy. The full-scale detector comprises 12 800 voxels and operates over a very short 6.3–8.9 m baseline from the reactor core. The detector segmentation and its three-dimensional imaging capabilities facilitate the extraction of the positron energy from the rest of the visible energy, allowing the latter to be utilized for signal-background discrimination. We present a result obtained from 280 reactor-on days (55 MW mean power) and 172 reactor-off days, respectively, of live data taking. A total of 29 479 ± 603 (stat) antineutrino candidates have been selected, corresponding to an average rate of 105 events per day and a signal-to-background ratio of 0.27. A search for disappearance of antineutrinos to a sterile state has been conducted using complementary model-dependent frequentist and Bayesian fits, providing constraints on the allowed region of the reactor antineutrino anomaly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Spectral density reconstruction with Chebyshev polynomials

Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semiexclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. Here, in the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to obtain the local density of states for graphene in a magnetic field as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Structure Factors of Neutron Matter at Finite Temperature

Here, we compute continuum and infinite volume limit extrapolations of the structure factors of neutron matter at finite temperature and density. Using a lattice formulation of leading-order pionless effective field theory, we compute the momentum dependence of the structure factors at finite temperature and at densities beyond the reach of the virial expansion. The Tan contact parameter is computed and the result agrees with the high momentum tail of the vector structure factor. All errors, statistical and systematic, are controlled for. This calculation is a first step towards a model-independent understanding of the linear response of neutron matter at finite temperature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Casimir Light in Dispersive Nanophotonics

Time-varying optical media, whose dielectric properties are actively modulated in time, introduce a host of novel effects in the classical propagation of light, and are of intense current interest. In the quantum domain, time-dependent media can be used to convert vacuum fluctuations (virtual photons) into pairs of real photons. Here, we refer to these processes broadly as “dynamical vacuum effects” (DVEs). Despite interest for their potential applications as sources of quantum light, DVEs are generally very weak, presenting many opportunities for enhancement through modern techniques in nanophotonics, such as using media which support excitations such as plasmon and phonon polaritons. Here, we present a theory of weakly modulated DVEs in arbitrary nanostructured, dispersive, and dissipative systems. A key element of our framework is the simultaneous incorporation of time-modulation and “dispersion” through time-translation-breaking linear response theory. As an example, we use our approach to propose a highly efficient scheme for generating entangled surface polaritons based on time-modulation of the optical phonon frequency of a polar insulator.

42 ENGINEERING↗