Search NASA⌕ Search

SEARCH · Search NASA

Results for “relativistic”

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 OSTI2025

Novel Relativistic Electronic Structure Theories for Actinide-Containing Compounds

Actinides of importance to basic energy sciences contain electrons moving at speed comparable to the speed of light. Reliable computational simulation of these electrons and hence actinide chemistry requires accurate description of relativistic effects. The present project advances computational actinide chemistry with development of new methodologies, algorithms, and computer programs in relativistic quantum chemistry, as well as applications to actinide chemistry and spectroscopy. A new “electrons-only” exact two-component approach has been developed to provide efficient treatments of relativistic effects, while maintaining chemical accuracy. New computational algorithms developed here extend the applicability of relativistic electron-correlation methods to larger molecules. The method-development work in this project also features the first implementation of analytic gradient technique for relativistic electron-correlation methods, which provides significantly enhanced ability to compute properties for molecules containing actinides. The applicability and usefulness of these new methods and computer programs have been demonstrated in calculations of actinide-containing molecules to facilitate understanding of actinide chemistry and spectroscopy.

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

Temporal evolution and polarization dependence of relativistic transparency in the ultra-short pulse regime

An ultra-intense laser pulse interacting with a solid target can heat electrons to relativistic energies, driving the plasma to transparency before it expands into the classically underdense regime. This effect, called relativistic transparency (RT), has wide ranging significance across many lines of inquiry in relativistic laser–plasma interactions. Here, we show the temporal evolution of the relativistically induced transparency of a laser heated target as measured by a lower intensity probe beam, providing the first time-resolved measurement of the return to opacity in a target undergoing RT. We also measure a shift in the ellipticity angle of the probe polarization by up to 7.8°. Supporting 3D particle-in-cell simulations corroborate these measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Nonparametric extensions of nuclear equations of state: Probing the breakdown scale of relativistic mean-field theory

Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. Here, we further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.

Equations of state of nuclear matter↗
DOE OSTI2025

Energy Spectrum and Mass Composition of Ultra-high-energy Cosmic Rays Originating from Relativistic Jets of Nearby Radio Galaxies

Relativistic jets of radio galaxies (RGs) are possible sources of ultra-high-energy cosmic rays (UHECRs). Recent studies combining relativistic hydrodynamic simulations with Monte Carlo particle transport have demonstrated that UHECRs can be accelerated to energies beyond 10 20 eV through shocks, turbulence, and relativistic shear in jet-induced flows of Fanaroff–Riley type RGs. The resulting time-asymptotic UHECR spectrum is well modeled by a double power law with an “extended” exponential cutoff, primarily shaped by relativistic shear acceleration. In this study, we adopt this novel source spectrum and simulate the propagation of UHECRs from nearby RGs using the CRPropa code. We focus on Virgo A (Vir A), Centaurus A (Cen A), Fornax A (For A), and Cygnus A (Cyg A), expected to be the most prominent UHECR sources among RGs. We then analyze the energy spectrum and mass composition of UHECRs arriving at Earth. We find that, due to the extended high-energy tail in the source spectrum, UHECRs from Vir A, which has a higher Lorentz factor, exhibit a higher flux at the highest energies and a lighter mass composition at Earth compared to those from Cen A and For A with lower Lorentz factors. Despite Cyg A having an even higher Lorentz factor, the large distance limits its contribution. With a small number of nearby prominent RGs, our findings suggest that if RGs are the major sources of UHECRs, the energy spectrum and mass composition of observed UHECRs would exhibit hemispheric differences between the Northern and Southern skies at the highest energies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Ultra high-energy cosmic rays from relativistic outflows in accretion induced collapse of white dwarfs

When a rapidly-rotating, highly magnetized white dwarf (WD) approaches the Chandrashekhar limit through mass accretion, it can undergo an accretion-induced collapse (AIC) to form a proto-neutron star or protomagnetar. The protomagnetar can drive a magnetically-dominated relativistic outflow, whose low entropy can lead to efficient formation of heavy nuclei. In this work, we propose that such relativistic outflows from AIC of WDs can contribute as sources of ultra high-energy cosmic rays (UHECRs). We model the acceleration of heavy nuclei in these relativistic outflows, and show that AICs can dominantly power the observed UHECRs, if a majority of them host relativistic outflows. Accounting for uncertainties in the acceleration mechanisms and AIC rates, AICs can contribute $\sim$ a few $10^{43} - 10^{45}\ {\rm erg \ Mpc}^{-3} {\rm yr}^{-1}$ in UHECR energy generation rate density, assuming iron-like nuclei.

Mukhopadhyay, Mainak [Fermilab; Chicago U., KICP] ↗
DOE OSTI2025

Relativistic core–valence-separated equation-of-motion coupled-cluster singles and doubles method: Efficient implementation and benchmark calculations

An efficient implementation for the relativistic exact two-component core–valence-separated equation-of-motion coupled-cluster singles and doubles (X2C-CVS-EOM-CCSD) method is reported. The explicit exclusion of pure valence excitations in the EOM-CCSD excited-state eigenvalue equations significantly improves the efficiency for calculations of core-excited states. Benchmark relativistic CVS-EOM-CC calculations with systematic inclusion of relativistic, correlation, and basis-set effects are shown to provide highly accurate results for core ionized and excited states involving heavy atoms.

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

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2026

Definitive Assessment of the Accuracy, Variationality, and Convergence of Relativistic Coupled Cluster and Density Matrix Renormalization Group in 100-Orbital Space

Accuracy, variationality, and convergence underpin the reliability of modern electronic structure methods, yet definitive benchmarks in the relativistic regime remain elusive due to the absence of numerically exact full configuration interaction (CI) references. Recent algorithmic advances in the CI framework, enabled by the small-tensor-product (STP) decomposition approach, have dramatically extended the tractable size of the configuration space, making numerically exact CI calculations feasible in large active spaces previously beyond reach. In this paper, we employ the recently developed STP-CI framework to perform large-scale numerically exact CI calculations and directly benchmark relativistic coupled cluster and density matrix renormalization group methods. Definitive benchmarking of approximate relativistic electronic structure methods is ensured through the application of the gap theorem, which provides rigorous error bounds on the CI reference and establishes a controlled standard for assessing accuracy, variationality, and convergence.

Chemical calculations↗
DOE OSTI2024

Collective excitations and low-energy ionization signatures of relativistic particles in silicon detectors

Abstract Solid-state detectors with a low energy threshold have several applications, including searches of non-relativistic halo dark-matter particles with sub-GeV masses. When searching for relativistic, beyond-the-Standard-Model particles with enhanced cross sections for small energy transfers, a small detector with a low energy threshold may have better sensitivity than a larger detector with a higher energy threshold. In this paper, we calculate the low-energy ionization spectrum from high-velocity particles scattering in a dielectric material. We consider the full material response including the excitation of bulk plasmons. We generalize the energy-loss function to relativistic kinematics, and benchmark existing tools used for halo dark-matter scattering against electron energy-loss spectroscopy data. Compared to calculations commonly used in the literature, such as the Photo-Absorption-Ionization model or the free-electron model, including collective effects shifts the recoil ionization spectrum towards higher energies, typically peaking around 4–6 electron-hole pairs. We apply our results to the three benchmark examples: millicharged particles produced in a beam, neutrinos with a magnetic dipole moment produced in a reactor, and upscattered dark-matter particles. Our results show that the proper inclusion of collective effects typically enhances a detector’s sensitivity to these particles, since detector backgrounds, such as dark counts, peak at lower energies.

Physics↗
DOE OSTI2025

The role of electric dominance for particle injection in relativistic reconnection

ABSTRACT Magnetic reconnection in relativistic plasmas – where the magnetization $\sigma \gg 1$ – is regarded as an efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. We employ two-dimensional (2D) particle-in-cell simulations of relativistic pair-plasma reconnection with vanishing guide field and outflow boundaries to quantify the impact of the energy gain occurring in regions of electric dominance ($E\gt B$) for the early stages of particle acceleration (i.e. the ‘injection’ stage). Given an injection threshold energy $\epsilon ^\ast =\sigma /4$ for the particles that eventually attain energy $\epsilon _{\rm T}$ by time T, we calculate the mean fractional contribution $\zeta (\epsilon ^\ast ,\epsilon _{\rm T})$ by $E\gt B$ fields to particle energization at the time when the threshold $\epsilon ^\ast$ is reached. We find that $\zeta$ monotonically increases with $\sigma$ and $\epsilon _{\rm T}$; for $\sigma \gtrsim 50$ and $\epsilon _{\rm T}/\sigma \gtrsim 8$, we find that $\gtrsim 80~{{\ \rm per\ cent}}$ of the energy gain obtained before reaching $\epsilon ^\ast =\sigma /4$ occurs in $E\gt B$ regions. We find that $\zeta$ is independent of simulation box size $L_x$, as long as $\epsilon _{\rm T}$ is normalized to the maximum particle energy, which scales as $\epsilon _{\rm max}\propto L_{\rm x}^{1/2}$ in 2D. The distribution of energy gains $\epsilon _{\chi }$ acquired in $E\gt B$ regions can be modelled as $\mathrm{ d}N/\mathrm{ d}\epsilon _{\chi }\propto \epsilon _{\chi }^{-0.35}\exp [-(\epsilon _{\chi }/0.06\, \sigma)^{0.5}]$. Our results help assess the role of electric dominance in relativistic reconnection with vanishing guide fields, which is realized in the magnetospheres of black holes and neutron stars.

Gupta, Sanya (ORCID:0000000151944384)↗
DOE OSTI2026

Relativistic Exact Two-Component Theory in the Generalized Pseudospectral Representation

We present a formulation and implementation of exact two-component (X2C) relativistic theory in the generalized pseudospectral representation. When combined with the Hartree-Fock-Slater framework, this approach enables efficient and accurate treatments of scalar-relativistic and spin-orbit effects in atomic electronic structures without the computational overhead of four-component methods. Benchmark calculations across light and heavy elements demonstrate that the our X2C scheme yields substantially more accurate relativistic corrections to core-electron binding energies than perturbational Breit-Pauli treatments while converging more rapidly with respect to basis size. The method provides an improved computational framework for modeling ultrafast x-ray-induced processes in heavy-element systems.

Wang, Xubo↗
DOE OSTI2026

Relativistic collapse of the classical triple bond in the CBi − molecular ion

The conventional framework for chemical bonding between main-group elements involves separate σ and π orbitals to describe multiple bonds. However, relativistic effects mix these orbitals in molecules containing heavy elements through spin–orbit coupling, leaving the total angular-momentum projection (ω) as the only good quantum number. Direct experimental evidence that relativistic effects change the σ-π bonding framework has remained elusive. Here, we probe the carbon-bismuth triple bond in the CBi − anion using high-resolution cryogenic photoelectron spectroscopy, coupled with relativistic four-component Dirac-Coulomb coupled-cluster calculations. Furthermore, even though the CBi − anion is isovalent to the well-known CN − species, we demonstrate that the traditional σ + 2π triple-bond picture collapses into a pure π-like |ω| = 3/2 and two |ω| = 1/2 Kramers pairs containing substantial σ/π mixing.

Kahraman, Deniz [Brown University, Providence, RI ↗
DOE OSTI2025

Tearing-mediated Alfvénic Turbulence in a Relativistic Plasma

In astrophysical relativistic plasmas, Alfvénic turbulence exists across a wide range of scales. Similar to a nonrelativistic case, this turbulence can be influenced by the tearing instability. We argue that in an ultrarelativistic pair plasma the tearing effects become significant at the critical scale given by ${a}_{t}\sim {d}_{\mathrm{rel}}{\left({L}_{\perp }/{d}_{\mathrm{rel}}\right)}^{n/(4+5n)}$, where L ⊥ represents the outer scale of the turbulence and d rel the relativistic electron inertial scale. Here n = 1 corresponds to the Harris magnetic configurations, while n = 2 applies to localized (non-Harris) magnetic profiles of turbulence-generated current sheets. Below the critical scale, the energy spectrum of magnetic fluctuations steepens to W(k)dk ∝ k -(7n+2)/3n dk. However, this happens only when $\tilde{\sigma }\ll {\left({L}_{\perp }/{d}_{\mathrm{rel}}\right)}^{(2+4n)/(4+5n)}$, where $\tilde{\sigma }$ is the magnetization parameter based on the magnetic fluctuations. If the opposite inequality holds, the Alfvénic cascade is instead affected by charge starvation at scales below ${a}_{\mathrm{ch}}\sim {d}_{\mathrm{rel}}{\left({d}_{\mathrm{rel}}/{L}_{\perp }\right)}^{1/3}{\tilde{\sigma }}^{2/3}$, and the energy spectrum becomes steeper than W(k)dk ∝ k -3 dk. When d rel is replaced by its nonrelativistic equivalent, our results agree with the nonrelativistic analysis of N. F. Loureiro & S. Boldyrev (2018). In the Harris case n = 1, they agree with the recent relativistic studies by I. Demidov & Y. Lyubarsky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Bridging Scales in Black Hole Accretion and Feedback: Relativistic Jet Linking the Horizon to the Host Galaxy

Simulating black hole (BH) accretion and feedback from the BH horizon to galactic scales is extremely challenging, as it involves a vast range of scales. Recently, our multizone method has successfully achieved global dynamical steady states of hot accretion flows in 3D general relativistic magnetohydrodynamic simulations by tracking the bidirectional interaction between a nonspinning BH and its host galaxy. In this paper, we present technical improvements to the method and apply it to spin a * = 0.9 BHs, which power relativistic jets. We first test the new multizone setup with a smaller Bondi radius, R B ≈ 400 r g , where r g is the gravitational radius. The strongly magnetized accretion launches a relativistic jet with an intermediate feedback efficiency η ∼ 30%, in between that of a prograde (η ∼ 100%) and retrograde (η ∼ 10%) torus. Interestingly, both prograde and retrograde simulations also eventually converge to the same intermediate efficiency when evolved long enough, as accumulated magnetic fields remove gas rotation. We then extend strongly magnetized simulations to larger Bondi radii, R B ≈ 2 × 10 3 , 2 × 10 4 , 2 × 10 5 r g . We find that the BH accretion rate $\dot{M}$ is suppressed with respect to the Bondi rate as $\dot{M}_{\textrm{B}}$ as $\dot{M}/\dot{M}_{\textrm{B}} ∝ R_{\textrm{B}}^{-1/2}$. However, despite some variability, the time-averaged feedback efficiency remains at η ∼ 30%, independent of R B . This suggests that BH feedback efficiency in hot accretion flows is mainly governed by the BH spin (a * ) rather than by the galactic properties (R B ). From these first-principles simulations, we provide a feedback subgrid prescription for cosmological simulations: $\dot{E}_{\textrm{fb}} = 2$ x $10^{-3}[R_{\textrm{B}}/(2$ x $10^5 r_g)]^{-1/2}$ $\dot{M}_{\textrm{B}}c^2$ for BH spin a * = 0.9.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Particle Acceleration and Pitch-angle Evolution in Relativistic Turbulence

Synchrotron radiation detected from relativistic astrophysical objects such as pulsar-wind nebulae and jets from active galactic nuclei depends on the magnetic fields and the distribution functions of energetic electrons in these systems. Relativistic magnetically dominated turbulence has been recognized as an efficient mechanism for structure formation and nonthermal particle acceleration in these environments. Recent numerical simulations of relativistic turbulence have provided insights into the energy distribution functions of accelerated electrons. Much less is currently understood about their pitch angle distributions, which are crucial for accurately interpreting the spectra of synchrotron radiation. We perform a detailed case study of the pitch angle distributions formed during the process of turbulent acceleration for B 0 /δB 0 = 10 and $\tilde{σ}$ 0 ~40, where B0 is the uniform component of the magnetic field, δB0 is the fluctuating component, and $\tilde{σ}$ 0 is the plasma magnetization based on the magnetic fluctuations. We find that even minimal numerical noise can cause substantial pitch angle scattering, but we demonstrate techniques for overcoming the numerical challenges associated with the evolution of very small pitch angles. Our numerical results are consistent with the phenomenological considerations found in C. Vega et al. (2024a, 2025).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The Effect of Inverse Compton Losses on Particle Acceleration in Three-dimensional Relativistic Reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived “free phase” where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a “trapped phase,” where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection, including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor γrad exceeds the magnetization σ. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is dN free /dγ ∝ γ −1 , as in the uncooled case. The spectrum of trapped electrons above the cooling break γcool (in the range γ cool < γ < γ rad ) is dN/dγ ∝ γ −3 , steeper than the scaling dN/dγ ∝ γ −2 of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by Zhang et al. for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. II. General Relativistic Module

Black hole (BH) environments often host plasmas that are fully collisionless or contain intrinsically collisionless regions, including relativistic jets and coronae, where particle energization is ubiquitous. Capturing the physics of these systems requires numerical methods capable of modeling relativistic, magnetized, collisionless plasmas in strong gravitational fields. In this work, we introduce the general relativistic module for Entity—the first open-source, coordinate-agnostic, performance-portable, particle-in-cell code. The code enables fast axisymmetric simulations of collisionless plasmas around BHs on any modern high-performance computing architecture (both GPUs and CPUs).

Galishnikova, Alisa [Flatiron Institute, New York,↗
DOE OSTI2026

The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Relativistic magnetic reconnection is a key mechanism for dissipating magnetic energy and accelerating particles in astrophysics. In the absence of radiative cooling, recent particle-in-cell (PIC) simulations have shown that high-energy particles gain most of their energy in the upstream region, during a short-lived "free phase" where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a "trapped phase", where no significant energization occurs. Here, we perform a suite of 3D PIC simulations of relativistic reconnection including inverse Compton (IC) losses in the weakly cooled regime in which the radiation-reaction-limited Lorentz factor $γ_{\rm rad}$ exceeds the magnetization $σ$. We show that electron cooling losses do not appreciably alter the reconnection rate, the structure of the layer, and the physics of particle acceleration in the free phase, so the spectrum of free electrons is $dN_{\rm free}/dγ\propto γ^{-1}$, as in the uncooled case. The spectrum of trapped electrons above the cooling break $γ_{\rm cool}$ (in the range $γ_{\rm cool}<γ<γ_{\rm rad}$) is $dN/dγ\propto γ^{-3}$, steeper than the scaling $dN/dγ\propto γ^{-2}$ of uncooled simulations. This confirms that no significant particle energization occurs during the trapped phase. Our results validate the model by arXiv:2302.12269 for particle acceleration in 3D relativistic reconnection, and imply that radiative emission models of reconnection-powered astrophysical sources should employ a two-zone structure, that differentiates between free, rapidly accelerating particles and trapped, passively cooling particles.

FOS: Physical sciences↗