Source author record

L. O. Silva

L. O. Silva 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

64works
12topics
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

64 published item(s)

preprint2024arXiv

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell's equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly-rotating magnetised neutron stars taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. $R_s \sim 0.5 [R_*]$, we show that the new terms lead to a percent order correction in the magnetic field orientation and strength compared to the case with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

preprint2022arXiv

Laser-Driven, Ion-Scale Magnetospheres in Laboratory Plasmas. I. Experimental Platform and First Results

Magnetospheres are a ubiquitous feature of magnetized bodies embedded in a plasma flow. While large planetary magnetospheres have been studied for decades by spacecraft, ion-scale "mini" magnetospheres can provide a unique environment to study kinetic-scale, collisionless plasma physics in the laboratory to help validate models of larger systems. In this work, we present preliminary experiments of ion-scale magnetospheres performed on a unique high-repetition-rate platform developed for the Large Plasma Device (LAPD) at UCLA. The experiments utilize a high-repetition-rate laser to drive a fast plasma flow into a pulsed dipole magnetic field embedded in a uniform magnetized background plasma. 2D maps of magnetic field with high spatial and temporal resolution are measured with magnetic flux probes to examine the evolution of magnetosphere and current density structures for a range of dipole and upstream parameters. The results are further compared to 2D PIC simulations to identify key observational signatures of the kinetic-scale structures and dynamics of the laser-driven plasma. We find that distinct 2D kinetic-scale magnetopause and diamagnetic current structures are formed at higher dipole moments, and their locations are consistent with predictions based on pressure balances and energy conservation.

preprint2021arXiv

Analysis of Proton Bunch Parameters in the AWAKE Experiment

A precise characterization of the incoming proton bunch parameters is required to accurately simulate the self-modulation process in the Advanced Wakefield Experiment (AWAKE). This paper presents an analysis of the parameters of the incoming proton bunches used in the later stages of the AWAKE Run 1 data-taking period. The transverse structure of the bunch is observed at multiple positions along the beamline using scintillating or optical transition radiation screens. The parameters of a model that describes the bunch transverse dimensions and divergence are fitted to represent the observed data using Bayesian inference. The analysis is tested on simulated data and then applied to the experimental data.

preprint2021arXiv

Slowdown of interpenetration of two counterpropagating plasma slab due to collective effects

The nonlinear evolution of electromagnetic instabilities driven by the interpenetration of two $e^-\,e^+$ plasma clouds is explored using {\it ab initio} kinetic plasma simulations. We show that the plasma clouds slow down due to both oblique and Weibel generated electromagnetic fields, which deflect the particle trajectories, transferring bulk forward momentum into transverse momentum and thermal velocity spread. This process causes the flow velocity $v_{inst}$ to decrease approximately by a factor of $\sqrt{1/3}$ in a time interval $Δt_{αB} ω_p \sim c/(v_{fl}\sqrt{α_B})$, where $α_B$ is the magnetic equipartition parameter determined by the non-linear saturation of the instabilities, $v_{fl}$ is the initial flow speed, and $ω_p$ is the plasma frequency. For the $α_B$ measured in our simulations, $Δt_{αB}$ is close to $10 \times$ the instability growth time. We show that as long as the plasma slab length $L > v_{fl} Δt_{αB}$, the plasma flow is expected to slow down by a factor close to $1/\sqrt{3}$.

preprint2020arXiv

Experimental study of extended timescale dynamics of a plasma wakefield driven by a self-modulated proton bunch

Plasma wakefield dynamics over timescales up to 800 ps, approximately 100 plasma periods, are studied experimentally at the Advanced Wakefield Experiment (AWAKE). The development of the longitudinal wakefield amplitude driven by a self-modulated proton bunch is measured using the external injection of witness electrons that sample the fields. In simulation, resonant excitation of the wakefield causes plasma electron trajectory crossing, resulting in the development of a potential outside the plasma boundary as electrons are transversely ejected. Trends consistent with the presence of this potential are experimentally measured and their dependence on wakefield amplitude are studied via seed laser timing scans and electron injection delay scans.

preprint2020arXiv

High-order harmonic generation in an electron-positron-ion plasma

The laser interaction with an electron-positron-ion mixed plasma is studied, from the perspective of the associated high-order harmonic generation. For an idealized mixed plasma which is assumed with a sharp plasma-vacuum interface and uniform density distribution, when it is irradiated by a weakly relativistic laser pulse, well-defined signals at harmonics of the plasma frequency in the harmonic spectrum are observed. These characteristic signals are attributed to the inverse two-plasmon decay of the counterpropagating monochromatic plasma waves which are excited by the energetic electrons and the positron beam accelerated by the laser. Particle-in-cell simulations show the signal at twice the plasma frequency can be observed for a pair density as low as $\sim 10^{-5}$ of the plasma density. In the self-consistent scenario of pair production by an ultraintense laser striking a solid target, particle-in-cell simulations, which account for quantum electrodynamic effects (photon emission and pair production), show that dense (greater than the relativistically-corrected critical density) and hot pair plasmas can be created. The harmonic spectrum shows weak low order harmonics, indicating a high laser absorption due to quantum electrodynamic effects. The characteristic signals at harmonics of the plasma frequency are absent, because broadband plasma waves are excited due to the high plasma inhomogeneity introduced by the interaction. However, the high frequency harmonics are enhanced due to the high-frequency modulations from the direct laser coupling with created pair plasmas.

preprint2020arXiv

Multidimensional Plasma Wake Excitation in the Non-linear Blowout Regime

Plasma accelerators can sustain very high acceleration gradients. They are promising candidates for future generations of particle accelerators for several scientific, medical and technological applications. Current plasma based acceleration experiments operate in the relativistic regime, where the plasma response is strongly non-linear. We outline some of the key properties of wakefield excitation in these regimes. We outline a multidimensional theory for the excitation of plasma wakefields in connection with current experiments. We then use these results and provide design guidelines for the choice of laser and plasma parameters ensuring a stable laser wakefield accelerator that maximizes the quality of the accelerated electrons. We also mention some of the future challenges associated with this technology.

preprint2020arXiv

Proton beam defocusing in AWAKE: comparison of simulations and measurements

In 2017, AWAKE demonstrated the seeded self-modulation (SSM) of a 400 GeV proton beam from the Super Proton Synchrotron (SPS) at CERN. The angular distribution of the protons deflected due to SSM is a quantitative measure of the process, which agrees with simulations by the two-dimensional (axisymmetric) particle-in-cell code LCODE. Agreement is achieved for beam populations between $10^{11}$ and $3 \times 10^{11}$ particles, various plasma density gradients ($-20 ÷20\%$) and two plasma densities ($2\times 10^{14} \text{cm}^{-3}$ and $7 \times 10^{14} \text{cm}^{-3}$). The agreement is reached only in the case of a wide enough simulation box (at least five plasma wavelengths).

preprint2020arXiv

The effects of collisions on the generation and suppression of temperature anisotropies and the Weibel instability

The expansion of plasma with non-parallel temperature and density gradients, and the generation of magnetic field via the Biermann battery is modeled using particle-in-cell simulations that include collisional effects via Monte Carlo methods. A scaling of the degree of collisionality shows that an anisotropy can be produced, and drive the Weibel instability, for gradient scales shorter than the mean free path. For larger collision rates, the Biermann battery dominates as the cause of magnetic field generation. When the most energetic particles remain collisionless, the Nernst effect causes the Biermann field to be dragged with the heat flux, piled up, and enhanced.

preprint2019arXiv

Interplay between the Weibel instability and the Biermann battery in realistic laser-solid interactions

A novel setup allows the Weibel instability and its interplay with the Biermann battery to be probed in laser-driven collisionless plasmas. Ab initio particle-in-cell (PIC) simulations of the interaction of short ($ \le ps$) intense $(a_0 \ge 1)$ laser-pulses with overdense plasma targets show observable Weibel generated magnetic fields. This field strength surpasses that of the Biermann battery, usually dominant in experiments, as long as the gradient scale length is much larger than the local electron inertial length; this is achievable by carefully setting the appropriate gradients in the front of the target e.g. by tuning the delay between the main laser pulse and the pre-pulse.

preprint2018arXiv

Magnetized Current Filaments as a Source of Circularly Polarized Light

We show that the Weibel or currente filamentation instability can lead to the emission of circularly polarized radiation. Using particle-in-cell (PIC) simulations and a radiation post-processing numerical algorithm, we demonstrate that the level of circular polarization increases with the initial plasma magnetization, saturating at ~13% when the magnetization, given by the ratio of magnetic energy density to the electron kinetic energy density, is larger than 0.05. Furthermore, we show that this effect requires an ion-electron mass ratio greater than unity. These findings, which could also be tested in currently available laboratory conditions, show that the recent observation of circular polarization in gamma ray burst afterglows could be attributed to the presence of magnetized current filaments driven by the Weibel or the current filamentation instability.

preprint2016arXiv

Amplification and generation of ultra-intense twisted laser pulses via stimulated Raman scattering

Twisted Laguerre-Gaussian lasers, with orbital angular momentum and characterised by doughnut shaped intensity profiles, provide a transformative set of tools and research directions in a growing range of fields and applications, from super-resolution microcopy and ultra-fast optical communications to quantum computing and astrophysics. The impact of twisted light is widening as recent numerical calculations provided solutions to long-standing challenges in plasma-based acceleration by allowing for high gradient positron acceleration. The production of ultrahigh intensity twisted laser pulses could then also have a broad influence on relativistic laser-matter interactions. Here we show theoretically and with ab-initio three-dimensional particle-in-cell simulations, that stimulated Raman backscattering can generate and amplify twisted lasers to Petawatt intensities in plasmas. This work may open new research directions in non-linear optics and high energy density science, compact plasma based accelerators and light sources.

preprint2016arXiv

Essential criteria for efficient pulse amplification via Raman and Brillouin scattering

Raman and Brillouin amplification are two schemes for amplifying and compressing short laser pulses in plasma. Analytical models have already been derived for both schemes, but the full consequences of these models are little known or used. Here, we present new criteria that govern the evolution of the attractor solution for the seed pulse in Raman and Brillouin amplification, and show how the initial laser pulses need to be shaped to control the properties of the final amplified seed and improve the amplification efficiency.

preprint2016arXiv

High orbital angular momentum harmonic generation

We identify and explore a high orbital angular momentum (OAM) harmonics generation and amplification mechanism that manipulates the OAM independently of any other laser property, by preserving the initial laser wavelength, through stimulated Raman backscattering in a plasma. The high OAM harmonics spectra can extend at least up to the limiting value imposed by the paraxial approximation. We show with theory and particle-in-cell simulations that the orders of the OAM harmonics can be tuned according to a selection rule that depends on the initial OAM of the interacting waves. We illustrate the high OAM harmonics generation in a plasma using several examples including the generation of prime OAM harmonics. The process can also be realised in any nonlinear optical Kerr media supporting three-wave interactions.

preprint2016arXiv

Physical mechanism of the transverse instability in radiation pressure ion acceleration

The transverse stability of the target is crucial for obtaining high quality ion beams using the laser radiation pressure acceleration (RPA) mechanism. In this letter, a theoretical model and supporting two-dimensional (2D) Particle-in-Cell (PIC) simulations are presented to clarify the physical mechanism of the transverse instability observed in the RPA process. It is shown that the density ripples of the target foil are mainly induced by the coupling between the transverse oscillating electrons and the quasi-static ions, a mechanism similar to the transverse two stream instability in the inertial confinement fusion (ICF) research. The predictions of the mode structure and the growth rates from the theory agree well with the results obtained from the PIC simulations in various regimes, indicating the model contains the essence of the underlying physics of the transverse break-up of the target.

preprint2016arXiv

The generation of magnetic fields by the Biermann battery and the interplay with the Weibel instability

An investigation of magnetic fields generated in an expanding bubble of plasma with misaligned temperature and density gradients (driving the Biermann battery mechanism) is performed. With gradient scales $L$, large-scale magnetic fields are generated by the Biermann battery mechanism with plasma $β\sim 1$, as long as $L$ is comparable to the ion inertial length $d_i$. For larger system sizes, $L/d_e > 100$ (where $d_e$ is the electron inertial length), the Weibel instability generates magnetic fields of similar magnitude but with wavenumber $k d_e \sim 0.2$. In both cases, the growth and saturation of these fields have a weak dependence on mass ratio $m_i/m_e$, indicating electron mediated physics. A scan in system size is performed at $m_i/m_e = 2000$, showing agreement with previous results with $m_i/m_e = 25$. In addition, the instability found at large system sizes is quantitatively demonstrated to be the Weibel instability. Furthermore, magnetic and electric energy spectra at scales below the electron Larmor radius are found to exhibit power law behavior with spectral indices $-16/3$ and $-7/3$, respectively.

preprint2015arXiv

Modelling radiation emission in the transition from the classical to the quantum regime

An emissivity formula is derived using the generalised Fermi-Weizacker-Williams method of virtual photons which accounts for the recoil the charged particle experiences as it emits radiation. It is found that through this derivation the formula obtained by Sokolov et al using QED perturbation theory is recovered. The corrected emissivity formula is applied to nonlinear Thomson scattering scenarios in the transition from the classical to the quantum regime, for small values of the nonlinear quantum parameter χ. Good agreement is found between this method and a QED probabilistic approach for scenarios where both are valid. In addition, signatures of the quantum corrections are identified and explored.

preprint2015arXiv

Persistence of magnetic field driven by relativistic electrons in a plasma

The onset and evolution of magnetic fields in laboratory and astrophysical plasmas is determined by several mechanisms, including instabilities, dynamo effects and ultra-high energy particle flows through gas, plasma and interstellar-media. These processes are relevant over a wide range of conditions, from cosmic ray acceleration and gamma ray bursts to nuclear fusion in stars. The disparate temporal and spatial scales where each operates can be reconciled by scaling parameters that enable to recreate astrophysical conditions in the laboratory. Here we unveil a new mechanism by which the flow of ultra-energetic particles can strongly magnetize the boundary between the plasma and the non-ionized gas to magnetic fields up to 10-100 Tesla (micro Tesla in astrophysical conditions). The physics is observed from the first time-resolved large scale magnetic field measurements obtained in a laser wakefield accelerator. Particle-in-cell simulations capturing the global plasma and field dynamics over the full plasma length confirm the experimental measurements. These results open new paths for the exploration and modelling of ultra high energy particle driven magnetic field generation in the laboratory.

preprint2015arXiv

Physics of collisionless shocks - theory and simulation

Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical environments and they are mediated by the Weibel or filamentation instability. In such shocks, charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are characterized by a strong electrostatic field, which leads to electron trapping. Ions are accelerated by reflection from the electrostatic potential. Shock formation and particle acceleration will be discussed in theory and simulations.

preprint2015arXiv

Slow down of a globally neutral relativistic $e^-e^+$ beam shearing the vacuum

The microphysics of relativistic collisionless sheared flows is investigated in a configuration consisting of a globally neutral, relativistic $e^-e^+$ beam streaming through a hollow plasma/dielectric channel. We show through multidimensional PIC simulations that this scenario excites the Mushroom instability (MI), a transverse shear instability on the electron-scale, when there is no overlap (no contact) between the $e^-e^+$ beam and the walls of the hollow plasma channel. The onset of the MI leads to the conversion of the beam's kinetic energy into magnetic (and electric) field energy, effectively slowing down a globally neutral body in the absence of contact. The collisionless shear physics explored in this configuration may operate in astrophysical environments, particularly in highly relativistic and supersonic settings where macroscopic shear processes are stable.

preprint2015arXiv

Spatial-temporal evolution of the current filamentation instability

The spatial-temporal evolution of the purely transverse current filamentation instability is analyzed by deriving a single partial differential equation for the instability and obtaining the analytical solutions for the spatially and temporally growing current filament mode. When the beam front always encounters fresh plasma, our analysis shows that the instability grows spatially from the beam front to the back up to a certain critical beam length; then the instability acquires a purely temporal growth. This critical beam length increases linearly with time and in the non-relativistic regime it is proportional to the beam velocity. In the relativistic regime the critical length is inversely proportional to the cube of the beam Lorentz factor $γ_{0b}$. Thus, in the ultra-relativistic regime the instability immediately acquires a purely temporal growth all over the beam. The analytical results are in good agreement with multidimensional particle-in-cell simulations performed with OSIRIS. Relevance of current study to recent and future experiments on fireball beams is also addressed.

preprint2015arXiv

Study on Coulomb explosions of ion mixtures

The paper presents a theoretical work on the dynamics of Coulomb explosion for spherical nanoplasmas composed by two different ion species. Particular attention has been dedicated to study the energy spectra of the ions with the larger charge-to-mass ratio. The connection between the formation of shock shells and the energy spread of the ions has been the object of a detailed analysis, showing that under particular conditions the width of the asymptotic energy spectrum tends to become very narrow, which leads to a multi-valued ion phase-space. The conditions to generate a quasi mono-energetic ion spectrum have been rigorously demonstrated and verifed by numerical simulations, using a technique that, exploiting the spherical symmetry of the problem, allows one to obtain very accurate and precise results.

preprint2015arXiv

Transverse electron-scale instability in relativistic shear flows

Electron-scale surface waves are shown to be unstable in the transverse plane of a shear flow in an initially unmagnetized plasma, unlike in the (magneto)hydrodynamics case. It is found that these unstable modes have a higher growth rate than the closely related electron-scale Kelvin-Helmholtz instability in relativistic shears. Multidimensional particle-in-cell simulations verify the analytic results and further reveal the emergence of mushroom-like electron density structures in the nonlinear phase of the instability, similar to those observed in the Rayleigh Taylor instability despite the great disparity in scales and different underlying physics. Macroscopic ($\gg c/ω_{pe}$) fields are shown to be generated by these microscopic shear instabilities, which are relevant for particle acceleration, radiation emission and to seed MHD processes at long time-scales.

preprint2014arXiv

Beam quality requirements for the Ion-Channel Laser

In this paper, we determine the electron beam quality requirements to obtain exponential radiation amplification in the ion-channel laser, where a relativistic electron beam wiggles in a focusing ion-channel that can be created in a wakefield accelerator. The beam energy and wiggler parameter spreads should be limited. Those spread limits are functions of the Pierce parameter, which is calculated here without neglecting the radiation diffraction. Two dimensional and three dimensional simulations of the self-consistent ion-channel laser confirm our theoretical predictions.

preprint2014arXiv

Boosting the performance of Brillouin amplification at sub-quarter-critical densities via reduction of parasitic Raman scattering

Raman and Brillouin amplification of laser pulses in plasma have been shown to produce picosecond pulses of petawatt power. In previous studies, filamentation of the probe pulse has been identified as the biggest threat to the amplification process, especially for Brillouin amplification, which employs the highest plasma densities. Therefore it has been proposed to perform Brillouin scattering at densities below $n_{cr}/4$ to reduce the influence of filamentation. However, parastic Raman scattering can become a problem at such densities, contrary to densities above $n_{cr}/4$, where it is suppressed. In this paper, we investigate the influence of parasitic Raman scattering on Brillouin amplification at densities below $n_{cr}/4$. We expose the specific problems posed by both Raman backward and forward scattering, and how both types of scattering can be mitigated, leading to an increased performance of the Brillouin amplification process.

preprint2014arXiv

Electromagnetic field generation in the downstream of electrostatic shocks due to electron trapping

A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $ \tilde \, 10^4 \, ω_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instability. The generated magnetic field is well-confined to the downstream region of the electrostatic shock. The shock formation process is not modified and the features of the shock front responsible for ion acceleration, which are currently probed in laser-plasma laboratory experiments, are maintained. However, such a strong magnetic field determines the particle trajectories downstream and has the potential to modify the signatures of the collisionless shock.

preprint2014arXiv

Electron trapping and acceleration by the plasma wakefield of a self-modulating proton beam

It is shown that co-linear injection of electrons or positrons into the wakefield of the self-modulating particle beam is possible and ensures high energy gain. The witness beam must co-propagate with the tail part of the driver, since the plasma wave phase velocity there can exceed the light velocity, which is necessary for efficient acceleration. If the witness beam is many wakefield periods long, then the trapped charge is limited by beam loading effects. The initial trapping is better for positrons, but at the acceleration stage a considerable fraction of positrons is lost from the wave. For efficient trapping of electrons, the plasma boundary must be sharp, with the density transition region shorter than several centimeters. Positrons are not susceptible to the initial plasma density gradient.

preprint2014arXiv

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

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

preprint2014arXiv

Exploring the nature of collisionless shocks under laboratory conditions

Collisionless shocks are pervasive in astrophysics and they are critical to understand cosmic ray acceleration. Laboratory experiments with intense lasers are now opening the way to explore and characterise the underlying microphysics, which determine the acceleration process of collisionless shocks. We determine the shock character - electrostatic or electromagnetic - based on the stability of electrostatic shocks to transverse electromagnetic fluctuations as a function of the electron temperature and flow velocity of the plasma components, and we compare the analytical model with particle-in-cell simulations. By making the connection with the laser parameters driving the plasma flows, we demonstrate that shocks with different and distinct underlying microphysics can be explored in the laboratory with state-of-the-art laser systems.

preprint2014arXiv

Implementation of a hybrid particle code with a PIC description in r-z and a gridless description in $ϕ$ into OSIRIS

For many plasma physics problems, three-dimensional and kinetic effects are very important. However, such simulations are very computationally intensive. Fortunately, there is a class of problems for which there is nearly azimuthal symmetry and the dominant three-dimensional physics is captured by the inclusion of only a few azimuthal harmonics. Recently, it was proposed [A. Lifschitz et al., J. Comp. Phys. 228 (5) (2009) 1803-1814] to model one such problem, laser wakefield acceleration, by expanding the fields and currents in azimuthal harmonics and truncating the expansion after only the first harmonic. The complex amplitudes of the fundamental and first harmonic for the fields were solved on an r-z grid and a procedure for calculating the complex current amplitudes for each particle based on its motion in Cartesian geometry was presented using a Marder's correction to maintain the validity of Gauss's law. In this paper, we describe an implementation of this algorithm into OSIRIS using a rigorous charge conserving current deposition method to maintain the validity of Gauss's law. We show that this algorithm is a hybrid method which uses a particles-in-cell description in r-z and a gridless description in $ϕ$. We include the ability to keep an arbitrary number of harmonics and higher order particle shapes. Examples, for laser wakefield acceleration, plasma wakefield acceleration, and beam loading are also presented and directions for future work are discussed.

preprint2014arXiv

Ion motion in the wake driven by long particle bunches in plasmas

We explore the role of the background plasma ion motion in self-modulated plasma wakefield accelerators. We employ J. Dawson's plasma sheet model to derive expressions for the transverse plasma electric field and ponderomotive force in the narrow bunch limit. We use these results to determine the on-set of the ion dynamics, and demonstrate that the ion motion could occur in self-modulated plasma wakefield accelerators. Simulations show the motion of the plasma ions can lead to the early suppression of the self-modulation instability and of the accelerating fields. The background plasma ion motion can nevertheless be fully mitigated by using plasmas with heavier plasmas.

preprint2014arXiv

Magnetic field generation and amplification in an expanding plasma

Particle-in-cell simulations are used to investigate the formation of magnetic fields, B, in plasmas with perpendicular electron density and temperature gradients. For system sizes, L, comparable to the ion skin depth, d_i, it is shown that B ~ d_i/L, consistent with the Biermann battery effect. However, for large L/d_i, it is found that the Weibel instability (due to electron temperature anisotropy) supersedes the Biermann battery as the main producer of B. The Weibel-produced fields saturate at a finite amplitude (plasma β\approx 100), independent of L. The magnetic energy spectra below the electron Larmor radius scale is well fitted by power law with slope -16/3, as predicted in Schekochihin et al., Astrophys. J. Suppl. Ser 182, 310 (2009).

preprint2014arXiv

Plasma wakefields driven by intense, broadband, incoherent electromagnetic radiation

Non-linear wave-driven processes in plasmas are normally described by either a monochromatic pump wave that couples to other monochromatic waves, or as a random phase wave coupling to other random phase waves. An alternative approach involves an incoherent, random or broadband pump coupling to monochromatic and/or coherent structures in the plasma. This approach can be implemented through the wave kinetic model. In this model, the incoming pump wave is described by either a bunch (for coherent waves) or a sea (for random phase waves) of quasi-particles. A particle-in-cell type code has been developed to perform numerical simulations of such interactions using the quasi-particle approach. This code allows for a comparatively easy description of both random phase and coherent pump pulses coupling to slow electrostatic plasma waves, while providing an extended range of powerful diagnostics leading to a deeper physical insight into the dynamics of the fast waves. As an example, the propagation of short, intense laser pulses through a plasma has been simulated. A sample of the phenomena that can be studied for this case includes modulational instabilities, photon wave breaking and turbulence, and pulse compression, stretching, and chirping. The extensibility of the numerical implementation to other types of fast wave-slow wave interactions is also discussed.

preprint2014arXiv

Proton-driven plasma wakefield acceleration: a path to the future of high-energy particle physics

New acceleration technology is mandatory for the future elucidation of fundamental particles and their interactions. A promising approach is to exploit the properties of plasmas. Past research has focused on creating large-amplitude plasma waves by injecting an intense laser pulse or an electron bunch into the plasma. However, the maximum energy gain of electrons accelerated in a single plasma stage is limited by the energy of the driver. Proton bunches are the most promising drivers of wakefields to accelerate electrons to the TeV energy scale in a single stage. An experimental program at CERN -- the AWAKE experiment -- has been launched to study in detail the important physical processes and to demonstrate the power of proton-driven plasma wakefield acceleration. Here we review the physical principles and some experimental considerations for a future proton-driven plasma wakefield accelerator.

preprint2014arXiv

Self-modulation instability of ultra-relativistic particle bunches with finite rise times

We study the evolution of the self-modulation instability using bunches with finite rise times. Using particle-in-cell simulations we show that unlike long bunches with sharp rise times, there are large variations of the wake amplitudes and wake phase velocity when bunches with finite rise times are used. These results show that use of bunches with sharp rise times is important to seed the self-modulation instability and to ensure stable acceleration regimes.

preprint2014arXiv

SEP acceleration in CME driven shocks using a hybrid code

We preform hybrid simulations of super Alfvenic quasi-parallel shock, driven by a Coronal Mass Ejection (CME), propagating in the Outer Coronal or Solar Wind at distances of between 3 to 6 solar radii. The hybrid treatment of the problem enable the study of the shock propagation on the ion time scale, preserving ion kinetics and allowing for a self consistent treatment of the shock propagation and particle acceleration. The CME plasma drags the embedded magnetic field lines stretching from the sun, and propagates out into interplanetary space at a greater velocity than the in-situ solar wind, driving the shock, and producing very energetic particles. Our results show electromagnetic Alfven waves are generated at the shock front. The waves propagate upstream of the shock and are produced by the counter streaming ions of the solar wind plasma being reflected at the shock. A significant fraction of the particles are accelerated in two distinct phases: first, particles drift from the shock and are accelerated in the upstream region and, second, particles arriving at the shock get trapped, and are accelerated at the shock front. A fraction of the particles diffused back to the shock, which is consistent with the Fermi acceleration mechanism.

preprint2013arXiv

Collisionless shock formation, spontaneous electromagnetic fluctuations and streaming instabilities

Collisionless shocks are ubiquitous in astrophysics and in the lab. Recent numerical simulations and experiments have shown how they can arise from the encounter of two collisionless plasma shells. When the shells interpenetrate, the overlapping region turns unstable, triggering the shock formation. As a first step towards a microscopic understanding of the process, we analyze here in detail the initial instability phase. On the one hand, 2D relativistic PIC simulations are performed where two symmetric initially cold pair plasmas collide. On the other hand, the instabilities at work are analyzed, as well as the field at saturation and the seed field which gets amplified. For mildly relativistic motions and onward, Weibel modes govern the linear phase. We derive an expression for the duration of the linear phase in good agreement with the simulations. This saturation time constitutes indeed a lower-bound for the shock formation time.

preprint2013arXiv

Ion acceleration from laser-driven electrostatic shocks

Multi-dimensional particle-in-cell simulations are used to study the generation of electrostatic shocks in plasma and the reflection of background ions to produce high-quality and high-energy ion beams. Electrostatic shocks are driven by the interaction of two plasmas with different density and/or relative drift velocity. The energy and number of ions reflected by the shock increase with increasing density ratio and relative drift velocity between the two interacting plasmas. It is shown that the interaction of intense lasers with tailored near-critical density plasmas allows for the efficient heating of the plasma electrons and steepening of the plasma profile at the critical density interface, leading to the generation of high-velocity shock structures and high-energy ion beams. Our results indicate that high-quality 200 MeV shock-accelerated ion beams required for medical applications may be obtained with current laser systems.

preprint2013arXiv

Laser-plasma interactions for fast ignition

In the electron-driven fast-ignition approach to inertial confinement fusion, petawatt laser pulses are required to generate MeV electrons that deposit several tens of kilojoules in the compressed core of an imploded DT shell. We review recent progress in the understanding of intense laser plasma interactions (LPI) relevant to fast ignition. Increases in computational and modeling capabilities, as well as algorithmic developments have led to enhancement in our ability to perform multi-dimensional particle-in-cell (PIC) simulations of LPI at relevant scales. We discuss the physics of the interaction in terms of laser absorption fraction, the laser-generated electron spectra, divergence, and their temporal evolution. Scaling with irradiation conditions such as laser intensity are considered, as well as the dependence on plasma parameters. Different numerical modeling approaches and configurations are addressed, providing an overview of the modeling capabilities and limitations. In addition, we discuss the comparison of simulation results with experimental observables. In particular, we address the question of surrogacy of today's experiments for the full-scale fast ignition problem.

preprint2013arXiv

Relativistic generalization of formation and ion reflection conditions in electrostatic shocks

The theoretical model by Sorasio et al. (2006) for the steady state Mach number of electrostatic shocks formed in the interaction of two plasma slabs of arbitrary density and temperature is generalized for relativistic electron and non-relativistic ion temperatures. We find that the relativistic correction leads to lower Mach numbers, and as a consequence, ions are reflected with lower energies. The steady state bulk velocity of the downstream population is introduced as an additional parameter to describe the transition between the minimum and maximum Mach numbers in dependence of the initial density and temperature ratios. In order to transform the soliton-like solution in the upstream region into a shock, a population of reflected ions is considered and differences to a zero-ion temperature model are discussed.

preprint2012arXiv

Acceleration in perpendicular relativistic shocks for plasmas consisting of leptons and hadrons

We investigate the acceleration of light particles in perpendicular shocks for plasmas consisting of a mixture of leptonic and hadronic particles. Starting from the full set of conservation equations for the mixed plasma constituents, we generalize the magneto-hydrodynamical jump conditions for a multi-component plasma, including information about the specific adiabatic constants for the different species. The impact of deviations from the standard model of an ideal gas is compared in theory and particle-in-cell simulations, showing that the standard-MHD model is a good approximation. The simulations of shocks in electron-positron-ion plasmas are for the first time multi-dimensional, transverse effects are small in this configuration and 1D simulations are a good representation if the initial magnetization is chosen high. 1D runs with a mass ratio of 1836 are performed, which identify the Larmor frequency ω_{ci} as the dominant frequency that determines the shock physics in mixed component plasmas. The maximum energy in the non-thermal tail of the particle spectra evolves in time according to a power-law proportional to t^αwith αin the range 1/3 < α< 1, depending on the initial parameters. A connection is made with transport theoretical models by Drury (1983) and Gargate & Spitkovsky (2011), which predict an acceleration time proportional to γand the theory for small wavelength scattering by Kirk & Reville (2010), which predicts a behavior rather as proportional to γ^2. Furthermore, we compare different magnetic field orientations with B_0 inside and out of the plane, observing qualitatively different particle spectra than in pure electron-ion shocks.

preprint2012arXiv

DC magnetic field generation in unmagnetized shear flows

The generation of DC magnetic fields in unmagnetized plasmas with velocity shear is predicted for non relativistic and relativistic scenarios either due to thermal effects or due to the onset of the Kelvin-Helmholtz instability (KHI). A kinetic model describes the growth and the saturation of the DC field. The predictions of the theory are confirmed by multidimensional particle-in-cell simulations, demonstrating the formation of long lived magnetic fields ($t \sim 100s ω_{pi}^{-1}$) along the full longitudinal extent of the shear layer, with transverse width on the electron length scale ($\sqrt{γ_0}c/ω_{pe}$), reaching magnitudes $eB_{\mathrm{DC}}/m_ecω_{pe}\sim β_0\sqrt{γ_0}$.

preprint2012arXiv

Efficient modeling of laser plasma interactions in high energy density scenarios

We describe how a new framework for coupling a full-PIC algorithm with a reduced PIC algorithm has been implemented into the code OSIRIS. We show that OSIRIS with this new hybrid-PIC algorithm can efficiently and accurately model high energy density scenarios such as ion acceleration in laser-solid interactions and fast ignition of fusion targets. We model for the first time the full density range of a fast ignition target in a fully self-consistent hybrid-PIC simulation, illustrating the possibility of stopping the laser generated electron flux at the core region with relatively high efficiencies. Computational speedups greater than 1000 times are demonstrated, opening the way for full-scale multi-dimensional modeling of high energy density scenarios and for the guiding of future experiments.

preprint2012arXiv

Influence of realistic parameters on state-of-the-art LWFA experiments

We examine the influence of non-ideal plasma-density and non-Gaussian transverse laser-intensity profiles in the laser wakefield accelerator analytically and numerically. We find that the characteristic amplitude and scale length of longitudinal density fluctuations impacts on the final energies achieved by electron bunches. Conditions that minimize the role of the longitudinal plasma density fluctuations are found. The influence of higher order Laguerre-Gaussian laser pulses is also investigated. We find that higher order laser modes typically lead to lower energy gains. Certain combinations of higher order modes may, however, lead to higher electron energy gains.

preprint2012arXiv

Laser-driven shock acceleration of monoenergetic ion beams

We show that monoenergetic ion beams can be accelerated by moderate Mach number collisionless, electrostatic shocks propagating in a long scale-length exponentially decaying plasma profile. Strong plasma heating and density steepening produced by an intense laser pulse near the critical density can launch such shocks that propagate in the extended plasma at high velocities. The generation of a monoenergetic ion beam is possible due to the small and constant sheath electric field associated with the slowly decreasing density profile. The conditions for the acceleration of high-quality, energetic ion beams are identified through theory and multidimensional particle-in-cell simulations. The scaling of the ion energy with laser intensity shows that it is possible to generate $\sim 200$ MeV proton beams with state-of-the-art 100 TW class laser systems.

preprint2012arXiv

Magnetically assisted self-injection and radiation generation for plasma based acceleration

It is shown through analytical modeling and numerical simulations that external magnetic fields can relax the self-trapping thresholds in plasma based accelerators. In addition, the transverse location where self-trapping occurs can be selected by adequate choice of the spatial profile of the external magnetic field. We also find that magnetic-field assisted self-injection can lead to the emission of betatron radiation at well defined frequencies. This controlled injection technique could be explored using state-of-the-art magnetic fields in current/next generation plasma/laser wakefield accelerator experiments.

preprint2012arXiv

The impact of kinetic effects on the properties of relativistic electron-positron shocks

We assess the impact of non-thermally shock-accelerated particles on the magnetohydrodynamic (MHD) jump conditions of relativistic shocks. The adiabatic constant is calculated directly from first principle particle-in-cell simulation data, enabling a semi-kinetic approach to improve the standard fluid model and allowing for an identification of the key parameters that define the shock structure. We find that the evolving upstream parameters have a stronger impact than the corrections due to non-thermal particles. We find that the decrease of the upstream bulk speed yields deviations from the standard MHD model up to 10%. Furthermore, we obtain a quantitative definition of the shock transition region from our analysis. For Weibel-mediated shocks the inclusion of a magnetic field in the MHD conservation equations is addressed for the first time.

preprint2012arXiv

The ion motion in self-modulated plasma wakefield accelerators

The effects of plasma ion motion in self-modulated plasma based accelerators is examined. An analytical model describing ion motion in the narrow beam limit is developed, and confirmed through multi-dimensional particle-in-cell simulations. It is shown that the ion motion can lead to the early saturation of the self-modulation instability, and to the suppression of the accelerating gradients. This can reduce the total energy that can be transformed into kinetic energy of accelerated particles. For the parameters of future proton-driven plasma accelerator experiments, the ion dynamics can have a strong impact. Possible methods to mitigate the effects of the ion motion in future experiments are demonstrated.

preprint2012arXiv

Three-dimensional simulations of laser-plasma interactions at ultrahigh intensities

Three-dimensional (3D) particle-in-cell (PIC) simulations are used to investigate the interaction of ultrahigh intensity lasers ($> 10^{20}$ W/cm$^{-2}$) with matter at overcritical densities. Intense laser pulses are shown to penetrate up to relativistic critical density levels and to be strongly self-focused during this process. The heat flux of the accelerated electrons is observed to have an annular structure when the laser is tightly focused, showing that a large fraction of fast electrons is accelerated at an angle. These results shed light into the multi-dimensional effects present in laser-plasma interactions of relevance to fast ignition of fusion targets and laser-driven ion acceleration in plasmas.

preprint2012arXiv

Transverse self-modulation of ultra-relativistic lepton beams in the plasma wakefield accelerator

The transverse self-modulation of ultra-relativistic, long lepton bunches in high-density plasmas is explored through full-scale particle-in-cell simulations. We demonstrate that long SLAC-type electron and positron bunches can become strongly self-modulated over centimeter distances, leading to wake excitation in the blowout regime with accelerating fields in excess of 20 GV/m. We show that particles energy variations exceeding 10 GeV can occur in meter-long plasmas. We find that the self-modulation of positively and negatively charged bunches differ when the blowout is reached. Seeding the self-modulation instability suppresses the competing hosing instability. This work reveals that a proof-of-principle experiment to test the physics of bunch self-modulation can be performed with available lepton bunches and with existing experimental apparatus and diagnostics.

preprint2012arXiv

Weibel-instability-mediated collisionless shocks in laboratory with ultraintense lasers

The formation of non-relativistic collisionless shocks in laboratory with ultrahigh intensity lasers is studied via \emph{ab initio} multi-dimensional particle-in-cell simulations. The microphysics behind shock formation and dissipation, and the detailed shock structure are analyzed, illustrating that the Weibel instability plays a crucial role in the generation of strong subequipartition magnetic fields that isotropize the incoming flow and lead to the formation of a collisionless shock, similarly to what occurs in astrophysical scenarios. The possibility of generating such collisionless shocks in laboratory opens the way to the direct study of the physics associated with astrophysical shocks.

preprint2011arXiv

Effect of the frequency chirp on laser wakefield acceleration

The role of laser frequency chirps in the laser wakefield accelerator is examined. We show that in the linear regime, the evolution of the laser pulse length is affected by the frequency chirp, and that positive (negative) chirp compresses (stretches) the laser pulse, thereby increasing (decreasing) the peak vector potential and wakefield amplitude. In the blowout regime, the frequency chirp can be used to fine tune the localized etching rates at the front of the laser. In our simulations, chirped laser pulses can lead to 15% higher self-trapped electrons, and 10% higher peak energies as compare to the transform-limited pulse. Chirps may be used to control the phase velocity of the wake, and to relax the self-guiding conditions at the front of the laser. Our predictions are confirmed by multi-dimensional particle-in-cell simulations with OSIRIS.

preprint2011arXiv

Large-scale magnetic field generation via the kinetic Kelvin-Helmholtz instability in unmagnetized scenarios

Collisionless plasma instabilities are fundamental in magnetic field generation in astrophysical scenarios, but their role has been addressed in scenarios where velocity shear is absent. In this work we show that velocity shears must be considered when studying realistic astrophysical scenarios, since these trigger the collisionless Kelvin-Helmholtz instability (KHI). We present the first self-consistent three-dimensional (3D) particle-in-cell (PIC) simulations of the KHI in conditions relevant for unmagnetized relativistic outflows with velocity shear, such as active galactic nuclei (AGN) and gamma-ray bursts (GRBs). We show the generation of a strong large-scale DC magnetic field, which extends over the entire shear-surface, reaching thicknesses of a few tens of electron skin depths, and persisting on time-scales much longer than the electron time scale. This DC magnetic field is not captured by MHD models since it arises from intrinsically kinetic effects. Our results indicate that the KHI can generate intense magnetic fields yielding equipartition values up to ε_B/ε_p ~ 10^-3 in the electron time-scale. The KHI-induced magnetic fields have a characteristic structure that will lead to a distinct radiation signature, and can seed the turbulent dynamo amplification process. The dynamics of the KHI are relevant for non-thermal radiation modeling and can also have a strong impact on the formation of relativistic shocks in presence of velocity shears.

preprint2011arXiv

Magnetic control of particle-injection in plasma based accelerators

The use of an external transverse magnetic field to trigger and to control electron self-injection in laser- and particle-beam driven wakefield accelerators is examined analytically and through full-scale particle-in-cell simulations. A magnetic field can relax the injection threshold and can be used to control main output beam features such as charge, energy, and transverse dynamics in the ion channel associated with the plasma blowout. It is shown that this mechanism could be studied using state-of-the-art magnetic fields in next generation plasma accelerator experiments.

preprint2011arXiv

Non-invasive characterization of transverse beam emittance of electrons from a laser-plasma wakefield accelerator in the bubble regime using betatron x-ray radiation

We propose and use a technique to measure the transverse emittance of a laser-wakefield accelerated beam of relativistic electrons. The technique is based on the simultaneous measurements of the electron beam divergence given by $v_{\perp}/v_{\parallel}$, the measured longitudinal spectrum $γ_\parallel$ and the transverse electron bunch size in the bubble $r_{\perp}$. The latter is obtained via the measurement of the source size of the x-rays emitted by the accelerating electron bunch in the bubble. We measure a \textit{normalised} RMS beam transverse emittance $<0.5$ $π$ mm$\:$mrad as an upper limit for a spatially gaussian, spectrally quasi-monoenergetic electron beam with 230 MeV energy in agreement with numerical modeling and analytic theory in the bubble regime.

preprint2011arXiv

Polarized Beam Conditioning in Plasma Based Acceleration

The acceleration of polarized electron beams in the blowout regime of plasma-based acceleration is explored. An analytical model for the spin precession of single beam electrons, and depolarization rates of zero emittance electron beams, is derived. The role of finite emittance is examined numerically by solving the equations for the spin precession with a spin tracking algorithm. The analytical model is in very good agreement with the results from 3D particle-in-cell simulations in the limits of validity of our theory. Our work shows that the beam depolarization is lower for high-energy accelerator stages, and that under the appropriate conditions, the depolarization associated with the acceleration of 100-500 GeV electrons can be kept below 0.1-0.2%.

preprint2011arXiv

Production of picosecond, kilojoule, petawatt laser pulses via Raman amplification of nanosecond pulses

Raman amplification in plasma has been promoted as a means of compressing picosecond optical laser pulses to femtosecond duration to explore the intensity frontier. Here we show for the first time that it can be used, with equal success, to compress laser pulses from nanosecond to picosecond duration. Simulations show up to 60% energy transfer from pump to probe pulses, implying that multi-kiloJoule ultra-violet petawatt laser pulses can be produced using this scheme. This has important consequences for the demonstration of fast-ignition inertial confinement fusion.

preprint2010arXiv

High-brilliance synchrotron radiation induced by the plasma magnetostatic mode

Using multi-dimensional PIC simulations we show that the magnetic undulator-type field of the plasma magnetostatic mode is indeed produced by the interaction of a laser pulse with a relativistic ionization front, as predicted by linear theory for a cold plasma. When the front with this magnetostatic mode is followed by a relativistic electron beam, the interaction of the beam with this magnetic field, produces FEL-type synchrotron radiation, providing a direct signature of the magnetostatic mode. The possibility of generating readily detectable ultrashort wavelength radiation using this mode, by employing state-of-the-art laser systems, is demonstrated, thus opening the way towards experimental observation of the hitherto unseen magnetostatic mode and the use of this plasma FEL mechanism to provide a source of high-brilliance ultrashort wavelength radiation.

preprint2010arXiv

The nonlinear saturation of the non-resonant kinetically driven streaming instability

A non-resonant instability for the amplification of the interstellar magnetic field in young Supernova Remnant (SNR) shocks was predicted by Bell (2004), and is thought to be relevant for the acceleration of cosmic ray (CR) particles. For this instability, the CRs streaming ahead of SNR shock fronts drive electromagnetic waves with wavelengths much shorter than the typical CR Larmor radius, by inducing a current parallel to the background magnetic field. We explore the nonlinear regime of the non-resonant mode using Particle-in-Cell (PIC) hybrid simulations, with kinetic ions and fluid electrons, and analyze the saturation mechanism for realistic CR and background plasma parameters. In the linear regime, the observed growth rates and wavelengths match the theoretical predictions; the nonlinear stage of the instability shows a strong reaction of both the background plasma and the CR particles, with the saturation level of the magnetic field varying with the CR parameters. The simulations with CR-to-background density ratios of n_CR/n_b=10^(-5) reveal the highest magnetic field saturation levels, with energy also being transferred to the background plasma and to the perpendicular velocity components of the CR particles. The results show that amplification factors >10 for the magnetic field can be achieved, and suggest that this instability is important for the generation of magnetic field turbulence, and for the acceleration of CR particles.

preprint2009arXiv

Ion dynamics and acceleration in relativistic shocks

Ab-initio numerical study of collisionless shocks in electron-ion unmagnetized plasmas is performed with fully relativistic particle in cell simulations. The main properties of the shock are shown, focusing on the implications for particle acceleration. Results from previous works with a distinct numerical framework are recovered, including the shock structure and the overall acceleration features. Particle tracking is then used to analyze in detail the particle dynamics and the acceleration process. We observe an energy growth in time that can be reproduced by a Fermi-like mechanism with a reduced number of scatterings, in which the time between collisions increases as the particle gains energy, and the average acceleration efficiency is not ideal. The in depth analysis of the underlying physics is relevant to understand the generation of high energy cosmic rays, the impact on the astrophysical shock dynamics, and the consequent emission of radiation.

preprint2008arXiv

Beam loading in the nonlinear regime of plasma-based acceleration

A theory that describes how to load negative charge into a nonlinear, three-dimensional plasma wakefield is presented. In this regime, a laser or an electron beam blows out the plasma electrons and creates a nearly spherical ion channel, which is modified by the presence of the beam load. Analytical solutions for the fields and the shape of the ion channel are derived. It is shown that very high beam-loading efficiency can be achieved, while the energy spread of the bunch is conserved. The theoretical results are verified with the Particle-In-Cell code OSIRIS.

preprint2006arXiv

dHybrid: a massively parallel code for hybrid simulations of space plasmas

A massively parallel simulation code, called \textit{dHybrid}, has been developed to perform global scale studies of space plasma interactions. This code is based on an explicit hybrid model; the numerical stability and parallel scalability of the code are studied. A stabilization method for the explicit algorithm, for regions of near zero density, is proposed. Three-dimensional hybrid simulations of the interaction of the solar wind with unmagnetized artificial objects are presented, with a focus on the expansion of a plasma cloud into the solar wind, which creates a diamagnetic cavity and drives the Interplanetary Magnetic Field out of the expansion region. The dynamics of this system can provide insights into other similar scenarios, such as the interaction of the solar wind with unmagnetized planets.

preprint2006arXiv

Physical problems (microphysics) in relativistic plasma flows

Many problems in astrophysics involve relativistic outflows. The plasma dynamics in these scenarios is critical to determine the conditions for the self-consistent evolution of the fields and particle acceleration. Advances in computer power now allow for kinetic plasma simulations, based on the particle-in-cell (PIC) paradigm, capable of providing information about the role of plasma instabilities in relativistic outflows. A discussion of the key issues associated with PIC simulations is presented, along with some the most important results and open questions, with a particular emphasis on the long time evolution of the filamentation, or Weibel, instability, and on the possible collisionless mechanisms for particle acceleration arising in the collision of relativistic plasma shells.