Hostname: page-component-cd9895bd7-gvvz8 Total loading time: 0 Render date: 2024-12-24T16:19:16.708Z Has data issue: false hasContentIssue false

Bright High-Harmonic Generation through Coherent Synchrotron Emission Based on the Polarization Gating Scheme

Published online by Cambridge University Press:  01 January 2024

Chuliang Zhou
Affiliation:
State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai 201800, China Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
Ye Tian*
Affiliation:
State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai 201800, China Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
Yushan Zeng
Affiliation:
State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai 201800, China Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
Zhinan Zeng*
Affiliation:
State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai 201800, China Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
Ruxin Li
Affiliation:
State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Shanghai 201800, China Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China School of Physical Science and Technology, ShanghaiTech University, Shanghai 201210, China
*
Correspondence should be addressed to Ye Tian; [email protected]
Correspondence should be addressed to Ye Tian; [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Relativistic surface high harmonics, combined with the use of polarization gating, present a promising route towards intense single attosecond pulses. However, they impose stringent requirements on ultra-high laser contrast and are restricted by large intensity losses in real experiments. Here, we numerically demonstrate that by setting an optimal time delay in the polarization gating scheme, the intensity of the generated single attosecond pulses can become approximately 100 times stronger than that with nonoptimal time delay in the coherent synchrotron emission process. When a petawatt-class driving laser irradiates a solid target, an ultra-dense electron nanobunch and a strong space-charge sheath develop, and the accumulated electrostatic energy is only released in half of the laser cycle when this electron nanobunch moves backward. This process results in the emission of intense high harmonics. Our study provides a reliable method for developing bright attosecond extreme ultraviolet pulses.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © 2022 Chuliang Zhou et al.

1. Introduction

High-harmonic generation (HHG) from an ultra-intense laser solid-target interaction has been proven to be an important method for producing bright coherent extreme ultraviolet rays and X-rays [Reference Dromey, Zepf and Gopal1Reference Thaury, Quéré and Geindre3]. Generation of attosecond pulse trains on the solid surface has also recently been demonstrated experimentally [Reference Nomura, Hörlein and Tzallas4, Reference Hörlein, Nomura and Tzallas5]. Moreover, the coherent wake emission (CWE) [Reference Thaury, Quéré and Geindre3, Reference Nomura, Hörlein and Tzallas4, Reference Quéré, Thaury and Monot6, Reference Quéré, Thaury, Geindre, Bonnaud, Monot and Martin7] and relativistic oscillating mirror (ROM) [Reference Dromey, Zepf and Gopal1, Reference Tsakiris, Eidmann, Meyer-ter-Vehn and Krausz8Reference Lichters, Meyer‐ter‐Vehn and Pukhov14] models successfully describe the high-harmonics generation process at laser intensities of I < 1018W/cm2 and I > 1018W/cm2, respectively. In recent years, nanobunching of relativistic electrons has been identified to be responsible for coherent synchrotron emission (CSE) [Reference Dromey, Cousens and Rykovanov15Reference Brügge and Pukhov18]. An alternative theory of relativistic electronic springs (RESs) has also been proposed [Reference Gonoskov, Korzhimanov, Kim, Marklund and Sergeev19]. In the CSE model, dense nanometer-scale bunches of electrons accelerated by the laser and plasma fields in relativistic trajectories radiate high-frequency light. The most widely investigated scheme is the ROM mechanism, because of the relativistic Doppler upshift. Light reflected from a periodic oscillating mirror comprises a broad harmonic spectrum. For near-normal incidence, the plasma oscillations are driven by the Lorentz force of the incident laser pulse and charge-separation-induced electrostatic fields. For a laser pulse with linear polarization, the Lorentz force varies as F p a 2 t 1 + sin 2 ω L t , where a 2 t = I L t λ L 2 / 1.37 × 10 18 W μ m 2 cm 2 is the normalized vector potential. In the case of the circularly polarized pulse, this force is of the form F p a 2 t . Thus, it exhibits no fast oscillations, creating a smooth density depression, with which no harmonics is allowed to generate.

The critical requirement for effective harmonic generation in high-intensity laser systems is the effective control over the purity of the laser pulses to avoid any effect of the pedestal on the target. It has been verified that the most efficient HHG happens when an optimal preplasma scale length Ls (0.05λL < L s < 0.2 λL) is introduced [Reference Thaury, Quéré and Geindre3, Reference Kahaly, Monchocé and Vincenti20, Reference Chopineau, Leblanc and Blaclard21]. In this case, an ultra-high laser contrast, usually higher than 10−10, will be required, which poses formidable experimental challenges to current high power laser facilities [Reference Zhang, Wu and Hu22]. Hence, a feasible method should be found to relax the severe requirement on the preplasma scale length.

Conventionally, the existence of a long scale, nonstep-like preplasma in front of the target would significantly decrease the coherence and conversion efficiency of the harmonics. To avoid such detrimental effects and obtain single attosecond pulses, one can use few-cycle laser pulses with carrier-envelope phase stabilization or alternatively use longer pulses with polarization gating, a well-known technique commonly applied in gas-phase HHG. However, in the pioneering work on polarization gating for plasma HHG [Reference Rykovanov, Geissler, Meyer-ter-Vehn and Tsakiris23], the density gradient scale length at the plasma-vacuum interface was not taken into consideration, which would have a significant effect on the plasma HHG.

Different polarization gating techniques have been proposed for many-cycle laser pulses that suffer from large intensity losses [Reference Rykovanov, Geissler, Meyer-ter-Vehn and Tsakiris23Reference Yeung, Bierbach and Eckner25]. Recently, noncollinear polarization gating [Reference Yeung, Bierbach and Eckner25] has been demonstrated, which is simpler and more practical than the conventional polarization gating [Reference Rykovanov, Geissler, Meyer-ter-Vehn and Tsakiris23, Reference Baeva, Gordienko and Pukhov24]. For a given gate, the amplitude in the linear region of the pulse is smaller when using a longer pulse duration for polarization gating. However, a large fraction of the energy cannot be exploited in the circular polarization part.

In this study, we discuss the generation of intense high harmonics using multicycle driving pulses, in combination with the polarization gating scheme (Figure 1), and even under nonoptimal plasma-scale-length conditions. The efficiency of HHG from laser-plasma interaction can be boosted by exploiting a polarization gating scheme, which results in a recession of the electron in preplasma due to the pressure of the laser. Subsequently, an ultra-dense electron bunch and a strong electrostatic potential are formed, resulting in strong CSE. It is verified that the polarization gating technique driven by a multicycle laser can be used to develop an intense extreme ultraviolet source effectively in laser systems with a lower contrast. Our findings can also be applied to noncollinear polarization gating [Reference Yeung, Bierbach and Eckner25]. In the noncollinear gating scheme, another noncollinear gating term that arises from the wavefront rotation due to the noncollinearity of the two beams will be introduced. But if the angle γ that the wave front of each half beam makes with the center axis is small (γ << 1), the effect of wavefront rotation can be ignored.

Figure 1: (a) Schematic of intense extreme ultraviolet pulse generation with polarization gating scheme, which consists of two counter-rotating circularly polarized (CP) laser pulses normally incident on a solid target. These collinear left and right CP laser pulses create a linear gate in the middle portion. (b) Illustration of the electron density profiles when the front CP laser part irradiates the target. The preplasma is compressed and an ultra-dense electron nanobunch is formed. (c) When the linear polarization of the laser reaches the target surface, the dense electron nanobunch is accelerated by the charge-separation and laser electric field to relativistic velocities in a half laser cycle; thus, intense extreme ultraviolet pulses would be emitted.

2. Methods

Polarization gating [Reference Rykovanov, Geissler, Meyer-ter-Vehn and Tsakiris23, Reference Corkum, Burnett and Ivanov26, Reference Platonenko and Strelkov27] is created by combining two same-color, counter-rotating, and circularly polarized multicycle Gaussian laser pulses with a certain time delay between them, as shown in Figure 1(a). This pulse possesses a time-dependent ellipticity. A gate of nearly linear polarization is opened when the two pulses have similar intensities. The peak field amplitude E 0, carrier frequency ω, full-width at half-maximum (FWHM) pulse duration τ, and carrier-envelope phase φ are identical for the two pulses, where the time delay between them is T d . The electric fields of the left and right circularly polarized pulses propagating in the x-direction are, respectively,

(1) E l = E 0 exp 2 ln 2 t τ 2 y cos ω t + φ + z sin ω t + φ ,
(2) E r = E 0 exp 2 ln 2 t T d τ 2 y cos ω t + φ z sin ω t + φ .

The combined electric field is E t = E l + E r . The ellipticity ξ(t) is defined as the ratio of the minor and major axes of the electric field vector ellipse, ξ t = E y / E z , and the evolution of ellipticity with time is given by [Reference Rykovanov, Geissler, Meyer-ter-Vehn and Tsakiris23]

(3) ξ t = exp 4 ln 2 × | t T d 2 t 2 τ 2 .

For a threshold ellipticity value of 0.4, the gating time is given by

(4) δ gating = 0.2 τ 2 T d .

The ellipticity evolution in polarization gating with different relative time delays is shown in Figure 2(a). When T d is close to the pulse duration, the gating time is about one laser cycle, the isolated attosecond pulse can be obtained by filtering out the low frequency. The relative delay will affect the gate width, the intensity, and the degree of polarization of the linearly polarized laser in the gate width. We performed particle-in-cell (PIC) simulations with the kinetic code EPOCH [Reference Arber, Bennett and Brady28] to study harmonic generation in relativistic laser interaction with overdense plasma. Two counter-rotating CP Gaussian pulses (its normalized intensity is a = eE 0/m e 0) were used as the driving pulses with frequencies of ω0 (λL = 800 nm) and copropagated along the x-direction. The two pulses both had a FWHM duration τ of 6T 0, where T0 = λ L / c is the laser period and c is the velocity of light in vacuum. When the time delay between them is 0, the linearly polarized laser is formed. In the following simulation, the two circularly polarized lasers are each composed of two linearly polarized laser pulses with equal intensity but different polarization: one is polarized orthogonally to the other with a π/2 phase shift such that their superposition can form a circularly polarized pulse. Therefore, under different time delays, the total energy of the incident laser remains unchanged. The solid target is considered fully ionized, which consists of a 2λL thick slab with n max = 100 n c and an exponentially declining preplasma with a scale length of Ls, where n c = ε 0 m e ω 0 2 /e 2 represents the critical density, and m e and e are the electron rest mass and charge, respectively. The laser duration and plasma density were the same for both one-dimensional (1D) and two-dimensional (2D) simulations. For the 1D simulation, the grid size was dx = λ L /2000, where each cell contained 1000 macroparticles. The ions were treated as a neutralizing immobile background. While for the 2D simulation, the grid size was dx = dy = λ L /800, and each cell contained 100 macroparticles. The reflected emission spectra were determined from the Fourier transform of the time-resolved reflected electric field at the vacuum gap.

Figure 2: (a) Evolution of the ellipticity in polarization gating with different relative time delays. (b) Electron density (n e ) right before the generation of the most intense attosecond pulse. PIC simulation parameters: a = 10, time delay T d = 16 fs, and Ls = 0.2 λ L . (c) Obtained radiation pulse intensity by applying a spectral filter to select the 30th to 100th harmonic orders. PIC simulation parameters: a = 10, time delay T d = 12 fs, and Ls = 0.1 λ L . (d, g, j) Dependence of the shortest width and number density of the electron bunch on the preplasma scale length, time delay, and laser intensity. (e, h, k) Dependence of efficiency of harmonic generation on preplasma scale length, time delay, and laser intensity. (f, i, l) The dependence of number and the peak intensity of gated attosecond pulses (30th–100th) on the preplasma scale length, time delay, and laser intensity. Simulation parameters (d, e, f) a = 10, time delay T d = 6 fs; (g, h, i) a = 10, Ls = 0.1 λ L ; (j, k, l) time delay T d = 6 fs, Ls = 0.1 λ L .

3. Results and Discussion

We first clarify the electron density modification caused by the ponderomotive force when the rising edge of the laser interacts with the target. Subjected to light pressure, the electrons in the preplasma are pushed towards the target at first, forming an electron-depleted layer. As a result, an electron density sheath is formed at the position of the vacuum-plasma interface, and the electron density can be substantially compressed to several nanometers (Figure 2(b)). During this stage, the longitudinal electric field component grows until the electrostatic force starts to pull back the bunch. The electrons are piled up in a narrow region under the effect of light pressure, and a charge-separation field is formed between the thin layer of ions and electron nanobunches on the target surface. The dependence of the width and number density of the electron bunch on the preplasma scale length and time delay are present in Figures 2(d) and 2(g); as shown, with increasing the preplasma scale length or time delay, the shortest width of the electron nanobunch has been kept around 2 nm, while the electron number density tends to decrease, which is more obvious in Figure 2(d). Accordingly, Figures 2(e) and 2(h) show the variation of efficiency of harmonic generation with preplasma scale length and time delay, respectively. It can be seen that the evolution of harmonic efficiency is consistent with the variation of nanoelectron layer density in Figures 2(d) and 2(g). Figure 2(j) shows the dependences of width and number density of the electron nanobunch on the intensity of incident laser, while Figure 2(k) gives the relation between the efficiency of high-harmonic generation and incident laser intensity. As the laser intensity increases, the number density of electron nanobunch will increase as well, while its shortest width shows yet no significant changes. This is analogous to the scenario of changing the preplasma scale length and time delay, the harmonic generation efficiency will increase as the laser intensity strengthens, which verified again the obvious positive correlation between the harmonic generation efficiency and the number density of electron layer. Hence, we can reasonably conclude that the higher electron density is favorable for coherent and stronger radiation. Owing to the relativistic character of electron motion, the displaced electrons can reemit their energy when the electrons propagate back towards the boundary, thereby compressing the reemission into an attosecond burst. An intense extreme ultraviolet pulse is emitted when the linear polarization of the incident laser reaches the target surface.

Figure 2(c) presents the temporal shape of the pulse intensity after passing a bandpass (30th∼100th) filter. The number of gated attosecond pulses (only the pulse intensity that is higher than the 1/e 2 of peak intensity is counted, which can represent the isolation degree) and the corresponding peak intensity for varied preplasma scale length, delays, and incident laser intensity are respectively presented in Figures 2(f)2(l). From the results displayed in Figure 2(f), the number of attosecond pulses enhances with increasing preplasma scale length, while a shorter preplasma scale length is more beneficial for obtaining isolated attosecond pulses. In Figure 2(i), the number and intensity of attosecond pulses both show a gradual decrease at longer delays (and thus shorter gate times), and for laser intensity, it can be seen in Figure 2(l) that the intensity of gated attosecond pulse traces approximately exponentially with the driving laser intensity, whereas the degree of isolation appears to be slightly less sensitive to laser intensity variation.

The present working scheme we proposed has been checked under different laser-plasma parameters. In Figure 3, we calculated the influence of plasma scale length Ls and time delay T d on high-order harmonic emission from 1D PIC simulations. Results for harmonics within the range from 10th to 20th (Figure 3(a)) and from 20th to 100th (Figure 3(b)) are presented. From these results, we can also see that when the preplasma scale length is too long (Ls > 0.4 λ L ), the harmonic efficiency will stay low regardless of the relative time delay. But this efficiency would become sensitive to the delay when a mildly long (Ls < 0.4λ L ) preplasma scale length is assumed. Furthermore, at larger preplasma scale length (beyond 0.2 λ L ), the optimal delay also seems to increase. From Figures 2(f) and 2(i) and Figures 3(c) and 3(d), we can conclude that the shorter preplasma scale length and shorter gate width are both beneficial to the generation of isolated attosecond pulses.

Figure 3: (a, b) Harmonics generation efficiency as a function of both plasma scale length and delay of the two circular polarized pulses in polarization gating. (c, d) The number and corresponding peak intensity of gated attosecond pulses (30th–100th) as a function of preplasma scale length and delay between the two circular polarized pulses in the polarization gating scheme. The intensity of the incident field is a = 10, the pulse duration τ is 16 fs, and solid-target density is 100 n c . Results for harmonics within the range from 10th to 20th (a) and from 20th to 100th (b) are presented.

In Figure 4, we compare the PIC simulation results respectively of two distinct time delays: 16 fs (where gating width is about one laser cycle) and T d = 0 (linearly polarized laser), while the other parameters are set to angle of incidence θ = 0, a = 10, and L s = 0.2 λ L . The corresponding results at these two time delays are shown in the left (16 fs) and right columns (T d = 0) in Figure 4, respectively. Figures 4(a) and 4(d) depict the electric field evolutions of the laser pulses in these two cases, inside which the amplitude of the reflected field exceeds that of the incoming field (for E y , it is up to 30%, while for E z , it is 90%) for the 16 fs scenario in Figure 4(a). We know that under different laser intensities, the harmonic generation efficiency is very different, and there is no simple linear correlation between the efficiency of harmonics generation and the intensity of the incident laser. As a result, a significant difference exists in the behavior of E y and E z. This is only the result of a 1D simulation in which the effect of plasma denting cannot be observed in the electric field. A curved surface can also boost incident lasers [Reference Quéré and Vincenti29]. This will be more obvious under the longer preplasma scale. The reflected pulse will be spatially focused, resulting in stronger field strength, whose divergence will also increase by certain extent. The attosecond pulse divergence dependents on the denting of the plasma surface which is not invariant over time. However, in our case, we only concern a small number of consecutive pulses gated from the pulse train. Hence, the divergence should not vary significantly. Owing to the relativistic character of electron motion, the displaced electrons are spatially squeezed to a high-density ultra-thin layer (Figure 2(b)), which can reemit their energy when propagating back towards the boundary, and thereby compressing the reemission into an attosecond burst with larger electric field amplitude than the incident light.

Figure 4: PIC simulation results for two different time delay cases: T d = 16 fs (left panel, originally thought to be suitable for isolated attosecond pulse generation) and T d = 0 (right panel, corresponding to linear polarized laser), other parameters: a = 10, L s = 0.2 λ L , and incidence angle θ = 0. (a, d) The incoming and reflected electric fields at a fixed position before the target surface, respectively. (b, e) Fourier spectra of the reflected pulses. The dashed line marks the universal scaling law I ω 8 / 3 predicted based on the Baeva–Gordienko–Pukhov (BGP) theory. (c, f) The time-frequency analyses of the HHG plotted on a logarithmic scale.

The corresponding harmonic spectra are shown in Figures 4(b) and 4(e), where the dashed line marks the universal scaling I ω 8 / 3 predicted by Baeva–Gordienko–Pukhov (BGP) [Reference Baeva, Gordienko and Pukhov10]. In both cases, the harmonic spectra remain close to the I ω 8 / 3 scaling law, which also implies the existence of zeros in the transverse momenta of plasma surface electrons. There are some fuzzy harmonic spectral structures observable in linear polarization. For the polarization gating scheme, the harmonic spectral structures are more discernible, especially for higher orders. The time-frequency analyses of the HHG are shown in Figures 4(c) and 4(f), where all frequencies are emitted within a fraction of the driving laser period and trains of attosecond pulses are emitted in the polarization gating scheme. But, we failed to obtain isolated attosecond pulses as expected. Thus, an optimal time delay still awaits to be found when a long preplasma is considered. The superposition of phase-locking-attosecond pulses in different periods contribute to the integrated harmonics spectrum. However, the phase between attosecond pulses generated in different periods will inevitably be accompanied by chirp and phase fluctuations, which will degrade the time coherence of attosecond pulses generated in many cycles (Figure 4(f)) than those generated in fewer cycles (Figure 4(c)). Also, the resulting phase-locking spectrum is not as clear as the polarization gating scheme.

Figure 5(a) shows the electron density overlaid with the longitudinal J x (x, t) current density components. The simulation parameters are the same as Figure 4(a). Here, these currents exhibit a periodic characteristic, in the sense that they undergo longitudinal oscillations twice inside each laser period. This apparent regularity emerges from the electron motion that is aperiodic at the single-particle level. Figure 5(b) shows that the harmonic pulses are emitted in the form of an attosecond pulse train, with two attosecond pulses emitted per optical cycle of the laser field. It is efficiently emitted at the moments when the longitudinal momentum of the electrons reaches its maximum. These results suggest that the mechanism is in accordance with the γ-spike model of the ROM.

Figure 5: (a) Spatiotemporal density evolution of the electrons overlaid with longitudinal current density J x (x, t) for the polarization gating case. The overall motion of the front surface consists of longitudinal oscillations twice every laser period. (b) Electron density on a logarithmic scale and the emission of harmonic pulse. The harmonic fields are frequency-filtered by N > 10th order, showing that the harmonic pulses are emitted in the form of an attosecond pulse train. Radiation is efficiently emitted at the moments when the longitudinal momenta of the electrons reach a maximum, which is also when both transverse momenta of the plasma surface electrons reach zero simultaneously.

The results in Figure 6 are based on 2D PIC simulations, where a same total laser energy is assumed (using the same laser and plasma quantities as in the 1D simulation, the transverse size of the focal spot (FWHM) is 7.5 λ L ). In view of the preplasma scale length in the forthcoming experiments (for example, Ls = 0.2 λ L ), we have comprehensively compared the efficiency of harmonic generation and the isolation degree of attosecond pulses under different delays. In all the cases where single attoseconds can be obtained after frequency filtering by selecting the 30th to 100th harmonic orders, we observed a significant increase and modulation of the harmonic spectrum, as shown in Figure 6(a) (red line) with proper time delays (6 fs). This result reveals complex plasma dynamics that have not been observed in a similar study [Reference Rykovanov, Geissler, Meyer-ter-Vehn and Tsakiris23]. Specifically, we found that an optimal time delay exists in polarization gating when a relatively long preplasma is introduced. The spectrum declines more slowly and shows clearer harmonic structures when T d is chosen as 6 fs in our case, and the intensity of HHG with optimal time delay can be increased by approximately 100 times when compared with the other nonoptimal time delay (here, the HHG intensity is also presented at 12 fs delay, where the gate times are shorter and are generally expected to produce single attosecond pulse more efficiently). The spectrum is very different from the BGP model. Efficient HHG from a relativistic plasma mirror requires (i) the longitudinal electron velocity reaching the maximum and, (ii) meanwhile, the transverse electron momentum reaching the minimum. When the delay is adjusted, the gate width, intensity, and time-varying ellipticity of the laser in the gate width will also be changed. The difference in HHG efficiency between the different delay cases can be attributed to the phase difference and field intensity, which leads to different plasma dynamics. On a macroscopic level, the spectrum is no longer well characterized by a universal decay index p = 8/3, or others. There is now a well-defined break in the spectrum in Figure 6(a), centered at n = ∼50th. From previous publications [Reference Kahaly, Monchocé and Vincenti20] and the results of Figure 2(e), the optimal plasma scale length in our case is ∼0.1 λ L . Also, from Figure 6(a), we can see that there exists an optimal time delay, even at a longer plasma scale length (here, we choose 0.2 λ L as an example), at which the harmonic structure in the spectrum is still clear and its intensity is also enhanced, suggesting that this process works well at a longer scale length than the typical optimal scale length. Thus, it opens the possibility of relaxing the severe requirement on the preplasma scale length. In the case of a polarization-gated pulse, the surface oscillations are highly suppressed during the circular polarization. The plasma denting effect, as mentioned before, was also observed. In addition, in Figures 6(b) and 6(c), it is observed that a train of attosecond electron bunches is emitted during the interaction. It should originate from the aforementioned electron nanobunch and shows some similarity with the electron bunching process [Reference Yan, Huang, Deng, Liu, Wang and Zhao30] in free-electron lasers.

Figure 6: 2D PIC simulation. (a) Intensity boost and modulation of harmonic spectrum in the polarization gating scheme. The intensity of HHG with optimal time delay (red curve) can be increased by approximately 100 times when compared with the other nonoptimal time delay (blue curve). The other simulation parameters are a = 10, τ = 6 T 0, n e = 100 n c , Ls = 0.2 λ L , and normal incidence. (b, c) Snapshots of electron density at different times when the time delay is 6 fs. The other simulation parameters are a = 10, τ = 6 T 0, n e = 100 n c , Ls = 0.2 λ L , and normal incidence.

Figures 7(a) and 7(b) present different spatiotemporal evolutions of the electron density, corresponding to the two different delays in Figure 6(a), at Ls = 0.2 λ L . It is shown that a more energetic electron bunch is formed and moves toward the incident laser when the T d is 6 fs. This indicates that the dominant radiation generation is the strong synchrotron emission from the ultra-relativistic electron bunch, which occurs when the compressed ultra-dense nanobunch electrons move back and are accelerated to the maximum longitudinal velocity by the combined field of space-charge field and the laser field. This shows that the harmonic generation can be boosted at an optimal time delay. The efficient CSE can be restricted to occur only once near the cycle of the peak laser intensity. At L s = 0.2 λ L and T d = 6 fs, the attosecond pulse profile is shown in Figure 7(c), which is filtered from the reflected field in the interval 30 < ω/ω 0 < 100. Figure 7(d) shows the zoomed single attosecond pulse with an envelope fitting (red), where the FWHM pulse duration is approximately 160 attoseconds. Its normalized intensity is boosted to ∼80, corresponding to 1.7 × 1020 W/cm2. These results show the potential of multicycle laser pulses to produce isolated and bright attosecond extreme ultraviolet and soft X-ray pulses.

Figure 7: (a, b) The spatiotemporal evolution of electron density with delay T d = 6 fs and 12 fs, respectively. (c) Obtained radiation pulse intensity by applying a spectral filter to select the 30th to 100th harmonic orders at T d = 6 fs. (d) Zoomed view of a single pulse with envelope fitting (red). The other simulation parameters are a = 10, τ = 6 T 0, L s = 0.2 λ L , n e = 100 n c , and normal incidence.

The ability to generate intense high harmonics with polarization gating is particularly important because it shows that the relativistic surface high harmonics have not only higher conversion efficiency but also fewer strict requirements on the driving field, making the relativistic surface high harmonics combined with polarization gating a very promising source of bright attosecond extreme ultraviolet and soft X-ray pulses. The proposed scheme is robust and can soon be tested in experiments with petawatt-femtosecond lasers.

4. Conclusions

We have presented results from a large number of particle-in-cell simulations, demonstrating that the polarization gating concept applied in relativistic surface high-harmonic generation can be used to overcome the limitation typically associated with the relatively long preplasma scale length. An ultra-dense electron bunch and a strong electrostatic potential are formed owing to the pondermotive force of the incident pulse, which results in a strong coherent synchrotron emission. We have demonstrated numerically that the intensity of high-harmonic generation from laser-plasma interaction can be enhanced while maintaining the same total laser energy. Specifically, we have shown that the time delay in the polarization gating concept is crucial for high-harmonic generation. The intensity of high-harmonic generation with optimal time delay can be boosted approximately 100 times when compared with the other nonoptimal time delay. The scheme proposed here points to the bright attosecond extreme ultraviolet pulse generation driven by petawatt-class multicycle laser pulses.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors would like to thank the National Supercomputing Center in Guangzhou (NSCC-GZ) for the computing source. This work was supported by the National Natural Science Foundation of China (grant nos. 11922412, 11874372, and 91950203), Strategic Priority Research Program (B) (grant no. XDB16), Key Research Program of Frontier Sciences, CAS (grant no. ZDBS-LY-SLH018), Youth Innovation Promotion Association of Chinese Academy of Sciences, China Postdoctoral Science Foundation (grant no. 2021M703329), and Shanghai Pilot Program for Basic Research–Chinese Academy of Science, Shanghai Branch.

References

Dromey, B., Zepf, M., Gopal, A. et al., “High harmonic generation in the relativistic limit,Nature Physics, vol. 2, no. 7, pp. 456459, 2006.10.1038/nphys338CrossRefGoogle Scholar
Pukhov, A., “X-rays in a flash,Nature Physics, vol. 2, no. 7, pp. 439440, 2006.10.1038/nphys347CrossRefGoogle Scholar
Thaury, C., Quéré, F., Geindre, J.-P. et al., “Plasma mirrors for ultrahigh-intensity optics,Nature Physics, vol. 3, no. 6, pp. 424429, 2007.10.1038/nphys595CrossRefGoogle Scholar
Nomura, Y., Hörlein, R., Tzallas, P. et al., “Attosecond phase locking of harmonics emitted from laser-produced plasmas,Nature Physics, vol. 5, no. 2, pp. 124128, 2009.10.1038/nphys1155CrossRefGoogle Scholar
Hörlein, R., Nomura, Y., Tzallas, P. et al., “Temporal characterization of attosecond pulses emitted from solid-density plasmas,New Journal of Physics, vol. 12, Article ID 043020, 2010.10.1088/1367-2630/12/4/043020CrossRefGoogle Scholar
Quéré, F., Thaury, C., Monot, P. et al., “Coherent wake emission of high-order harmonics from overdense plasmas,Physical Review Letters, vol. 96, Article ID 125004, 2006.10.1103/PhysRevLett.96.125004CrossRefGoogle ScholarPubMed
Quéré, F., Thaury, C., Geindre, J. P., Bonnaud, G., Monot, P., and Martin, P., “Phase properties of laser high-order harmonics generated on plasma mirrors,Physical Review Letters, vol. 100, Article ID 095004, 2008.10.1103/PhysRevLett.100.095004CrossRefGoogle ScholarPubMed
Tsakiris, G. D., Eidmann, K., Meyer-ter-Vehn, J., and Krausz, F., “Route to intense single attosecond pulses,New Journal of Physics, vol. 8, 2006.10.1088/1367-2630/8/1/019CrossRefGoogle Scholar
Bulanov, S. V., Naumova, N. M., and Pegoraro, F., “Interaction of an ultrashort, relativistically strong laser pulse with an overdense plasma,Physics of Plasmas, vol. 1, no. 3, pp. 745757, 1994.10.1063/1.870766CrossRefGoogle Scholar
Baeva, T., Gordienko, S., and Pukhov, A., “Theory of high-order harmonic generation in relativistic laser interaction with overdense plasma,Physical review. E, Statistical, nonlinear, and soft matter physics, vol. 74, Article ID 046404, 2006.10.1103/PhysRevE.74.046404CrossRefGoogle ScholarPubMed
Dromey, B., Kar, S., Bellei, C. et al., “Bright multi-keV harmonic generation from relativistically oscillating plasma surfaces,Physical Review Letters, vol. 99, Article ID 085001, 2007.10.1103/PhysRevLett.99.085001CrossRefGoogle ScholarPubMed
Dromey, B., Adams, D., Hörlein, R. et al., “Diffraction-limited performance and focusing of high harmonics from relativistic plasmas,Nature Physics, vol. 5, no. 2, pp. 146152, 2009.10.1038/nphys1158CrossRefGoogle Scholar
Plaja, L., Roso, L., Rzazewski, K., and Lewenstein, M., “Generation of attosecond pulse trains during the reflection of a very intense laser on a solid surface,Journal of the Optical Society of America B, vol. 15, no. 7, pp. 19041911, 1998.10.1364/JOSAB.15.001904CrossRefGoogle Scholar
Lichters, R., Meyer‐ter‐Vehn, J., and Pukhov, A., “Short‐pulse laser harmonics from oscillating plasma surfaces driven at relativistic intensity,Physics of Plasmas, vol. 3, no. 9, pp. 34253437, 1996.10.1063/1.871619CrossRefGoogle Scholar
Dromey, B., Cousens, S., Rykovanov, S. et al., “Coherent synchrotron emission in transmission from ultrathin relativistic laser plasmas,New Journal of Physics, vol. 15, 2013.10.1088/1367-2630/15/1/015025CrossRefGoogle Scholar
Dromey, B., Rykovanov, S., Yeung, M. et al., “Coherent synchrotron emission from electron nanobunches formed in relativistic laser-plasma interactions,Nature Physics, vol. 8, no. 11, pp. 804808, 2012.10.1038/nphys2439CrossRefGoogle Scholar
Pukhov, A., An der Brugge, D., and Kostyukov, I., “Relativistic laser plasmas for electron acceleration and short wavelength radiation generation,Plasma Physics and Controlled Fusion, vol. 52, 2010.10.1088/0741-3335/52/12/124039CrossRefGoogle Scholar
Brügge, D. a. d. and Pukhov, A., “Enhanced relativistic harmonics by electron nanobunching,Physics of Plasmas, vol. 17, Article ID 033110, 2010.Google Scholar
Gonoskov, A. A., Korzhimanov, A. V., Kim, A. V., Marklund, M., and Sergeev, A. M., “Ultrarelativistic nanoplasmonics as a route towards extreme-intensity attosecond pulses,Physical review. E, Statistical, Nonlinear, and Soft Matter Physics, vol. 84, Article ID 046403, 2011.10.1103/PhysRevE.84.046403CrossRefGoogle ScholarPubMed
Kahaly, S., Monchocé, S., Vincenti, H. et al., “Direct observation of density-gradient effects in harmonic generation from plasma mirrors,Physical Review Letters, vol. 110, no. 17, Article ID 175001, 2013.10.1103/PhysRevLett.110.175001CrossRefGoogle ScholarPubMed
Chopineau, L., Leblanc, A., Blaclard, G. et al., “Identification of coupling mechanisms between ultraintense laser light and dense plasmas,Phys.Rev.X, vol. 9, Article ID 011050, 2019.Google Scholar
Zhang, Z., Wu, F., Hu, J. et al., “The laser beamline in SULF facility,High Power Laser Science and Engineering, vol. 8, p. e4, 2020.10.1017/hpl.2020.3CrossRefGoogle Scholar
Rykovanov, S. G., Geissler, M., Meyer-ter-Vehn, J., and Tsakiris, G. D., “Intense single attosecond pulses from surface harmonics using the polarization gating technique,New Journal of Physics, vol. 10, no. 2, Article ID 025025, 2008.10.1088/1367-2630/10/2/025025CrossRefGoogle Scholar
Baeva, T., Gordienko, S., and Pukhov, A., “Relativistic plasma control for single attosecond x-ray burst generation,Physical review. E, Statistical, nonlinear, and soft matter physics, vol. 74, Article ID 065401, 2006.10.1103/PhysRevE.74.065401CrossRefGoogle ScholarPubMed
Yeung, M., Bierbach, J., Eckner, E. et al., “Noncollinear polarization gating of attosecond pulse trains in the relativistic regime,Physical Review Letters, vol. 115, Article ID 193903, 2015.10.1103/PhysRevLett.115.193903CrossRefGoogle ScholarPubMed
Corkum, P. B., Burnett, N. H., and Ivanov, M. Y., “Subfemtosecond pulses,Optics Letters, vol. 19, no. 22, pp. 18701872, 1994.10.1364/OL.19.001870CrossRefGoogle ScholarPubMed
Platonenko, V. T. and Strelkov, V. V., “Single attosecond soft-x-ray pulse generated with a limited laser beam,Journal of the Optical Society of America B, vol. 16, no. 3, pp. 435440, 1999.10.1364/JOSAB.16.000435CrossRefGoogle Scholar
Arber, T. D., Bennett, K., Brady, C. S. et al., “Contemporary particle-in-cell approach to laser-plasma modelling,Plasma Physics and Controlled Fusion, vol. 57, no. 11, Article ID 113001, 2015.10.1088/0741-3335/57/11/113001CrossRefGoogle Scholar
Quéré, F. and Vincenti, H., “Reflecting petawatt lasers off relativistic plasma mirrors: a realistic path to the Schwinger limit,High Power Laser Sci. Eng., vol. 9, 2021.10.1017/hpl.2020.46CrossRefGoogle Scholar
Yan, J., Huang, N., Deng, H., Liu, B., Wang, D., and Zhao, Z., “First observation of laser–beam interaction in a dipole magnet,Advanced Photonics, vol. 3, Article ID 045003, 2021.10.1117/1.AP.3.4.045003CrossRefGoogle Scholar
Figure 0

Figure 1: (a) Schematic of intense extreme ultraviolet pulse generation with polarization gating scheme, which consists of two counter-rotating circularly polarized (CP) laser pulses normally incident on a solid target. These collinear left and right CP laser pulses create a linear gate in the middle portion. (b) Illustration of the electron density profiles when the front CP laser part irradiates the target. The preplasma is compressed and an ultra-dense electron nanobunch is formed. (c) When the linear polarization of the laser reaches the target surface, the dense electron nanobunch is accelerated by the charge-separation and laser electric field to relativistic velocities in a half laser cycle; thus, intense extreme ultraviolet pulses would be emitted.

Figure 1

Figure 2: (a) Evolution of the ellipticity in polarization gating with different relative time delays. (b) Electron density (ne) right before the generation of the most intense attosecond pulse. PIC simulation parameters: a = 10, time delay Td = 16 fs, and Ls = 0.2 λL. (c) Obtained radiation pulse intensity by applying a spectral filter to select the 30th to 100th harmonic orders. PIC simulation parameters: a = 10, time delay Td = 12 fs, and Ls = 0.1 λL. (d, g, j) Dependence of the shortest width and number density of the electron bunch on the preplasma scale length, time delay, and laser intensity. (e, h, k) Dependence of efficiency of harmonic generation on preplasma scale length, time delay, and laser intensity. (f, i, l) The dependence of number and the peak intensity of gated attosecond pulses (30th–100th) on the preplasma scale length, time delay, and laser intensity. Simulation parameters (d, e, f) a = 10, time delay Td = 6 fs; (g, h, i) a = 10, Ls = 0.1 λL; (j, k, l) time delay Td = 6 fs, Ls = 0.1 λL.

Figure 2

Figure 3: (a, b) Harmonics generation efficiency as a function of both plasma scale length and delay of the two circular polarized pulses in polarization gating. (c, d) The number and corresponding peak intensity of gated attosecond pulses (30th–100th) as a function of preplasma scale length and delay between the two circular polarized pulses in the polarization gating scheme. The intensity of the incident field is a = 10, the pulse duration τ is 16 fs, and solid-target density is 100 nc. Results for harmonics within the range from 10th to 20th (a) and from 20th to 100th (b) are presented.

Figure 3

Figure 4: PIC simulation results for two different time delay cases: Td = 16 fs (left panel, originally thought to be suitable for isolated attosecond pulse generation) and Td = 0 (right panel, corresponding to linear polarized laser), other parameters: a = 10, Ls = 0.2 λL, and incidence angle θ = 0. (a, d) The incoming and reflected electric fields at a fixed position before the target surface, respectively. (b, e) Fourier spectra of the reflected pulses. The dashed line marks the universal scaling law I∝ω−8/3 predicted based on the Baeva–Gordienko–Pukhov (BGP) theory. (c, f) The time-frequency analyses of the HHG plotted on a logarithmic scale.

Figure 4

Figure 5: (a) Spatiotemporal density evolution of the electrons overlaid with longitudinal current density Jx(x, t) for the polarization gating case. The overall motion of the front surface consists of longitudinal oscillations twice every laser period. (b) Electron density on a logarithmic scale and the emission of harmonic pulse. The harmonic fields are frequency-filtered by N > 10th order, showing that the harmonic pulses are emitted in the form of an attosecond pulse train. Radiation is efficiently emitted at the moments when the longitudinal momenta of the electrons reach a maximum, which is also when both transverse momenta of the plasma surface electrons reach zero simultaneously.

Figure 5

Figure 6: 2D PIC simulation. (a) Intensity boost and modulation of harmonic spectrum in the polarization gating scheme. The intensity of HHG with optimal time delay (red curve) can be increased by approximately 100 times when compared with the other nonoptimal time delay (blue curve). The other simulation parameters are a = 10, τ = 6 T0, ne = 100 nc, Ls = 0.2 λL, and normal incidence. (b, c) Snapshots of electron density at different times when the time delay is 6 fs. The other simulation parameters are a = 10, τ = 6 T0, ne = 100 nc, Ls = 0.2 λL, and normal incidence.

Figure 6

Figure 7: (a, b) The spatiotemporal evolution of electron density with delay Td = 6 fs and 12 fs, respectively. (c) Obtained radiation pulse intensity by applying a spectral filter to select the 30th to 100th harmonic orders at Td = 6 fs. (d) Zoomed view of a single pulse with envelope fitting (red). The other simulation parameters are a = 10, τ = 6 T0, Ls = 0.2 λL, ne = 100 nc, and normal incidence.