BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Lorenzo Sironi

Lorenzo Sironi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEphysics.plasm-phastro-ph.SRphysics.space-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

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 &#34;free phase&#34; where they meander between the two sides of the layer; when they get captured/trapped by the downstream flux ropes, they undergo a &#34;trapped phase&#34;, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Filamentation of Fast Radio Bursts in magnetar winds

Magnetars are the most promising progenitors of Fast Radio Bursts (FRBs). Strong radio waves propagating through the magnetar wind are subject to non-linear effects, including modulation/filamentation instabilities. We derive the dispersion relation for modulations of strong waves propagating in magnetically-dominated pair plasmas focusing on dimensionless strength parameters $a_0\lesssim 1$, and discuss implications for FRBs. As an effect of the instability, the FRB radiation intensity develops sheets perpendicular to the direction of the wind magnetic field. When the FRB front expands outside the radius where the instability ends, the radiation sheets are scattered due to diffraction. The FRB scattering timescale depends on the properties of the magnetar wind. In a cold wind, the typical scattering timescale is $τ_{\rm sc}\sim{\rm\; μs-ms}$ at the frequency $ν\sim 1{\rm\; GHz}$. The scattering timescale increases at low frequencies, with the scaling $τ_{\rm sc}\proptoν^{-2}$. The frequency-dependent broadening of the brightest pulse of FRB 181112 is consistent with this scaling. From the scattering timescale of the pulse, one can estimate that the wind Lorentz factor is larger than a few tens. In a warm wind, the scattering timescale can approach $τ_{\rm sc}\sim{\rm\; ns}$. Then scattering produces a frequency modulation of the observed intensity with a large bandwidth, $Δν\sim 1/τ_{\rm sc}\gtrsim 100{\rm\; MHz}$. Broadband frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Global Kinetic Modeling of the Intrabinary Shock in Spider Pulsars

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 lightcurves which are in good agreement with X-ray observations: (1) the synchrotron spectrum is nearly flat, $F_ν\propto {\rm 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 $\sim 0.8\, {\rm rad}$; (3) the peak flux exceeds the one at inferior conjunction by a factor of ten. We demonstrate that the double-peaked signature in the lightcurve is due to Doppler boosting in the post-shock flow.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ion and Electron Acceleration in Fully Kinetic Plasma Turbulence

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-$β$ plasmas ($β$ is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher $β$. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions stay roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations $σ\lesssim 10^{-4}$ and pair-loading factors $Z_\pm \lesssim 10$ are studied, where $Z_\pm$ 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 $σ_{\rm L}$ that decreases with $Z_\pm$. At $σ\lesssimσ_{\rm L}$ the shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease with $Z_\pm$, leading to lower $σ_{\rm 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 $\overline{E}_{\rm e}\propto (Z_\pm+1)^{-1}$. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as $σ\approx 5\times 10^{-6}$. The ions then become essentially thermal with mean energy $\overline{E}_{\rm i}$, while electrons form a nonthermal tail, extending from $E\sim (Z_\pm + 1)^{-1}\overline{E}_{\rm i}$ to $\overline{E}_{\rm 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-ideal fields solve the injection problem in relativistic reconnection

Magnetic reconnection in relativistic plasmas is well established as a fast and efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. With particle-in-cell simulations, we demonstrate the importance of non-ideal fields for the early stages (&#34;injection&#34;) of particle acceleration. Most of the particles ending up with high energies (near or above the mean magnetic energy per particle) must have passed through non-ideal regions where the assumptions of ideal magnetohydrodynamics are broken (i.e., regions with $E>B$ or nonzero $E_\parallel={\bf E}\cdot {\bf B}/B$), whereas particles that do not experience non-ideal fields end up with Lorentz factors of order unity. Thus, injection by non-ideal fields is a necessary prerequisite for further acceleration. Our results have important implications for the origin of nonthermal particles in high-energy astrophysical sources.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-resonant particle acceleration in strong turbulence: comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of non-thermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D PIC, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of non-resonant particle acceleration in strongly magnetized turbulence~\cite{2021PhRvD.104f3020L}, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this non-resonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. We also observe that the parallel shear contribution tends to dominate the physics of energization in the PIC simulations, while in the MHD incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Relativistic Asymmetric Reconnection

We derive basic scaling equations for relativistic magnetic reconnection in the general case of asymmetric inflow conditions and obtain predictions for the outflow Lorentz factor and the reconnection rate. Kinetic Particle-in-Cell simulations show that the outflow speeds as well as the nonthermal spectral index are constrained by the inflowing plasma with the weaker magnetic energy per particle, in agreement with the scaling predictions. These results are significant for understanding non-thermal emission from reconnection in magnetically-dominated, astrophysical systems, many of which may be asymmetric in nature. The results provide a quantitative approach for including asymmetry on reconnection in the relativistic regime.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A fully-kinetic model for orphan gamma-ray flares in blazars

Blazars emit a highly-variable non-thermal spectrum. It is usually assumed that the same non-thermal electrons are responsible for the IR-optical-UV emission (via synchrotron) and the gamma-ray emission (via inverse Compton). Hence, the light curves in the two bands should be correlated. Orphan gamma-ray flares (i.e., lacking a luminous low-frequency counterpart) challenge our theoretical understanding of blazars. By means of large-scale two-dimensional radiative particle-in-cell simulations, we show that orphan gamma-ray flares may be a self-consistent by-product of particle energization in turbulent magnetically-dominated pair plasmas. The energized particles produce the gamma-ray flare by inverse Compton scattering an external radiation field, while the synchrotron luminosity is heavily suppressed since the particles are accelerated nearly along the direction of the local magnetic field. The ratio of inverse Compton to synchrotron luminosity is sensitive to the initial strength of turbulent fluctuations (a larger degree of turbulent fluctuations weakens the anisotropy of the energized particles, thus increasing the synchrotron luminosity). Our results show that the anisotropy of the non-thermal particle population is key to modeling the blazar emission.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Toy Model for the Time-Frequency Structure of Fast Radio Bursts: Implications for the CHIME Burst Dichotomy

We introduce a toy model for the time-frequency structure of fast radio bursts (FRB), in which the observed emission is produced as a narrowly-peaked intrinsic spectral energy distribution sweeps down in frequency across the instrumental bandpass as a power-law in time. Though originally motivated by emission models which invoke a relativistic shock, the model could in principle apply to a wider range of emission scenarios. We quantify the burst&#39;s detectability using the frequency bandwidth over which most of its signal-to-noise ratio (SNR) is accumulated. We demonstrate that by varying just a single parameter of the toy model-the power-law index βof the frequency drift rate-one can transform a long (and hence preferentially time-resolved) burst with a narrow time-integrated spectrum into a shorter burst with a broad power-law time-integrated spectrum. We suggest that burst-to-burst diversity in the value of βcould generate the dichotomy between burst duration and frequency-width recently found by CHIME. In shock models, the value of βis related to the radial density profile of external medium, which in light of the preferentially longer duration of bursts from repeating sources may point to diversity in the external environments surrounding repeating versus one-off FRB sources.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Pitch Angle Anisotropy Controls Particle Acceleration and Cooling in Radiative Relativistic Plasma Turbulence

Nature&#39;s most powerful high-energy sources are capable of accelerating particles to high energy and radiate it away on extremely short timescales, even shorter than the light crossing time of the system. It is yet unclear what physical processes can produce such an efficient acceleration, despite the copious radiative losses. By means of radiative particle-in-cell simulations, we show that magnetically dominated turbulence in pair plasmas subject to strong synchrotron cooling generates a nonthermal particle spectrum with a hard power-law range (slope $p \sim 1$) within a few eddy turnover times. Low pitch-angle particles can significantly exceed the nominal radiation-reaction limit, before abruptly cooling down. The particle spectrum becomes even harder ($p < 1$) over time owing to particle cooling with an energy-dependent pitch-angle anisotropy. The resulting synchrotron spectrum is hard ($νF_ν\propto ν^s$ with $s \sim 1$). Our findings have important implications for understanding the nonthermal emission from high-energy astrophysical sources, most notably the prompt phase of gamma-ray bursts and gamma-ray flares from the Crab nebula.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Shock-powered radio precursors of neutron star mergers from accelerating relativistic binary winds

During the final stages of a compact object merger, if at least one of the binary components is a magnetized neutron star (NS), then its orbital motion substantially expands the NS&#39;s open magnetic flux -- and hence increases its wind luminosity -- relative to that of an isolated pulsar. As the binary orbit shrinks due to gravitational radiation, the power and speed of this binary-induced inspiral wind may (depending on pair loading) secularly increase, leading to self-interaction and internal shocks in the outflow beyond the binary orbit. The magnetized forward shock can generate coherent radio emission via the synchrotron maser process, resulting in an observable radio precursor to binary NS merger. We perform 1D relativistic hydrodynamical simulations of shock interaction in the accelerating binary NS wind, assuming that the inspiral wind efficiently converts its Poynting flux into bulk kinetic energy prior to the shock radius. This is combined with the shock maser spectrum from particle-in-cell simulations, to generate synthetic radio light curves. The precursor burst with a fluence of $\sim1$ Jy$\cdot$ms at $\sim$GHz frequencies lasts $\sim 1-500$ ms following the merger for a source at $\sim3$ Gpc ($B_{\rm d}/10^{12}$ G)$^{8/9}$, where $B_{\rm d}$ is the dipole field strength of the more strongly-magnetized star. Given an outflow geometry concentrated along the binary equatorial, the signal may be preferentially observable for high-inclination systems, i.e. those least likely to produce a detectable gamma-ray burst.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stochastic Electron Acceleration by Temperature Anisotropy Instabilities Under Solar Flare Plasma Conditions

Using 2D particle-in-cell (PIC) plasma simulations we study electron acceleration by temperature anisotropy instabilities, assuming conditions typical of above-the-loop-top (ALT) sources in solar flares. We focus on the long-term effect of $T_{e,\perp} > T_{e,\parallel}$ instabilities by driving the anisotropy growth during the entire simulation time, through imposing a shearing or a compressing plasma velocity ($T_{e,\perp}$ and $T_{e,\parallel}$ are the temperatures perpendicular and parallel to the magnetic field). This magnetic growth makes $T_{e,\perp}/T_{e,\parallel}$ grow due to electron magnetic moment conservation, and amplifies the ratio $ω_{ce}/ω_{pe}$ from $\sim 0.53$ to $\sim 2$ ($ω_{ce}$ and $ω_{pe}$ are the electron cyclotron and plasma frequencies, respectively). In the regime $ω_{ce}/ω_{pe}\lesssim 1.2-1.7$ the instability is dominated by oblique, quasi-electrostatic (OQES) modes, and the acceleration is inefficient. When $ω_{ce}/ω_{pe}$ has grown to $ω_{ce}/ω_{pe}\gtrsim 1.2-1.7$, electrons are efficiently accelerated by the inelastic scattering provided by unstable parallel, electromagnetic z (PEMZ) modes. After $ω_{ce}/ω_{pe}$ reaches $\sim 2$, the electron energy spectra show nonthermal tails that differ between the shearing and compressing cases. In the shearing case, the tail resembles a power-law of index $α_s \sim$ 2.9 plus a high-energy bump reaching $\sim 300$ keV. In the compressing runs, $α_s \sim$ 3.7 with a spectral break above $\sim 500$ keV. This difference can be explained by the different temperature evolutions in these two types of simulations, suggesting a critical role played by the type of anisotropy driving, $ω_{ce}/ω_{pe}$ and the electron temperature in the efficiency of the acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A First-Principle Model for Polarization Swings during Reconnection-Powered Flares

We show that magnetic reconnection in a magnetically-dominated fast-cooling plasma can naturally produce bright flares accompanied by rotations in the synchrotron polarization vector. With particle-in-cell simulations of reconnection, we find that flares are powered by efficient particle acceleration at the interface of merging magnetic flux ropes, or &#34;plasmoids&#34;. The accelerated particles stream through the post-merger plasmoid towards the observer, thus progressively illuminating regions with varying plane-of-sky field direction, and so leading to a rotation in the observed polarization vector. Our results provide evidence for magnetic reconnection as the physical cause of high-energy flares from the relativistic jets of blazars (which recent observations have shown to be frequently associated with polarization rotations), and provide a first-principle physical mechanism for such flares.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraints on the Engines of Fast Radio Bursts

We model the sample of fast radio bursts (FRB), including the newly discovered CHIME repeaters, using the synchrotron blast wave model of Metzger, Margalit & Sironi (2019). This model postulates that FRBs are precursor radiation from ultra-relativistic magnetized shocks generated as flare ejecta from a central engine collide with an effectively stationary external medium. Downward drifting of the burst frequency structure naturally arises from deceleration of the blast-wave. The data are consistent with FRBs being produced by flares of energy $E_{\rm flare} \sim 10^{43}-10^{46}(f_ξ/10^{-3})^{-4/5}$ erg, where $f_ξ$ is the maser efficiency, and minimum bulk Lorentz factors $Γ\approx 10^2-10^3$, which generate the observed FRBs at shock radii $r_{\rm sh} \sim 10^{12}-10^{13}$ cm. We infer upstream densities $n_{\rm ext}(r_{\rm sh}) \sim 10^{2}-10^{4}$ cm$^{-3}$ and radial profiles $n_{\rm ext} \propto r^{-k}$ showing a range of slopes $k \approx [-2,1]$ (which are seen to evolve between bursts), broadly consistent with the upstream medium being the inner edge of an ion-loaded shell released by a recent energetic flare. The burst timescales, energetics, rates, and external medium properties are consistent with repeating FRBs arising from young, hyper-active flaring magnetars, but the methodology presented is generally applicable to any central engine which injects energy impulsively into a dense magnetized medium. Uncertainties and variations of the model are discussed, including the effects of the strong electric field of the FRB wave (strength parameter $a \gg 1$) on the upstream medium. One-dimensional particle-in-cell simulations of magnetized shocks into a pair plasma are presented which demonstrate that high maser efficiency can be preserved, even in the limit $a \gg 1$ in which the FRB wave accelerates the upstream electrons to ultra-relativistic speeds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electron Heating in Perpendicular Low-Beta Shocks

Collisionless shocks heat electrons in the solar wind, interstellar blast waves, and hot gas permeating galaxy clusters. How much shock heating goes to electrons instead of ions, and what plasma physics controls electron heating? We simulate 2-D perpendicular shocks with a fully kinetic particle-in-cell code. For magnetosonic Mach number $\mathcal{M}_\mathrm{ms} \sim 1$-$10$ and plasma beta $β_\mathrm{p} \lesssim 4$, the post-shock electron/ion temperature ratio $T_\mathrm{e}/T_\mathrm{i}$ decreases from $1$ to $0.1$ with increasing $\mathcal{M}_\mathrm{ms}$. In a representative $\mathcal{M}_\mathrm{ms}=3.1$, $β_\mathrm{p}=0.25$ shock, electrons heat above adiabatic compression in two steps: ion-scale $E_\parallel = \vec{E} \cdot \hat{b}$ accelerates electrons into streams along $\vec{B}$, which then relax via two-stream-like instability. The $\vec{B}$-parallel heating is mostly induced by waves; $\vec{B}$-perpendicular heating is mostly adiabatic compression by quasi-static fields.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hard Synchrotron Spectra from Magnetically Dominated Plasma Turbulence

Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically-dominated ($σ_0\gg1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim σ_0 γ_{th0}$ (here, $γ_{th0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are roughly isotropic. This energy-dependent anisotropy leads to a synchrotron spectral flux $νF_ν\propto ν^s$ that is much harder than for isotropic particles. Remarkably, for $σ_0\gg1$ we find that the solid-angle-averaged spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim δB_{\rm rms0}^2/B_0^2\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Our findings may help explaining the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of Pulsar Wind Nebulae.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kinetic Simulations of Radiative Magnetic Reconnection in the Coronae of Accreting Black Holes

We perform two- and three-dimensional particle-in-cell simulations of reconnection in magnetically-dominated pair plasmas subject to strong Compton cooling. Reconnection under such conditions operates in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses and locked in magnetically-dominated plasmoids with small thermal pressure. Compton drag clears cavities inside plasmoids and affects their bulk motions. These effects, however, weakly change the reconnection rate and the plasmoid size distribution from those in non-radiative reconnection. This demonstrates that the reconnection dynamics is governed by similar magnetic stresses in both cases and weakly affected by thermal pressure. We examine the energy distribution of particles energized by radiative reconnection and observe two components. (1) A mildly-relativistic peak, which results from bulk motions of cooled plasmoids. This component receives most of the dissipated reconnection power and dominates the output X-ray emission. The peak has a quasi-Maxwellian shape with an effective temperature of 100 keV. Thus, it mimics thermal Comptonization used previously to fit hard-state spectra of accreting black holes. (2) A high-energy tail, which receives 20% of dissipated reconnection power. It is populated by particles accelerated impulsively at X-points or &#34;picked up&#34; by fast outflows from X-points. The high-energy particles immediately cool, and their inverse Compton emission explains the MeV spectral tail detected in the hard state of Cyg X-1. Our first-principle simulations support reconnection as a mechanism powering hard X-ray emission from accreting black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reconnection-driven particle acceleration in relativistic shear flows

Particle energization in shear flows is invoked to explain non-thermal emission from the boundaries of relativistic astrophysical jets. Yet, the physics of particle injection, i.e., the mechanism that allows thermal particles to participate in shear-driven acceleration, remains unknown. With particle-in-cell simulations, we study the development of Kelvin-Helmholtz (KH) instabilities seeded by the velocity shear between a relativistic magnetically-dominated electron-positron jet and a weakly magnetized electron-ion ambient plasma. We show that, in their nonlinear stages, KH vortices generate kinetic-scale reconnection layers, which efficiently energize the jet particles, thus providing a first-principles mechanism for particle injection into shear-driven acceleration. Our work lends support to spine-sheath models of jet emission - with a fast core/spine surrounded by a slower sheath - and can explain the origin of radio-emitting electrons at the boundaries of relativistic jets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The synchrotron maser emission from relativistic magnetized shocks: Dependence on the pre-shock temperature

Electromagnetic precursor waves generated by the synchrotron maser instability at relativistic magnetized shocks have been recently invoked to explain the coherent radio emission of Fast Radio Bursts. By means of two-dimensional particle-in-cell simulations, we explore the properties of the precursor waves in relativistic electron-positron perpendicular shocks as a function of the pre-shock magnetization $σ\gtrsim 1$ (i.e., the ratio of incoming Poynting flux to particle energy flux) and thermal spread $Δγ\equiv kT/mc^2=10^{-5}-10^{-1}$. We measure the fraction $f_ξ$ of total incoming energy that is converted into precursor waves, as computed in the post-shock frame. At fixed magnetization, we find that $f_ξ$ is nearly independent of temperature as long as $Δγ\lesssim 10^{-1.5}$ (with only a modest decrease of a factor of three from $Δγ=10^{-5}$ to $Δγ=10^{-1.5}$), but it drops by nearly two orders of magnitude for $Δγ\gtrsim 10^{-1}$. At fixed temperature, the scaling with magnetization $f_ξ\sim 10^{-3}\,σ^{-1}$ is consistent with our earlier one-dimensional results. For our reference $σ=1$, the power spectrum of precursor waves is relatively broad (fractional width $\sim 1-3$) for cold temperatures, whereas it shows pronounced line-like features with fractional width $\sim 0.2$ for $10^{-3} \lesssim Δγ\lesssim 10^{-1.5} $. For $σ\gtrsim 1$, the precursor waves are beamed within an angle $\simeq σ^{-1/2}$ from the shock normal (as measured in the post-shock frame), as required so they can outrun the shock. Our results can provide physically-grounded inputs for FRB emission models based on maser emission from relativistic shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tilting Instability of Magnetically Confined Spheromaks

We consider the tilting instability of a magnetically confined spheromak using 3D MHD and relativistic PIC calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter alignment of the spheromak&#39;s intrinsic magnetic dipole with the external magnetic field. Initially the spheromak rotates - tilts - trying to lower its magnetic potential energy. As a result a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on few Alfvén time scales. In the case of higher order (second order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfvén waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetars&#39; giant flares, fast development of tilting instabilities, and no stabilization of the higher order spheromaks, make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Numerical methods for General Relativistic particles

We present recent developments on numerical algorithms for computing photon and particle trajectories in the surrounding of compact objects. Strong gravity around neutron stars or black holes causes relativistic effects on the motion of massive particles and distorts light rays due to gravitational lensing. Efficient numerical methods are required for solving the equations of motion and compute i) the black hole shadow obtained by tracing light rays from the object to a distant observer, and ii) obtain information on the dynamics of the plasma at the microscopic scale. Here, we present generalized algorithms capable of simulating ensembles of photons or massive particles in any spacetime, with the option of including external forces. The coupling of these tools with GRMHD simulations is the key point for obtaining insight on the complex dynamics of accretion disks and jets and for comparing simulations with upcoming observational results from the Event Horizon Telescope.

0 citations0 reviewsNo code linkedOpen access