Engineering topics

Sironi, Lorenzo (ORCID:0000000212272754)

Publications and source records attributed to Sironi, Lorenzo (ORCID:0000000212272754).

9 records

DOE OSTI2025

Bulk Motions in the Black Hole Jet Sheath as a Candidate for the Comptonizing Corona

Abstract Using two-dimensional general relativistic resistive magnetohydrodynamic simulations, we investigate the properties of the sheath separating the black hole jet from the surrounding medium. We find that the electromagnetic power flowing through the jet sheath is comparable to the overall accretion power of the black hole. The sheath is an important site of energy dissipation as revealed by the copious appearance of reconnection layers and plasmoid chains. About 20% of the sheath power is dissipated between 2 and 10 gravitational radii. The plasma in the dissipative sheath moves along a nearly paraboloidal surface with transrelativistic bulk motions dominated by the radial component, whose dimensionless 4-velocity is ∼1.2 ± 0.5. In the frame moving with the mean (radially dependent) velocity, the distribution of stochastic bulk motions resembles a Maxwellian with an “effective bulk temperature” of ∼100 keV. Scaling the global simulation to Cygnus X-1 parameters gives a rough estimate of the Thomson optical depth across the jet sheath, ∼0.01–0.1, and it may increase in future magnetohydrodynamic simulations with self-consistent radiative losses. These properties suggest that the dissipative jet sheath may be a viable “coronal” region, capable of upscattering seed soft photons into a hard, nonthermal tail, as seen during the hard states of X-ray binaries and active galactic nuclei.

Sridhar, Navin (ORCID:0000000255199550)↗

DOE OSTI2024

A Magnetized Strongly Turbulent Corona as the Source of Neutrinos from NGC 1068

Abstract The cores of active galactic nuclei are potential accelerators of 10–100 TeV cosmic rays, in turn producing high-energy neutrinos. This picture was confirmed by the compelling evidence of a TeV neutrino signal from the nearby active galaxy NGC 1068, leaving open the question of what is the site and mechanism of cosmic-ray acceleration. One candidate is the magnetized turbulence surrounding the central supermassive black hole. Recent particle-in-cell simulations of magnetized turbulence indicate that stochastic cosmic-ray acceleration is nonresonant, in contrast to the assumptions of previous studies. We show that this has important consequences on a self-consistent theory of neutrino production in the corona, leading to a more rapid cosmic-ray acceleration than previously considered. The turbulent magnetic-field fluctuations needed to explain the neutrino signal are consistent with a magnetically powered corona. We find that strong turbulence, with turbulent magnetic energy density higher than 1% of the rest-mass energy density, naturally explains the normalization of the IceCube neutrino flux, in addition to the neutrino spectral shape. Only a fraction of the protons in the corona, which can be directly inferred from the neutrino signal, are accelerated to high energies. Thus, in this framework, the neutrino signal from NGC 1068 provides a testbed for particle acceleration in magnetized turbulence.

Fiorillo, Damiano F. G. (ORCID:0000000349279850)↗

DOE OSTI2024

Electron Heating in the Transrelativistic Perpendicular Shocks of Tilted Accretion Flows

Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is β pi1 ∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed is M s ∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperature T e2 exceeds the adiabatic expectation T e2,ad by an amount T e 2 / T e 2 , ad − 1 ≃ 0.0016 M s 3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratio T e1 / T i1 , which we vary from 0.1 to unity. We investigate the heating physics for M s ∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping at T e1 / T i1 = 1, whereas heating by B -parallel electric fields (i.e., parallel to the local magnetic field) dominates at T e1 / T i1 = 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.

Astronomy & Astrophysics↗

DOE OSTI2024

Electron Heating in 2D Particle-in-cell Simulations of Quasi-perpendicular Low-beta Shocks

Abstract We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta ( β p = 0.25) collisionless ion–electron shocks with mass ratio m i / m e = 200, fast Mach number M ms = 1 –4, and upstream magnetic field angle θ Bn = 55°–85° from the shock normal n ˆ . It is known that shock electron heating is described by an ambipolar, B -parallel electric potential jump, Δ ϕ ∥ , that scales roughly linearly with the electron temperature jump. Our simulations have Δ ϕ ∥ / ( 0.5 m i u sh 2 ) ∼ 0.1 –0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measure ϕ ∥ , including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate of ϕ ∥ in our low- β p shocks. We further focus on two θ Bn = 65° shocks: a M s = 4 ( M A = 1.8 ) case with a long, 30 d i precursor of whistler waves along n ˆ , and a M s = 7 ( M A = 3.2 ) case with a shorter, 5 d i precursor of whistlers oblique to both n ˆ and B ; d i is the ion skin depth. Within the precursors, ϕ ∥ has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the M s = 4 , θ Bn = 65° case, ϕ ∥ shows a weak dependence on the electron plasma-to-cyclotron frequency ratio ω pe /Ω ce , and ϕ ∥ decreases by a factor of 2 as m i / m e is raised to the true proton–electron value of 1836.

Tran, Aaron (ORCID:0000000334834890)↗

DOE OSTI2024

Long-term Evolution of Relativistic Unmagnetized Collisionless Shocks

Abstract We study a relativistic collisionless electron–positron shock propagating into an unmagnetized ambient medium using 2D particle-in-cell simulations of unprecedented duration and size. The shock generates intermittent magnetic structures of increasingly larger size as the simulation progresses. Toward the end of our simulation, at around 26,000 plasma times, the magnetic coherence scale approaches λ ∼ 100 plasma skin depths, both ahead and behind the shock front. We anticipate a continued growth of λ beyond the time span of our simulation, as long as the shock accelerates particles to increasingly higher energies. The post-shock field is concentrated in localized patches, which maintain a local magnetic energy fraction ε B ∼ 0.1. Particles randomly sampling the downstream fields spend most of their time in low field regions ( ε B ≪ 0.1) but emit a large fraction of the synchrotron power in the localized patches with strong fields ( ε B ∼ 0.1). Our results have important implications for models of gamma-ray burst afterglows.

Astronomy & Astrophysics↗

DOE OSTI2024

TeV Neutrinos and Hard X-Rays from Relativistic Reconnection in the Corona of NGC 1068

Abstract The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of nonthermal protons within a compact “coronal” region close to the central black hole. The acceleration mechanism of these nonthermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and nonthermal protons is established in the reconnection region. Motivated by recent 3D particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is dn p / dE p ∝ E p − 1 below the break and dn p / dE p ∝ E p − s above the break, with IceCube neutrino observations suggesting s ≃ 3. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron–positron pairs that makes the corona moderately Compton thick.

Astronomy & Astrophysics↗

DOE OSTI2023

The Origin of Power-law Spectra in Relativistic Magnetic Reconnection

Abstract Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factorγ≳ 3σ(here,σis the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is t acc ∼ γ η rec − 1 ω c − 1 ≃ 20 γ ω c − 1 , whereη rec ≃ 0.06 is the inflow speed in units of the speed of light andω c =eB 0 /mcis the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftert esc , when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuret esc in our simulations and find thatt esc ∼t acc forσ≳ few. This leads to a universal (i.e.,σ-independent) power-law spectrum dN free / d γ ∝ γ − 1 for the particles undergoing active acceleration, and dN / d γ ∝ γ − 2 for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

Astronomy & Astrophysics↗

DOE OSTI2022

Global Kinetic Modeling of the Intrabinary Shock in Spider Pulsars

Abstract Spider pulsars are compact binary systems composed of a millisecond pulsar and a low-mass companion. The relativistic magnetically dominated pulsar wind impacts onto the companion, ablating it and slowly consuming its atmosphere. The interaction forms an intrabinary shock, a proposed site of particle acceleration. We perform global fully kinetic particle-in-cell simulations of the intrabinary shock, assuming that the pulsar wind consists of plane-parallel stripes of alternating polarity and that the shock wraps around the companion. We find that particles are efficiently accelerated via shock-driven reconnection. We extract first-principles synchrotron spectra and light curves, which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat, F ν ∝ const; (2) when the pulsar spin axis is nearly aligned with the orbital angular momentum, the light curve displays two peaks, just before and after the pulsar eclipse (pulsar superior conjunction), separated in phase by ∼0.8 rad; (3) the peak flux exceeds the one at the inferior conjunction by a factor of 10.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

Abstract We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations σ ≲ 10 −4 and pair-loading factors Z ± ≲ 10 are studied, where Z ± is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field when σ exceeds a critical value σ L that decreases with Z ± . At σ ≲ σ L the shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease with Z ± , leading to lower σ L . (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as E ¯ e ∝ Z ± + 1 − 1 . (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as σ ≈ 5 × 10 −6 . The ions then become essentially thermal with mean energy E ¯ i , while electrons form a nonthermal tail, extending from E ∼ Z ± + 1 − 1 E ¯ i to E ¯ i . When σ = 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗