Abstract
Vortex beams with ultra-high brilliance can greatly enrich the light and matter interaction process and even shed light on the unexpected information in relativistic nonlinear optics. Thus, we propose a scheme for relativistic intense vortex harmonic radiation by use of bi-circular Laguerre–Gaussian lasers irradiating relativistic plasmas. Three-dimensional particle-in-cell simulation results show that the emitted harmonics own controllable spin and orbital angular momentum simultaneously, which can be attributed to the vortex mirror model and the related conservation laws. Meanwhile, the conversion efficiency of harmonic generation can be tuned through adjusting the intensity ratio of the two driving field components.
Introduction
Light waves can carry two distinct forms of angular momenta in the direction of propagation: spin and orbital angular momentum (SAM and OAM). The SAM is manifest as right or left circular polarization of light ( per photon), whereas the OAM is associated with its characteristic helical phase profile which can be described by an angular-dependent phase , where is the azimuthal coordinate and is the topological charge (known as vortex beams). Different values of represent different vortex modes, indicating that the OAM carried by each photon is equal to . In addition, vortex beams do not have a well-defined phase at the center, resulting in the typical zero intensity center (donut-shaped intensity profile) [1, 2]. These remarkable properties have led to a wide range of applications going from optical communications, bio-photonics, optical micromanipulations to ion acceleration [3], etc. In recent years, high-order vortex beams in the extreme-ultraviolet (XUV) spectral range have attracted more and more attention for their novel potential applications. However, the traditional techniques using spiral phase plates or forked diffraction gratings are not effective due to the limitation of etching resolution. Thus, many novel schemes have been attempted: the free-electron lasers (FELs) can be used to generate such vortex beams, but the huge size and high cost greatly limit the scope of their application. As a table-top and alternative to large-scale instruments, high-harmonic generation (HHG) based on the frequency up-conversion has attracted an increasing attention. This nonlinear process includes laser-gas HHG [4, 5] and laser-plasma HHG [6–22]. These two HHG mechanisms are fundamentally different, and the intensity of emitted harmonics in the latter case can be increased by several orders of magnitude compared with that of the former, in which the limitations of the driving laser intensity (often below ) exist to avoid strong ionization. Among the existing researches of HHG, two-color laser scheme has been extensively investigated in both gas and plasma for its many favorable properties of the harmonic radiation, including high efficiency, spectra tunability and so on [4, 15, 17]. Initially, this scheme was to use the circularly polarized (CP) fundamental driving laser with a Gaussian spatial profile mixed with its counter-rotating second harmonic in collinear illumination, in which the spin angular momentum (SAM) carried by the driving lasers can be successfully transferred to the harmonics. Comparing with the non-collinear one, higher HHG efficiency can be acquired via collinear scheme, which has been proved in Ref. [16]. Based on this, if the driving fields are replaced by two-color twisted CP driving pulses, the OAM of the helical pulses can also be transferred to the harmonics. This point has already been proved in laser-gas HHG process, in which high harmonics with both SAM and OAM were generated and a form of simultaneous SAM and OAM conservation laws were uncovered through theoretical analysis and experimental investigation [23]. However, with the urgent requirement for intense vortex beams up to the relativistic regime, whether this scheme should work in relativistic plasma HHG process needs to be validated.
In this paper, we use full three-dimensional particle-in-cell (PIC) simulations to verify the effectiveness of two-color twisted CP laser scheme in the relativistic plasma HHG regime and find that the similar SAM and OAM conservation laws uncovered in laser-gas HHG are applicable to laser-plasma HHG. As a result, intense harmonic radiation with controllable SAM and OAM can be acquired by use of two-color Laguerre–Gaussian (LG) driving pulses irradiating relativistic plasmas. Its underlying mechanism can be attributed to the “vortex oscillating mirror” (VOM) model [2] and the conservation laws of photon’s energy, SAM and OAM. Meanwhile, the vortex mode and conversion efficiency of radiated harmonics can be flexibly tuned. Therefore, this work provides an efficient and practical approach to produce bright, spectral tunable harmonic radiation with designable SAM and OAM, which may own the potential of extending the existing application areas of vortex beams to ultra-high intensity regime and nanometer spatial, sub-femtosecond temporal scales.
Particle-in-cell simulation results
We perform full three-dimensional PIC simulations to investigate the two-color counter-rotating twisted laser scheme with the open-source EPOCH code [24]. The total electric field of the dichromatic twisted pulses can be expressed as:Where , . The fundamental laser has a frequency of corresponding to the wavelength of and the topological charge of , the second one owns frequency of corresponding to the wavelength of and the topological charge of . For the fundamental pulse, the normalized peak amplitude is , whereis the electric field amplitude, are the unit electron charge, the electron mass and the light speed in vacuum, respectively. This corresponds to a laser peak intensity of about. The temporal intensity profile owns the full-width at half-maximum pulse duration of . The laser spot radius is. The intensity, pulse duration and spot radiusof the second driving laser component are kept the same with those of the fundamental one. The two driving pulses are spatiotemporally overlapped and normally irradiate the solid target along the X direction to drive the laser-plasma HHG process. The target is fully ionized with a peak electron density of , whereis the plasma critical density with respect to the fundamental laser frequency andrepresents the vacuum dielectric constant. The target thickness isand has pre-plasma with an exponential density distribution of at the front surface of the plasma target. The pre-plasma scale length is set as. The ions are assumed to be immobile. The grid size of the simulation box is. The grid cell size is and each cell is filled with 8 macro-particles. Absorption boundary conditions are used for both fields and particles in the whole simulation process. The sketch of our scheme is shown in Figure 1 (the inset displays the enlarged view of the vortex modes of the fundamental and second-harmonic driving pulses):
FIGURE 1
First of all, the spectra feature of produced harmonics including their intensity and ellipticity distribution versusand(are the spatial frequencies along X and Y directions normalized by the fundamental laser wave number , respectively) are shown in Figure 2. Here, the ellipticity of emitted harmonicscan be acquired as follows [19]: First, we perform Fourier transform of the two electric filed components of generated harmonics: , , whereis the frequency of harmonics, are the strength of the harmonic field along the direction, respectively. Then, the ellipticity of generated harmonics can be calculated by. In this equation,,. From Figure 2A, one can see only the harmonics with orders of are produced, while every third one is missing. On the other hand, Figure 2B shows that all of the radiated harmonics are quasi-CP (the SAM per photon of harmonics are close to). In addition, the SAM carried by photons of harmonics () are the same with the fundamental driving pulse (), while those carried by photons of harmonics () are consistent with the double-frequency one (). These conclusions in frequency domain can be verified by the two panels of the inset, which display the reconstructed three-dimensional image of the electric field vectors of 4th and 5th harmonics (red) and the waveforms of the two orthogonal electric field components (black), as well as the projection of (blue). From the inset, one can observe that the amplitude ratios between the two electric components () are both close to unity, while the phase shift of the 4th harmonic and that of the 5th areandrespectively. These results indicate that the selection rules for the allowed orders and SAM of generated harmonics in traditional bi-circular HHG are still upheld in the twisted laser scheme [15].
FIGURE 2
Moreover, it can be noticed here that only low order harmonics () are displayed in Figure 2, which can be explained as follows: according to the Shannon criterion and the resolution of the full three-dimensional PIC simulation, the highest resolvable harmonic order can be estimated as around 10. Thus, we only choose the harmonic with () for display with sufficiently high resolution. If the resolution is enhanced, much higher order harmonics can certainly be discernible.
In the following, the vortex characteristics of generated harmonics are investigated, which are displayed in Figures 3–5. The six panels of Figure 3 and Figure 4 show the electric field distributions in Z-Y plane atand the corresponding electric field isosurface of the fourth, fifth, seventh, eighth, tenth and eleventh harmonics atrespectively. According to the number of intertwined helices, one can see that the corresponding topological charges of the harmonics emitted are also four, five, seven, eight, ten and eleven, which are consistent with the harmonic orders (). These results indicate that the SAM and OAM conservation rules () which was analytically derived in Ref. [23] are also valid in our scheme (if we substitute our simulation parameters () into the SAM and OAM conservation rules, the topological charges of high harmonics can be easily acquired as). Moreover, the transverse intensity distribution (normalized to) and its corresponding profile with Y of the fourth harmonic (at the same time with Figures 3, 4) are shown in Figures 5A,B, from which the typical donut-shaped intensity distribution of optical vortices can be clearly observed. These results demonstrate that the radiated field by two-color LG laser scheme owns a complicated structure which contains both quasi-circular polarization and high-order LG modes scaling with the harmonic orders (the n-order harmonic exhibitsmode) [2]. In other word, besides the SAM, the OAM of driving lasers can also be transferred to the high-harmonics in this HHG scheme.
FIGURE 3
FIGURE 4
FIGURE 5
Up to now, it has been confirmed that the two-color LG laser-gas HHG scheme is also applicable in relativistic plasma HHG process, and the intensity of produced harmonics can be obviously enhanced.
Theoretical analysis
In this part, the physical mechanism of vortex harmonic radiation will be discussed. For the convenience of theoretical derivation, the expression of two-color counter-rotating twisted laser pulses can be simplified as follows [18]:whereis the propagating factor, is the duration of laser pulse, are the amplitudes of the normalized electric fields of laser pulses with the frequency of andare their corresponding topological charges.
By squaring the superposed laser electric field, the normalized laser intensity can be derived:where . From Eq. 4, we can see that the plasma electrons driven by dichromatic LG lasers can be described as a “vortex oscillating mirror” (VOM) [2]. On one side, it oscillates at the frequency of in the longitudinal direction, on the other side, it is helical in the azimuthal direction. As a result, structured high-harmonics will be efficiently produced due to the VOM effect.
The analysis above preliminarily explains the underlying mechanism of HHG through two-color LG laser scheme. However, they are not able to accurately predict the spectra features of produced harmonics (including the allowed harmonic orders, the polarization features and the vortex modes). To predict these, its characteristic selection rules can be acquired based on the conservation laws for energy, parity, SAM and OAM of photos. Considering the HHG process of our scheme in a simple photon-exchange picture, the conservation of energy, SAM and OAM can be expressed as [2, 15]:
In Eqs 5–7, andare integers associated with the number of driving photos annihilated in the HHG process of the driving pulses at angular frequenciesand. The parity constraints require thatmust be odd (). In addition, represent the SAM values (in units of ) of the emitted harmonic photon and the two driving lasers withand, and are the OAM values(also in units of ) of the radiated harmonics, the fundamental laser driver and the second one. According to our simulation parameters, we can set, and. Thus, Eqs 5–7 can be expressed as,and, respectively. Here, only two cases are allowed: 1) , resulting in, 2) , leading to(is equivalent to). These analytical results are consistent with the SAM and OAM of emitted harmonic radiation which are numerically shown in Figures 2–4.
Discussions
According to the conservation laws for energy, parity, SAM and OAM of photos, it can be easily seen that the vortex modes of emitted harmonics can be tuned just through changing the topological charges of driving pulses. For example, if we set (with other parameter unchanged), Eqs 5–7 can be expressed as,and . Similarly, if, we can get. Whilecan be acquired when. Thus, some new vortex modes which are forbidden withmay be acquired. Specifically, the vortex modes of the fourth and thirteenth harmonics will berespectively, as shown in Figure 6. This demonstrates that the OAM of generated harmonics can be flexibly designed through the bi-circular twisted laser scheme.
FIGURE 6
Finally, the tunability of HHG efficiency (which acts as the key parameter of HHG process) through adjusting the intensity ratio of the two driving field components are also discussed, as displayed in Figure 7. Here the total intensity of the two driving pulses is kept constant (, whereare the intensity of the driving pulses with frequency of () respectively, as introduced in the part of particle-in-cell simulation results) while the intensity ratiois varying from 0 to 0.1, 0.33 and 1. The remaining parameters are kept the same as those in Figures 2–5. The Lissajous figures of the superposed electric fields are shown in Figures 7A–D, and the corresponding harmonic spectra are displayed in Figures 7E–H. When, almost no high harmonics appear (as shown in Figure 7E) due to the quasi-CP driving field (as shown in Figure 7A). As the intensity ratio is increased, the superposed driving fields exhibit three-fold symmetry and gradually show linearly-polarized (LP) patterns (as Figures 7B–D display), leading to gradually enhanced HHG efficiencies (as Figures 7F–H display). For the case of , the efficiency is close to thescaling, as the typical value of HHG efficiency via relativistically oscillating mirror (ROM) effect [10]. This demonstrates that the HHG efficiency of our bi-circular twisted laser scheme can also be sufficiently high. In addition, the tunability of HHG efficiency can be attributed to the different driving field vector pattern induced by different laser intensity ratio [15].
FIGURE 7
Conclusion
With three-dimensional PIC simulations, we have demonstrated a tunable, table-top, intense harmonic radiation source with designable SAM and OAM driven by bi-circular LG lasers. The HHG efficiency can be adjusted through changing the intensity ratio of the two driving laser components, which is easy to implement in experiments. The underlying mechanism can be understood by the VOM model and related conservation laws. Such bright, circularly polarized vortex beams in XUV or soft X-ray region may open a wide range of applications including spatially resolved circular dichroism spectroscopies and the SAM and OAM dichroism X-ray absorption measurements, etc.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
DX and HZ conducted the simulations and drafted the manuscript. YY and TY supervised the work. All authors discussed the results and reviewed the manuscript.
Funding
This work was supported by the National Key R&D Program of China (Grant No. 2018YFA0404802), Scientific Research Foundation of Hunan Provincial Education Department (No. 20A042), National Natural Science Foundation of China (Grant Nos.. 12135009, 11774430 and 11875319) and Science and Technology Innovation Program of Hunan Province (Grant No. 2020RC4020).
Acknowledgments
We would like to thank the National Supercomputing Center in Guangzhou (NSCC-GZ) for providing their computing facilities.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
AllenLBeijersbergenMWSpreeuwRJCWoerdmanJP. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Phys Rev A (Coll Park) (1992) 45:8185–9. 10.1103/physreva.45.8185
2.
ZhangX-MShenB-FShiYWangX-FZhangL-GWangW-Pet alGeneration of intense high-order vortex harmonics. Phys Rev Lett (2015) 114:173901. 10.1103/physrevlett.114.173901
3.
JuLBHuangTWLiRJiangKWuCNZhangHet alTopological control of laser-driven acceleration structure for producing extremely bright ion beams. Nucl Fusion (2021) 61:066006. 10.1088/1741-4326/abeed6
4.
GauthierDRebernikPAdhikaryGCamperAChappuisCCuciniRet alTunable orbital angular momentum in high-harmonic generation. Nat Commun (2017) 8:14971. 10.1038/ncomms14971
5.
Hernandez-GarciaCPiconASan RomanJPlajaL. Attosecond extreme ultraviolet vortices from high-order harmonic generation. Phys Rev Lett (2013) 111:083602. 10.1103/physrevlett.111.083602
6.
QuereFThauryCMonotPDoboszSMarinPGeindreJPet alCoherent wake emission of high-order harmonics from overdense plasmas. Phys Rev Lett (2006) 96:125004. 10.1103/physrevlett.96.125004
7.
NomuraYHorleinRTzallasPDromeyBRykovanovSMajorZet alAttosecond phase locking of harmonics emitted from laser-produced plasmas. Nat Phys (2009) 5:124–8. 10.1038/nphys1155
8.
BulanovSVNaumovaNMPegoraroF. Interaction of an ultrashort, relativistically strong laser pulse with an overdense plasma. Phys Plasmas (1994) 1:745–57. 10.1063/1.870766
9.
LichtersRMeyer-ter-VehnJPukhovA. Short-pulse laser harmonics from oscillating plasma surfaces driven at relativistic intensity. Phys Plasmas (1996) 3:3425–37. 10.1063/1.871619
10.
von der LindeDRzazewskiK. High-order optical harmonic generation from solid surfaces. Appl Phys B (1996) 63:499–506. 10.1007/bf01828947
11.
BaevaTGordienkoSPukhovA. Theory of high-order harmonic generation in relativistic laser interaction with overdense plasma. Phys Rev E (2006) 74:046404. 10.1103/physreve.74.046404
12.
ThauryCQuéréFGeindreJ-PLevyACeccottiTMonotPet alPlasma mirrors for ultrahigh-intensity optics. Nat Phys (2007) 3:424–9. 10.1038/nphys595
13.
DromeyBKarSBelleiCCarrollDCClarkeRJGreenJSet alBright multi-kev harmonic generation from relativistically oscillating plasma surfaces. Phys Rev Lett (2007) 99:085001. 10.1103/physrevlett.99.085001
14.
ChenZ-YPukhovA. Bright high-order harmonic generation with controllable polarization from a relativistic plasma mirror. Nat Commun (2016) 7:12515. 10.1038/ncomms12515
15.
ChenZ-Y. Spectral control of high harmonics from relativistic plasmas using bicircular fields. Phys Rev E (2018) 97:043202. 10.1103/physreve.97.043202
16.
XieDZhuo*H-BJiaoJ-LZhangS-JZhaoNZhouH-Y. 3D particle simulation of harmonic radiation excited by double circularly-polarized laser pulses irradiating over-dense plasma surface. Appl Phys B (2019) 7:148. 10.1007/s00340-019-7263-3
17.
XieDYinYYuT-PZhouH-YChenZ-YZhuoH-B. High-harmonic generation driven by two-color relativistic circularly polarized laser pulses at various frequency ratios. Plasma Sci Technol (2021) 23:045502. 10.1088/2058-6272/abe848
18.
LiQ-NXuX-RWuY-BYinYZouD-BYuT-P. Efficient high-order harmonics generation from overdense plasma irradiated by a two-color co-rotating circularly polarized laser pulse. Opt Express (2022) 30:15470–81. 10.1364/oe.459866
19.
ChenZ-YLiX-YLiB-YChenMLiuF. Isolated elliptically polarized attosecond soft X-ray with high-brilliance using polarization gating of harmonics from relativistic plasmas at oblique incidence. Opt Express (2018) 26:4572–80. 10.1364/oe.26.004572
20.
PukhovAan der BruggeDKostyukovI. Relativistic laser plasmas for electron acceleration and short wavelength radiation generation. Plasma Phys Control Fusion (2010) 52:124039–10. 10.1088/0741-3335/52/12/124039
21.
an der BruggeDPukhovA. Enhanced relativistic harmonics by electron nanobunching. Phys Plasmas (2010) 17:033110. 10.1063/1.3353050
22.
DromeyBRykovanovSYeungMHorleinRJungDGautierDCet alCoherent synchrotron emission from electron nanobunches formed in relativistic laser-plasma interactions. Nat Phys (2012) 8:804–8. 10.1038/nphys2439
23.
DorneyK-MRegoLBrooksN-JRomanJ-SLiaoC-TEllisJ-Let alControlling the polarization and vortex charge of attosecond high-harmonic beams via simultaneous spin-orbit momentum conservation. Nat Photon (2019) 13:123–30. 10.1038/s41566-018-0304-3
24.
ArberTDBennettKBradyCSLawrence-DouglasARamsayMGSircombeNJet alContemporary particle-in-cell approach to laser-plasma modelling. Plasma Phys Control Fusion (2015) 57:113001. 10.1088/0741-3335/57/11/113001
Summary
Keywords
vortex beams, Laguerre-Gaussian lasers, relativistic plasmas, orbital angular momentum, vortex mirror model, conservation laws
Citation
Xie D, Yin Y, Yu T, Zhang H and Zhou H (2022) High-order vortex harmonics generation by bi-circular Laguerre-Gaussian laser fields with relativistic plasmas. Front. Phys. 10:962956. doi: 10.3389/fphy.2022.962956
Received
07 June 2022
Accepted
01 August 2022
Published
30 August 2022
Volume
10 - 2022
Edited by
Jinqing Yu, Hunan University, China
Reviewed by
Xiaofei Shen, Heinrich Heine University of Düsseldorf, Germany
Willem Boutu, CEA Saclay, France
Updates
Copyright
© 2022 Xie, Yin, Yu, Zhang and Zhou.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Duan Xie, 627065391@qq.com
This article was submitted to Fusion Plasma Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.