-
The Eel Pulsar Wind Nebula: a PeVatron-Candidate Origin for HAWC J1826-128 and HESS J1826-130
Authors:
Daniel A. Burgess,
Kaya Mori,
Joseph D. Gelfand,
Charles J. Hailey,
Yarone M. Tokayer,
Jooyun Woo,
Hongjun An,
Kelly Malone,
Stephen P. Reynolds,
Samar Safi-Harb,
Tea Temim
Abstract:
HAWC J1826-128 is one of the brightest Galactic TeV gamma-ray sources detected by the High Altitude Water Cherenkov (HAWC) Observatory, with photon energies extending up to nearly $\sim$100 TeV. This HAWC source spatially coincides with the H.E.S.S. TeV source HESS J1826-130 and the "Eel" pulsar wind nebula (PWN), which is associated with the GeV pulsar PSR J1826-1256. In the X-ray band, Chandra a…
▽ More
HAWC J1826-128 is one of the brightest Galactic TeV gamma-ray sources detected by the High Altitude Water Cherenkov (HAWC) Observatory, with photon energies extending up to nearly $\sim$100 TeV. This HAWC source spatially coincides with the H.E.S.S. TeV source HESS J1826-130 and the "Eel" pulsar wind nebula (PWN), which is associated with the GeV pulsar PSR J1826-1256. In the X-ray band, Chandra and XMM-Newton revealed that the Eel PWN is composed of both a compact nebula ($\sim$15") and diffuse X-ray emission ($\sim$6'$\times$2') extending away from the pulsar. Our NuSTAR observation detected hard X-ray emission from the compact PWN up to $\sim$20 keV and evidence of the synchrotron burn-off effect. In addition to the spatial coincidence between HESS J1826-130 and the diffuse X-ray PWN, our multi-wavelength spectral energy distribution (SED) analysis using X-ray and gamma-ray data establishes a leptonic origin of the TeV emission associated with the Eel PWN. Furthermore, our evolutionary PWN SED model suggests (1) a low PWN B-field of $\sim$1 $μ$G, (2) a significantly younger pulsar age ($t \sim5.7$ kyr) than the characteristic age ($τ= 14.4$ kyr) and (3) a maximum electron energy of $E_{max} = 2$ PeV. The low B-field as well as the putative supersonic motion of the pulsar may account for the asymmetric morphology of the diffuse X-ray emission. Our results suggest that the Eel PWN may be a leptonic PeVatron particle accelerator powered by the $\sim$6-kyr-old pulsar PSR J1826-1256 with a spin-down power of $3.6 \times 10^{36}$ erg s$^{-1}$.
△ Less
Submitted 6 April, 2022;
originally announced April 2022.
-
ULF wave transmission across collisionless shocks: 2.5D local hybrid simulations
Authors:
Primoz Kajdic,
Yann Pfau-Kempf,
Lucile Turc,
Andrew P Dimmock,
Minna Palmroth,
Kazue Takahashi,
Eemilia Kilpua,
Jan Soucek,
Naoko Takahashi,
Luis Preisser,
Xochitl Blanco-Cano,
Domenico Trotta,
David Burgess
Abstract:
We study the interaction of upstream ultra-low frequency (ULF) waves with collisionless shocks by analyzing the outputs of eleven 2D local hybrid simulation runs. Our simulated shocks have Alfvénic Mach numbers between 4.29-7.42 and their $θ_{BN}$ angles are 15$^\circ$, 30$^\circ$, 45$^\circ$ and 50$^\circ$. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude…
▽ More
We study the interaction of upstream ultra-low frequency (ULF) waves with collisionless shocks by analyzing the outputs of eleven 2D local hybrid simulation runs. Our simulated shocks have Alfvénic Mach numbers between 4.29-7.42 and their $θ_{BN}$ angles are 15$^\circ$, 30$^\circ$, 45$^\circ$ and 50$^\circ$. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude of the upstream waves exhibit a complex dependence on the shock's M$_A$ and $θ_{BN}$. The wavelength positively correlates with both parameters, with the dependence on $θ_{BN}$ being much stronger. The amplitude of the ULF waves is proportional to the product of the reflected beam velocity and density, which also depend on M$_A$ and $θ_{BN}$. The interaction of the ULF waves with the shock causes large-scale (several tens of upstream ion inertial lengths) shock rippling. The properties of the shock ripples are related to the ULF wave properties, namely thier wavelength and amplitude. In turn, the ripples have a large impact on the ULF wave transmission across the shock because they change local shock properties ($θ_{BN}$, strength), so that different sections of the same ULF wave front encounter shock with different characteristics. Downstream fluctuations do not resemble the upstream waves in terms the wavefront extension, orientation or their wavelength. However some features are conserved in the Fourier spectra of downstream compressive waves that present a bump or flattening at wavelengths approximately corresponding to those of the upstream ULF waves. In the transverse downstream spectra these features are weaker.
△ Less
Submitted 25 January, 2022;
originally announced January 2022.
-
NuSTAR broad-band X-ray observational campaign of energetic pulsar wind nebulae in synergy with VERITAS, HAWC and Fermi gamma-ray telescopes
Authors:
Kaya Mori,
Hongjun An,
Daniel Burgess,
Massimo Capasso,
Brenda Dingus,
Joseph Gelfand,
Charles Hailey,
Brian Humensky,
Kelly Malone,
Reshmi Mukherjee,
Nahee Park,
Isaac Pope,
Stephen Reynolds,
Samar Safi-Harb,
Jooyun Woo,
Galactic TeV source collaboration
Abstract:
We report recent progress on the on-going NuSTAR observational campaign of 8 TeV-detected pulsar wind nebulae (PWNe). This campaign constitutes a major part of our NuSTAR study of some of the most energetic TeV sources in our Galaxy detected by VERITAS and HAWC. NuSTAR is the first focusing X-ray telescope operating above 10 keV in space with sub-arcminute angular resolution. Broad-band X-ray imag…
▽ More
We report recent progress on the on-going NuSTAR observational campaign of 8 TeV-detected pulsar wind nebulae (PWNe). This campaign constitutes a major part of our NuSTAR study of some of the most energetic TeV sources in our Galaxy detected by VERITAS and HAWC. NuSTAR is the first focusing X-ray telescope operating above 10 keV in space with sub-arcminute angular resolution. Broad-band X-ray imaging and spectroscopy data, obtained by NuSTAR, allow us to probe sub-PeV electron populations through detecting synchrotron X-ray radiation. Our targets include PeVatron candidates detected by HAWC, the Boomerang nebula, PWNe crushed by supernova remnant shocks (or relic PWNe) and G0.9+0.1 in the Galactic Center. Using Fermi-LAT data and available TeV data, we aim to provide a complete, multi-wavelength view of a diverse class of middle-aged (~10-100 kyrs old) PWNe. Our NuSTAR analysis detected hard X-ray emission from the Eel and Boomerang PWNe and characterized their broad-band X-ray spectra most accurately. We plan to apply both time-evolution and multi-zone PWN models to multi-wavelength spectral energy distribution (SED) data over the radio, X-ray, GeV and TeV bands. In this proceeding, we will review our observational campaign and discuss the preliminary results for some PWNe.
△ Less
Submitted 1 August, 2021;
originally announced August 2021.
-
Direct Evidence for Magnetic Reflection of Heavy Ions from High Mach Number Collisionless Shocks
Authors:
Hadi Madanian,
Steven J. Schwartz,
Stephen A. Fuselier,
David Burgess,
Drew L. Turner,
Li-Jen Chen,
Mihir I. Desai,
Michael J. Starkey
Abstract:
Strong shocks in collisionless plasmas, such as supernovae shocks and shocks driven by coronal mass ejections, are known to be a primary source of energetic particles. Due to their different mass per charge ratio, the interaction of heavy ions with the shock layer differs from that of protons, and injection of these ions into acceleration processes is a challenge. Here we show the first direct obs…
▽ More
Strong shocks in collisionless plasmas, such as supernovae shocks and shocks driven by coronal mass ejections, are known to be a primary source of energetic particles. Due to their different mass per charge ratio, the interaction of heavy ions with the shock layer differs from that of protons, and injection of these ions into acceleration processes is a challenge. Here we show the first direct observational evidence of magnetic reflection of alpha particles from a high Mach number quasi-perpendicular shock using in-situ spacecraft measurements. The intense magnetic amplification at the shock front associated with nonstationarity modulates the trajectory of alpha particles, some of which travel back upstream as they gyrate in the enhanced magnetic field and experience further acceleration in the upstream region. Our results in particular highlight the important role of high magnetic amplification in seeding heavy ions into the energization processes at nonstationary reforming shocks.
△ Less
Submitted 20 August, 2021; v1 submitted 18 June, 2021;
originally announced June 2021.
-
The interaction between shocks and plasma turbulence: phase space transport
Authors:
Domenico Trotta,
Francesco Valentini,
David Burgess,
Sergio Servidio
Abstract:
The interaction of collisionless shocks with fully developed plasma turbulence is numerically investigated. Hybrid kinetic simulations, where a turbulent jet is slammed against an oblique shock, are employed to address the role of upstream turbulence on plasma transport. A novel technique, using coarse-graining of the Vlasov equation, is proposed, showing that the transport properties strongly dep…
▽ More
The interaction of collisionless shocks with fully developed plasma turbulence is numerically investigated. Hybrid kinetic simulations, where a turbulent jet is slammed against an oblique shock, are employed to address the role of upstream turbulence on plasma transport. A novel technique, using coarse-graining of the Vlasov equation, is proposed, showing that the transport properties strongly depend on upstream turbulence strength, with turbulent structures-modulated patterns. These results might be relevant for the understanding of acceleration and heating processes in space plasmas.
△ Less
Submitted 20 November, 2020;
originally announced November 2020.
-
The Evolution and Role of Solar Wind Turbulence in the Inner Heliosphere
Authors:
C. H. K. Chen,
S. D. Bale,
J. W. Bonnell,
D. Borovikov,
T. A. Bowen,
D. Burgess,
A. W. Case,
B. D. G. Chandran,
T. Dudok de Wit,
K. Goetz,
P. R. Harvey,
J. C. Kasper,
K. G. Klein,
K. E. Korreck,
D. Larson,
R. Livi,
R. J. MacDowall,
D. M. Malaspina,
A. Mallet,
M. D. McManus,
M. Moncuquet,
M. Pulupa,
M. Stevens,
P. Whittlesey
Abstract:
The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy leve…
▽ More
The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy levels by more than an order of magnitude, a magnetic field spectral index of -3/2 matching that of the velocity and both Elsasser fields, a lower magnetic compressibility consistent with a smaller slow-mode kinetic energy fraction, and a much smaller outer scale that has had time for substantial nonlinear processing. There is also an overall increase in the dominance of outward-propagating Alfvénic fluctuations compared to inward-propagating ones, and the radial variation of the inward component is consistent with its generation by reflection from the large-scale gradient in Alfvén speed. The energy flux in this turbulence at 0.17 au was found to be ~10% of that in the bulk solar wind kinetic energy, becoming ~40% when extrapolated to the Alfvén point, and both the fraction and rate of increase of this flux towards the Sun is consistent with turbulence-driven models in which the solar wind is powered by this flux.
△ Less
Submitted 4 December, 2019;
originally announced December 2019.
-
Modeling Kelvin-Helmholtz instability-driven turbulence with hybrid simulations of Alfvénic turbulence
Authors:
Luca Franci,
Julia E. Stawarz,
Emanuele Papini,
Petr Hellinger,
Takuma Nakamura,
David Burgess,
Simone Landi,
Andrea Verdini,
Lorenzo Matteini,
Robert Ergun,
Olivier Le Contel,
Per-Arne Lindqvist
Abstract:
Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and…
▽ More
Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfvénic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfvén wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath.
△ Less
Submitted 17 November, 2019;
originally announced November 2019.
-
Fast acceleration of transrelativistic electrons in astrophysical turbulence
Authors:
Domenico Trotta,
Luca Franci,
David Burgess,
Petr Hellinger
Abstract:
Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the ene…
▽ More
Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the energization of transrelativistic electrons in turbulence using hybrid particle-in-cell, which provide a realistic model of Alfvénic turbulence from MHD to sub-ion scales, and test particle simulations for electrons. We find that, depending on the electron initial energy and turbulence strength, electrons may undergo a fast and efficient phase of energization due to the magnetic curvature drift during the time they are trapped in dynamic magnetic structures. In addition, electrons are accelerated stochastically which is a slower process that yields lower maximum energies. The combined effect of these two processes determines the overall electron acceleration. With appropriate turbulence parameters, we find that superthermal electrons can be accelerated up to relativistic energies. For example, with heliospheric parameters and a relatively high turbulence level, rapid energization to MeV energies is possible.
△ Less
Submitted 17 April, 2020; v1 submitted 21 October, 2019;
originally announced October 2019.
-
Statistical analysis of ions in two-dimensional plasma turbulence
Authors:
Francesco Pecora,
Francesco Pucci,
Giovanni Lapenta,
David Burgess,
Sergio Servidio
Abstract:
The statistical properties of ions in two-dimensional fully developed turbulence have been compared between two different numerical algorithms. In particular, we compare Hybrid Particle In Cell (hybrid PIC with fluid electrons) and full PIC simulations, focusing on particle diffusion and acceleration phenomena. To investigate several heliospheric plasma conditions, a series of numerical simulation…
▽ More
The statistical properties of ions in two-dimensional fully developed turbulence have been compared between two different numerical algorithms. In particular, we compare Hybrid Particle In Cell (hybrid PIC with fluid electrons) and full PIC simulations, focusing on particle diffusion and acceleration phenomena. To investigate several heliospheric plasma conditions, a series of numerical simulations has been performed by varying the plasma $β$ - the ratio between kinetic and magnetic pressure. These numerical studies allow the exploration of different scenarios, going from the solar corona (low $β$) to the solar wind ($β\sim 1$), as well as the Earth's magnetosheath (high $β$). It has been found that the two approaches compare pretty well, especially for the spectral properties of the magnetic field and the ion diffusion statistics. Small differences among the models have been found regarding the electric field behaviour at sub-ion scales and the acceleration statistics, due evidently to the more consistent treatment of the plasma in the full PIC approach.
△ Less
Submitted 7 August, 2019;
originally announced August 2019.
-
A Case for Electron-Astrophysics
Authors:
Daniel Verscharen,
Robert T. Wicks,
Olga Alexandrova,
Roberto Bruno,
David Burgess,
Christopher H. K. Chen,
Raffaella D'Amicis,
Johan De Keyser,
Thierry Dudok de Wit,
Luca Franci,
Jiansen He,
Pierre Henri,
Satoshi Kasahara,
Yuri Khotyaintsev,
Kristopher G. Klein,
Benoit Lavraud,
Bennett A. Maruca,
Milan Maksimovic,
Ferdinand Plaschke,
Stefaan Poedts,
Chirstopher S. Reynolds,
Owen Roberts,
Fouad Sahraoui,
Shinji Saito,
Chadi S. Salem
, et al. (5 additional authors not shown)
Abstract:
A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour,…
▽ More
A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.
△ Less
Submitted 6 August, 2019;
originally announced August 2019.
-
First observations of irregular surface of interplanetary shocks at ion scales by Cluster
Authors:
Primož Kajdič,
Luis Preisser,
Xóchitl Blanco-Cano,
David Burgess,
Domenico Trotta
Abstract:
We present the first observational evidence of the irregular surface of interplanetary (IP) shocks by using multi-spacecraft observations of the Cluster mission. In total we discuss observations of four IP shocks that exhibit moderate Alfvénic Mach numbers (M$_A\leq$6.5). Three of them are high-$β$ shocks with upstream $β$ = 2.2--3.7. During the times when these shocks were observed, the Cluster s…
▽ More
We present the first observational evidence of the irregular surface of interplanetary (IP) shocks by using multi-spacecraft observations of the Cluster mission. In total we discuss observations of four IP shocks that exhibit moderate Alfvénic Mach numbers (M$_A\leq$6.5). Three of them are high-$β$ shocks with upstream $β$ = 2.2--3.7. During the times when these shocks were observed, the Cluster spacecraft formed constellations with inter-spacecraft separations ranging from less than one upstream ion inertial length (d$_i$) up to 100~d$_i$. Expressed in kilometers, the distances ranged between 38~km and $\sim$10$^4$~km. We show that magnetic field profiles and the local shock normals of observed shocks are very similar when the spacecraft are of the order of one d$_i$ apart, but are strikingly different when the distances increase to ten or more d$_i$. We interpret these differences to be due to the irregular surface of IP shocks and discuss possible causes for such irregularity. We strengthen our interpretation by comparing observed shock profiles with profiles of simulated shocks. The latter had similar characteristics (M$_A$, $θ_{BN}$, upstream ion $β$) as observed shocks and the profiles were obtained at separations across the simulation domain equivalent to the Cluster inter-spacecraft distances.
△ Less
Submitted 27 March, 2019;
originally announced March 2019.
-
Electron acceleration at quasi-perpendicular shocks in sub- and supercritical regimes: 2D and 3D simulations
Authors:
D. Trotta,
D. Burgess
Abstract:
Shock accelerated electrons are found in many astrophysical environments, and the mechanisms by which they are accelerated to high energies are still not completely clear. For relatively high Mach numbers, the shock is supercritical, and its front exhibit broadband fluctuations, or ripples. Shock surface fluctuations have been object of many observational and theoretical studies, and are known to…
▽ More
Shock accelerated electrons are found in many astrophysical environments, and the mechanisms by which they are accelerated to high energies are still not completely clear. For relatively high Mach numbers, the shock is supercritical, and its front exhibit broadband fluctuations, or ripples. Shock surface fluctuations have been object of many observational and theoretical studies, and are known to be important for electron acceleration. We employ a combination of hybrid Particle-In-Cell and test-particle methods to study how shock surface fluctuations influence the acceleration of suprathermal electrons in fully three dimensional simulations, and we give a complete comparison for the 2D and 3D cases. A range of different quasi-perpendicular shocks in 2D and 3D is examined, over a range of parameters compatible with the ones observed in the solar wind. Initial electron velocity distributions are taken as kappa functions, consistent with solar wind \emph{in-situ} measurements. Electron acceleration is found to be enhanced in the supercritical regime compared to subcritical. When the fully three-dimensional structure of the shock front is resolved, slightly larger energisation for the electrons is observed, and we suggest that this is due to the possibility for the electrons to interact with more than one surface fluctuation per interaction. In the supecritical regime, efficient electron energisation is found also at shock geometries departing from $θ_{Bn}$ very close to 90$^\circ$. Two dimensional simulations show indications of unrealistic electron trapping, leading to slightly higher energisation in the subcritical cases.
△ Less
Submitted 2 August, 2018;
originally announced August 2018.
-
An in situ comparison of electron acceleration at collisionless shocks under differing upstream magnetic field orientations
Authors:
A. Masters,
A. H. Sulaiman,
Ł. Stawarz,
B. Reville,
N. Sergis,
M. Fujimoto,
D. Burgess,
A. J. Coates,
M. K. Dougherty
Abstract:
A leading explanation for the origin of Galactic cosmic rays is acceleration at high-Mach number shock waves in the collisionless plasma surrounding young supernova remnants. Evidence for this is provided by multi-wavelength non-thermal emission thought to be associated with ultrarelativistic electrons at these shocks. However, the dependence of the electron acceleration process on the orientation…
▽ More
A leading explanation for the origin of Galactic cosmic rays is acceleration at high-Mach number shock waves in the collisionless plasma surrounding young supernova remnants. Evidence for this is provided by multi-wavelength non-thermal emission thought to be associated with ultrarelativistic electrons at these shocks. However, the dependence of the electron acceleration process on the orientation of the upstream magnetic field with respect to the local normal to the shock front (quasi-parallel/quasi-perpendicular) is debated. Cassini spacecraft observations at Saturn's bow shock has revealed examples of electron acceleration under quasi-perpendicular conditions, and the first in situ evidence of electron acceleration at a quasi-parallel shock. Here we use Cassini data to make the first comparison between energy spectra of locally accelerated electrons under these differing upstream magnetic field regimes. We present data taken during a quasi-perpendicular shock crossing on 2008 March 8 and during a quasi-parallel shock crossing on 2007 February 3, highlighting that both were associated with electron acceleration to at least MeV energies. The magnetic signature of the quasi-perpendicular crossing has a relatively sharp upstream-downstream transition, and energetic electrons were detected close to the transition and immediately downstream. The magnetic transition at the quasi-parallel crossing is less clear, energetic electrons were encountered upstream and downstream, and the electron energy spectrum is harder above ~100 keV. We discuss whether the acceleration is consistent with diffusive shock acceleration theory in each case, and suggest that the quasi-parallel spectral break is due to an energy-dependent interaction between the electrons and short, large-amplitude magnetic structures.
△ Less
Submitted 31 May, 2017;
originally announced May 2017.
-
Multiple current sheet systems in the outer heliosphere: Energy release and turbulence
Authors:
David Burgess,
P. W. Gingell,
Lorenzo Matteini
Abstract:
In the outer heliosphere, beyond the solar wind termination shock, it is expected that the warped heliospheric current sheet forms a region of closely-packed, multiple, thin current sheets. Such a system may be subject to the ion-kinetic tearing instability, and hence generate magnetic islands and hot populations of ions associated with magnetic reconnection. Reconnection processes in this environ…
▽ More
In the outer heliosphere, beyond the solar wind termination shock, it is expected that the warped heliospheric current sheet forms a region of closely-packed, multiple, thin current sheets. Such a system may be subject to the ion-kinetic tearing instability, and hence generate magnetic islands and hot populations of ions associated with magnetic reconnection. Reconnection processes in this environment have important implications for local particle transport, and for particle acceleration at reconnection sites and in turbulence. We study this complex environment by means of three-dimensional hybrid simulations over long time scales, in order to capture the evolution from linear growth of the tearing instability to a fully developed turbulent state at late times. The final state develops from the highly ordered initial state via both forward and inverse cascades. Component and spectral anisotropy in the magnetic fluctuations is present when a guide field is included. The inclusion of a population of new-born interstellar pickup protons does not strongly affect these results. Finally, we conclude that reconnection between multiple current sheets can act as an important source of turbulence in the outer heliosphere, with implications for energetic particle acceleration and propagation.
△ Less
Submitted 22 March, 2016;
originally announced March 2016.
-
Microstructure in two- and three-dimensional hybrid simulations of perpendicular collisionless shocks
Authors:
David Burgess,
Petr Hellinger,
Peter W. Gingell,
Pavel M. Trávníček
Abstract:
Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been…
▽ More
Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been carried out using the hybrid method, wherein the ions are treated as particles and the electron response is modelled as a massless fluid. We confirm the results of earlier 2-D simulations showing both field-parallel aligned propagating fluctuations and fluctuations carried by the reflected-gyrating ions. In addition, it is shown that, for 2-D simulations of the shock coplanarity plane, the presence of short-wavelength fluctuations in all magnetic components is associated with the ion Weibel instability driven at the upstream edge of the foot by the reflected-gyrating ions. In 3-D simulations we show for the first time that the dominant microstructure is due to a coupling between field-parallel propagating fluctuations in the ramp and the motion of the reflected ions. This results in a pattern of fluctuations counter-propagating across the surface of the shock at an angle inclined to the magnetic field direction, due to a combination of field-parallel motion at the Alfvén speed of the ramp, and motion in the sense of gyration of the reflected ions.
△ Less
Submitted 7 January, 2016;
originally announced January 2016.
-
Ion Acceleration at the Quasi-Parallel Bow Shock: Decoding the Signature of Injection
Authors:
Torbjörn Sundberg,
Christopher T. Haynes,
David Burgess,
Christian X. Mazelle
Abstract:
Collisionless shocks are efficient particle accelerators. At Earth, ions with energies exceeding 100 keV are seen upstream of the bow shock when the magnetic geometry is quasi-parallel, and large-scale supernova remnant shocks can accelerate ions into cosmic rays energies. This energization is attributed to diffusive shock acceleration, however, for this process to become active the ions must firs…
▽ More
Collisionless shocks are efficient particle accelerators. At Earth, ions with energies exceeding 100 keV are seen upstream of the bow shock when the magnetic geometry is quasi-parallel, and large-scale supernova remnant shocks can accelerate ions into cosmic rays energies. This energization is attributed to diffusive shock acceleration, however, for this process to become active the ions must first be sufficiently energized. How and where this initial acceleration takes place has been one of the key unresolved issues in shock acceleration theory. Using Cluster spacecraft observations, we study the signatures of ion reflection events in the turbulent transition layer upstream of the terrestrial bow shock, and with the support of a hybrid simulation of the shock, we show that these reflection signatures are characteristic of the first step in the ion injection process. These reflection events develop in particular in the region where the trailing edge of large-amplitude upstream waves intercept the local shock ramp and the upstream magnetic field changes from quasi-perpendicular to quasi-parallel. The dispersed ion velocity signature observed can be attributed to a rapid succession of ion reflections at this wave boundary. After the ions' initial interaction with the shock, they flow upstream along the quasi-parallel magnetic field. Each subsequent wave front in the upstream region will sweep the ions back toward the shock, where they gain energy with each transition between the upstream and the shock wave frames. Within three to five gyroperiods, some ions have gained enough parallel velocity to escape upstream, thus completing the injection process.
△ Less
Submitted 20 January, 2016; v1 submitted 6 November, 2015;
originally announced November 2015.
-
Quasiperpendicular high Mach number Shocks
Authors:
A. H. Sulaiman,
A. Masters,
M. K. Dougherty,
D. Burgess,
M. Fujimoto,
G. B. Hospodarsky
Abstract:
Shock waves exist throughout the universe and are fundamental to understanding the nature of collisionless plasmas. Reformation is a process, driven by microphysics, which typically occurs at high Mach number supercritical shocks. While ongoing studies have investigated this process extensively both theoretically and via simulations, their observations remain few and far between. In this letter we…
▽ More
Shock waves exist throughout the universe and are fundamental to understanding the nature of collisionless plasmas. Reformation is a process, driven by microphysics, which typically occurs at high Mach number supercritical shocks. While ongoing studies have investigated this process extensively both theoretically and via simulations, their observations remain few and far between. In this letter we present a study of very high Mach number shocks in a parameter space that has been poorly explored and we identify reformation using in situ magnetic field observations from the Cassini spacecraft at 10 AU. This has given us an insight into quasi-perpendicular shocks across two orders of magnitude in Alfven Mach number (MA) which could potentially bridge the gap between modest terrestrial shocks and more exotic astrophysical shocks. For the first time, we show evidence for cyclic reformation controlled by specular ion reflection occurring at the predicted timescale of ~0.3 τc, where τc is the ion gyroperiod. In addition, we experimentally reveal the relationship between reformation and MA and focus on the magnetic structure of such shocks to further show that for the same MA, a reforming shock exhibits stronger magnetic field amplification than a shock that is not reforming.
△ Less
Submitted 17 September, 2015;
originally announced September 2015.
-
Magnetic Field Rotations in the Solar Wind at Kinetic Scales
Authors:
C. H. K. Chen,
L. Matteini,
D. Burgess,
T. S. Horbury
Abstract:
The solar wind magnetic field contains rotations at a broad range of scales, which have been extensively studied in the MHD range. Here we present an extension of this analysis to the range between ion and electron kinetic scales. The distribution of rotation angles was found to be approximately log-normal, shifting to smaller angles at smaller scales almost self-similarly, but with small, statist…
▽ More
The solar wind magnetic field contains rotations at a broad range of scales, which have been extensively studied in the MHD range. Here we present an extension of this analysis to the range between ion and electron kinetic scales. The distribution of rotation angles was found to be approximately log-normal, shifting to smaller angles at smaller scales almost self-similarly, but with small, statistically significant changes of shape. The fraction of energy in fluctuations with angles larger than $α$ was found to drop approximately exponentially with $α$, with e-folding angle $9.8^\circ$ at ion scales and $0.66^\circ$ at electron scales, showing that large angles ($α> 30^\circ$) do not contain a significant amount of energy at kinetic scales. Implications for kinetic turbulence theory and the dissipation of solar wind turbulence are discussed.
△ Less
Submitted 15 September, 2015; v1 submitted 28 July, 2015;
originally announced July 2015.
-
Electron temperature anisotropy in an expanding plasma: Particle-in-Cell simulations
Authors:
Enrico Camporeale,
David Burgess
Abstract:
We perform fully-kinetic particle-in-cell simulations of an hot plasma that expands radially in a cylindrical geometry. The aim of the paper is to study the consequent development of the electron temperature anisotropy in an expanding plasma flow as found in a collisionless stellar wind. Kinetic plasma theory and simulations have shown that the electron temperature anisotropy is controlled by fl…
▽ More
We perform fully-kinetic particle-in-cell simulations of an hot plasma that expands radially in a cylindrical geometry. The aim of the paper is to study the consequent development of the electron temperature anisotropy in an expanding plasma flow as found in a collisionless stellar wind. Kinetic plasma theory and simulations have shown that the electron temperature anisotropy is controlled by fluctuations driven by electromagnetic kinetic instabilities. In this study the temperature anisotropy is driven self-consistently by the expansion. While the expansion favors an increase of parallel anisotropy ($T_\parallel>T_\perp$), the onset of the firehose instability will tend to decrease it. We show the results for a supersonic, subsonic, and static expansion flows, and suggest possible applications of the results for the solar wind and other stellar winds.
△ Less
Submitted 11 February, 2010;
originally announced February 2010.
-
Simulations of Electron Acceleration at Collisionless Shocks: The Effects of Surface Fluctuations
Authors:
D. Burgess
Abstract:
Energetic electrons are a common feature of interplanetary shocks and planetary bow shocks, and they are invoked as a key component of models of nonthermal radio emission, such as solar radio bursts. A simulation study is carried out of electron acceleration for high Mach number, quasi-perpendicular shocks, typical of the shocks in the solar wind. Two dimensional self-consistent hybrid shock sim…
▽ More
Energetic electrons are a common feature of interplanetary shocks and planetary bow shocks, and they are invoked as a key component of models of nonthermal radio emission, such as solar radio bursts. A simulation study is carried out of electron acceleration for high Mach number, quasi-perpendicular shocks, typical of the shocks in the solar wind. Two dimensional self-consistent hybrid shock simulations provide the electric and magnetic fields in which test particle electrons are followed. A range of different shock types, shock normal angles, and injection energies are studied. When the Mach number is low, or the simulation configuration suppresses fluctuations along the magnetic field direction, the results agree with theory assuming magnetic moment conserving reflection (or Fast Fermi acceleration), with electron energy gains of a factor only 2 - 3. For high Mach number, with a realistic simulation configuration, the shock front has a dynamic rippled character. The corresponding electron energization is radically different: Energy spectra display: (1) considerably higher maximum energies than Fast Fermi acceleration; (2) a plateau, or shallow sloped region, at intermediate energies 2 - 5 times the injection energy; (3) power law fall off with increasing energy, for both upstream and downstream particles, with a slope decreasing as the shock normal angle approaches perpendicular; (4) sustained flux levels over a broader region of shock normal angle than for adiabatic reflection. All these features are in good qualitative agreement with observations, and show that dynamic structure in the shock surface at ion scales produces effective scattering and can be responsible for making high Mach number shocks effective sites for electron acceleration.
△ Less
Submitted 24 October, 2006;
originally announced October 2006.