A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm

A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different k-space operators in Faraday’s and Ampere’s law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. The new algorithm was implemented into Osiris by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, eA/mc2, exceeding 104 for typical cell sizes and time steps.info:eu-repo/semantics/publishedVersio

A Particle Module for the PLUTO code: II - Hybrid Framework for Modeling Non-thermal emission from Relativistic Magnetized flows

We describe a new hybrid framework to model non-thermal spectral signatures from highly energetic particles embedded in a large-scale classical or relativistic MHD flow. Our method makes use of \textit{Lagrangian} particles moving through an Eulerian grid where the (relativistic) MHD equations are solved concurrently. Lagrangian particles follow fluid streamlines and represent ensembles of (real) relativistic particles with a finite energy distribution. The spectral distribution of each particle is updated in time by solving the relativistic cosmic ray transport equation based on local fluid conditions. This enables us to account for a number of physical processes, such as adiabatic expansion, synchrotron and inverse Compton emission. An accurate semi-analytically numerical scheme that combines the method of characteristics with a Lagrangian discretization in the energy coordinate is described. In presence of (relativistic) magnetized shocks, a novel approach to consistently model particle energization due to diffusive shock acceleration has been presented. Our approach relies on a refined shock-detection algorithm and updates the particle energy distribution based on the shock compression ratio, magnetic field orientation and amount of (parameterized) turbulence. The evolved distribution from each \textit{Lagrangian} particle is further used to produce observational signatures like emission maps and polarization signals accounting for proper relativistic corrections. We further demonstrate the validity of this hybrid framework using standard numerical benchmarks and evaluate the applicability of such a tool to study high energy emission from extra-galactic jets.Comment: 23 pages, 14 figures, Accepted for publication in The Astrophysical Journa

Synthetic radiation diagnostics as a pathway for studying plasma dynamics from advanced accelerators to astrophysical observations

In this thesis, two novel diagnostic techniques for the identi1cation of plasma dynamics and thequanti cation of essential parameters of the dynamics by means of electromagnetic plasmaradiation are presented. Based on particle-in-cell simulations, both the radiation signatures of micrometer-sized laser plasma accelerators and light-year-sized plasma jets are simulated with the same highly parallel radiation simulation framework, in-situ to the plasma simulation. The basics and limits of classical radiation calculation, as well as the theoretical and technical foundation of modern plasma simulation using the particle-in-cell method, are brie2y introduced. The combination of previously independent methods in an in-situ analysis code as well as its validation and extension with newly developed algorithms for the simultaneous quantitative prediction of both coherent and incoherent radiation and the prevention of numerical artifacts is outlined in the initial chapters. For laser wake1eld acceleration, a hitherto unknown off-axis beam signature is observed,which can be used to identify the so-called blowout regime during laser defocusing. Since signi cant radiation is emitted only after the minimum spot size is reached, this signature is ideally suited to determine the laser focus position itself in the plasma to below 100 _m and thus to quantify the in2uence of relativistic self-focusing. A simple semi-analytical scattering model was developed to explain the blowout radiation signature. The spectral signature predicted by the model is veri1ed using both a large-scale explorative simulation and a simulation parameter study, based on an experiment conducted at the HZDR. Identi1ed by the simulations, a temporal asymmetry in the scattered laser light, which cannot be described by state of the art quasi-static models of the blowout regime, makes it possible to determine the focus position precisely by using this radiation signature

Study of transport of laser-driven relativistic electrons in solid materials

With the ultra intense lasers available today, it is possible to generate very hot electron beams in solid density materials. These intense laser-matter interactions result in many applications which include the generation of ultrashort secondary sources of particles and radiation such as ions, neutrons, positrons, x-rays, or even laser-driven hadron therapy. For these applications to become reality, a comprehensive understanding of laser-driven energy transport including hot electron generation through the various mechanisms of ionization, and their subsequent transport in solid density media is required. This study will focus on the characterization of electron transport effects in solid density targets using the state-of- the-art particle-in-cell code PICLS. A number of simulation results will be presented on the topics of ionization propagation in insulator glass targets, non-equilibrium ionization mod- eling featuring electron impact ionization, and electron beam guiding by the self-generated resistive magnetic field. An empirically derived scaling relation for the resistive magnetic in terms of the laser parameters and material properties is presented and used to derive a guiding condition. This condition may prove useful for the design of future laser-matter interaction experiments