Source author record

R. A. Fonseca

R. A. Fonseca 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

45works
11topics
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

45 published item(s)

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.

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

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

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

Formation of Lunar Swirls

In this paper we show a plausible mechanism that could lead to the formation of the Dark Lanes in Lunar Swirls, and the electromagnetic shielding of the lunar surface that results in the preservation of the white colour of the lunar regolith. We present the results of a fully self-consistent 2 and 3 dimensional particle-in-cell simulations of mini-magnetospheres that form above the lunar surface and show that they are consistent with the formation of `lunar swirls' such as the archetypal formation Reiner Gamma. The simulations show how the microphysics of the deflection/shielding of plasma operates from a kinetic-scale cavity, and show that this interaction leads to a footprint with sharp features that could be the mechanism behind the generation of `dark lanes'. The physics of mini-magnetospheres is described and shown to be controlled by space-charge fields arising due to the magnetized electrons and unmagnetized ions. A comparison between model and observation is shown for a number of key plasma parameters.

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

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

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

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

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

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.

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

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

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

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

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

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.