BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Lorenzo Sironi

Lorenzo Sironi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEphysics.plasm-phastro-ph.GAastro-ph.SRastro-ph.COphysics.space-phastro-ph.IMphysics.comp-ph

Catalog footprint

What is connected

47works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

47 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 "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.

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 ("injection") 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'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'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'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 "plasmoids". 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 "picked up" 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'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' 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
preprint2016arXiv

Blazar flares powered by plasmoids in relativistic reconnection

Powerful flares from blazars with short ($\sim$ min) variability timescales are challenging for current models of blazar emission. Here, we present a physically motivated ab initio model for blazar flares based on the results of recent particle-in-cell (PIC) simulations of relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids filled with high-energy particles and magnetic fields are a self-consistent by-product of the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration profile, particle and magnetic content) with a kinetic equation for the evolution of the electron distribution function we demonstrate that relativistic reconnection in blazar jets can produce powerful flares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron self-Compton flares of duration of several hours--days powered by the largest and slowest moving plasmoids that form in the reconnection layer. Smaller and faster plasmoids produce flares of sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed fluence in both types of flares is similar. Multiple flares with a range of flux-doubling timescales (minutes to several hours) observed over a longer period of flaring activity (days or longer) may be used as a probe of the reconnection layer's orientation and the jet's magnetization. Our model shows that blazar flares are naturally expected as a result of magnetic reconnection in a magnetically-dominated jet.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Radio synchrotron emission from secondary electrons in interaction-powered supernovae

Several supernovae (SNe) with an unusually dense circumstellar medium (CSM) have been recently observed at radio frequencies. Their radio emission is powered by relativistic electrons that can be either accelerated at the SN shock (primaries) or injected as a by-product (secondaries) of inelastic proton-proton collisions. We investigate the radio signatures from secondary electrons, by detailing a semi-analytical model to calculate the temporal evolution of the distributions of protons, primary and secondary electrons. With our formalism, we track the cooling history of all the particles that have been injected into the emission region up to a given time, and calculate the resulting radio spectra and light curves. For a SN shock propagating through the progenitor wind, we find that secondary electrons control the early radio signatures, but their contribution decays faster than that of primary electrons. This results in a flattening of the light curve at a given radio frequency that depends only upon the radial profiles of the CSM density and of the shock velocity, $v_0$. The relevant transition time at the peak frequency is $\sim 190 \, K_{\rm ep,-3}^{-1} A_{\rm w, 16} β_{0, -1.5}^{-2}\,{\rm d}$, where $A_{\rm w}$ is the wind mass-loading parameter, $β_0=v_0/c$ and $K_{\rm ep}$ is the electron-to-proton ratio of accelerated particles. We explicitly show that late peak times at 5 GHz (i.e., $t_{\rm pk} \gtrsim 300-1000$ d) suggest a shock wave propagating in a dense wind ($A_{\rm w} \gtrsim 10^{16}-10^{17}$ gr cm$^{-1}$), where secondary electrons are likely to power the observed peak emission.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks

We formulate a magnetohydrodynamic-particle-in-cell (MHD-PIC) method for describing the interaction between collisionless cosmic ray (CR) particles and a thermal plasma. The thermal plasma is treated as a fluid, obeying equations of ideal MHD, while CRs are treated as relativistic Lagrangian particles subject to the Lorentz force. Backreaction from CRs to the gas is included in the form of momentum and energy feedback. In addition, we include the electromagnetic feedback due to CR-induced Hall effect that becomes important when the electron-ion drift velocity of the background plasma induced by CRs approaches the Alfvén velocity. Our method is applicable on scales much larger than the ion inertial length, bypassing the microscopic scales that must be resolved in conventional PIC methods, while retaining the full kinetic nature of the CRs. We have implemented and tested this method in the Athena MHD code, where the overall scheme is second-order accurate and fully conservative. As a first application, we describe a numerical experiment to study particle acceleration in non-relativistic shocks. Using a simplified prescription for ion injection, we reproduce the shock structure and the CR energy spectra obtained with more self-consistent hybrid-PIC simulations, but at substantially reduced computational cost. We also show that the CR-induced Hall effect reduces the growth rate of the Bell instability and affects the gas dynamics in the vicinity of the shock front. As a step forward, we are able to capture the transition of particle acceleration from non relativistic to relativistic regimes, with momentum spectrum $f(p)\sim p^{-4}$ connecting smoothly through the transition, as expected from the theory of Fermi acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic Jets Shine through Shocks or Magnetic Reconnection?

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative efficiency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are efficient at dissipating energy, they typically do not accelerate particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger than the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating particles, which naturally accounts for a commonly observed property of blazar jets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic Shocks: Particle Acceleration and Magnetization

We review the physics of relativistic shocks, which are often invoked as the sources of non-thermal particles in pulsar wind nebulae (PWNe), gamma-ray bursts (GRBs), and active galactic nuclei (AGN) jets, and as possible sources of ultra-high energy cosmic-rays. We focus on particle acceleration and magnetic field generation, and describe the recent progress in the field driven by theory advances and by the rapid development of particle-in-cell (PIC) simulations. In weakly magnetized or quasi parallel-shocks (where the magnetic field is nearly aligned with the flow), particle acceleration is efficient. The accelerated particles stream ahead of the shock, where they generate strong magnetic waves which in turn scatter the particles back and forth across the shock, mediating their acceleration. In contrast, in strongly magnetized quasi-perpendicular shocks, the efficiencies of both particle acceleration and magnetic field generation are suppressed. Particle acceleration, when efficient, modifies the turbulence around the shock on a long time scale, and the accelerated particles have a characteristic energy spectral index of ~ 2.2 in the ultra-relativistic limit. We discuss how this novel understanding of particle acceleration and magnetic field generation in relativistic shocks can be applied to high-energy astrophysical phenomena, with an emphasis on PWNe and GRB afterglows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electron Heating by the Ion Cyclotron Instability in Collisionless Accretion Flows. I. Compression-Driven Instabilities and the Electron Heating Mechanism

In systems accreting well below the Eddington rate, the plasma in the innermost regions of the disk is collisionless and two-temperature, with the ions hotter than the electrons. Yet, whether a collisionless faster-than-Coulomb energy transfer mechanism exists in two-temperature accretion flows is still an open question. We study the physics of electron heating during the growth of ion velocity-space instabilities, by means of multi-dimensional particle-in-cell (PIC) simulations. A large-scale compression - embedded in a novel form of the PIC equations - continuously amplifies the field. This constantly drives a pressure anisotropy P_perp > P_parallel, due to the adiabatic invariance of the particle magnetic moments. We find that, for ion plasma beta values beta_i ~ 5-30 appropriate for the midplane of low-luminosity accretion flows, mirror modes dominate if the electron-to-proton temperature ratio is > 0.2, whereas if it is < 0.2 the ion cyclotron instability triggers the growth of strong Alfven-like waves, that pitch-angle scatter the ions to maintain marginal stability. We develop an analytical model of electron heating during the growth of the ion cyclotron instability, which we validate with PIC simulations. We find that for cold electrons (beta_e < m_e/m_i), the electron energy gain is controlled by the magnitude of the E-cross-B velocity induced by the ion cyclotron waves. This term is independent of the initial electron temperature, so it provides a solid energy floor even for electrons starting with extremely low temperatures. On the other hand, the electron energy gain for beta_e > m_e/m_i - governed by the conservation of the magnetic moment in the growing fields of the instability - is proportional to the initial electron temperature. Our results have implications for two-temperature accretion flows as well as the solar wind and intracluster plasmas. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electron Heating by the Ion Cyclotron Instability in Collisionless Accretion Flows. II. Electron Heating Efficiency as a Function of Flow Conditions

In the innermost regions of low-luminosity accretion flows, including Sgr A* at the center of our Galaxy, the frequency of Coulomb collisions is so low that the plasma is two-temperature, with the ions substantially hotter than the electrons. This paradigm assumes that Coulomb collisions are the only channel for transferring the ion energy to the electrons. In this work, the second of a series, we assess the efficiency of electron heating by ion velocity-space instabilities in collisionless accretion flows. The instabilities are seeded by the pressure anisotropy induced by magnetic field amplification, coupled to the adiabatic invariance of the particle magnetic moments. Using two-dimensional (2D) particle-in-cell (PIC) simulations, we showed in Paper I that if the electron-to-ion temperature ratio is < 0.2, the ion cyclotron instability is the dominant mode for values of ion beta_i ~ 5-30 (here, beta_i is the ratio of ion thermal pressure to magnetic pressure), as appropriate for the midplane of low-luminosity accretion flows. In this work, we employ analytical theory and 1D PIC simulations (with the box aligned with the fastest growing wavevector of the ion cyclotron mode) to fully characterize how the electron heating efficiency during the growth of the ion cyclotron instability depends on the electron-to-proton temperature ratio, the plasma beta, the Alfven speed, the amplification rate of the mean field (in units of the ion Larmor frequency) and the proton-to-electron mass ratio. Our findings can be incorporated as a physically-grounded sub-grid model into global fluid simulations of low-luminosity accretion flows, thus helping to assess the validity of the two-temperature assumption.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Thermal Electron Acceleration in Low Mach Number Collisionless Shocks. I. Particle Energy Spectra and Acceleration Mechanism

Electron acceleration to non-thermal energies in low Mach number (M<5) shocks is revealed by radio and X-ray observations of galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Diffusive shock acceleration, also known as first-order Fermi acceleration, cannot be directly invoked to explain the acceleration of electrons. Rather, an additional mechanism is required to pre-accelerate the electrons from thermal to supra-thermal energies, so they can then participate in the Fermi process. In this work, we use two- and three-dimensional particle-in-cell plasma simulations to study electron acceleration in low Mach number shocks. We focus on the particle energy spectra and the acceleration mechanism in a reference run with M=3 and a quasi-perpendicular pre-shock magnetic field. We find that about 15 percent of the electrons can be efficiently accelerated, forming a non-thermal power-law tail in the energy spectrum with a slope of p~2.4. Initially, thermal electrons are energized at the shock front via shock drift acceleration. The accelerated electrons are then reflected back upstream, where their interaction with the incoming flow generates magnetic waves. In turn, the waves scatter the electrons propagating upstream back toward the shock, for further energization via shock drift acceleration. In summary, the self-generated waves allow for repeated cycles of shock drift acceleration, similarly to a sustained Fermi-like process. This mechanism offers a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-Thermal Electron Acceleration in Low Mach Number Collisionless Shocks. II. Firehose-Mediated Fermi Acceleration and its Dependence on Pre-Shock Conditions

Electron acceleration to non-thermal energies is known to occur in low Mach number (M<5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M=3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves, that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (>20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic magnetic reconnection in pair plasmas and its astrophysical applications

This review discusses the physics of magnetic reconnection, a process in which the magnetic field topology changes and magnetic energy is converted to kinetic energy, in pair plasmas in the relativistic regime. We focus on recent progress in the field driven by theory advances and the maturity of particle-in-cell codes. This work shows that fragmentation instabilities at the current sheet can play a critical role in setting the reconnection speed and affect the resulting particle acceleration, anisotropy, bulk flows, and radiation. Then, we discuss how this novel understanding of relativistic reconnection can be applied to high-energy astrophysical phenomena, with an emphasis on pulsars, pulsar wind nebulae, and active galactic nucleus jets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relativistic Reconnection: an Efficient Source of Non-Thermal Particles

In magnetized astrophysical outflows, the dissipation of field energy into particle energy via magnetic reconnection is often invoked to explain the observed non-thermal signatures. By means of two- and three-dimensional particle-in-cell simulations, we investigate anti-parallel reconnection in magnetically-dominated electron-positron plasmas. Our simulations extend to unprecedentedly long temporal and spatial scales, so we can capture the asymptotic state of the system beyond the initial transients, and without any artificial limitation by the boundary conditions. At late times, the reconnection layer is organized into a chain of large magnetic islands connected by thin X-lines. The plasmoid instability further fragments each X-line into a series of smaller islands, separated by X-points. At the X-points, the particles become unmagnetized and they get accelerated along the reconnection electric field. We provide definitive evidence that the late-time particle spectrum integrated over the whole reconnection region is a power-law, whose slope is harder than -2 for magnetizations sigma>10. Efficient particle acceleration to non-thermal energies is a generic by-product of the long-term evolution of relativistic reconnection in both two and three dimensions. In three dimensions, the drift-kink mode corrugates the reconnection layer at early times, but the long-term evolution is controlled by the plasmoid instability, that facilitates efficient particle acceleration, in analogy to the two-dimensional physics. Our findings have important implications for the generation of hard photon spectra in pulsar winds and relativistic astrophysical jets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the effect of the Sgr A* accretion flow on the appearance of G2 after pericenter

We study the dynamical interaction of the G2 cloud with the accretion flow around Sgr A* by means of three-dimensional, hydrodynamic simulations. We show the effects of the rotating accretion flow on the evolution of G2 by projecting the cloud density onto the plane of the sky, and extracting position-velocity diagrams. We study a number of possible orientations of the cloud orbit with respect to the disk. We find that once the center of mass of the cloud has crossed the pericenter, the differences between models becomes significant. Models with the cloud counter-rotating with respect to the disk are expected to reach higher blue-shifted line of sight velocities. The spatial extent of the emission depends strongly on the cloud-to-disk inclination angle. Future imaging and spectroscopy of G2 emission will shed light both on the structure of the Sgr A* disk and on the properties of the cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Late-Time Flattening of Light Curves in Gamma-Ray Burst Afterglows

The afterglow emission from Gamma-Ray Bursts (GRBs) is usually interpreted as synchrotron radiation from relativistic electrons accelerated at the GRB external shock, that decelerates from ultra-relativistic to non-relativistic speeds as it sweeps up the surrounding medium. We investigate the temporal decay of the emission from GRB afterglows at late times, when the bulk of the shock-accelerated electrons are non-relativistic. For a uniform circumburst medium, we show that such "deep Newtonian phase" begins at t_{DN} ~ 3 epsilon_{e,-1}^{5/6} t_{ST}, where t_{ST} marks the transition of the blast wave to the non-relativistic spherically-symmetric Sedov-Taylor solution, and epsilon_e = 0.1 epsilon_{e,-1} quantifies the amount of shock energy transferred to the post-shock electrons. For typical parameters, the deep Newtonian stage starts ~ 0.5-several years after the GRB. The radio flux in this phase decays as F_nu ~ t^{-3(p+1)/10}, with decay slope between -0.9 and -1.2, for a power-law distribution of shock-accelerated electrons with index 2 < p < 3. This is shallower than the commonly assumed scaling F_nu ~ t^{-3(5p-7)/10} (with decay slope between -0.9 and -2.4) derived by Frail et al. (2000), which only applies if the GRB shock is non-relativistic, but the electron distribution still peaks at ultra-relativistic energies (a regime that is relevant for a narrow time interval, and only if t_{DN} > t_{ST}, or epsilon_e > 0.03). We discuss how the deep Newtonian phase can be reliably used for GRB calorimetry, and we comment on the good detection prospects of trans-relativistic blast waves at 0.1-10 GHz with EVLA and LOFAR.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Photon Feedback: Screening and the Eddington Limit

Bright star-forming galaxies radiate well below their Eddington Limit. The value of the flux-mean opacity that mediates the radiation force onto matter is orders of magnitude smaller than the UV or optical dust opacity. On empirical grounds, it is shown that high-redshift ULIRGs radiate at two orders of magnitude below their Eddington Limit, while the local starbursters M82 and Arp 220 radiate at a few percent of their Eddington Limit. A simple model for the radiative transfer of UV and optical light in dust-rich environments is considered. Radiation pressure on dust does not greatly affect the large-scale gas dynamics of star-forming galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic Pair Beams from TeV Blazars: A Source of Reprocessed GeV Emission rather than IGM Heating

The interaction of TeV photons from blazars with the extragalactic background light produces a relativistic beam of electron-positron pairs streaming through the intergalactic medium (IGM). The fate of the beam energy is uncertain. By means of two- and three-dimensional particle-in-cell simulations, we study the non-linear evolution of dilute ultra-relativistic pair beams propagating through the IGM. We explore a wide range of beam Lorentz factors gamma_b>>1 and beam-to-plasma density ratios alpha<<1, so that our results can be extrapolated to the extreme parameters of blazar-induced beams (gamma_b~10^6 and alpha~10^(-15), for the most powerful blazars). For cold beams, we show that the oblique instability governs the early stages of evolution, but its exponential growth terminates - due to self-heating of the beam in the transverse direction - when only a negligible fraction ~(alpha/gamma_b)^(1/3)~10^(-7) of the beam energy has been transferred to the IGM plasma. Further relaxation of the beam proceeds through quasi-longitudinal modes, until the momentum dispersion in the direction of propagation saturates at DeltaP_{b,parallel}/gamma_b m_e c ~ 0.2. This corresponds to a fraction ~10% of the beam energy being ultimately transferred to the IGM plasma, irrespective of gamma_b or alpha. If the initial dispersion in beam momentum satisfies DeltaP_{b,parallel}/gamma_b m_e c > 0.2 (as typically expected for blazar-induced beams), the fraction of beam energy deposited into the IGM is much smaller than ~10%. It follows that at least ~90% of the beam energy is still available to power the GeV emission produced by inverse Compton up-scattering of the Cosmic Microwave Background by the beam pairs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Maximum Energy of Accelerated Particles in Relativistic Collisionless Shocks

The afterglow emission from gamma-ray bursts (GRBs) is usually interpreted as synchrotron radiation from electrons accelerated at the GRB external shock, that propagates with relativistic velocities into the magnetized interstellar medium. By means of multi-dimensional particle-in-cell simulations, we investigate the acceleration performance of weakly magnetized relativistic shocks, in the magnetization range 0<sigma<1e-1. The pre-shock magnetic field is orthogonal to the flow, as generically expected for relativistic shocks. We find that relativistic perpendicular shocks propagating in electron-positron plasmas are efficient particle accelerators if the magnetization is sigma<1e-3. For electron-ion plasmas, the transition to efficient acceleration occurs for sigma<3e-5. Here, the acceleration process proceeds similarly for the two species, since the electrons enter the shock nearly in equipartition with the ions, as a result of strong pre-heating in the self-generated upstream turbulence. In both electron-positron and electron-ion shocks, we find that the maximum energy of the accelerated particles scales in time as t^(1/2). This scaling is shallower than the so-called (and commonly assumed) Bohm limit, and it naturally results from the small-scale nature of the Weibel turbulence generated in the shock layer. In magnetized plasmas, the energy of the accelerated particles increases until it reaches a saturation value that scales with the magnetization as sigma^(-1/4). Further energization is prevented by the fact that the self-generated turbulence is confined within a finite region of thickness proportional to sigma^(-1/2) around the shock. Our results can provide physically-grounded inputs for models of non-thermal emission from a variety of astrophysical sources, with particular relevance to GRB afterglows.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The S2 star as a probe of the accretion disk of Sgr A*

How accretion proceeds around the massive black hole in the Galactic center and other highly sub-Eddington accretors remains poorly understood. The orbit of the S2 star in the Galactic center passes through the accretion disk of the massive black hole and any observational signature from such interaction may be used as an accretion probe. Because of its early stellar type, S2 is expected to possess a fairly powerful wind. We show here that the ram pressure of the accretion disk shocks the stellar wind fairly close to the star. The shocked fluid reaches a temperature of ~ 1 keV and cools efficiently through optically thin, thermal bremsstrahlung emission. The radiation from the shocked wind peaks around the epoch of the pericenter passage of the star at a luminosity potentially comparable to the quiescent emission detected from Sgr A*. Detection of shocked wind radiation can constrain the density of the accretion disk at a distance of several thousands of gravitational radii from the black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Particle-in-cell simulations of shock-driven reconnection in relativistic striped winds

By means of two- and three-dimensional particle-in-cell simulations, we investigate the process of driven magnetic reconnection at the termination shock of relativistic striped flows. In pulsar winds and in magnetar-powered relativistic jets, the flow consists of stripes of alternating magnetic field polarity, separated by current sheets of hot plasma. At the wind termination shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength "lambda" or the wind magnetization "sigma" (in the regime sigma>>1 of magnetically-dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor increases by a factor of sigma with respect to the pre-shock value. In the limit lambda/(r_L*sigma)>>1, where r_L is the relativistic Larmor radius in the wind, the post-shock particle spectrum approaches a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. The presence of a current-aligned "guide" magnetic field suppresses the acceleration of particles only when the guide field is stronger than the alternating component. Our findings place important constraints on the models of non-thermal radiation from Pulsar Wind Nebulae and relativistic jets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radio Synchrotron Emission from a Bow Shock around the Gas Cloud G2 Heading toward the Galactic Center

A dense ionized cloud of gas has been recently discovered to be moving directly toward the supermassive black hole, Sgr A*, at the Galactic Center. In June 2013, at the pericenter of its highly eccentric orbit, the cloud will be approximately 3100 Schwarzschild radii from the black hole and will move supersonically through the ambient hot gas with a velocity of v_p ~ 5400 km/s. A bow shock is likely to form in front of the cloud and could accelerate electrons to relativistic energies. We estimate via particle-in-cell simulations the energy distribution of the accelerated electrons and show that the non-thermal synchrotron emission from these electrons might exceed the quiescent radio emission from Sgr A* by a factor of several. The enhanced radio emission should be detectable at GHz and higher frequencies around the time of pericentric passage and in the following months. The bow shock emission is expected to be displaced from the quiescent radio emission of Sgr A* by ~33 mas. Interferometric observations could resolve potential changes in the radio image of Sgr A* at wavelengths < 6 cm.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Acceleration of Particles at the Termination Shock of a Relativistic Striped Wind

The relativistic wind of obliquely-rotating pulsars consists of toroidal stripes of opposite magnetic field polarity, separated by current sheets of hot plasma. By means of two- and three-dimensional particle-in-cell simulations, we investigate particle acceleration and magnetic field dissipation at the termination shock of a relativistic striped wind. At the shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength "lambda" or the wind magnetization "sigma" (in the regime sigma>>1 of magnetically-dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor increases by a factor of sigma with respect to the pre-shock value. The shape of the post-shock spectrum depends primarily on the ratio lambda/(r_L*sigma), where r_L is the relativistic Larmor radius in the wind. The spectrum becomes broader as the value of lambda/(r_L*sigma) increases, passing from a relativistic Maxwellian to a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. Close to the equatorial plane of the wind, where the stripes are symmetric, the highest energy particles resulting from magnetic reconnection can escape ahead of the shock, and be injected into a Fermi-like acceleration process. In the post-shock spectrum, they populate a power-law tail with slope around -2.5, that extends beyond the flat component produced by reconnection. Our study suggests that the spectral break between the radio and the optical band in Pulsar Wind Nebulae can be a natural consequence of particle acceleration at the termination shock of striped pulsar winds.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Particle Acceleration in Relativistic Magnetized Collisionless Electron-Ion Shocks

We investigate shock structure and particle acceleration in relativistic magnetized collisionless electron-ion shocks by means of 2.5D particle-in-cell simulations with ion-to-electron mass ratios (m_i/m_e) ranging from 16 to 1000. We explore a range of inclination angles between the pre-shock magnetic field and the shock normal. In "subluminal" shocks, where relativistic particles can escape ahead of the shock along the magnetic field lines, ions are efficiently accelerated via a Fermi-like mechanism. The downstream ion spectrum consists of a relativistic Maxwellian and a high-energy power-law tail, which contains ~5% of ions and ~30% of ion energy. Its slope is -2.1. Upstream electrons enter the shock with lower energy than ions, so they are more strongly tied to the field. As a result, only ~1% of the incoming electrons are Fermi-accelerated at the shock before being advected downstream, where they populate a steep power-law tail (with slope -3.5). For "superluminal" shocks, where relativistic particles cannot outrun the shock along the field, the self-generated turbulence is not strong enough to permit efficient Fermi acceleration, and the ion and electron downstream spectra are consistent with thermal distributions. The incoming electrons are heated up to equipartition with ions, due to strong electromagnetic waves emitted by the shock into the upstream. Thus, efficient electron heating (>15% of the upstream ion energy) is the universal property of relativistic electron-ion shocks, but significant nonthermal acceleration of electrons (>2% by number, >10% by energy, with slope flatter than -2.5) is hard to achieve in magnetized flows and requires weakly magnetized shocks (magnetization <1e-3). These findings place important constraints on the models of AGN jets and Gamma Ray Bursts that invoke particle acceleration in relativistic magnetized electron-ion shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Particle Acceleration in Relativistic Magnetized Collisionless Pair Shocks: Dependence of Shock Acceleration on Magnetic Obliquity

We investigate shock structure and particle acceleration in relativistic magnetized collisionless pair shocks by means of 2.5D and 3D particle-in-cell simulations. We explore a range of inclination angles between the pre-shock magnetic field and the shock normal. We find that only magnetic inclinations corresponding to "subluminal" shocks, where relativistic particles following the magnetic field can escape ahead of the shock, lead to particle acceleration. The downstream spectrum in such shocks consists of a relativistic Maxwellian and a high-energy power-law tail with exponential cutoff. For increasing magnetic inclination in the subluminal range, the high-energy tail accounts for an increasing fraction of particles (from ~1% to ~2%) and energy (from ~4% to ~12%). The spectral index of the power law increases with angle from -2.8+-0.1 to -2.3+-0.1. Particle energization is driven by the Diffusive Shock Acceleration process for nearly parallel shocks, and switches to Shock-Drift Acceleration for larger subluminal inclinations. For "superluminal" shocks, the downstream particle spectrum does not show any significant suprathermal tail. As seen from the upstream frame, efficient acceleration in relativistic (Lorentz factor gamma0 > 5) magnetized (sigma > 0.03) flows exists only for a very small range of magnetic inclination angles (< 34/gamma0 degrees), so relativistic astrophysical pair shocks have to be either nearly parallel or weakly magnetized to generate nonthermal particles. These findings place constraints on the models of AGN jets, Pulsar Wind Nebulae and Gamma Ray Bursts that invoke particle acceleration in relativistic magnetized shocks. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Polarized Infrared Emission by Polycyclic Aromatic Hydrocarbons resulting from Anisotropic Illumination

We study the polarized infrared emission by Polycyclic Aromatic Hydrocarbons (PAHs), when anisotropically illuminated by UV photons. PAH molecules are modeled as planar disks with in-plane and out-of-plane vibrational dipoles. As first pointed out by Leger (1988), infrared emission features resulting from in-plane and out-of-plane modes should have orthogonal polarization directions. We show analytically how the degree of polarization depends on the viewing geometry and the molecule's internal alignment between principal axis of inertia and angular momentum, which gets worse after photon absorption. Longer wavelength features, emitted after better internal alignment is recovered, should be more strongly polarized. The degree of polarization for uni-directional illumination (e.g., by a star) is larger than for diffuse illumination (e.g., by a disk galaxy), all else being equal. For PAHs in the Cold Neutral Medium, the predicted polarization is probably too small to distinguish from the contribution of linear dichroism by aligned foreground dust. The level of polarization predicted for PAH emission from the Orion Bar is only ~0.06% at 3.3 microns; Sellgren et al. (1988) report a much larger value, 0.86+-0.28%, which suggests that the smallest PAHs may have moderately suprathermal rotation rates. Future observations of (or upper limits on) the degree of polarization for the Orion Bar or for dust above edge-on galaxies (e.g., NGC 891 or M82) may constrain the internal alignment of emitting PAHs, thus providing clues to their rotational dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Synthetic Spectra from PIC Simulations of Relativistic Collisionless Shocks

We extract synthetic photon spectra from first-principles particle-in-cell simulations of relativistic shocks propagating in unmagnetized pair plasmas. The two basic ingredients for the radiation, namely accelerated particles and magnetic fields, are produced self-consistently as part of the shock evolution. We use the method of Hededal & Nordlund (2005) and compute the photon spectrum via Fourier transform of the electric far-field from a large number of particles, sampled directly from the simulation. We find that the spectrum from relativistic collisionless shocks is entirely consistent with synchrotron radiation in the magnetic fields generated by Weibel instability. We can recover the so-called "jitter'' regime only if we artificially reduce the strength of the electromagnetic fields, such that the wiggler parameter K = qB lambda/mc^2 becomes much smaller than unity ("B" and "lambda" are the strength and scale of the magnetic turbulence, respectively). These findings may place constraints on the origin of non-thermal emission in astrophysics, especially for the interpretation of the hard (harder than synchrotron) low-frequency spectrum of Gamma-Ray Bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Eddington Limit in Cosmic Rays: An Explanation for the Observed Lack of Low-Mass Radio-Loud Quasars and the M_{BH}-M_{Bulge} Relation

We present a feedback mechanism for supermassive black holes and their host bulges that operates during epochs of radio-loud quasar activity. In the radio cores of relativistic quasar jets, internal shocks convert a fraction of ordered bulk kinetic energy into randomized relativistic ions, or in other words cosmic rays. By employing a phenomenologically-motivated jet model, we show that enough 1-100 GeV cosmic rays escape the radio core into the host galaxy to break the Eddington limit in cosmic rays. As a result, hydrostatic balance is lost and a cosmic ray momentum-driven wind develops, expelling gas from the host galaxy and thus self-limiting the black hole and bulge growth. Although the interstellar cosmic ray power is much smaller than the quasar photon luminosity, cosmic rays provide a stronger feedback than UV photons, since they exchange momentum with the galactic gas much more efficiently. The amount of energy released into the host galaxy as cosmic rays, per unit of black hole rest mass energy, is independent of black hole mass. It follows that radio-loud jets should be more prevalent in relatively massive systems since they sit in galaxies with relatively deep gravitational potentials. Therefore, jet-powered cosmic ray feedback not only self-regulates the black hole and bulge growth, but also provides an explanation for the lack of radio-loud activity in relatively small galaxies. By employing basic known facts regarding the physical conditions in radio cores, we approximately reproduce both the slope and the normalization of the M_{BH}-M_{Bulge} relation.

0 citations0 reviewsNo code linkedOpen access