Source author record

Patrick Kilian

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

14 published item(s)

preprint2022arXiv

Particle Acceleration in Magnetic Reconnection with Ad hoc Pitch-angle Scattering

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.

preprint2020arXiv

Exploring the acceleration mechanisms for particle injection and power-law formation during trans-relativistic magnetic reconnection

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.

preprint2020arXiv

Fermi-type particle acceleration from magnetic reconnection at the termination shock of a relativistic striped wind

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).

preprint2020arXiv

Ion acceleration in non-relativistic quasi-parallel shocks using fully kinetic simulations

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.

preprint2020arXiv

Recent Progress on Particle Acceleration and Reconnection Physics during Magnetic Reconnectionin the Magnetically-dominated Relativistic Regime

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.

preprint2020arXiv

Studying particle acceleration from driven magnetic reconnection at the termination shock of a relativistic striped wind using particle-in-cell simulations

A rotating pulsar creates a surrounding pulsar wind nebula (PWN) by steadily releasing an energetic wind into the interior of the expanding shockwave of supernova remnant or interstellar medium. At the termination shock of a PWN, the Poynting-flux-dominated relativistic striped wind is compressed. Magnetic reconnection is driven by the compression and converts magnetic energy into particle kinetic energy and accelerating particles to high energies. We carrying out particle-in-cell (PIC) simulations to study the shock structure as well as the energy conversion and particle acceleration mechanism. By analyzing particle trajectories, we find that many particles are accelerated by Fermi-type mechanism. The maximum energy for electrons and positrons can reach hundreds of TeV.

preprint2020arXiv

Time-step dependent force interpolation scheme for suppressing numerical Cherenkov instability in relativistic particle-in-cell simulations

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.

preprint2015arXiv

Energy loss of intergalactic pair beams: Particle-in-Cell simulation

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.

preprint2015arXiv

PICPANTHER: A simple, concise implementation of the relativistic moment implicit Particle-in-Cell method

A three-dimensional, parallelized implementation of the electromagnetic relativistic moment implicit particle-in-cell method in Cartesian geometry (Noguchi et. al., 2007) is presented. Particular care was taken to keep the C++11 codebase simple, concise, and approachable. GMRES is used as a field solver and during the Newton-Krylov iteration of the particle pusher. Drifting Maxwellian problem setups are available while more complex simulations can be implemented easily. Several test runs are described and the code's numerical and computational performance is examined. Weak scaling on the SuperMUC system is discussed and found suitable for large-scale production runs.

preprint2014arXiv

Fundamental and harmonic plasma emission in different plasma environments

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.

preprint2012arXiv

Emission of Type II Radio Bursts - Single-Beam versus Two-Beam Scenario

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.

preprint2012arXiv

Nonlinear Wave Interactions as Emission Process of Type II Radio Bursts

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.

preprint2012arXiv

Note on the use of Yee-lattices in (semi-) implicit Particle-in-cell codes

A modification of the implicit algorithm for particle-in-cell simulations proposed by Petrov and Davis [2011] is presented. The original lattice arrangement is not inherently divergence-free, possibly leading to unphysical results. This arrangement is replaced by a staggered mesh resulting in a reduction of the divergence of the magnetic field by several orders of magnitude.

preprint2012arXiv

Numerical challenges in kinetic simulations of three-wave interactions

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.