Laser-driven Rayleigh-Taylor instability: plasmonic effects and three-dimensional structures

Improved energy spread in the radiation pressure acceleration of protons with a linearly polarized laser

Degradation in the energy spread of accelerated protons due to the transverse instability induced transparency is one of the critical issues in the laser-driven radiation pressure acceleration (RPA) scheme. This issue is more severe for linearly polarized lasers due to enhanced heating of electrons. Therefore, in spite of being experimentally challenging, most of the numerical studies are performed with circularly polarized lasers. In this work, through particle-in-cell simulations, we demonstrate a significant improvement in the energy spread of the accelerated protons when a multilayered target is irradiated by a linearly polarized laser. This multilayered target consists of a near-critical-density (NCD) layer, sandwiched between a thick metallic foil and a thin RPA target. The role of the NCD target is to suppress the laser transparency to increase the coupling of laser momentum to the RPA protons. On the other hand, the metallic foil utilizes the kinetic energy of the escaping fast electrons to form an electrostatic sheath to filter the low-energy RPA protons. This results in significant improvement in the accelerated proton spectrum, even with a linearly polarized laser.

Stabilized Radiation Pressure Acceleration and Neutron Generation in Ultrathin Deuterated Foils

Premature relativistic transparency of ultrathin, laser-irradiated targets is recognized as an obstacle to achieving a stable radiation pressure acceleration in the "light sail" (LS) mode. Experimental data, corroborated by 2D PIC simulations, show that a few-nm thick overcoat surface layer of high Z material significantly improves ion bunching at high energies during the acceleration. This is diagnosed by simultaneous ion and neutron spectroscopy following irradiation of deuterated plastic targets. In particular, copious and directional neutron production (significantly larger than for other in-target schemes) arises, under optimal parameters, as a signature of plasma layer integrity during the acceleration.

Rayleigh-Taylor instability in strongly coupled plasma

Rayleigh–Taylor instability (RTI) is the prominent energy mixing mechanism when heavy fluid lies on top of light fluid under the gravity. In this work, the RTI is studied in strongly coupled plasmas using two-dimensional molecular dynamics simulations. The motivation is to understand the evolution of the instability with the increasing correlation (Coulomb coupling) that happens when the average Coulombic potential energy becomes comparable to the average thermal energy. We report the suppression of the RTI due to a decrease in growth rate with increasing coupling strength. The caging effect is expected a physical mechanism for the growth suppression observed in both the exponential and the quadratic growth regimes. We also report that the increase in shielding due to background charges increases the growth rate of the instability. Moreover, the increase in the Atwood number, an entity to quantify the density gradient, shows the enhancement of the growth of the instability. The dispersion relation obtained from the molecular dynamics simulation of strongly coupled plasma shows a slight growth enhancement compared to the hydrodynamic viscous fluid. The RTI and its eventual impact on turbulent mixing can be significant in energy dumping mechanisms in inertial confinement fusion where, during the compressed phases, the coupling strength approaches unity.

Dimensional effects in Rayleigh-Taylor mixing

We study the effects of dimensional confinement on the evolution of incompressible Rayleigh-Taylor mixing both in a bulk flow and in porous media by means of numerical simulations of the transport equations. In both cases, the confinement to two-dimensional flow accelerates the mixing process and increases the speed of the mixing layer. Dimensional confinement also produces stronger correlations between the density and the velocity fields affecting the efficiency of the mass transfer, quantified by the dependence of the Nusselt number on the Rayleigh number. This article is part of the theme issue 'Scaling the turbulence edifice (part 2)'.

Incompressible Rayleigh-Taylor mixing in circular and spherical geometries

We present a numerical investigation of the turbulent evolution of the mixing layer developing from the Rayleigh-Taylor instability for miscible incompressible fluids in circular (in two dimensions) and in spherical (in three dimensions) geometries in the Boussinesq approximation. We show that the main difference caused by the convergent geometry with respect to the planar case is that the center of the mixing layer drifts toward the center of the domain during the evolution of the mixing layer. A similar effect is observed for the radial profile of the density flux. We derive a simple geometrical relation for this inward drift based on mass conservation. In the late stage of the evolution we observe also the appearance of an inward-outward asymmetry in the radial profiles.

PIC methods in astrophysics: simulations of relativistic jets and kinetic physics in astrophysical systems

The Particle-In-Cell (PIC) method has been developed by Oscar Buneman, Charles Birdsall, Roger W. Hockney, and John Dawson in the 1950s and, with the advances of computing power, has been further developed for several fields such as astrophysical, magnetospheric as well as solar plasmas and recently also for atmospheric and laser-plasma physics. Currently more than 15 semi-public PIC codes are available which we discuss in this review. Its applications have grown extensively with increasing computing power available on high performance computing facilities around the world. These systems allow the study of various topics of astrophysical plasmas, such as magnetic reconnection, pulsars and black hole magnetosphere, non-relativistic and relativistic shocks, relativistic jets, and laser-plasma physics. We review a plethora of astrophysical phenomena such as relativistic jets, instabilities, magnetic reconnection, pulsars, as well as PIC simulations of laser-plasma physics (until 2021) emphasizing the physics involved in the simulations. Finally, we give an outlook of the future simulations of jets associated to neutron stars, black holes and their merging and discuss the future of PIC simulations in the light of petascale and exascale computing.

Scaling laws for laser-driven ion acceleration from nanometer-scale ultrathin foils

Laser-driven ion acceleration has attracted global interest for its potential towards the development of a new generation of compact, low-cost accelerators. Remarkable advances have been seen in recent years with a substantial proton energy increase in experiments, when nanometer-scale ultrathin foil targets and high-contrast intense lasers are applied. However, the exact acceleration dynamics and particularly the ion energy scaling laws in this novel regime are complex and still unclear. Here, we derive a scaling law for the attainable maximum ion energy from such laser-irradiated nanometer-scale foils based on analytical theory and multidimensional particle-in-cell simulations, and further show that this scaling law can be used to accurately describe experimental data over a large range of laser and target parameters on different facilities. This provides crucial references for parameter design and experimentation of the future laser devices towards various potential applications.

Effects of the Transverse Instability and Wave Breaking on the Laser-Driven Thin Foil Acceleration

Acceleration of ultrathin foils by the laser radiation pressure promises a compact alternative to the conventional ion sources. Among the challenges on the way to practical realization, one fundamental is a strong transverse plasma instability, which develops density perturbations and breaks the acceleration. In this Letter, we develop a theoretical model supported by three-dimensional numerical simulations to explain the transverse instability growth from noise to wave breaking and its crucial effect on stopping the acceleration. The wave-broken nonlinear mode triggers rapid stochastic heating that finally explodes the target. Possible paths to mitigate this problem for getting efficient ion acceleration are discussed.

Proton sheet crossing in thin relativistic plasma irradiated by a femtosecond petawatt laser pulse

Leveraging on analyses of Hamiltonian dynamics to examine the ion motion, we explicitly demonstrate that the proton sheet crossing and plateau-type energy spectrum are two intrinsic features of the effectively accelerated proton beams driven by a drift quasistatic longitudinal electric field. Via two-dimensional particle-in-cell simulations, we show the emergence of proton sheet crossing in a relativistically transparent plasma foil irradiated by a linearly polarized short pulse with the power of one petawatt. Instead of successively blowing the whole foil forward, the incident laser pulse readily penetrates through the plasma bulk, where the proton sheet crossing takes place and the merged self-generated longitudinal electric field traps and reflects the protons to yield a group of protons with plateau-type energy spectrum.

Physical mechanism of the electron-ion coupled transverse instability in laser pressure ion acceleration for different regimes

In radiation pressure ion acceleration (RPA) research, the transverse stability within laser plasma interaction has been a long-standing, crucial problem over the past decades. In this paper, we present a one-dimensional two-fluid theory extended from a recent work Wan et al. Phys. Rev. Lett. 117, 234801 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.234801 to clearly clarify the origin of the intrinsic transverse instability in the RPA process. It is demonstrated that the purely growing density fluctuations are more likely induced due to the strong coupling between the fast oscillating electrons and quasistatic ions via the ponderomotive force with spatial variations. The theory contains a full analysis of both electrostatic (ES) and electromagnetic modes and confirms that the ES mode actually dominates the whole RPA process at the early linear stage. By using this theory one can predict the mode structure and growth rate of the transverse instability in terms of a wide range of laser plasma parameters. Two-dimensional particle-in-cell simulations are systematically carried out to verify the theory and formulas in different regimes, and good agreements have been obtained, indicating that the electron-ion coupled instability is the major factor that contributes the transverse breakup of the target in RPA process.

Radiation-pressure-driven ion Weibel instability and collisionless shocks

The Weibel instability from counterstreaming plasma flows is a basic process highly relevant for collisionless shock formation in astrophysics. In this paper we investigate, via two- and three-dimensional simulations, suitable configurations for laboratory investigations of the ion Weibel instability (IWI) driven by a fast quasineutral plasma flow launched into the target via the radiation pressure of an ultra-high-intensity laser pulse ("hole-boring" process). The use of S-polarized light at oblique incidence is found to be an optimal configuration for driving IWI, as it prevents the development of surface rippling observed at normal incidence that would lead to strong electron heating and would favor competing instabilities. Conditions for the evolution of IWI into a collisionless shock are also investigated.

Relativistic Electron Streaming Instabilities Modulate Proton Beams Accelerated in Laser-Plasma Interactions

We report experimental evidence that multi-MeV protons accelerated in relativistic laser-plasma interactions are modulated by strong filamentary electromagnetic fields. Modulations are observed when a preplasma is developed on the rear side of a μm-scale solid-density hydrogen target. Under such conditions, electromagnetic fields are amplified by the relativistic electron Weibel instability and are maximized at the critical density region of the target. The analysis of the spatial profile of the protons indicates the generation of B>10  MG and E>0.1  MV/μm fields with a μm-scale wavelength. These results are in good agreement with three-dimensional particle-in-cell simulations and analytical estimates, which further confirm that this process is dominant for different target materials provided that a preplasma is formed on the rear side with scale length ≳0.13λ_{0}sqrt[a_{0}]. These findings impose important constraints on the preplasma levels required for high-quality proton acceleration for multipurpose applications.

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 quasistatic ions, a mechanism similar to the oscillating two stream instability in the inertial confinement fusion 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 breakup of the target.

Surface modulation and back reflection from foil targets irradiated by a Petawatt femtosecond laser pulse at oblique incidence

A significant level of back reflected laser energy was measured during the interaction of ultra-short, high contrast PW laser pulses with solid targets at 30° incidence. 2D PIC simulations carried out for the experimental conditions show that at the laser-target interface a dynamic regular structure is generated during the interaction, which acts as a grating (quasi-grating) and reflects back a significant amount of incident laser energy. With increasing laser intensity above 10¹⁸ W/cm² the back reflected fraction increases due to the growth of the surface modulation to larger amplitudes. Above 10²⁰ W/cm² this increase results in the partial destruction of the quasi-grating structure and, hence, in the saturation of the back reflection efficiency. The PIC simulation results are in good agreement with the experimental findings, and, additionally, demonstrate that in presence of a small amount of pre-plasma this regular structure will be smeared out and the back reflection reduced.

Towards optical polarization control of laser-driven proton acceleration in foils undergoing relativistic transparency

Control of the collective response of plasma particles to intense laser light is intrinsic to relativistic optics, the development of compact laser-driven particle and radiation sources, as well as investigations of some laboratory astrophysics phenomena. We recently demonstrated that a relativistic plasma aperture produced in an ultra-thin foil at the focus of intense laser radiation can induce diffraction, enabling polarization-based control of the collective motion of plasma electrons. Here we show that under these conditions the electron dynamics are mapped into the beam of protons accelerated via strong charge-separation-induced electrostatic fields. It is demonstrated experimentally and numerically via 3D particle-in-cell simulations that the degree of ellipticity of the laser polarization strongly influences the spatial-intensity distribution of the beam of multi-MeV protons. The influence on both sheath-accelerated and radiation pressure-accelerated protons is investigated. This approach opens up a potential new route to control laser-driven ion sources.

Intense laser pulse interaction with ultra-thin foils constitutes a promising approach for proton acceleration. Here the authors show that the degree of ellipticity in the laser beam polarization can be used to control the proton beam profile.

Chirped-Standing-Wave Acceleration of Ions with Intense Lasers

We propose a novel mechanism for ion acceleration based on the guided motion of electrons from a thin layer. The electron motion is locked to the moving nodes of a standing wave formed by a chirped laser pulse reflected from a mirror behind the layer. This provides a stable longitudinal field of charge separation, thus giving rise to chirped-standing-wave acceleration of the residual ions of the layer. We demonstrate, both analytically and numerically, that stable proton beams, with energy spectra peaked around 100 MeV, are feasible for pulse energies at the level of 10 J. Moreover, a scaling law for higher laser intensities and layer densities is presented, indicating stable GeV-level energy gains of dense ion bunches, for soon-to-be-available laser intensities.

Electron heating in radiation-pressure-driven proton acceleration with a circularly polarized laser

Dynamics of electron heating in the radiation-pressure-driven acceleration through self-induced transparency (SIT) is investigated with the help of particle-in-cell simulations. The SIT is achieved through laser filamentation which is seeded by the transverse density modulations due to the Rayleigh-Taylor-like instability. We observe stronger SIT induced electron heating for the longer duration laser pulses leading to deterioration of accelerated ion beam quality (mainly energy spread). Such heating can be controlled to obtain a quasimonoenergetic beam by cascaded foils targets where a second foil behind the main accelerating foil acts as a laser reflector to suppress the SIT.

Numerical investigation of the transverse instability on the radiation-pressure-driven foil

The development of transverse instability in the radiation-pressure-acceleration dominant laser-foil interaction is numerically examined by two-dimensional particle-in-cell simulations. When a plane laser impinges on a foil with modulated surface, the transverse instability is incited, and periodic perturbations of the proton density develop. The growth rate of the transverse instability is numerically diagnosed. It is found that the linear growth of the transverse instability lasts only a few laser periods, then the instability gets saturated. In order to optimize the modulation wavelength of the target, a method of information entropy is put forward to describe the chaos degree of the transverse instability. With appropriate modulation, the transverse instability shows a low chaos degree, and a quasi-monoenergetic proton beam is produced.