Hollow screw-like drill in plasma using an intense Laguerre–Gaussian laser

With the development of ultra-intense laser technology, MeV ions can be obtained from laser–foil interactions in the laboratory. These energetic ion beams can be applied in fast ignition for inertial confinement fusion, medical therapy, and proton imaging. However, these ions are mainly accelerated in the laser propagation direction. Ion acceleration in an azimuthal orientation was scarcely studied. In this research, a doughnut Laguerre–Gaussian (LG) laser is used for the first time to examine laser–plasma interaction in the relativistic intensity regime in three-dimensional particle-in-cell simulations. Studies have shown that a novel rotation of the plasma is produced from the hollow screw-like drill of an mode laser. The angular momentum of particles in the longitudinal direction produced by the LG laser is enhanced compared with that produced by the usual laser pulses, such as linearly and circularly polarized Gaussian pulses. Moreover, the particles (including electrons and ions) can be trapped and uniformly compressed in the dark central minimum of the doughnut LG pulse. The hollow-structured LG laser has potential applications in the generation of x-rays with orbital angular momentum, plasma accelerators, fast ignition for inertial confinement fusion, and pulsars in the astrophysical environment.

With the development of ultra-intense laser technology, MeV ions can be obtained from laser-foil interactions in the laboratory. These energetic ion beams can be applied in fast ignition for inertial confinement fusion, medical therapy, and proton imaging. However, these ions are mainly accelerated in the laser propagation direction. Ion acceleration in an azimuthal orientation was scarcely studied. In this research, a doughnut Laguerre-Gaussian (LG) laser is used for the first time to examine laser-plasma interaction in the relativistic intensity regime in three-dimensional particle-in-cell simulations. Studies have shown that a novel rotation of the plasma is produced from the hollow screw-like drill of an LG 1 0 mode laser. The angular momentum of particles in the longitudinal direction produced by the LG laser is enhanced compared with that produced by the usual laser pulses, such as linearly and circularly polarized Gaussian pulses. Moreover, the particles (including electrons and ions) can be trapped and uniformly compressed in the dark central minimum of the doughnut LG pulse. The hollow-structured LG laser has potential applications in the generation of x-rays with orbital angular momentum, plasma accelerators, fast ignition for inertial confinement fusion, and pulsars in the astrophysical environment. N owadays, laser intensity can increase up to 10 22 W/cm 2 1,2 . Energetic protons have been obtained through different mechanisms, such as target normal sheath acceleration [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17] , radiation pressure acceleration [18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37] , collisionless shock acceleration [38][39][40][41][42] , breakout afterburner 43,44 , and a combination of different mechanisms [43][44][45][46] . However, these ions are mainly accelerated in the laser propagation direction. Ion acceleration in azimuthal orientation is scarcely mentioned. A circularly polarized (CP) light may carry the angular momentum 47 . The main reason for this phenomenon is that a CP light carries an orbital angular momentum s z h (s z 561) per photon. More than 70 years ago, the mechanical torque created by the transfer of angular momentum of a CP light was first observed in Beth's 47 experiments. However, the small quantities of the optical angular momentum are difficult to detect in the CP light experiments.
A laser with a Gaussian mode, such as Laguerre-Gaussian (LG) mode, also possesses an orbital angular momentum 48 . A linearly polarized (LP) LG laser with a helical wave-front structure has a central phase singularity 49 . The angular momentum produced by such structure is sometimes referred to an orbital angular momentum, which is different from the spin angular momentum produced by the CP laser pulse 47 . LG laser pulse is circularly symmetric in the cross-section with respect to the optical axis [the direction of light propagation, Figs. 1(b) to 1(e)]. The mode of the LG laser pulse (LG l p ) is described by integer indices l and p, where l denotes the number of 2p phase cycles around the circumference and (p 1 1) denotes the number of radial nodes in the mode profile. This study discusses the LG 1 0 mode, where l ? 0 indicates the presence of an azimuthal phase term exp(-il w) in the laser mode. LG 1 0 laser carries an orbital angular momentum l h per photon. The nature of the orbital angular momentum of different LG modes has been investigated in optical trapping experiments. Allen et al. 48 showed that an LG mode has a well-defined orbital angular momentum. They also observed the torque on suspended cylindrical lenses arising from the reversal helicity of an LG mode. He et al. 50 demonstrated that the absorptive particles trapped in the dark central minimum of a doughnut laser pulse are set into rotation. Furthermore, the rotation particles are controlled using both the spin and the orbital angular momentum of light. In such case, the LG light is beneficial because it reduces the ac stark shift and the broadening of transitions at the trap center. Kuga et al. 51 proposed to trap atoms along the beam center using an LG light. They exploited the spatial profile of LG modes with p 5 0, which has the form of a ring of light. This feature is important in laser cooling and trapping experiments because the repulsive optical dipole force for blue detuned laser light restricts the atoms to the inner dark region of the laser beam, where photon scattering and the associated heating are minimized. Such hollow-structured LG laser can be used to investigate some difficult problems, such as generation of x-rays with orbital angular momentum 52-54 , plasma accelerators 55,56 , fast ignition for inertial confinement fusion [57][58][59] , and pulsars in the astrophysical environment 60 .

Results
In this letter, the doughnut LG laser is used for the first time in relativistic intensity laser plasma interaction. The LG laser rotates electrons and protons in the azimuthal orientation. Unlike conventional laser pulses, such as LP and CP, enhancement of the proton angular momentum along the longitudinal direction is obtained when an intense LP LG laser pulse irradiates on a thin foil. Threedimensional (3D) PIC simulations are performed to investigate the LG laser interaction on a foil. In the simulation, the field amplitude E (LG l p ) of an LG laser with mode (l, p) is given by where E 0 is the peak amplitude of the electric field, r is the radius, W H is the radius at which the Gaussian term falls to 1/e of its on-axis value, W x is the pulse length in the x direction, 0 , t , 2 W x /c, L l P is the generalized Laguerre polynomial, k is the wave number, v is the laser frequency, x is the distance from the beam waist, x r is the Rayleigh range, x/x r is the Guoy phase of the mode, and w is the azimuthal angle 61 . This study mainly discusses the mode of LG 1 0 , and thus, p 5 0 and l 5 1 are used in Eq. (1). Fig. 2 shows the total angular momentum of the particles (electrons and protons) in 3D PIC simulations. The detailed simulation parameters are shown in the Methods section. To describe the rotation effects of LG 1 0 mode on the plasma, the angular momentum of the particles in the x direction (the longitudinal direction) m e (yp z 2 zp y ) 1 m p (yp z 2 zp y ) is calculated, where p y, z 5 cv y, z is the velocity and c 5 (1 2 v 2 y,z ) 21/2 is the relativistic factor. Fig. 2 further shows the simulation results at different times. The LG laser irradiates on the foil at t 5 30 T, and is totally reflected by the foil at t , 54 T. The angular momentum of the electrons and protons increases up to 21.55 3 10 217 kg?m 2 /s until the laser pulse is totally reflected. Thus, the trapped particles rotate by the helicity of the LG 1 0 laser. Fig. 3 shows the angular momentums of particles at 60 T for different laser amplitudes a 0 . The radius W H and duration t L of the laser pulse remain constant for different a 0 . The angular momentums increase to 21.32 3 10 216 kg?m 2 /s for a 0 5 14. From Fig. 3, a critical condition for the rotation of the particles is observed at a 0 < 5. Detailed simulations have proved this condition. The LG laser is almost reflected for a 0 5 1 at t 5 50 T (Fig. 4). The laser transmits through the foil at a 0 5 5, and more portions of the laser pulse transmit at a 0 5 10. Thus, the critical condition for the proton rotation may be related to the transmission of the foil 18,35,36,[62][63][64][65] , and may be expressed by a 0 < pnd [35][36][37] , where the foil density n is normalized by n c and the foil thickness d is normalized by laser wavelength l. The theory estimates that the laser starts to transmit through the foil at a 0 , 6.28 for n 5 2 and d 5 1 according to a 0 , pnd, which is larger than the simulation results in Figs. 4(d) to 4(f), where the laser begins to transmit through the foil at a 0 5 5. This phenomenon is attributed to the enhancement of the transmission of the laser pulse caused by the self-focusing of the laser pulse and multi-dimensional instabilities, as demonstrated in our previous simulation study 35,36 . Fig. 3 shows that the angular momentums of particles increase when the LG laser transmits through the foil. The main reason is that more laser energy is absorbed by the foil when the laser transmits through the foil. In realistic interaction, the laser transmission may be realized by LG laser pulses with larger intensity and a larger angular momentum of the foil can be generated, similar to the case in Fig. 3.  Fig. 6, which clearly show the process of trapping and compression with time. Fig. 6(a) shows that the electrons are first dragged along the tangential direction at t 5 35 T. The protons remains almost at rest due to   , which may realize the screw-like drilling in the inertial confinement fusion and laser-driven particle accelerations.

Discussion
The particles are rotated by the LG laser (Fig. 2). The transfer of angular momentum from the laser to the particles is then theoretically estimated in realistic cases. In terms of quantum mechanics, the rotation can be caused by the angular momentum of photons. The LG mode can be seen as the eigenmode of the angular momentum operator of L z 48 and carries an orbital angular momentum of l h per photon. The angular momentum carried by a photon of a polarized LG mode laser is (l 1 s z ) h, where s z is 61 for the CP laser and 0 for the LP laser. The total angular momentum absorbed from the laser can then be approximately expressed as where g is the absorbing ratio from the laser pulse during the interaction, S , p W 2 \ and W H are the radius at which the Gaussian term falls to 1/e of its on-axis value, h 5 6.63 3 10 234 J?s is Plank constant, n 5 1/T is the frequency of light, hn is the energy of single photon, and t L is the pulse duration. P angular 5 8.22 3 10 216 kg?m 2 /s is then obtained for an LP LG 1 0 laser (l 5 1 and s z 5 0) with a 0 5 10, W H 5 4 mm, and t L 5 40 fs. A ring of the target (the radius of the inner and outer ring is 2 and 4 mm, respectively) is assumed to be rotated by the hollow structure of the LG 1 0 mode laser. Assuming that the angular momentum of the laser is totally transferred to the foil (g 5 where the foil density is 2 3 10 3 kg/m 3 , foil thickness is d 5 1 mm, c is the relativistic factor, v is the angular velocity of the foil, and r is the radius of the foil. Afterward, v < 1.09 3 10 12 rad/s is obtained and the velocity of the particle at r 5 4 mm is approximately 4.36 3 10 6 m/s. The angular momentum is proportional to the laser amplitude a 2 0 and laser duration t L , indicating that the angular momentum absorbed from the LG laser can be increased to a certain extent using a high-intensity long pulse based on Eq. (2). Notably, the increase of the angular momentums of particles mainly depends on the laser energy and the absorbing ratio from the laser pulse during the interaction g, as shown in Eq. (2). Clearly, a more accurate measurement is taken into account, which states that the total angular momentum is divided by the laser energy. P angular /P laser / g is obtained when the laser energy is considered as P laser 5 a 2 0 S t L . Previous studies have shown that hole boring is deeper for a higher a 0 , and the time of hole boring is longer for a larger t L 36,37 . The particles can absorb more energy from the laser pulse (corresponding to a larger g) with the enhancement of the hole boring. Then, a larger P angular is obtained according to Eq. (2). It should be noted that the total angular momentum transferred from laser to particles can be calculated from Eq. (2) when the laser pulse does not transmit the foil. In this case, the angular momentum P angular is proportional to the laser energy absorbed by the foil (g a 2 0 S t L ). Less angular momentum may be generated when the foil is totally destroyed. In this case, a smaller g is obtained and less energy of the LG laser is contributed to total angular momentum of protons.
To show the difference of LG laser on the rotation of the particles, LP and CP laser pulses are also considered. The total angular momentum density per unit power has been defined by Allen et al., where the cases of LP (s z 5 0) and CP (s z 5 61) are    Fig. 7 shows that the angular momentums of the particles in the longitudinal direction m e (yp z 2 zp y ) 1 m p (yp z 2 zp y ) for LG, LP, and CP laser pulses. The amplitude of the LP pulse is sin v L t ð Þê y (a 0 5 7.4), and the amplitude of the CP pulse is a~a 0 exp 2), where 0 , t , 2 W x /c. The amplitude of the LG pulse is expressed by Eq. (1), where a 0 5 10, W H 5 4 mm, and W x 5 12 mm; d 5 1 mm, and n 0 5 4 n c are used in three cases. The values of laser energy in three cases are similar. Compared with LP and CP lasers, the LG laser can generate larger angular momentum of particles (Fig. 7), indicating that the LG laser is beneficial for the rotation of the particles. Compared with case of CP, the doughnut-structured LG pulse (see Fig. 1) has an S-shape potential well just as shown in Fig. 1. More electrons and ions can be trapped in the potential well as the LG laser screw-like drill into the plasma. The LG laser also works as a fan to blow the plasmas forward, which may enhance the absorption of the laser pulse g (Fig. 5). Thus the total angular momentums of protons can be raised, just as shown in Fig. 7.
In conclusion, the particles rotation alignment in the tangential direction is realized with the use of an intense LG laser pulse. Compared with LP and CP lasers, the enhancement of the proton angular momentum in the longitudinal direction is obtained when an intense LP LG laser pulse irradiating on a thin foil. The PIC simulations show that the LG pulse can drill into the plasma like a screw. The angular momentum of the particles produced by the LG laser is enhanced compared with that produced by the LP and CP pulses with similar pulse energies. It is also found that electrons and protons are trapped and uniformly compressed in the dark central minimum of the doughnut LG pulse. LG laser has been generated by several techniques. For example, a high-order Hermite-Gaussian (HG) mode can be generated by inserting an intra-cavity cross-wire into a laser cavity. An LG laser can then be obtained using a mode converter on this HG laser 66 . A spiral phase plate 67 and a computergenerated hologram 49 may be used to generate the LG modes from a fundamental Gaussian mode (TEM 00 ). The intense LG laser exhibits potential applications in field of relativistic intensity, such as laserdriven plasma accelerators. Such LG laser can be applied in the TNSA or RPA experiments to generate a particle beam with angular momentum in the future.

Methods
The 3D simulations are performed with VORPAL, which is a relativistic, arbitrary, and dimensional hybrid plasma and beam simulation code. It includes the utilities for data analysis and scripts for data visualization. Particle-in-cell (PIC) algorithm is used in VORPAL to describe the kinetic plasma model. A charge-conserving current deposition algorithm is applied in electromagnetic limit to enable the integration of Maxwell's equations without any additional divergence correction. The instantaneous charge distribution is used to calculate Poisson's equation at every time step in the electrostatic limit.
In our simulation, a 40 fs p-polarized LG 1 0 laser pulse is incident from the left side on the foil. The dimensionless peak amplitude of the incident laser pulse is a 0 5 eE 0 / m e v L c 5 5 (the intensity is I 5 3.4 3 10 19 W/cm 2 ), where v L is the laser frequency, l 5 1 mm is the laser wavelength, c is the light speed in a vacuum, and m e and e are the electron mass and charge at rest. The radius of the laser is W H 5 4 mm and W x 5 12 mm. The laser front reaches the front surface of the foil at t 5 30 T (T 5 l/c). The foil thickness is d 5 1 mm and the front surface of the foil is at x 5 30 mm. The transverse range of the foil is 214 mm , y , 14 mm and 214 mm , z ,14 mm. The foil is assumed to be fully ionized into protons and electrons before the arrival of the main pulse. The foil density is n 0 5 2 n c , where n c~v 2 L m e =4pe 2 is the critical density. A low-density step-like density profile is used to simplify the model and reduce the 3D PIC simulation time. The size of the simulation box is (60 3 60 3 60) mm, and the cell number is 600 3 600 3 600. Each cell is filled with 10 protons and 10 electrons.