Sites Of Particle Acceleration Research Articles

Future X-ray Polarimetry of Relativistic Accelerators: Pulsar Wind Nebulae and Supernova Remnants

Supernova remnants (SNRs) and pulsar wind nebulae (PWNs) are among the most significant sources of non-thermal X-rays in the sky, and the best means by which relativistic plasma dynamics and particle acceleration can be investigated. Being strong synchrotron emitters, they are ideal candidates for X-ray polarimetry, and indeed the Crab nebula is up to present the only object where X-ray polarization has been detected with a high level of significance. Future polarimetric measures will likely provide us with crucial information on the level of turbulence that is expected at particle acceleration sites, together with the spatial and temporal coherence of magnetic field geometry, enabling us to set stronger constraints on our acceleration models. PWNs will also allow us to estimate the level of internal dissipation. I will briefly review the current knowledge on the polarization signatures in SNRs and PWNs, and I will illustrate what we can hope to achieve with future missions such as IXPE/XIPE.

Ultrahigh-energy cosmic rays from tidally-ignited white dwarfs

Ultra-high-energy cosmic rays (UHECRs) can be accelerated by tidal disruption events of stars by black holes. We suggest a novel mechanism for UHECR acceleration wherein white dwarfs (WDs) are tidally compressed by intermediate-mass black holes (IMBHs), leading to their ignition and subsequent explosion as a supernova. Cosmic rays accelerated by the supernova may receive an energy boost when crossing the accretion-powered jet. The rate of encounters between WDs and IMBHs can be relatively high, as the number of IMBHs may be substantially augmented once account is taken of their likely presence in dwarf galaxies. Here we show that this kind of tidal disruption event naturally provides an intermediate composition for the observed UHECRs, and suggest that dwarf galaxies and globular clusters are suitable sites for particle acceleration to ultra-high energies.

Pulsar Emission Physics: The First Fifty Years

AbstractOver the last fifty years since the discovery of pulsars, our understanding of where and how pulsars emit the radiation we observe has undergone significant revision. The location and mechanisms of high-energy radiation are intimately tied to the sites of particle acceleration. The evolution of emission models has paralleled the development of increasingly more sensitive telescopes, especially at high energies. I will review the history of pulsar emission modeling, from the early days of gaps at the polar caps, to outer gaps and slot gaps in the outer magnetosphere, to the present era of global magnetosphere simulations that locate most acceleration and high-energy emission in the current sheets.

Energetic ion leakage from foreshock transient cores

AbstractEarth's foreshock is filled with backstreaming particles that can interact with the ambient solar wind and its discontinuities to form foreshock transients. Many foreshock transients have a core with low dynamic pressure that can significantly perturb the bow shock and the magnetosphere‐ionosphere system. Foreshock transients have also been recently recognized as sites of particle acceleration, which may be important for seeding the parent shock with energetic particles. A relevant step of this seeding would be energetic ion leakage into the surrounding foreshock environment. On the other hand, such leakage would also suppress the energetic particle flux contrast across foreshock transients' boundaries masking their perceived contribution to ion energization. To further examine this hypothesis of ion leakage, we report on multipoint case studies of three foreshock transient events selected from a large database. The cases were selected to exemplify, in increasing complexity, the nature and consequences of energetic ion leakage. Ion energy dispersion, observed upstream and/or downstream of the foreshock transients, is explained with a simple, ballistic model of ions leaking from the foreshock transients. Larger energies are required for leaked ions to reach the spacecraft as the distance between the transient and spacecraft increases. Our model, which explains well the observed ion energy dispersion and velocity distributions, can also be used to reveal the shape of the foreshock transients in three dimensions. Our results suggest that ion leakage from foreshock transient cores needs to be accounted for both in statistical studies and in global models of ion acceleration under quasi‐parallel foreshock conditions.

General relativistic study of astrophysical jets with internal shocks

We explore the possibility of formation of steady internal shocks in jets around black holes. We consider a fluid described by a relativistic equation of state, flowing about the axis of symmetry ($\theta=0$) in a Schwarzschild metric. We use two models for the jet geometry, (i) a conical geometry and (ii) a geometry with non-conical cross-section. Jet with conical geometry is smooth flow. While the jet with non-conical cross section undergoes multiple sonic point and even standing shock. The jet shock becomes stronger, as the shock location is situated further from the central black hole. Jets with very high energy and very low energy do not harbour shocks, but jets with intermediate energies do harbour shocks. One advantage of these shocks, as opposed to shocks mediated by external medium is that, these shocks have no effect on the jet terminal speed, but may act as possible sites for particle acceleration. Typically, a jet with energy $1.8~c^2$, will achieve a terminal speed of $v_\infty=0.813c$ for jet with any geometry. But for a jet of non-conical cross-section for which the length scale of the inner torus of the accretion disc is $40\rg$, then in addition, a steady shock will form at $\rsh \sim 7.5\rg$ and compression ratio of $R\sim 2.7$. Moreover, electron-proton jet seems to harbour the strongest shock. We discuss possible consequences of such a scenario.

Particle Acceleration Due to Coronal Non-null Magnetic Reconnection

Various topological features, for example magnetic null points and separators, have been inferred as likely sites of magnetic reconnection and particle acceleration in the solar atmosphere. In fact, magnetic reconnection is not constrained to solely take place at or near such topological features and may also take place in the absence of such features. Studies of particle acceleration using non-topological reconnection experiments embedded in the solar atmosphere are uncommon. We aim to investigate and characterise particle behaviour in a model of magnetic reconnection which causes an arcade of solar coronal magnetic field to twist and form an erupting flux rope, crucially in the absence of any common topological features where reconnection is often thought to occur. We use a numerical scheme that evolves the gyro-averaged orbit equations of single electrons and protons in time and space, and simulate the gyromotion of particles in a fully analytical global field model. We observe and discuss how the magnetic and electric fields of the model and the initial conditions of each orbit may lead to acceleration of protons and electrons up to 2 MeV in energy (depending on model parameters). We describe the morphology of time-dependent acceleration and impact sites for each particle species and compare our findings to those recovered by topologically based studies of three-dimensional (3D) reconnection and particle acceleration. We also broadly compare aspects of our findings to general observational features typically seen during two-ribbon flare events.

RADIO DIAGNOSTICS OF ELECTRON ACCELERATION SITES DURING THE ERUPTION OF A FLUX ROPE IN THE SOLAR CORONA

ABSTRACT Electron acceleration in the solar corona is often associated with flares and the eruption of twisted magnetic structures known as flux ropes. However, the locations and mechanisms of such particle acceleration during the flare and eruption are still subject to much investigation. Observing the exact sites of particle acceleration can help confirm how the flare and eruption are initiated and how they evolve. Here we use the Atmospheric Imaging Assembly to analyze a flare and erupting flux rope on 2014 April 18, while observations from the Nançay Radio Astronomy Facility allow us to diagnose the sites of electron acceleration during the eruption. Our analysis shows evidence of a pre-formed flux rope that slowly rises and becomes destabilized at the time of a C-class flare, plasma jet, and the escape of ≳75 keV electrons from the rope center into the corona. As the eruption proceeds, continued acceleration of electrons with energies of ∼5 keV occurs above the flux rope for a period over 5 minutes. At the flare peak, one site of electron acceleration is located close to the flare site, while another is driven by the erupting flux rope into the corona at speeds of up to 400 km s−1. Energetic electrons then fill the erupting volume, eventually allowing the flux rope legs to be clearly imaged from radio sources at 150–445 MHz. Following the analysis of Joshi et al. (2015), we conclude that the sites of energetic electrons are consistent with flux rope eruption via a tether cutting or flux cancellation scenario inside a magnetic fan-spine structure. In total, our radio observations allow us to better understand the evolution of a flux rope eruption and its associated electron acceleration sites, from eruption initiation to propagation into the corona.

SHEAR ACCELERATION IN EXPANDING FLOWS

ABSTRACT Shear flows are naturally expected to occur in astrophysical environments and potential sites of continuous non-thermal Fermi-type particle acceleration. Here we investigate the efficiency of expanding relativistic outflows to facilitate the acceleration of energetic charged particles to higher energies. To this end, the gradual shear acceleration coefficient is derived based on an analytical treatment. The results are applied to the context of the relativistic jets from active galactic nuclei. The inferred acceleration timescale is investigated for a variety of conical flow profiles (i.e., power law, Gaussian, Fermi–Dirac) and compared to the relevant radiative and non-radiative loss timescales. The results exemplify that relativistic shear flows are capable of boosting cosmic-rays to extreme energies. Efficient electron acceleration, on the other hand, requires weak magnetic fields and may thus be accompanied by a delayed onset of particle energization and affect the overall jet appearance (e.g., core, ridge line, and limb-brightening).

Hard X-ray morphology of the X1.3 April 25, 2014 partially occulted limb solar flare

At hard X-ray energies, the bright footpoint emission from solar flare loops often prevents a detailed analysis of the weaker loop-top source morphology due to the limited dynamic range available for X-ray imaging. Here, we study the X1.3 April 25, 2014 flare with the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). This partially occulted limb flare allows the analysis of the loop-top emission in isolation. We present results on the flare light curve at different energies, the source morphology from X-ray imaging and a detailed spectral analysis of the different source components by imaging spectroscopy. The loop-top source, a likely site of particle acceleration, shows a clear composition of different emission components. The results indicate the opportunities that detailed imaging of hard X-rays can provide to learn about particle acceleration, transport and heating processes in solar flares.

Multi-D magnetohydrodynamic modelling of pulsar wind nebulae: recent progress and open questions

In the last decade, the relativistic magnetohydrodynamic (MHD) modelling of pulsar wind nebulae, and of the Crab nebula in particular, has been highly successful, with many of the observed dynamical and emission properties reproduced down to the finest detail. Here, we critically discuss the results of some of the most recent studies: namely the investigation of the origin of the radio emitting particles and the quest for the acceleration sites of particles of different energies along the termination shock, by using wisp motions as a diagnostic tool; the study of the magnetic dissipation process in high magnetization nebulae by means of new long-term three-dimensional simulations of the pulsar wind nebula evolution; the investigation of the relativistic tearing instability in thinning current sheets, leading to fast reconnection events that might be at the origin of the Crab nebula gamma-ray flares.

CRITICAL DIFFERENCES OF ASYMMETRIC MAGNETIC RECONNECTION FROM STANDARD MODELS

ABSTRACT We have clarified the structure of asymmetric magnetic reconnection in detail as the result of the spontaneous evolutionary process. The asymmetry is imposed as ratio k of the magnetic field strength in both sides of the initial current sheet (CS) in the isothermal equilibrium. The MHD simulation is carried out by the HLLD code for the long-term temporal evolution with very high spatial resolution. The resultant structure is drastically different from the symmetric case (e.g., the Petschek model) even for slight asymmetry k = 2. (1) The velocity distribution in the reconnection jet clearly shows a two-layered structure, i.e., the high-speed sub-layer in which the flow is almost field aligned and the acceleration sub-layer. (2) Higher beta side (HBS) plasma is caught in a lower beta side plasmoid. This suggests a new plasma mixing process in the reconnection events. (3) A new large strong fast shock in front of the plasmoid forms in the HBS. This can be a new particle acceleration site in the reconnection system. These critical properties that have not been reported in previous works suggest that we contribute to a better and more detailed knowledge of the reconnection of the standard model for the symmetric magnetic reconnection system.

STUDYING THE SGR 1806-20/Cl* 1806-20 REGION USING THE FERMI LARGE AREA TELESCOPE

ABSTRACT The region around SGR 1806-20 and its host stellar cluster Cl* 1806-20 is a potentially important site of particle acceleration. The soft γ-ray repeater and Cl* 1806-20, which also contains several very massive stars including a luminous blue variable hypergiant LBV 1806-20, are capable of depositing a large amount of energy to the surroundings. Using the data taken with the Fermi Large Area Telescope (LAT), we identified an extended LAT source to the southwest of Cl* 1806-20. The centroid of the 1–50 GeV emission is consistent with that of HESS J1808-204 (until now unidentified). The LAT spectrum is best-fit by a broken power law with the break energy . The index above E b is 2.60 ± 0.04 and is consistent with the flux and spectral index above 100 GeV for HESS J1808-204, suggesting an association between the two sources. Meanwhile, the interacting supernova remnant SNR G9.7-0.0 is also a potential contributor to the LAT flux. A tentative flux enhancement at the MeV band during a 45 day interval (2011 January 21–March 7) is also reported. We discuss possible origins of the extended LAT source in the context of both leptonic and hadronic scenarios.

Spectral lags of flaring events in LS I+61°303 from RXTE Observations

This work reports the first discovery of (negative) spectral lags in X-ray emission below 10keV from the gamma ray binary LS I+61° 303 during large flaring episodes using Rossi X-ray Timing Explorer (RXTE) observations. It is found from the RXTE data that during the flares, low energy (3–5 keV) variations lead the higher energy (8–10 keV) variations by a few tens of seconds whereas no significant time lag is observed during the non-flaring states. The observed spectral lag features for flaring events suggest that inverse Compton scattering may be operating, at least in some part of the system. Another possibility is that the sites of particle acceleration may be different for flaring and non-flaring events such as in the microquasar model, in which the flaring radiation may come from hot spots sitting above the black hole while steady state emissions are due to the jets.

Modeling nonthermal emission from stellar bow shocks

Runaway O- and early B-type stars passing throughout the interstellar medium at supersonic velocities and characterized by strong stellar winds may produce bow shocks that can serve as particle acceleration sites. Previous theoretical models predict the production of high energy photons by non-thermal radiative processes, but their efficiency is still debated. We aim to test and explain the possibility of emission from the bow shocks formed by runaway stars traveling through the interstellar medium by using previous theoretical models. We apply our model to AE Aurigae, the first reported star with an X-ray detected bow shock, to BD+43 3654, in which the observations failed in detecting high energy emission, and to the transition phase of a supergiant star in the late stages of its life.From our analysis, we confirm that the X-ray emission from the bow shock produced by AE Aurigae can be explained by inverse Compton processes involving the infrared photons of the heated dust. We also predict low high energy flux emission from the bow shock produced by BD+43 3654, and the possibility of high energy emission from the bow shock formed by a supergiant star during the transition phase from blue to red supergiant.Bow shock formed by different type of runaway stars are revealed as a new possible source of high energy photons in our neighbourhood.

ON THE DISTRIBUTION OF PARTICLE ACCELERATION SITES IN PLASMOID-DOMINATED RELATIVISTIC MAGNETIC RECONNECTION

We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection.

The inner knot of the Crab nebula

We model the inner knot of the Crab Nebula as a synchrotron emission coming from the non-spherical MHD termination shock of relativistic pulsar wind. The post-shock flow is mildly relativistic; as a result the Doppler-beaming has a strong impact on the shock appearance. The model can reproduce the knot location, size, elongation, brightness distribution, luminosity and polarization provided the effective magnetization of the section of the pulsar wind producing the knot is low, $\sigma \leq 1$. In the striped wind model, this implies that the striped zone is rather wide, with the magnetic inclination angle of the Crab pulsar $\ge 45^\circ$; this agrees with the previous model-dependent estimate based on the gamma-ray emission of the pulsar. We conclude that the tiny knot is indeed a bright spot on the surface of a quasi-stationary magnetic relativistic shock and that this shock is a site of efficient particle acceleration. On the other hand, the deduced low magnetization of the knot plasma implies that this is an unlikely site for the Crab's gamma-ray flares, if they are related to the fast relativistic magnetic reconnection events.

HORIZON-SCALE LEPTON ACCELERATION IN JETS: EXPLAINING THE COMPACT RADIO EMISSION IN M87

It has now become clear that the radio jet in the giant elliptical galaxy M87 must turn on very close to the black hole. This implies the efficient acceleration of leptons within the jet at scales much smaller than feasible by the typical dissipative events usually invoked to explain jet synchrotron emission. Here we show that the stagnation surface, the separatrix between material that falls back into the black hole and material that is accelerated outward forming the jet, is a natural site of pair formation and particle acceleration. This occurs via an inverse-Compton pair catastrophe driven by unscreened electric fields within the charge-starved region about the stagnation surface and substantially amplified by a post-gap cascade. For typical estimates of the jet properties in M87, we find excellent quantitive agreement between the predicted relativistic lepton densities and those required by recent high-frequency radio observations of M87. This mechanism fails to adequately fill a putative jet from Sagittarius A* with relativistic leptons, which may explain the lack of an obvious radio jet in the Galactic center. Finally, this process implies a relationship between the kinetic jet power and the gamma-ray luminosity of blazars, produced during the post-gap cascade.

Constraints on particle acceleration sites in the Crab nebula from relativistic magnetohydrodynamic simulations

The Crab Nebula is one of the most efficient accelerators in the Galaxy and the only galactic source showing direct evidence of PeV particles. In spite of this, the physical process behind such effective acceleration is still a deep mystery. While particle acceleration, at least at the highest energies, is commonly thought to occur at the pulsar wind termination shock, the properties of the upstream flow are thought to be non-uniform along the shock surface, and important constraints on the mechanism at work come from exact knowledge of where along this surface particles are being accelerated. Here we use axisymmetric relativistic MHD simulations to obtain constraints on the acceleration site(s) of particles of different energies in the Crab Nebula. Various scenarios are considered for the injection of particles responsible for synchrotron radiation in the different frequency bands, radio, optical and X-rays. The resulting emission properties are compared with available data on the multi wavelength time variability of the inner nebula. Our main result is that the X-ray emitting particles are accelerated in the equatorial region of the pulsar wind. Possible implications on the nature of the acceleration mechanism are discussed.

Neutrinos from colliding wind binaries: future prospects for PINGU and ORCA

Abstract. Massive stars play an important role in explaining the cosmic ray spectrum below the knee, possibly even up to the ankle, i.e. up to energies of 1015 or 1018.5 eV, respectively. In particular, Supernova Remnants are discussed as one of the main candidates to explain the cosmic ray spectrum. Even before their violent deaths, during the stars' regular life times, cosmic rays can be accelerated in wind environments. High-energy gamma-ray measurements indicate hadronic acceleration binary systems, leading to both periodic gamma-ray emission from binaries like LSI + 60 303 and continuous emission from colliding wind environments like η-Carinae. The detection of neutrinos and photons from hadronic interactions are one of the most promising methods to identify particle acceleration sites. In this paper, future prospects to detect neutrinos from colliding wind environments in massive stars are investigated. In particular, the seven most promising candidates for emission from colliding wind binaries are investigated to provide an estimate of the signal strength. The expected signal of a single source is about a factor of 5–10 below the current IceCube sensitivity and it is therefore not accessible at the moment. What is discussed in addition is future the possibility to measure low-energy neutrino sources with detectors like PINGU and ORCA: the minimum of the atmospheric neutrino flux at around 25 GeV from neutrino oscillations provides an opportunity to reduce the background and increase the significance to searches for GeV–TeV neutrino sources. This paper presents the first idea, detailed studies including the detector's effective areas will be necessary in the future to test the feasibility of such an approach.

Electron-scale shear instabilities: magnetic field generation and particle acceleration in astrophysical jets

Strong shear flow regions found in astrophysical jets are shown to be important dissipation regions, where the shear flow kinetic energy flow is converted into electric and magnetic field energy via shear instabilities. The emergence of these self-consistent fields makes shear flows significant sites for radiation emission and particle acceleration. We focus on electron-scale instabilities, namely the collisionless, unmagnetized electron-scale Kelvin–Helmholtz instability (ESKHI) and a large-scale DC magnetic field generation mechanism on the electron scales. We show that these processes are important candidates to generate magnetic fields in the presence of strong velocity shears, which may naturally originate in energetic matter outbursts of active galactic nuclei and gamma-ray bursters. We show that the ESKHI is robust to density jumps between shearing flows, thus operating in various scenarios with different density contrasts. Multidimensional particle-in-cell (PIC) simulations of the ESKHI, performed with OSIRIS, reveal the emergence of a strong and large-scale DC magnetic field component, which is not captured by the standard linear fluid theory. This DC component arises from kinetic effects associated with the thermal expansion of electrons of one flow into the other across the shear layer, whilst ions remain unperturbed due to their inertia. The electron expansion forms DC current sheets, which induce a DC magnetic field. Our results indicate that most of the electromagnetic energy developed in the ESKHI is stored in the DC component, reaching values of equipartition on the order of in the electron time-scale, and persists longer than the proton time-scale. Particle scattering/acceleration in the self-generated fields of these shear flow instabilities is also analyzed.