Observation of postsoliton expansion following laser propagation through an underdense plasma

Compression and acceleration of ions by ultrashort, ultraintense azimuthally polarized light

An efficient plasma compression scheme using azimuthally polarized light is proposed. Azimuthally polarized light possesses a donutlike intensity pattern, enabling it to compress and accelerate ions toward the optical axis across a wide range of parameters. When the light intensity reaches the relativistic regime of 10^{18}W/cm^{2}, and the plasma density is below the critical density, protons can be compressed and accelerated by the toroidal soliton formed by the light. The expansion process of the soliton can be well described by the snowplow model. Three-dimensional particle-in-cell simulations show that within the soliton regime, despite the ion density exceeding ten times the critical density, the ions' energy is insufficient for efficient neutron production. When the light intensity increases to 10^{22}W/cm^{2}, and the plasma density reaches several tens of times the critical density, deuterium ions can be compressed to thousands of times the critical density and simultaneously accelerated to the MeV level by tightly focused azimuthally polarized light during the hole-boring process. This process is far more dramatic compared to the soliton regime and can produce up to 10^{4} neutrons in a few light cycles. Moreover, in the subsequent beam-target stage, neutron yield is estimated to exceed 10^{8}. Finally, we present a comparison with the results obtained using a radially polarized light to examine the influence of light polarization.

Generation of electromagnetic solitons with angular momentum

The optical vortex has been widely studied owing to its specific characteristics such as the orbital angular momentum, hollow intensity distribution, and topological charge. We report the generation of electromagnetic solitons with angular momentum and the conversion of angular momentum via a circularly polarized (CP) laser and underdense plasma interactions on the basis of three-dimensional particle-in-cell simulations. We find that when a CP laser is incident into the underdense plasma, a longitudinal current will be induced off the laser axis, which is critical for the angular momentum conversion. This novel, to the best of our knowledge, regime will allow potential applications such as optical control and electron manipulation.

Electrostatic shock acceleration of ions in near-critical-density plasma driven by a femtosecond petawatt laser

With the recent advances in ultrahigh intensity lasers, exotic astrophysical phenomena can be investigated in laboratory environments. Collisionless shock in a plasma, prevalent in astrophysical events, is produced when a strong electric or electromagnetic force induces a shock structure in a time scale shorter than the collision time of charged particles. A near-critical-density (NCD) plasma, generated with an intense femtosecond laser, can be utilized to excite a collisionless shock due to its efficient and rapid energy absorption. We present electrostatic shock acceleration (ESA) in experiments performed with a high-density helium gas jet, containing a small fraction of hydrogen, irradiated with a 30 fs, petawatt laser. The onset of ESA exhibited a strong dependence on plasma density, consistent with the result of particle-in-cell simulations on relativistic plasma dynamics. The mass-dependent ESA in the NCD plasma, confirmed by the preferential reflection of only protons with two times the shock velocity, opens a new possibility of selective acceleration of ions by electrostatic shock.

Energy absorption in the laser-QED regime

A theoretical and numerical investigation of non-ponderomotive absorption at laser intensities relevant to quantum electrodynamics is presented. It is predicted that there is a regime change in the dependence of fast electron energy on incident laser energy that coincides with the onset of pair production via the Breit-Wheeler process. This prediction is numerically verified via an extensive campaign of QED-inclusive particle-in-cell simulations. The dramatic nature of the power law shift leads to the conclusion that this process is a candidate for an unambiguous signature that future experiments on multi-petawatt laser facilities have truly entered the QED regime.

Channel optimization of high-intensity laser beams in millimeter-scale plasmas

Channeling experiments were performed at the OMEGA EP facility using relativistic intensity (>10^{18}W/cm^{2}) kilojoule laser pulses through large density scale length (∼390-570 μm) laser-produced plasmas, demonstrating the effects of the pulse's focal location and intensity as well as the plasma's temperature on the resulting channel formation. The results show deeper channeling when focused into hot plasmas and at lower densities, as expected. However, contrary to previous large-scale particle-in-cell studies, the results also indicate deeper penetration by short (10 ps), intense pulses compared to their longer-duration equivalents. This new observation has many implications for future laser-plasma research in the relativistic regime.

Ion-beam-plasma interaction effects on electrostatic solitary wave propagation in ultradense relativistic quantum plasmas

Understanding the transport properties of charged particle beams is important not only from a fundamental point of view but also due to its relevance in a variety of applications. A theoretical model is established in this article, to model the interaction of a tenuous positively charged ion beam with an ultradense quantum electron-ion plasma, by employing a rigorous relativistic quantum-hydrodynamic (fluid plasma) electrostatic model proposed in McKerr et al. [M. McKerr, F. Haas, and I. Kourakis, Phys. Rev. E 90, 033112 (2014)PLEEE81539-375510.1103/PhysRevE.90.033112]. A nonlinear analysis is carried out to elucidate the propagation characteristics and the existence conditions of large amplitude electrostatic solitary waves propagating in the plasma in the presence of the beam. Anticipating stationary profile excitations, a pseudomechanical energy balance formalism is adopted to reduce the fluid evolution equation to an ordinary differential equation. Exact solutions are thus obtained numerically, predicting localized excitations (pulses) for all of the plasma state variables, in response to an electrostatic potential disturbance. An ambipolar electric field form is also obtained. Thorough analysis of the reality conditions for all variables is undertaken in order to determine the range of allowed values for the solitonic pulse speed and how it varies as a function of the beam characteristics (beam velocity and density).

Relativistic breather-type solitary waves with linear polarization in cold plasmas

Linearly polarized solitary waves, arising from the interaction of an intense laser pulse with a plasma, are investigated. Localized structures, in the form of exact numerical nonlinear solutions of the one-dimensional Maxwell-fluid model for a cold plasma with fixed ions, are presented. Unlike stationary circularly polarized solitary waves, the linear polarization gives rise to a breather-type behavior and a periodic exchange of electromagnetic energy and electron kinetic energy at twice the frequency of the wave. A numerical method based on a finite-differences scheme allows us to compute a branch of solutions within the frequency range Ωmin<Ω<ωpe, where ωpe and Ωmin are the electron plasma frequency and the frequency value for which the plasma density vanishes locally, respectively. A detailed description of the spatiotemporal structure of the waves and their main properties as a function of Ω is presented. Small-amplitude oscillations appearing in the tail of the solitary waves, a consequence of the linear polarization and harmonic excitation, are explained with the aid of the Akhiezer-Polovin system. Direct numerical simulations of the Maxwell-fluid model show that these solitary waves propagate without change for a long time.

Improving the self-guiding of an ultraintense laser by tailoring its longitudinal profile

Self-guiding of an ultraintense laser requires the refractive index to build up rapidly to a sufficient value before the main body of the pulse passes by. We show that placing a low-intensity precursor in front of the main pulse mitigates the diffraction of its leading edge and facilitates reaching a self-guided state that remains stable for more than 10 Rayleigh lengths. Furthermore, this precursor slows the phase slippage between the trapped electrons and the wakefield and leads to an accelerating structure that is more stable, contains more energy, and is sustained longer. Examples from three-dimensional particle-in-cell simulations show that the conversion efficiency from the laser to the self-trapped electrons increases by an order of magnitude when using the precursor.

Modeling relativistic soliton interactions in overdense plasmas: a perturbed nonlinear Schrödinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrödinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrödinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude.

Invited article: Relation between electric and magnetic field structures and their proton-beam images

Proton imaging is commonly used to reveal the electric and magnetic fields that are found in high energy density plasmas. Presented here is an analysis of this technique that is directed towards developing additional insight into the underlying physics. This approach considers: formation of images in the limits of weak and strong intensity variations; caustic formation and structure; image inversion to obtain line-integrated field characteristics; direct relations between images and electric or magnetic field structures in a plasma; imaging of sharp features such as Debye sheaths and shocks. Limitations on spatial and temporal resolution are assessed, and similarities with optical shadowgraphy are noted. Synthetic proton images are presented to illustrate the analysis. These results will be useful for quantitatively analyzing experimental proton imaging data and verifying numerical codes.

Anticorrelation between ion acceleration and nonlinear coherent structures from laser-underdense plasma interaction

In laser-plasma experiments, we observed that ion acceleration from the Coulomb explosion of the plasma channel bored by the laser is prevented when multiple plasma instabilities, such as filamentation and hosing, and nonlinear coherent structures (vortices or postsolitons) appear in the wake of an ultrashort laser pulse. The tailoring of the longitudinal plasma density ramp allows us to control the onset of these instabilities. We deduced that the laser pulse is depleted into these structures in our conditions, when a plasma at about 10% of the critical density exhibits a gradient on the order of 250 μm (Gaussian fit), thus hindering the acceleration. A promising experimental setup with a long pulse is demonstrated enabling the excitation of an isolated coherent structure for polarimetric measurements and, in further perspectives, parametric studies of ion plasma acceleration efficiency.

Two-dimensional s-polarized solitary waves in plasmas. II. Stability, collisions, electromagnetic bursts, and post-soliton evolution

The dynamics of two-dimensional s-polarized solitary waves is investigated with the aid of particle-in-cell (PIC) simulations. Instead of the usual excitation of the waves with a laser pulse, the PIC code was directly initialized with the numerical solutions from the fluid plasma model. This technique allows the analysis of different scenarios including the theoretical problems of the solitary wave stability and their collision as well as features already measured during laser-plasma experiments such as the emission of electromagnetic bursts when the waves reach the plasma-vacuum interface, or their expansion on the ion time scale, usually named post-soliton evolution. Waves with a single density depression are stable whereas multihump solutions decay to several waves. Contrary to solitons, two waves always interact through a force that depends on their relative phases, their amplitudes, and the distance between them. On the other hand, the radiation pattern at the plasma-vacuum interface was characterized, and the evolution of the diameter of different waves was computed and compared with the "snow plow" model.

Two-dimensional s-polarized solitary waves in relativistic plasmas. I. The fluid plasma model

The properties of two-dimensional linearly s-polarized solitary waves are investigated by fluid-Maxwell equations and particle-in-cell (PIC) simulations. These self-trapped electromagnetic waves appear during laser-plasma interactions, and they have a dominant electric field component E(z), normal to the plane of the wave, that oscillates at a frequency below the electron plasma frequency ω(pe). A set of equations that describe the waves are derived from the plasma fluid model in the case of cold or warm plasma and then solved numerically. The main features, including the maximum value of the vector potential amplitude, the total energy, the width, and the cavitation radius are presented as a function of the frequency. The amplitude of the vector potential increases monotonically as the frequency of the wave decreases, whereas the width reaches a minimum value at a frequency of the order of 0.82 ω(pe). The results are compared with a set of PIC simulations where the solitary waves are excited by a high-intensity laser pulse.