Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RELATIVISTIC PARTICLE”

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 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 OSTI2024

Larmor power limit for cyclotron radiation of relativistic particles in a waveguide

Cyclotron radiation emission spectroscopy (CRES) is a modern technique for high-precision energy spectroscopy, in which the energy of a charged particle in a magnetic field is measured via the frequency of the emitted cyclotron radiation. The He6-CRES collaboration aims to use CRES to probe beyond the standard model physics at the TeV scale by performing high-resolution and low-background beta-decay spectroscopy of 6 He and 19 Ne. Having demonstrated the first observation of individual, high-energy (0.1–2.5 MeV) positrons and electrons via their cyclotron radiation, the experiment provides a novel window into the radiation of relativistic charged particles in a waveguide via the time-derivative (slope) of the cyclotron radiation frequency, df c /dt. We show that analytic predictions for the total cyclotron radiation power emitted by a charged particle in circular and rectangular waveguides are approximately consistent with the Larmor formula, each scaling with the Lorentz factor of the underlying e ± as γ 4 . This hypothesis is corroborated with experimental CRES slope data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Comparison of integral equations used to study ${T}_{cc}^{+}$ for a stable D *

We perform a detailed comparison between three formalisms used in recent studies of DD* scattering at heavier-than-physical pion masses, which aim to understand the properties of the doubly-charmed tetraquark, ${T}_{cc}^{+}$ (3875). These methods are the three-particle relativistic field theory (RFT) formalism, the two-body Lippmann-Schwinger (LS) equation with chiral effective field theory potentials, and the two-particle relativistic framework proposed by Baião Raposo and Hansen (BRH approach). In a simplified single-channel setting, we derive the conditions under which the infinite-volume integral equations from the RFT and BRH approaches reduce to the LS form. We present numerical examples showing that differences between these methods can be largely removed by adjusting short-range couplings. We also address a number of technical issues in the RFT approach.

Hadronic Spectroscopy↗
DOE OSTI2025

Fluorescence time profile measurement of LAB based liquid scintillator in response to medium relativistic ion particles

Liquid scintillator is widely used in particle physics experiments due to its high light yield, good timing resolution, scalability and low cost. Certain liquid scintillators exhibit pulse shape discrimination capabilities because of difference in fluorescence timing properties induced by different particles. Its fluoresence timing properties have been measured mostly for radioactive decay sources at MeV energies. Here, we present a novel measurement of fluorescence time properties of Linear Alkyl Benzene (LAB) based liquid scintillator in response to high-energy ions of hydrogen (Z = 1), helium (Z = 2) and krypton at around 200–300 MeV/u for the first time. We compared the results to those from radioactive sources and observed a distinct dE / dX dependence, regardless of the particle type. These findings are essential for physics searches such as the diffuse supernova neutrino background in large liquid scintillator detectors like JUNO, and are also critical towards understanding the underlying scintillation timing mechanism.

Ion identification systems↗
DOE OSTI2025

Fermionic sub-GeV dark matter from evaporating primordial black holes at DarkSide-50

We present a search for boosted dark matter from primordial black hole (PBH) evaporation using the DarkSide-50 ionization-signal-only dataset corresponding to the experiment’s ( 12202 ± 180 ) kg d exposure. We focus on evaporation of PBHs with masses in the range [ 10 14 , 10 16 ] g producing Dirac fermionic dark matter particles with sub-GeV kinetic energy. These relativistic particles, with energies up to hundreds of MeV, can generate detectable signals for masses below O ( 100 ) MeV . The absence of a signal enables setting complementary limits to those derived from cosmological observations and direct detection searches for cosmic-ray-boosted dark matter.

Agnes, P. [Gran Sasso; GSSI, Aquila]↗
DOE OSTI2025

Effect of interparticle fields and radiation reaction on beam dynamics

The dynamics of relativistic particles in an intense electromagnetic field can be described by the Landau-Lifshitz (LL) equation, where the radiation reaction (RR) is accounted for via a self-force, and interparticle fields are often neglected as an approximation. However, the inclusion of interparticle fields is necessary to ensure energy-momentum conservation, particularly during coherent emission. Here we present (i) an analytical proof showing that the energy-momentum conservation law of the Hamilton-Rohrlich-Dirac action, which is divergence free and describes a generic system of interacting charges, respects causality and provides physically sensible results; (ii) a simple generalization of the LL equation for many particles evaluated as a function of the total field, i.e., the sum of the external and interparticle fields. By performing first-principles numerical simulations of a neutral, relativistic bunch of electrons and positrons (e − /e + ) colliding with a laser pulse, this theory is shown to satisfy energy-momentum conservation when interparticle fields and RR are simultaneously taken into account; and (iii) the combined effect of interparticle fields and RR primarily affects the tail of the particle energy distribution. Additionally, our first-principles simulations show that the effect of interparticle fields on beam energy loss becomes smaller when most of the radiated energy is incoherent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Cosmology of dark energy radiation

In this work, we quantify the cosmological signatures of dark energy radiation—a novel description of dark energy, which proposes that the dynamical component of dark energy is comprised of a thermal bath of relativistic particles sourced by thermal friction from a slowly rolling scalar field. For a minimal model with particle production emerging from first principles, we find that the abundance of radiation sourced by dark energy can be as large as Ω DER = 0.03, exceeding the bounds on relic dark radiation by three orders of magnitude. Although the background and perturbative evolution of dark energy radiation are distinct from Quintessence, we find that current and near-future cosmic microwave background and supernova data will not distinguish these models of dark energy. We also find that our constraints on all models are dominated by their impact on the expansion rate of the Universe. Considering extensions that allow the dark radiation to populate neutrinos, axions, and dark photons, we evaluate the direct detection prospects of a thermal background comprised of these candidates consistent with cosmological constraints on dark energy radiation. Our study indicates that a resolution of ~6 meV is required to achieve sensitivity to relativistic neutrinos compatible with dark energy radiation in a neutrino capture experiment on tritium. We also find that dark matter axion experiments lack sensitivity to a relativistic thermal axion background, even if enhanced by dark energy radiation, and dedicated search strategies are required to probe new parameter space. We derive constraints arising from a dark photon background from oscillations into visible photons, and find that viable parameter space can be explored with the late dark energy radiation experiment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Particle Acceleration in Relativistic Alfvénic Turbulence

Abstract Strong magnetically dominated Alfvénic turbulence is an efficient engine of nonthermal particle acceleration in a relativistic collisionless plasma. We argue that in the limit of strong magnetization, the type of energy distribution attained by accelerated particles depends on the relative strengths of turbulent fluctuations δ B 0 and the guide field B 0 . If δ B 0 ≪ B 0 , the particle magnetic moments are conserved, and the acceleration is provided by magnetic curvature drifts. Curvature acceleration energizes particles in the direction parallel to the magnetic field lines, resulting in log-normal tails of particle energy distribution functions. Conversely, if δ B 0 ≳ B 0 , interactions of energetic particles with intense turbulent structures can scatter particles, creating a population with large pitch angles. In this case, magnetic mirror effects become important, and turbulent acceleration leads to power-law tails of the energy distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

hymera

Hymera is a performance portable multiscale simulation framework based on parthenon and Kokkos, that designed to push relativistic particles, governed by slowly varying, quasi static background fields

Edelmann, Philipp [@LANL]↗
DOE OSTI2025

Pair Cascades in Magnetar Magnetospheres

Resonant inverse Compton scattering (RICS) of soft thermal photons by relativistic particles on closed magnetic field loops has been proposed to explain the hard emission observed up to and beyond 200 keV from magnetars. If particles injected at the base of the loops have Lorentz factors ≳10 2 , the RICS spectra will be attenuated by both one-photon pair production and photon splitting in the ultrastrong magnetar fields, producing additional spectral components from pair synchrotron radiation and split photons that produce further generations of pairs and split photons. We investigate such cascades initiated by the primary injected electrons through a Monte Carlo simulation, and study the cascade spectra and pair distributions. For most observer angles, the pair synchrotron and split-photon spectra dominate the RICS primary spectra and produce complex polarization signals. In particular, the synchrotron spectra are highly polarized with degree 40%–80%, are softer than the RICS spectra, and may account for the high polarization of some magnetar spectra observed by the Imaging X-ray Polarimetry Explorer above 3 keV.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2024

Measuring differential particle correlations in relativistic nuclear collisions

This study explores the transverse momentum ( p T ) dependencies of symmetric and asymmetric correlations (SC and ASC) with one and two particles of interest in Au + Au collisions at 200 GeV. Leveraging a multiphase transport (AMPT) model, the investigation delves into the sensitivity of these correlations to the final state effects, providing valuable insights into their potential for constraining the final state effects' p T dependencies. The HIJING model is employed as a benchmark for nonflow correlations, shedding light on their impact on interpreting SC and ASC data. Moreover, the study points out that differential SC and ASC with one and two particles of interest (POIs) typically incorporate contributions from event-plane angle fluctuations. Consequently, this work highlights the significance of SC and ASC with one and two POIs as valuable tools for investigating the p T nature of the final state effects and advocates for comprehensive experimental measurements across various beam energies and system sizes to enhance understanding and provide additional constraints for theoretical models. Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Acceleration of the particle-in-cell code Osiris with graphics processing units

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Particle Acceleration in Collisionless Magnetically Arrested Disks

We present the first collisionless realization of two-dimensional axisymmetric black hole accretion consistent with a persistent magnetically arrested disk state. The accretion flow, consisting of an ion-electron disk plasma combined with magnetospheric pair creation effects, is simulated using first-principles general-relativistic particle-in-cell methods. The simulation is evolved over significant dynamical timescales during which a quasisteady accretion state is reached with several magnetic flux eruption cycles. We include a realistic treatment of inverse Compton scattering and pair production, which allows for studying the interaction between the collisionless accretion flow and pair-loaded jet. Our findings indicate that magnetic flux eruptions associated with equatorial magnetic reconnection within the black hole magnetosphere and the formation of spark gaps are locations of maximal particle acceleration. Flux eruptions, starting near the central black hole, can trigger Kelvin-Helmholtz-like vortices at the jet-disk interface that facilitate efficient mixing between disk and jet plasma in this region. Transient periods of increased pair production following magnetic flux eruptions and reconnection events are responsible for most of the highly accelerated particles.

Accretion disk & black-hole plasma↗
DOE OSTI2025

Transition to Petschek Reconnection in Subrelativistic Pair Plasmas: Implications for Particle Acceleration

While relativistic magnetic reconnection in pair plasmas has emerged in recent years as a candidate for the origin of radiation from extreme astrophysical environments, the corresponding subrelativistic pair-plasma regime has remained less explored, leaving open the question of how relativistic physics affects reconnection. In this paper, we investigate the differences between these regimes by contrasting two-dimensional particle-in-cell simulations of reconnection in pair plasmas with relativistic magnetization (σ ≫ 1) and subrelativistic magnetization (σ < 1). By utilizing unprecedentedly large domain sizes and outflow boundary conditions, we demonstrate that lowering the magnetization results in a change in the reconnection geometry from a plasmoid chain to a Petschek geometry, where laminar exhausts bounded by slow-mode shocks emanate from a single diffusion region. We attribute this change to the reduced plasmoid production rate in the low-σ case: When the secondary tearing rate is sufficiently low, plasmoids are too few in number to prevent the system from relaxing into a stable Petschek configuration. This geometric change also affects particle energization: We show that while high-σ plasmoid chains generate power-law energy spectra, low-σ Petschek exhausts merely heat incoming plasma and yield negligible nonthermal acceleration. These results have implications for predicting the global current sheet geometry and the resulting energy spectra in a variety of systems.

Plasma astrophysics↗
DOE OSTI2025

Effect of magnetic field inclination on black hole jet power and particle acceleration

Rotating black holes are known to launch relativistic jets and accelerate particles if they accrete a magnetized plasma. It remains unclear, however, how the global magnetic field orientation affects the jet powering efficiency. We propose the first kinetic study of a collisionless plasma around a Kerr black hole embedded in a magnetic field that is inclined with respect to the black hole spin axis. Using three-dimensional general relativistic particle-in-cell simulations, we show that while oblique magnetic field configurations significantly reduce the jet power, particle acceleration still remains highly efficient. This suggests that black holes producing a weak jet might still be bright sources of nonthermal radiation and cosmic rays.

acceleration of particles↗