Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kernel”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗
DOE OSTI2021

Predicting impurity spectral functions using machine learning

The Anderson Impurity Model (AIM) is a canonical model of quantum many-body physics. Here we investigate whether machine learning models, both neural networks (NN) and kernel ridge regression (KRR), can accurately predict the AIM spectral function in all of its regimes, from empty orbital, to mixed valence, to Kondo. To tackle this question, we construct two large spectral databases containing approximately 410 000 and 600 000 spectral functions of the single-channel impurity problem. We show that the NN models can accurately predict the AIM spectral function in all of its regimes, with pointwise mean absolute errors down to 0.003 in normalized units. We find that the trained NN models outperform models based on KRR and enjoy a speedup on the order of 10 5 over traditional AIM solvers. Finally, the required size of the training set of our model can be significantly reduced using farthest point sampling in the AIM parameter space, which is important for generalizing our method to more complicated multichannel impurity problems of relevance to predicting the properties of real materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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

Dark-exciton energy splitting in monolayer W Se 2 : Insights from time-dependent density functional theory

Here, we present here a formalism based on time-dependent density functional theory to describe characteristics of both intra- and intervalley excitons in semiconductors, the latter of which had remained a challenge. Through the usage of an appropriate exchange-correlation kernel (nanoquanta), we trace the energy difference between the intra- and intervalley dark excitons in monolayer (1L) W⁢ Se 2 to the domination of the exchange part in the exchange-correlation energies of these states. Furthermore, our calculated transition contribution maps establish the momentum resolved weights of the electron-hole excitations in both bright and dark excitons thereby providing a comprehensive understanding of excitonic properties of 1L W⁢ Se 2 . We find that the states consist of hybridized excitations around the corresponding valleys which leads to brightening of the dark excitons, i.e., significantly decreasing their lifetime which is reflected in the PL spectrum. Using many-body perturbation theory, we calculate the phonon contribution to the energy band gap and the linewidths of the excited electrons, holes, and (bright) exciton to find that as the temperature increases the band gap significantly decreases, while the linewidths increase. Our work paves for describing the ultrafast charge dynamics of transition metal dichalcogenide within an ab initio framework.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2024

Absence of quantization in the circular photogalvanic effect in disordered chiral Weyl semimetals

The circularly polarized photogalvanic effect (CPGE) is studied in chiral Weyl semimetals with short-range quenched disorder. Without disorder, the topological properties of chiral Weyl semimetals lead to quantization of the CPGE, which is a second-order optical response. Furthermore, using a combination of diagrammatic perturbation theory in the continuum and exact numerical calculations via the kernel polynomial method on a lattice model, we show that disorder perturbatively destabilizes the quantization of the CPGE.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Euclidean formulation of relativistic quantum mechanics of N particles

A Euclidean formulation of relativistic quantum mechanics for systems of a finite number of degrees of freedom is discussed. Relativistic treatments of quantum theory are needed to study hadronic systems at subhadronic distance scales. While direct interaction approaches to relativistic quantum mechanics have proved to be useful, they have two disadvantages. One is that cluster properties are difficult to realize for systems of more than two particles. The second is that the relation to quantum field theories is indirect. Euclidean formulations of relativistic quantum mechanics provide an alternative representation that does not have these difficulties. More surprising, the theory can be formulated entirely in the Euclidean representation without the need for analytic continuation. In this work a Euclidean representation of a relativistic N-particle system is discussed. Kernels for systems of N free particles of any spin are given and shown to be reflection positive. Explicit formulas for generators of the Poincaré group for any spin are constructed and shown to be self-adjoint on the Euclidean representation of the Hilbert space. The structure of correlations that preserve both the Euclidean covariance and reflection positivity is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Non-Gaussian estimates of tensions in cosmological parameters

We discuss how to efficiently and reliably estimate the level of agreement and disagreement on parameter determinations from different experiments, fully taking into account non-Gaussianities in the parameter posteriors. We develop two families of scalable algorithms that allow us to perform this type of calculations in increasing number of dimensions and for different levels of tensions. One family of algorithms rely on kernel density estimates of posterior distributions while the other relies on machine learning modeling of the posterior distribution with normalizing flows. We showcase their effectiveness and accuracy with a set of benchmark examples and find both methods agree with each other and the true tension within 0.5σ in difficult cases and generally to 0.2σ or better. This allows us to study the level of internal agreement between different measurements of the clustering of cosmological structures from the Dark Energy Survey and their agreement with measurements of the Cosmic Microwave Background from the Planck satellite.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Reionization effective likelihood from Planck 2018 data

We release relike (reionization effective likelihood), a fast and accurate effective likelihood code based on the latest Planck 2018 data that allows one to constrain any model for reionization between 6 < z < 30 using five constraints from the CMB reionization principal components (PC). We tested the code on two example models which showed excellent agreement with sampling the exact Planck likelihoods using either a simple Gaussian PC likelihood or its full kernel density estimate. This code enables a fast and consistent means for combining Planck constraints with other reionization data sets, such as kinetic Sunyaev-Zeldovich effects, line-intensity mapping, luminosity function, star formation history, quasar spectra, etc., where the redshift dependence of the ionization history is important. Since the PC technique tests any reionization history in the given range, we also derive model-independent constraints for the total Thomson optical depth τ PC = $0.0619$$^{+0.0056}_{–0.0068}$ and its 15 ≤ z ≤ 30 high redshift component τ PC (15,30) < 0.020 (95% C.L.). Furthermore, the upper limits on the high-redshift optical depth is a factor of ~3 larger than those reported in the Planck 2018 cosmological parameter paper using the FlexKnot method and we validate our results with a direct analysis of a two-step model which permits this small high-z component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Energy loss of a heavy fermion in a collisional QED plasma

We compute the energy loss of heavy fermions moving in a plasma, taking into account the modification of the photon collective modes induced by collisions using a Bhatnagar-Gross–Krook collisional kernel. We include contributions from both hard and soft scatterings of the heavy fermion using a collisionally modified hard-thermal-loop resummed propagator. Using this method, one does not need to introduce a separation scale between hard- and soft-momentum exchanges. To place our calculation in context, we review other theoretical approaches to computing the collisional energy loss of fermions and discuss the systematics and results obtained in each approach compared to using a resummed propagator for both hard and soft momentum exchanges. Our final results indicate that self-consistently including the effect of collisions in the self-energies of the resummed propagator results in an increased energy loss compared to using collisionless hard-thermal-loop propagators. The effect becomes larger as the magnitude of the coupling constant and the velocity of the fermion increase. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Unpolarized transverse momentum dependent parton distributions of the nucleon from lattice QCD

We present a first lattice QCD calculation of the unpolarized nucleon’s isovector transverse-momentum-dependent parton distribution functions (TMDPDFs), which are essential to predict observables of multiscale, semi-inclusive processes in the standard model. We use a N f = 2 + 1 + 1 MILC ensemble with valence clover fermions on a highly improved staggered quark (HISQ) sea to compute the quark momentum distributions in a large-momentum nucleon on the lattice. The state-of-the-art techniques in renormalization and extrapolation in the correlation distance on the lattice are adopted. The perturbative kernel up to next-to-next-to-leading order is taken into account, and the dependence on the pion mass and the hadron momentum is explored. Our results are qualitatively comparable with phenomenological TMDPDFs, which provide an opportunity to predict high energy scatterings from first principles. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Revisiting shear stress tensor evolution: Nonresistive magnetohydrodynamics with momentum-dependent relaxation time

This study aims to develop second-order relativistic viscous magnetohydrodynamics (MHD) derived from kinetic theory within an extended relaxation time approximation (momentum/energy dependent) for the collision kernel. The investigation involves a detailed examination of shear stress tensor evolution equations and associated transport coefficients. The Boltzmann equation is solved using a Chapman-Enskog-like gradient expansion for a charge-conserved conformal system, incorporating a momentum-dependent relaxation time. The derived relativistic nonresistive, viscous second-order MHD equations for the shear stress tensor reveal significant modifications in the coupling with dissipative charge current and magnetic field due to the momentum dependence of the relaxation time. By utilizing a power law parametrization to quantify the momentum dependence of the relaxation time, the anisotropic magnetic field-dependent shear coefficients in the Navier-Stokes limit have been investigated. The resulting viscous coefficients are seen to be sensitive to the momentum dependence of the relaxation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Modeling stochastic fluctuations in relativistic kinetic theory

Using the information current, we develop a Lorentz-covariant framework for modeling equilibrium fluctuations in relativistic kinetic theory in the grand-canonical ensemble. The resulting stochastic theory is proven to be causal and covariantly stable, and its predictions do not depend on the choice of spacetime foliation used to define the grand-canonical probabilities. As expected, in a box containing N > 5 particles, Boltzmann’s molecular chaos postulate is broken with (almost exact) probability N -1/2 , leading to a breakdown of the Boltzmann equation in small systems. Here, we also verify that, in ultrarelativistic gases, transient hydrodynamics already accounts for at least 80% of the equilibrium fluctuations of the stress-energy tensor at a given time. Finally, we compute the correlators at nonequal times for two selected collision kernels: that of a chemically active diluted solution, and that of ultrarelativistic scalar particles self-interacting via a quartic potential. For the former, we compute the density-density correlators analytically in real space, and dehydrodynamization of the stochastic theory is proven to occur whenever the mean free path diverges at high energy.

Astronomy & Astrophysics↗
DOE OSTI2024

Phenomenology of TMD parton distributions in Drell-Yan and 𝑍 0 boson production in a hadron structure oriented approach

We present a first practical implementation of a recently proposed hadron structure oriented (HSO) approach to transverse momentum dependent (TMD) phenomenology applied to Drell-Yan-like processes, including lepton pair production at moderate 𝑄 2 and 𝑍 0 boson production. We compare and contrast general features of our methodology with other common practices and emphasize aspects that we view as improvements. Separately, we discuss the importance of treatments that preserve various theoretical properties of the nonperturbative parts of TMD parametrizations for future applications that have the extraction of hadron structure as a primary goal. These include the HSO’s preservation of a basic TMD parton-model-like framework even while accounting for full TMD factorization and evolution, explicit preservation of the integral relationship between TMD and collinear parton density functions, and the ability to meaningfully compare different theoretical models of nonperturbative TMD parton distributions. In our examples, we show that there is significant sensitivity at moderate 𝑄 2 to both the form of the nonperturbative transverse momentum dependence and the parametrization of collinear parton densities. However, we also find that evolving to 𝑄 2 =𝑀$^2_𝑍$, without fitting, results in a satisfactory postdiction of existing data for 𝑍0 production, nearly independently of the modeling of nonperturbative transverse momentum behavior. We argue that this demonstrates that moderate 𝑄 measurements should be given greater weight than high 𝑄 measurements in extractions of nonperturbative transverse momentum dependence. We also obtain new extractions of the nonperturbative Collins-Soper kernel within the HSO approach. We discuss its features and compare with some earlier extractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Influence of neutrino-electron scattering and neutrino-pair annihilation on hypermassive neutron star

Here, we investigate the influence of inelastic neutrino microphysics in general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. In particular, we include species/energy groups coupled neutrino-matter interactions, such as inelastic neutrino-electron scattering and electron-positron annihilation kernels, into simulations up to 50 ms. Neutrino-electron inelastic scattering is known to have effective neutrino-matter energy exchange. We show that, with neutrino-electron inelastic scattering, simulations predict 75% higher disc mass with slightly different mass-averaged compositions, and 18% more ejected mass with similar distributions. The enhancement of the mass of the disc and the ejecta results in stronger baryon pollution, leading to less favorable jet launching environments. Furthermore, neutrino luminosities are about 50%, 40%, and 30% higher for electron neutrino, electron antineutrino, and heavy-lepton neutrinos. In contrast, we do not see any significant impacts due to electron-positron annihilation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Effective field theory for radiative corrections to charged-current processes. II. Axial-vector coupling

Here we discuss the hadronic structure-dependent radiative corrections to the axial-vector coupling that controls single-nucleon weak charged-current processes—commonly denoted by 𝑔 𝐴 . We match the Standard Model at the GeV scale onto chiral perturbation theory at next-to-leading order in the one-nucleon sector, in the presence of electromagnetic and weak interactions. As a result, we provide a representation for the corrections to 𝑔 𝐴 in terms of infrared finite convolutions of simple kernels with the single-nucleon matrix elements of time-ordered products of two and three quark bilinears (vector, axial-vector, and pseudoscalar). We discuss strategies to determine the required nonperturbative input from data, lattice-QCD (+QED), and possibly hadronic models. This work paves the way for a precise comparison of the values of the ratio 𝑔 𝐴 ⁡/𝑔 𝑉 extracted from experiment and from lattice QCD, which constrain physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

QCD factorization with multihadron fragmentation functions

Important aspects of quantum chromodynamics (QCD) factorization theorems are the properties of the objects involved that can be identified as universal. One example is that the definitions of parton densities and fragmentation functions for different types of hadrons differ only in the identity of the nonperturbative states that form the matrix elements, but are otherwise the same. This leads to independence of perturbative calculations on nonperturbative details of external states. It also lends support to interpretations of correlation functions as encapsulations of intrinsic nonperturbative properties. These characteristics have usually been presumed to still hold true in fragmentation functions even when the observed nonperturbative state is a small-mass cluster of n hadrons rather than simply a single isolated hadron. However, the multidifferential aspect of cross sections that rely on these latter types of fragmentation functions complicates the treatment of kinematical approximations in factorization derivations. That has led to recent claims that the operator definitions for fragmentation functions need to be modified from the single hadron case with nonuniversal prefactors. With such concerns as our motivation, we retrace the steps for factorizing the unpolarized semi-inclusive e + e − annihilation cross section and confirm that they do apply without modification to the case of a small-mass multihadron observed in the final state. In particular, we verify that the standard operator definition from single hadron fragmentation, with its usual prefactor, remains equally valid for the small-mass n -hadron case with the same hard parts and evolution kernels, whereas the more recently proposed definitions with nonuniversal prefactors do not. Our results reaffirm the reliability of most past phenomenological applications of dihadron fragmentation functions. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗