-
Particle In Cell Simulations of Mildly Relativistic Outflows in Kilonova Emissions
Authors:
Mohira Rassel,
Patrick Kilian,
Vito Aberham,
Felix Spanier,
Nicole Lloyd-Ronning,
Chris L. Fryer
Abstract:
The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet as well as thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical…
▽ More
The electromagnetic emission from neutron star mergers is comprised of multiple components. Synchrotron emission from the disk-powered jet as well as thermal emission from the merger ejecta (powered by a variety of sources) are among the most studied sources. The low masses and high velocities of the merger ejecta quickly develop conditions where emission from collisionless shocks becomes critical and synchrotron emission from the merger ejecta constitutes a third component to the observed signal. The aim of this project is to examine shock development, magnetic field generation and particle acceleration in the case of mildly relativistic shocks, which are expected when the tidal ejecta of neutron star mergers drive a shock into the external medium. Using LANL's VPIC (vector particle-in-cell) code, we have run simulations of such mildly-relativistic, collisionless, weakly-magnetized plasmas and compute the resultant magnetic fields and particle energy spectra. We show the effects of varying plasma conditions, as well as explore the validity of using different proton to electron mass ratios in VPIC. Our results have implications for observing late-time electromagnetic counterparts to gravitational wave detections of neutron star mergers.
△ Less
Submitted 19 May, 2023;
originally announced May 2023.
-
Comment on "Nonideal Fields Solve the Injection Problem in Relativistic Reconnection"
Authors:
Fan Guo,
Xiaocan Li,
Omar French,
William Daughton,
William Matthaeus,
Qile Zhang,
Yi-Hsin Liu,
Patrick Kilian,
Grant Johnson,
Hui Li
Abstract:
Recently, Sironi (PRL, 128, 145102; S22) reported the correlation between particles accelerated into high energy and their crossings of regions with electric field larger than magnetic field (E>B regions) in kinetic simulations of relativistic magnetic reconnection. They claim that electric fields in E>B regions (for a vanishing guide field) dominate in accelerating particles to the injection ener…
▽ More
Recently, Sironi (PRL, 128, 145102; S22) reported the correlation between particles accelerated into high energy and their crossings of regions with electric field larger than magnetic field (E>B regions) in kinetic simulations of relativistic magnetic reconnection. They claim that electric fields in E>B regions (for a vanishing guide field) dominate in accelerating particles to the injection energy. S22 presented test-particle simulations showing that if particle energies are reset to low energies in E>B regions, efficient injection is suppressed. This Comment re-examines these claims by analyzing a simulation resembling the reference case in S22. We show that during crossings E>B acceleration only contributes a small fraction to the injection energy as E>B regions only host particles for a short duration. The energization before any E>B crossings has a comparable contribution, indicating E>B regions are not unique in pre-accelerating particles. A new test-particle simulation shows that zero-outing electric fields in E>B regions does not strongly influence the injection. We suggest that the procedure used in S22 to exclude E>B acceleration partly removes acceleration outside E>B regions, leading to a false conclusion.
△ Less
Submitted 5 August, 2022;
originally announced August 2022.
-
Single Pulse Dispersion Measure of the Crab Pulsar
Authors:
N. Lewandowska,
P. B. Demorest,
M. A. McLaughlin,
P. Kilian,
T. H. Hankins
Abstract:
We investigate the use of bright single pulses from the Crab pulsar to determine separately the dispersion measure (DM) for the Main Pulse and Interpulse components. We develop two approaches using cross correlation functions (CCFs). The first method computes the CCF of the total intensity of each of 64 frequency channels with a reference channel and converts the time lag of maximum correlation in…
▽ More
We investigate the use of bright single pulses from the Crab pulsar to determine separately the dispersion measure (DM) for the Main Pulse and Interpulse components. We develop two approaches using cross correlation functions (CCFs). The first method computes the CCF of the total intensity of each of 64 frequency channels with a reference channel and converts the time lag of maximum correlation into a DM. The second method separately computes the CCF between every pair of channels for each individual bright pulse and extracts an average DM from the distribution of all channel-pair DMs. Both methods allow the determination of the DM with a relative uncertainty of better than 10^-5 and provide robust estimates for the uncertainty of the best-fit value. We find differences in DM between the Main Pulse, the Low Frequency Interpulse, and the High Frequency Interpulse using both methods in a frequency range from 4 to 6 GHz. Earlier observations of the High Frequency Interpulse carried out by Hankins et al. (2016) resulted in DM_HFIP-DM_MP of 0.010 +- 0.016 pc cm^-3. Our results indicate a DM_HFIP-DM_MP of 0.0127 +- 0.0011 pc cm^-3 (with DM_comp being the DM value of the respective emission component), confirming earlier results with an independent method. During our studies we also find a relation between the brightness of single pulses in the High Frequency Interpulse and their DM. We also discuss the application of the developed methods on the identification of substructures in the case of Fast Radio Bursts.
△ Less
Submitted 9 July, 2022;
originally announced July 2022.
-
Particle Acceleration in Magnetic Reconnection with Ad hoc Pitch-angle Scattering
Authors:
Grant Johnson,
Patrick Kilian,
Fan Guo,
Xiaocan Li
Abstract:
Particle acceleration during magnetic reconnection is a long-standing topic in space, solar and astrophysical plasmas. Recent 3D particle-in-cell simulations of magnetic reconnection show that particles can leave flux ropes due to 3D field-line chaos, allowing particles to access additional acceleration sites, gain more energy through Fermi acceleration, and develop a power-law energy distribution…
▽ More
Particle acceleration during magnetic reconnection is a long-standing topic in space, solar and astrophysical plasmas. Recent 3D particle-in-cell simulations of magnetic reconnection show that particles can leave flux ropes due to 3D field-line chaos, allowing particles to access additional acceleration sites, gain more energy through Fermi acceleration, and develop a power-law energy distribution. This 3D effect does not exist in traditional 2D simulations, where particles are artificially confined to magnetic islands due to their restricted motions across field lines. Full 3D simulations, however, are prohibitively expensive for most studies. Here, we attempt to reproduce 3D results in 2D simulations by introducing ad hoc pitch-angle scattering to a small fraction of the particles. We show that scattered particles are able to transport out of 2D islands and achieve more efficient Fermi acceleration, leading to a significant increase of energetic particle flux. We also study how the scattering frequency influences the nonthermal particle spectra. This study helps achieve a complete picture of particle acceleration in magnetic reconnection.
△ Less
Submitted 17 May, 2022;
originally announced May 2022.
-
Fermi-type particle acceleration from magnetic reconnection at the termination shock of a relativistic striped wind
Authors:
Yingchao Lu,
Fan Guo,
Patrick Kilian,
Hui Li,
Chengkun Huang,
Edison Liang
Abstract:
An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell (PIC) simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechani…
▽ More
An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell (PIC) simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultra-relativistic regime with bulk Lorentz factor up to γ_{0}=10^{6}. Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. We find that the resulting energy spectra crucially depend on λ/d_{e}. When λ/d_{e} is large (λ\gtrsim40d_{e}) , the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime with upstream magnetization parameter σ_{0}=10. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as γ_{0} and σ_{0} increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor γ_{0}\approx10^{6} and magnetization parameter σ_{0}=10, which can explain the recent observations of high-energy gamma-rays from pulsar wind nebulae (PWNe).
△ Less
Submitted 27 October, 2020;
originally announced October 2020.
-
Magnetic Energy Release, Plasma Dynamics, and Particle Acceleration during Relativistic Turbulent Magnetic Reconnection
Authors:
Fan Guo,
Xiaocan Li,
William Daughton,
Hui Li,
Patrick Kilian,
Yi-Hsin Liu,
Qile Zhang,
Haocheng Zhang
Abstract:
In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be a primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored…
▽ More
In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be a primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored issue is the consequence of 3D dynamics, where turbulent structures are naturally generated as various types of instabilities develop. This paper presents a series of 3D, fully-kinetic simulations of relativistic turbulent magnetic reconnection (RTMR) in positron-electron plasmas with system domains much larger than kinetic scales. Our simulations start from a force-free current sheet with several different modes of long wavelength magnetic field perturbations, which drive additional turbulence in the reconnection region. Because of this, the current layer breaks up and the reconnection region quickly evolves into a turbulent layer filled with coherent structures such as flux ropes and current sheets. We find that plasma dynamics in RTMR is vastly different from their 2D counterparts in many aspects. The flux ropes evolve rapidly after their generation, and can be completely disrupted due to the secondary kink instability. This turbulent evolution leads to superdiffusion behavior of magnetic field lines as seen in MHD studies of turbulent reconnection. Meanwhile, nonthermal particle acceleration and energy-release time scale can be very fast and do not strongly depend on the turbulence amplitude. The main acceleration mechanism is a Fermi-like acceleration process supported by the motional electric field, whereas the non-ideal electric field acceleration plays a subdominant role. We discuss possible observational implications of 3D RTMR in high-energy astrophysics.
△ Less
Submitted 29 July, 2021; v1 submitted 6 August, 2020;
originally announced August 2020.
-
Recent Progress on Particle Acceleration and Reconnection Physics during Magnetic Reconnectionin the Magnetically-dominated Relativistic Regime
Authors:
Fan Guo,
Yi-Hsin Liu,
Xiaocan Li,
Hui Li,
William Daughton,
Patrick Kilian
Abstract:
Magnetic reconnection in strongly magnetized astrophysical plasma environments is believed to be the primary process for fast energy release and particle energization. Currently there is strong interest in relativistic magnetic reconnection, in that it may provide a new explanation for high-energy particle acceleration and radiation in strongly magnetized astrophysical systems. We review recent ad…
▽ More
Magnetic reconnection in strongly magnetized astrophysical plasma environments is believed to be the primary process for fast energy release and particle energization. Currently there is strong interest in relativistic magnetic reconnection, in that it may provide a new explanation for high-energy particle acceleration and radiation in strongly magnetized astrophysical systems. We review recent advances in particle acceleration and reconnection physics in the magnetically-dominated regime. More discussion is focused on the physics of particle acceleration, power-law formation as well as the reconnection rate problem. In addition, we provide an outlook for studying reconnection acceleration mechanisms and kinetic physics in the next step.
△ Less
Submitted 27 June, 2020;
originally announced June 2020.
-
Ion acceleration in non-relativistic quasi-parallel shocks using fully kinetic simulations
Authors:
Cedric Schreiner,
Patrick Kilian,
Felix Spanier,
Patricio A. Muñoz,
Jörg Büchner
Abstract:
The formation of collisionless shock fronts is an ubiquitous phenomenon in space plasma environments. In the solar wind shocks might accompany coronal mass ejections, while even more violent events, such as supernovae, produce shock fronts traveling at relativistic speeds. While the basic concepts of shock formation and particle acceleration in their vicinity are known, many details on a micro-phy…
▽ More
The formation of collisionless shock fronts is an ubiquitous phenomenon in space plasma environments. In the solar wind shocks might accompany coronal mass ejections, while even more violent events, such as supernovae, produce shock fronts traveling at relativistic speeds. While the basic concepts of shock formation and particle acceleration in their vicinity are known, many details on a micro-physical scope are still under discussion. In recent years the hybrid kinetic simulation approach has allowed to study the dynamics and acceleration of protons and heavier ions in great detail. However, Particle-in-Cell codes allow to study the process including also electron dynamics and the radiation pressure. Additionally a further numerical method allows for crosschecking results. We therefore investigate shock formation and particle acceleration with a fully kinetic particle-in-cell code. Besides electrons and protons we also include helium and carbon ions in our simulations of a quasi-parallel shock. We are able to reproduce characteristic features of the energy spectra of the particles, such as the temperature ratios of the different ion species in the downstream which scale with the ratio of particle mass to charge. We also find that approximately 12-15% of the energy of the unperturbed upstream is transferred to the accelerated particles escaping the shock.
△ Less
Submitted 16 March, 2020;
originally announced March 2020.
-
Exploring the acceleration mechanisms for particle injection and power-law formation during trans-relativistic magnetic reconnection
Authors:
Patrick Kilian,
Xiaocan Li,
Fan Guo,
Hui Li
Abstract:
Magnetic reconnection in the relativistic and trans-relativistic regimes is able to accelerate particles to hard power law energy spectra $f \propto γ^{-p}$ (approaching $p=1$). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during mag…
▽ More
Magnetic reconnection in the relativistic and trans-relativistic regimes is able to accelerate particles to hard power law energy spectra $f \propto γ^{-p}$ (approaching $p=1$). The underlying acceleration mechanism that determines the spectral shape is currently a topic of intense investigation. By means of fully kinetic plasma simulations, we carry out a study of particle acceleration during magnetic reconnection in the trans-relativistic regime of a proton-electron plasma. While earlier work in this parameter regime has focused on the effects of electric field parallel to the local magnetic field on the particle injection (from thermal energy to the lower energy bound of the power-law spectrum), here we examine the roles of both parallel and perpendicular electric fields to gain a more complete understanding on the injection process and further development of a power-law spectrum. We show that the parallel electric field does contribute significantly to particle injection, and is more important in the initial phase of magnetic reconnection. However, as the simulation proceeds, the acceleration by the perpendicular electric field becomes more important for particle injection and completely dominates the acceleration responsible for the high-energy power-law spectrum. This holds robustly, in particular for longer reconnection times and larger systems, i.e. in simulations that are more indicative of the processes in astrophysical sources.
△ Less
Submitted 2 July, 2020; v1 submitted 8 January, 2020;
originally announced January 2020.
-
Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations
Authors:
Yingchao Lu,
Patrick Kilian,
Fan Guo,
Hui Li,
Edison Liang
Abstract:
The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show…
▽ More
The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.
△ Less
Submitted 11 April, 2020; v1 submitted 20 September, 2019;
originally announced September 2019.
-
Formation of Power-law Electron Energy Spectra in Three-dimensional Low-$β$ Magnetic Reconnection
Authors:
Xiaocan Li,
Fan Guo,
Hui Li,
Adam Stanier,
Patrick Kilian
Abstract:
While observations have suggested that power-law electron energy spectra are a common outcome of strong energy release during magnetic reconnection, e.g., in solar flares, kinetic simulations have not been able to provide definite evidence of power-laws in energy spectra of non-relativistic reconnection. By means of 3D large-scale fully kinetic simulations, we study the formation of power-law elec…
▽ More
While observations have suggested that power-law electron energy spectra are a common outcome of strong energy release during magnetic reconnection, e.g., in solar flares, kinetic simulations have not been able to provide definite evidence of power-laws in energy spectra of non-relativistic reconnection. By means of 3D large-scale fully kinetic simulations, we study the formation of power-law electron energy spectra in non-relativistic low-$β$ reconnection. We find that both the global spectrum integrated over the entire domain and local spectra within individual regions of the reconnection layer have power-law tails with a spectral index $p \sim 4$ in the 3D simulation, which persist throughout the non-linear reconnection phase until saturation. In contrast, the spectrum in the 2D simulation rapidly evolves and quickly becomes soft. We show that 3D effects such as self-generated turbulence and chaotic magnetic field lines enable the transport of high-energy electrons across the reconnection layer and allow them to access several main acceleration regions. This leads to a sustained and nearly constant acceleration rate for electrons at different energies. We construct a model that explains the observed power-law spectral index in terms of the dynamical balance between particle acceleration and escape from main acceleration regions, which are defined based upon a threshold for the curvature drift acceleration term. This result could be important for explaining the formation of power-law energy spectrum in solar flares.
△ Less
Submitted 4 September, 2019;
originally announced September 2019.
-
High-Energy Polarimetry - a new window to probe extreme physics in AGN jets
Authors:
B. Rani,
H. Zhang,
S. D. Hunter,
F. Kislat,
M. Böttcher,
J. E. McEnery,
D. J. Thompson,
D. Giannios,
F. Guo,
H. Li,
M. Baring,
I. Agudo,
S. Buson,
M. Petropoulou,
V. Pavlidou,
E. Angelakis,
I. Myserlis,
Z. Wadiasingh,
R. M. Curado da Silva,
P. Kilian,
S. Guiriec,
V. V. Bozhilov,
J. Hodgson,
S. Antón,
D. Kazanas
, et al. (9 additional authors not shown)
Abstract:
The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is m…
▽ More
The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Determining the Dominant Acceleration Mechanism during Relativistic Magnetic Reconnection in Large-scale Systems
Authors:
Fan Guo,
Xiaocan Li,
William Daughton,
Patrick Kilian,
Hui Li,
Yi-Hsin Liu,
Wangcheng Yan,
Dylan Ma
Abstract:
While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, we demonstrate that Fermi-type acceleration within the large-scale motional electric fields domin…
▽ More
While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, since it opens up the possiblity of modeling the energetic spectra without resolving microscopic diffusion regions.
△ Less
Submitted 20 June, 2019; v1 submitted 24 January, 2019;
originally announced January 2019.
-
Pressure Balance and Intrabinary Shock Stability in Rotation-Powered State Redback and Transitional Millisecond Pulsar Binary Systems
Authors:
Zorawar Wadiasingh,
Christo Venter,
Alice K. Harding,
Markus Böttcher,
Patrick Kilian
Abstract:
A number of low-mass millisecond pulsar (MSP) binaries in their rotation-powered state exhibit double-peaked X-ray orbital modulation centered at inferior pulsar conjunction. This state, which has been known to persist for years, has recently been interpreted as emission from a shock that enshrouds the pulsar. However, the pressure balance for such a configuration is a crucial unresolved issue. We…
▽ More
A number of low-mass millisecond pulsar (MSP) binaries in their rotation-powered state exhibit double-peaked X-ray orbital modulation centered at inferior pulsar conjunction. This state, which has been known to persist for years, has recently been interpreted as emission from a shock that enshrouds the pulsar. However, the pressure balance for such a configuration is a crucial unresolved issue. We consider two scenarios for pressure balance: a companion magnetosphere and stellar mass loss with gas dominance. It is found that the magnetospheric scenario requires several kilogauss poloidal fields for isobaric surfaces to enshroud the MSP as well as for the magnetosphere to remain stable if there is significant mass loss. For the gas-dominated scenario, it is necessary that the companion wind loses angular momentum prolifically as an advection or heating-dominated flow. Thermal bremsstrahlung cooling in the flow may be observable as a UV to soft X-ray component independent of orbital phase if the mass rate is high. We formulate the general requirements for shock stability against gravitational influences in the pulsar rotation-powered state for the gas-dominated scenario. We explore stabilizing mechanisms, principally irradiation feedback, which anticipates correlated shock emission and companion variability and predicts $ F_γ/F_{\rm X} \lesssim 14$ for the ratio of pulsar magnetospheric $γ$-ray to total shock soft-to-hard X-ray fluxes. This stability criterion implies an unbroken extension of X-ray power-law emission to hundreds of keV for some systems. We explore observational discriminants between the gas-dominated and magnetospheric scenarios, motivating contemporaneous radio through $γ$-ray monitoring of these systems.
△ Less
Submitted 30 October, 2018;
originally announced October 2018.
-
A comment on "Importance of resolving the spectral support of beam-plasma instabilities in simulations" by M. Shalaby et al
Authors:
F. Spanier,
P. Kilian,
R. Schlickeiser
Abstract:
Understanding the physics of beam-plasma systems is of great importance for a large number of astrophysical systems. Recently the validity of particle-in-cell simulations of these systems has been questioned by Shalaby et al. (2017) We analyse one specific beam-plasma and its numerical simulations Kempf et al. (2016) that has been critized and prove that in fact the simulations performed fulfill t…
▽ More
Understanding the physics of beam-plasma systems is of great importance for a large number of astrophysical systems. Recently the validity of particle-in-cell simulations of these systems has been questioned by Shalaby et al. (2017) We analyse one specific beam-plasma and its numerical simulations Kempf et al. (2016) that has been critized and prove that in fact the simulations performed fulfill the new validity criteria introduced by Shalaby et al. (2016).
△ Less
Submitted 25 July, 2017; v1 submitted 28 April, 2017;
originally announced April 2017.
-
Particle Scattering off of Right-Handed Dispersive Waves
Authors:
Cedric Schreiner,
Patrick Kilian,
Felix Spanier
Abstract:
Resonant scattering of fast particles off low frequency plasma waves is a major process determining transport characteristics of energetic particles in the heliosphere and contributing to their acceleration. Usually, only Alfvén waves are considered for this process, although dispersive waves are also present throughout the heliosphere. We investigate resonant interaction of energetic electrons wi…
▽ More
Resonant scattering of fast particles off low frequency plasma waves is a major process determining transport characteristics of energetic particles in the heliosphere and contributing to their acceleration. Usually, only Alfvén waves are considered for this process, although dispersive waves are also present throughout the heliosphere. We investigate resonant interaction of energetic electrons with dispersive, right-handed waves. For the interaction of particles and a single wave a variable transformation into the rest frame of the wave can be performed. Here, wellestablished analytic models derived in the framework of magnetostatic quasi-linear theory (QLT) can be used as a reference to validate simulation results. However, this approach fails as soon as several dispersive waves are involved. Based on analytic solutions modeling the scattering amplitude in the magnetostatic limit, we present an approach to modify these equations for the use in the plasma frame. Thereby we aim at a description of particle scattering in the presence of several waves. A Particle-in-Cell (PiC) code is employed to study wave-particle scattering on a micro-physically correct level and to test the modified model equations. We investigate the interactions of electrons at different energies (from 1 keV to 1 MeV) and right-handed waves with various amplitudes. Differences between model and simulation arise in the case of high amplitudes or several waves. Analyzing the trajectories of single particles we find no microscopic diffusion in the case of a single plasma wave, although a broadening of the particle distribution can be observed.
△ Less
Submitted 1 November, 2016;
originally announced November 2016.
-
Energy loss of intergalactic pair beams: Particle-in-Cell simulation
Authors:
Andreas Kempf,
Patrick Kilian,
Felix Spanier
Abstract:
The change of the distribution function of electron-positron pair beams determines whether GeV photons can be produced as secondary radiation from TeV photons. We will discuss the instabilities driven by pair beams. The system of a thermal proton-electron plasma and the electron-positron beam is collision free. We have, therefore, used the Particle-in-Cell simulation approach. It was necessary to…
▽ More
The change of the distribution function of electron-positron pair beams determines whether GeV photons can be produced as secondary radiation from TeV photons. We will discuss the instabilities driven by pair beams. The system of a thermal proton-electron plasma and the electron-positron beam is collision free. We have, therefore, used the Particle-in-Cell simulation approach. It was necessary to alter the physical parameters, but the ordering of growth rates has been retained. We were able to show that plasma instabilities can be recovered in particle-in-cell simulations, but their effect on the pair distribution function is negligible for beam-background energy density ratios typically found in blazars.
△ Less
Submitted 2 December, 2015;
originally announced December 2015.
-
Fundamental and harmonic plasma emission in different plasma environments
Authors:
Urs Ganse,
Patrick Kilian,
Felix Spanier,
Rami Vainio
Abstract:
Aims: Emission of radio waves from plasmas through plasma emission with fundamental and harmonic frequencies is a familiar process known from solar type II radio bursts. Current models assume the existence of counterstreaming electron beam populations excited at shocks as sources for these emission features, which limits the plasma parameters to reasonable heliospheric shock conditions. However, s…
▽ More
Aims: Emission of radio waves from plasmas through plasma emission with fundamental and harmonic frequencies is a familiar process known from solar type II radio bursts. Current models assume the existence of counterstreaming electron beam populations excited at shocks as sources for these emission features, which limits the plasma parameters to reasonable heliospheric shock conditions. However, situations in which counterstreaming electron beams are present can also occur with different plasma parameters, such as higher magnetisation, including but not limited to our Sun. Similar radio emissions might also occur from these situations.
Methods: We used particle-in-cell simulations, to compare plasma microphysics of radio emission processes from counterstreaming beams in different plasma environments that differed in density and magnetization.
Results: Although large differences in wave populations are evident, the emission process of type II bursts appears to be qualitatively unaffected and shows the same behaviour in all environments.
△ Less
Submitted 10 March, 2014;
originally announced March 2014.
-
Numerical challenges in kinetic simulations of three-wave interactions
Authors:
Urs Ganse,
Patrick Kilian,
Stefan Siegel,
Felix Spanier
Abstract:
Generation of radio bursts in CME foreshock regions and turbulent cascades in the solar wind are assumed to be results of three-wave interaction processes of dispersive plasma modes. Using our Particle in Cell code ACRONYM, we have studied the behaviour of kinetic wavemodes in the presence of beamed electron populations, with a focus on type II radio burst emission processes. We discuss the numeri…
▽ More
Generation of radio bursts in CME foreshock regions and turbulent cascades in the solar wind are assumed to be results of three-wave interaction processes of dispersive plasma modes. Using our Particle in Cell code ACRONYM, we have studied the behaviour of kinetic wavemodes in the presence of beamed electron populations, with a focus on type II radio burst emission processes. We discuss the numerical challenges in generating and analyzing self-consistently evolving wave coupling processes with a PiC-Code and present preliminary results of said project.
△ Less
Submitted 30 July, 2012;
originally announced July 2012.
-
Emission of Type II Radio Bursts - Single-Beam versus Two-Beam Scenario
Authors:
Urs Ganse,
Patrick Kilian,
Rami Vainio,
Felix Spanier
Abstract:
The foreshock region of a CME shock front, where shock accelerated electrons form a beam population in the otherwise quiescent plasma is generally assumed to be the source region of type II radio bursts. Nonlinear wave interaction of electrostatic waves excited by the beamed electrons are the prime candidates for the radio waves' emission.
To address the question whether a single, or two counter…
▽ More
The foreshock region of a CME shock front, where shock accelerated electrons form a beam population in the otherwise quiescent plasma is generally assumed to be the source region of type II radio bursts. Nonlinear wave interaction of electrostatic waves excited by the beamed electrons are the prime candidates for the radio waves' emission.
To address the question whether a single, or two counterpropagating beam populations are a requirement for this process, we have conducted 2.5D particle in cell simulations using the fully relativistic ACRONYM code.
Results show indications of three wave interaction leading to electromagnetic emission at the fundamental and harmonic frequency for the two-beam case. For the single-beam case, no such signatures were detectable.
△ Less
Submitted 25 June, 2012;
originally announced June 2012.
-
Nonlinear Wave Interactions as Emission Process of Type II Radio Bursts
Authors:
Urs Ganse,
Patrick Kilian,
Felix Spanier,
Rami Vainio
Abstract:
The emission of fundamental and harmonic frequency radio waves of type II radio bursts are assumed to be products of three-wave interaction processes of beam-excited Langmuir waves. Using a particle-in-cell code, we have performed simulations of the assumed emission region, a CME foreshock with two counterstreaming electron beams. Analysis of wavemodes within the simulation shows self-consistent e…
▽ More
The emission of fundamental and harmonic frequency radio waves of type II radio bursts are assumed to be products of three-wave interaction processes of beam-excited Langmuir waves. Using a particle-in-cell code, we have performed simulations of the assumed emission region, a CME foreshock with two counterstreaming electron beams. Analysis of wavemodes within the simulation shows self-consistent excitation of beam driven modes, which yield interaction products at both fundamental and harmonic emission frequencies. Through variation of the beam strength, we have investigated the dependence of energy transfer into electrostatic and electromagnetic modes, confirming the quadratic dependence of electromagnetic emission on electron beam strength.
△ Less
Submitted 2 April, 2012;
originally announced April 2012.