Enhanced XUV Harmonics Generation with an Intense Laser Field in the Overdriven Regime

: High-order harmonic generation with high photon ﬂux has been a challenging task in strong-ﬁeld physics. According to the high-order harmonic generation process, the essential requirements for achieving efﬁcient harmonic radiations inside a gas medium are the improvement of the induced atomic dipole moment amplitude of the single-atom response in the microscopic and the phase matching of the high harmonics in the macroscopic medium. In this work, we demonstrated a feasible approach to enhance the extreme-ultraviolet harmonics in the plateau region by increasing the intensity of the driving laser while keeping the laser energy constant. The simulation results showed that by increasing the laser intensity to the overdriven regime, the average extreme-ultraviolet harmonics yield in the plateau region is approximately twice as high as that obtained optimally in the conventional loose focusing geometry scheme by utilizing a relatively low-intensity driving laser with the same laser energy. The quantitative analysis of the harmonics generation process in the macroscopic medium and the phase matching revealed that the observed enhancement in harmonics can be attributed to the ampliﬁcation of the induced atomic dipole moment amplitude of the single-atom response in the high-intensity driving laser and the favorable transient phase matching in the overdriven regime. Furthermore, the investigation of the driving laser indicated that the favorable transient phase matching is caused by the spatiotemporal reshaping of the driving laser in the overdriven regime.


Introduction
High-order harmonic generation (HHG) resulting from the interaction of intense laser pulses with a gaseous medium has become a promising tabletop coherent attosecond light source spanning from the extreme-ultraviolet (XUV) to water window X-ray [1][2][3][4][5][6][7][8][9].With such light sources, many ultrafast processes, such as ultrafast electronic dynamics inside atoms, time-resolved X-ray absorption fine structure spectroscopy, and nanoscale transient gratings excitation can be observed [10][11][12][13].Usually, high photon flux HHG is essential for these applications.However, due to the low conversion efficiency in gas HHG, the photon flux of the high harmonics is still low to nJ-level in the water window [14].According to the semiclassical three-step model [15], the cutoff photon energy of HHG is given by E cuto f f = I p + 3.17U p , where I p is the ionization potential of the gas medium and U p ∝ Iλ 2 denotes the pondermotive potential of the driving laser field.Here I and λ correspond to the intensity and wavelength of the driving laser, respectively.The feasible approaches to extend the HHG cutoff photon energy are to use long wavelength lasers or increase the laser intensity under certain conditions.However, the conversion efficiency of HHG from a single atom drops significantly with long wavelengths, e.g., λ −(4∼6) [16] and will be further reduced when taking the macroscopic propagation of the harmonics in the gas medium into consideration [17,18].Although this unfavorable scaling law could Photonics 2023, 10, 964 2 of 16 be partially compensated by increasing the gas pressure of the gas medium [14,19,20], shorter wavelength lasers are still preferable for generating stronger HHG with high photon flux [21].
Therefore, increasing the driving laser intensity moderately is another feasible way to extend the cutoff photon energy of the HHG.Usually, in a weakly ionized gas medium, the driving laser intensity is low enough that the phase mismatch caused by the plasma dispersion can be compensated by neutral atom dispersion.This mechanism is referred to as the conventional phase matching and the corresponding ionization level is below a "critical" ionization fraction [19,22].However, as the laser intensity increases, the gas medium would be ionized excessively and the plasma dispersion cannot be balanced by the neutral atom dispersion, which will cause a phase mismatch between the harmonics and the driving laser during the macroscopic propagation in the medium and prevent efficient buildup of the harmonic field [19,23].Several quasi-phase matching methods [24], including the use of multiple gas jets [25], counterpropagating lasers [26], or a periodically modulating waveguide [27,28], have been proposed to overcome this conventional phase matching limit and enhance the HHG yield.However, these quasi-phase matching methods are complicated to implement experimentally and limited to relatively low harmonic order [24].Recently, an alternative method was proposed to extend the phase matching condition in a high-intensity laser field [5,[29][30][31].This method takes advantage of laser defocusing in a strongly ionized gas medium and is commonly referred to as the overdriven regime [23].In the ordinary regime, the presence of ionization-induced plasma always hinders the scaling of HHG towards higher cutoff photon energy, as the plasma defocusing leads to a decrease in intensity and a blueshift in the frequency of the driving laser pulses.However, in the overdriven regime, this large plasma dispersion can also cause a significant spatiotemporal reshaping of the driving laser pulses, resulting in a transient buildup of harmonics over a distance of a few hundred micrometers (µm) [14,32].This propagation effect of the driving laser in the overdriven regime and the transient phase matching of the HHG can be exploited for the generation of isolated attosecond pulses [33,34] and for the optimization of HHG in the XUV region [29] up to the water window [5,31].For instance, Johnson et al. obtained sub-pJ soft X-ray emissions experimentally in the water window region in the overdriven regime [5], where the ionization fraction in the phase matching region is around 1%, significantly higher than the "critical" level.Schötz et al. demonstrated in experiments the generation of isolated attosecond pulses with a pulse width of 323-as at 80 eV in the overdriven regime, where the ionization fraction was up to 50% [33].In their studies, the optimized harmonics were found to be located near the cutoff region.However, as the driving laser intensity passes through the gas medium, it is strongly clamped due to the plasma defocusing, resulting in a transient phase matching that is restricted only to the entrance of the gas medium.As a result, the harmonic yield buildup near the cutoff region occurs rapidly at the entrance of the medium but then descends gradually along the propagation axis [32].Meanwhile, as the gas pressure increases, the plasma defocusing in the overdriven regime also leads to a decrease in laser intensity and cutoff photon energy along the propagation axis [35].
In addition, both theoretical and experimental findings have demonstrated that by increasing the intensity of the driving laser, it is possible to not only extend the cutoff photon energy but also significantly enhance the harmonic yield in the plateau region [29,31,35,36].However, the intensity of the driving laser in these studies was adjusted by attenuating laser energy using polarizers or by truncating the laser beam with an iris.Although the driving laser intensity can be precisely adjusted by a polarizer or iris in experiments, this implies that different laser intensities correspond to varying laser energy.Hence, the effects of the different laser intensities on the harmonic yield of the plateau region under the same laser energy have not been verified.
Therefore, different from the previous studies that focus on harmonic generation in the cutoff region, in this paper we present a theoretical study on the generation of XUV harmonics in the plateau region with different laser intensities while keeping the laser energy constant.In experiments, different laser intensities with the same laser energy can be achieved by using a focal mirror with different focal lengths.The optimal harmonic emission is extremely sensitive to the experimental parameters, such as laser field intensity and gas pressure [29].Therefore, in this work, we investigated the HHG driven by different laser intensities at different gas pressures by changing the laser focal spot size while keeping the peak power constant and obtained the optimal XUV harmonics emission within the plateau region in the overdriven regime.The single-atom response in our numerical study is calculated by the strong field approximation (SFA) model, and the macroscopic propagation of the driving laser and harmonic fields is obtained by solving the three-dimensional Maxwell's wave equations with the full electric field.Unlike the phase matching of the cutoff harmonics that occurred at the entrance of the gas medium [5,33], the phase matching of harmonics in the plateau region, as observed in our study, occurs in the second half of the gas medium.Furthermore, our findings revealed a gradual decrease in the phase mismatching of the XUV harmonics in the plateau region along the propagation axis in the overdriven regime, leading to a harmonic buildup over the entire gas medium.Therefore, by increasing the laser intensity to the overdriven regime, the average XUV harmonics yield in the plateau region is approximately twice as high as that obtained optimally in the conventional loose focusing geometry scheme by utilizing a relatively low-intensity driving laser with the same laser energy.This enhanced XUV harmonics in our study can be attributed to the enhancement of the single-atom response at an intense laser field and the transient phase matching of the harmonics caused by the spatiotemporal reshaping of the driving laser pulses in the overdriven regime.

The Macroscopic Propagation of the Driving Laser and Harmonic Fields
In our numerical model, the single-atom response of the harmonic radiation is calculated using the SFA, or the so-called Lewenstein model [37].This approach has been shown to provide good qualitative agreement with the numerical solution of the Schrödinger equation.The time-dependent induced nonlinear dipole moment due to a linearly polarized laser field is described as In Equation ( 1), E(t) is the electric field of the laser pulse, A(t) is its associated vector potential, is a positive regularization constant, p st and S are the stationary values of the momentum and quasi-classical action with the vector potential of the laser field, g(t) is the ground state amplitude calculated with the ADK theory [38], and d x is the dipole matrix element for bound-free transitions.It can be seen from Equation (1) that the dipole moment of the single atom is positively correlated with the laser electric field E(t), so the induced atomic dipole moment amplitude of the single-atom response can be enhanced by increasing the laser intensity, which contributes to the harmonic yield buildup.
In order to simulate the macroscopic effects of the harmonics generated in the gas medium, the propagation of the fundamental laser and harmonic fields should be considered, which can be described by Maxwell's equations.The single-atom response is inserted as a source term into the wave equations of the fundamental and the harmonic fields, which are solved in cylindrical coordinates and assuming radial symmetry.The propagation of the fundamental driving laser and the harmonic field in a gas medium is described by [39][40][41] ∇ where the source term P ion (r, z, t) in Equation ( 2) denotes the effective polarization of the driving field under the strong ionization, which includes neutral atomic dispersion, the Kerr effect, ionization, absorption of the ionized medium, and plasma defocusing.The source term P(r, z, t) of the harmonic fields in Equation ( 3) is related to the single-atom induced dipole moment in Equation ( 1) and depends upon the fundamental driving laser field E 1 (r, z, t).In the harmonic propagation equation, we have neglected the free-electron dispersion because the plasma frequency is much lower than the harmonics frequencies.Equations ( 2) and ( 3) can be solved by the Crank-Nicolson method in a moving coordinate frame (z = z, t = t − z/c) [40].

Transient Phase Matching in the Overdriven Regime
The above theoretical model is sufficient to simulate the HHG in a macroscopic medium, while for a deeper insight into the underlying physical mechanism, the phase matching of the HHG should be investigated.Phase matching refers to the constructive accumulation of radiation from coherent sources along the beam propagation direction, which is not only composed of the wavevector mismatching between the fundamental field and harmonics but also involves the mismatching induced by the intrinsic phase of the HHG dipole moment.Since the HHG is an extremely nonlinear process, both the single-atom response and the macroscopic effects could cause the phase mismatching of the harmonics.Usually, the phase mismatching degree of the qth harmonic can be expressed as [42,43] ∆k q = ∆k d + ∆k n + ∆k p + ∆k g .
Here ∆k n (z, t) = 2πqP(z)[1 − p(z, t)]∆n/λ and ∆k p (z, t) = −P(z)N atm p(z, t)r e λ q 2 − 1 /q are the phase mismatching due to the dispersion of the neutral atomic and free electrons, respectively, in which P(z) is the gas pressure, p(z, t) is the ionization probability, ∆n is the refractive index difference per atmosphere between the qth harmonic and the drive laser with wavelength λ, and N atm is the initial atomic number density per atmosphere.∆k g (z, t) = −q∂Φ g (z, t)/∂z is the phase mismatching arising from the Gouy phase shift Φ g (z, t), resulting from the geometrical focusing of the drive laser beam.The Gouy focal phase variation in space for a specific time moment t i can be expressed by the argument of the driving laser field as Φ g (r, z, t i ) = arg[E 1 (r, z, t i )], which can be directly extracted from the simulation.Lastly, in the overdriven regime, the intrinsic harmonic dipole phase mismatching induced by the intensity-dependent trajectories in single-atom response can be given by [32,33] which contains the effect of the blueshift and intensity variation on the dipole contribution, where U p (z, t) is the ponderomotive potential of the driving laser and α i is a photon energy- dependent coefficient that differs for long and short trajectories and can be obtained by the semiclassical three-step model or SFA model [33,43].

Effects of the Laser Intensity and Gas Pressure on the HHG
In our simulations, we adopt an 800 nm laser field with a pulse energy of about 0.36 mJ and a pulse duration of 15 fs.The peak laser intensity was changed by using different focal beam waists in simulations.The peak laser intensities used in our simulations are approximately 7.5 × 10 14 W/cm 2 (7.5 I 0 ), 1.2 × 10 15 W/cm 2 (12 I 0 ), 1.5 × 10 15 W/cm 2 (15 I 0 ), and 2.0 × 10 15 W/cm 2 (20 I 0 ), which corresponds to the laser beam waists of 45 µm, 36 µm, 32 µm, and 28 µm, respectively.The laser pulse is focused at the center of the helium (He) gas target with a Gaussian gas pressure distribution, and the medium length in each kind of laser intensity is chosen to be equal to the Rayleigh length (z 0 ) of each beam waist, which is about 8 mm, 5 mm, 4 mm, and 3 mm, respectively.To optimize harmonic yields in the overdriven regime, the gas pressure was varied.Figure 1a-d showed the simulated macroscopic HHG spectra of He atoms driven by the four kinds of laser intensities at different gas pressures.The harmonic yields are collected in the whole exit plane of the gas medium.For each laser intensity, one can see that with an increase in gas pressure, the harmonic yields in the plateau are progressively increased and the harmonic cutoff energy is slightly decreased [35,44].Therefore, an optimal pressure for each kind of laser intensity can be achieved to reach the best balance between the harmonic yield and plateau spectral range.The optimal pressures for the four driving laser intensities are approximately 0.3 bar, 0.45 bar, 0.6 bar, and 0.9 bar, respectively.approximately 7.5 × 10 14 W/cm 2 (7.5 I0), 1.2 × 10 15 W/cm 2 (12 I0), 1.5 × 10 15 W/cm 2 (15 I0), and 2.0 × 10 15 W/cm 2 (20 I0), which corresponds to the laser beam waists of 45 μm, 36 μm, 32 μm, and 28 μm, respectively.The laser pulse is focused at the center of the helium (He) gas target with a Gaussian gas pressure distribution, and the medium length in each kind of laser intensity is chosen to be equal to the Rayleigh length (z0) of each beam waist, which is about 8 mm, 5 mm, 4 mm, and 3 mm, respectively.To optimize harmonic yields in the overdriven regime, the gas pressure was varied.Figure 1a-d showed the simulated macroscopic HHG spectra of He atoms driven by the four kinds of laser intensities at different gas pressures.The harmonic yields are collected in the whole exit plane of the gas medium.For each laser intensity, one can see that with an increase in gas pressure, the harmonic yields in the plateau are progressively increased and the harmonic cutoff energy is slightly decreased [35,44].Therefore, an optimal pressure for each kind of laser intensity can be achieved to reach the best balance between the harmonic yield and plateau spectral range.The optimal pressures for the four driving laser intensities are approximately 0.3 bar, 0.45 bar, 0.6 bar, and 0.9 bar, respectively.In order to investigate the effect of the laser intensity on the HHG, we summarized the optimal HHG spectrum in each kind of laser intensity, as shown in Figure 2a.On the one hand, it is evident that as the laser intensity increases gradually to 1.5 × 10 15 W/cm 2 , not only does the harmonic cutoff energy is greatly extended, but there is also a notable enhancement in the harmonic yields within the plateau region.On the other hand, when the laser intensity increases from 1.5 × 10 15 W/cm 2 to 2.0 × 10 15 W/cm 2 , despite the continued extension of the harmonic cutoff energy, there is a reduction in the harmonic yields within the plateau region.It can be seen from Figure 2a that by moderately increasing the driving laser intensity and gas pressure, the harmonic yields within the whole plateau region of the low-intensity case (7.5 × 10 14 W/cm 2 ) can be significantly enhanced.In addition, the optimal laser intensity and gas pressure, which correspond to the highest harmonic yield, vary across different harmonic orders.By striking a balance between yield In order to investigate the effect of the laser intensity on the HHG, we summarized the optimal HHG spectrum in each kind of laser intensity, as shown in Figure 2a.On the one hand, it is evident that as the laser intensity increases gradually to 1.5 × 10 15 W/cm 2 , not only does the harmonic cutoff energy is greatly extended, but there is also a notable enhancement in the harmonic yields within the plateau region.On the other hand, when the laser intensity increases from 1.5 × 10 15 W/cm 2 to 2.0 × 10 15 W/cm 2 , despite the continued extension of the harmonic cutoff energy, there is a reduction in the harmonic yields within the plateau region.It can be seen from Figure 2a that by moderately increasing the driving laser intensity and gas pressure, the harmonic yields within the whole plateau region of the low-intensity case (7.5 × 10 14 W/cm 2 ) can be significantly enhanced.In addition, the optimal laser intensity and gas pressure, which correspond to the highest harmonic yield, vary across different harmonic orders.By striking a balance between yield enhancement and harmonic orders, the optimal laser intensity and gas pressure in our simulations are determined to be 1.5 × 10 15 W/cm 2 and gas pressure is 0.6 bar, respectively, which corresponds to the highest harmonic yield in the plateau region ranging from H65 to H95, as shown in Figure 2a. Figure 2b,c showed the detailed harmonic spectrum and corresponding harmonic yields of H65 to H95 driven by the four different laser intensities.With the increase in harmonic order, the growth factor of harmonic yield gradually diminishes.The average harmonic yield between H65 and H95, driven by the laser intensity of 1.5 × 10 15 W/cm 2 , is approximately twice as high as that driven by the laser intensity of 7.5 × 10 14 W/cm 2 with the same laser energy.In order to quantitatively investigate the effects of different laser intensities at different gas pressures on the harmonics within the plateau region, the yield of a single harmonic order was analyzed in our simulations.The selected harmonic orders driven by the four laser intensities are shown by the arrow in Figure 2b, which was located around H80.These specific harmonic orders were selected because they are situated at the center of H65 to H95, and their corresponding harmonic yields closely approximate the average value.Additionally, the selected harmonic orders were not exactly the same due to the blueshift of the spectrum during macroscopic propagation [45][46][47].Hence, for the convenience of description, we uniformly refer to them as H80 in this paper.The yield of H80 in the plateau region driven by different laser intensities at different gas pressures is shown in Figure 2d.One can see that for each laser intensity, there is an optimal gas pressure that maximizes the harmonic yield, and the optimal gas pressure increases with the laser intensity.In addition, it was observed that increasing the laser intensity can also significantly enhance the harmonic yield, even at the same gas pressure, except for the low gas pressures (<0.4 bar).The maximum yield of H80 is approximately 1.5 × 10 8 with the laser intensity of 1.5 × 10 15 W/cm 2 , which is approximately twice as high as that obtained optimally in the laser intensity of 7.5 × 10 14 W/cm 2 with the same laser energy.Therefore, by moderately increasing the laser intensity, not only can the harmonic cutoff energy be extended, but there will also be a notable enhancement in the harmonic yields within the plateau region.To verify the accuracy of our model, we have compared the single-atom HHG at each driving laser intensity calculated based on the SFA theory and numerical solution of the time-dependent Schrödinger equation (TDSE), respectively (see Supplementary Materials).
enhancement and harmonic orders, the optimal laser intensity and gas pressure in our simulations are determined to be 1.5 × 10 15 W/cm 2 and gas pressure is 0.6 bar, respectively, which corresponds to the highest harmonic yield in the plateau region ranging from H65 to H95, as shown in Figure 2a. Figure 2b,c showed the detailed harmonic spectrum and corresponding harmonic yields of H65 to H95 driven by the four different laser intensities.With the increase in harmonic order, the growth factor of harmonic yield gradually diminishes.The average harmonic yield between H65 and H95, driven by the laser intensity of 1.5 × 10 15 W/cm 2 , is approximately twice as high as that driven by the laser intensity of 7.5 × 10 14 W/cm 2 with the same laser energy.In order to quantitatively investigate the effects of different laser intensities at different gas pressures on the harmonics within the plateau region, the yield of a single harmonic order was analyzed in our simulations.The selected harmonic orders driven by the four laser intensities are shown by the arrow in Figure 2b, which was located around H80.These specific harmonic orders were selected because they are situated at the center of H65 to H95, and their corresponding harmonic yields closely approximate the average value.Additionally, the selected harmonic orders were not exactly the same due to the blueshift of the spectrum during macroscopic propagation [45][46][47].Hence, for the convenience of description, we uniformly refer to them as H80 in this paper.The yield of H80 in the plateau region driven by different laser intensities at different gas pressures is shown in Figure 2d.One can see that for each laser intensity, there is an optimal gas pressure that maximizes the harmonic yield, and the optimal gas pressure increases with the laser intensity.In addition, it was observed that increasing the laser intensity can also significantly enhance the harmonic yield, even at the same gas pressure, except for the low gas pressures (<0.4 bar).The maximum yield of H80 is approximately 1.5 × 10 8 with the laser intensity of 1.5 × 10 15 W/cm 2 , which is approximately twice as high as that obtained optimally in the laser intensity of 7.5 × 10 14 W/cm 2 with the same laser energy.Therefore, by moderately increasing the laser intensity, not only can the harmonic cutoff energy be extended, but there will also be a notable enhancement in the harmonic yields within the plateau region.To verify the accuracy of our model, we have compared the single-atom HHG at each driving laser intensity calculated based on the SFA theory and numerical solution of the time-dependent Schrödinger equation (TDSE), respectively (see Supplementary Materials).In order to understand the laser intensity and gas pressure dependence obtained from the simulation shown in Figure 2, we revisited the harmonics generation process in the macroscopic medium with a one-dimensional (1D) model.For the qth harmonic, the instantaneous intensity generated at the exit of the medium, I q , can be expressed as [48,49] Here, L med is the length of the medium, z is the propagation coordinate, η(z) is the ionization probability, ρ(z) is the initial atomic number density.L abs (z) = 1/σ q ρ(z) is the absorption length, in which σ q is the absorption cross-section of the qth harmonic [50].ϕ q (z) is the accumulated phase difference between the driving laser and the generated qth harmonic field at the exit of the medium, and D q is the induced atomic dipole moment amplitude of the single-atom response in Equation ( 1), which is positively correlated with the laser intensity I(z).It can be seen from Equation (6) that, the qth harmonic intensity is positively correlated with the gas pressure and driving laser intensity, but negatively correlated with the ionization probability.In addition, the accumulated phase difference ϕ q (z) of the qth harmonic is also dependent on the gas pressure and laser intensity.Therefore, the HHG intensity can be optimized by changing the laser intensity and the gas pressure.
In order to investigate the dependence of the HHG intensity on the laser intensity and the gas pressure, we numerically integrated Equation (1) over the propagation distance for each laser intensity.The results are shown in Figure 3. Here, we have neglected the ϕ q (z) dependence (ϕ q (z) = 0) on the harmonic intensity for simplicity.Figure 3a,b showed the typical evolution of the on-axis laser intensity and ionization probability in the gas medium.One can see that the laser intensity decreased rapidly along the propagation distance due to the nonlinear propagation effects, such as spatial diffraction, material dispersion, nonlinear self-focusing, ionization, plasma defocusing, and so on, which lead to the spatial reshaping of the driving laser.Moreover, as the initial laser intensity increases, there is a faster decrease in laser intensity, resulting in the more pronounced spatial reshaping of the driving laser.The maximum ionization probabilities at the four laser intensities are about 0.75%, 5.5%, 13%, and 40%, respectively, which is much higher than the "critical" ionization probability in the conventional phase matching.The envelope of the shadow in Figure 3a (dark blue line) denotes the Gaussian gas pressure distribution used in the simulations.By inserting the laser intensity and ionization probability evolution in Figure 3a,b into Equation (6), the integrated intensity of H80 for different gas pressures at the four laser intensities is shown in Figure 3c.It can be seen that the evolution of the harmonic intensity calculated by Equation (6) for different gas pressures at the four laser intensities is similar to the results obtained in Figure 2d.Therefore, the harmonic intensity dependence on the gas pressure and laser intensity in Figures 2d and 3c can be explained by the balance between the induced atomic dipole moment amplitude, ionization probability, gas pressure, and absorption length.On the one hand, as the laser intensity increases, the induced atomic dipole moment amplitude D q of the single-atom response can be enhanced, which contributed to the increase in the HHG intensity.However, the ionization probability η(z) also increases with the laser intensity, which hinders the increase of the HHG intensity.Therefore, there exists an optimal intensity that makes the HHG intensity maximum for a certain gas pressure, as shown in Figures 2d and 3c.On the other hand, the atomic number density ρ(z) increases with the gas pressure, a high gas pressure contributes to generating the HHG with high intensity.However, the absorption length L abs is inversely proportional to the gas pressure, a high gas pressure corresponds to a short absorption length, which can also impede the increase of the HHG intensity.Therefore, there also exists an optimal gas pressure that makes the HHG intensity maximum for a certain laser intensity, as shown in Figures 2d and 3c.This dependence of the gas pressure and laser intensity on the HHG intensity shown in Figures 2d and 3c exhibits similarity to the findings reported in previous studies [48,51].
intensity.Therefore, there also exists an optimal gas pressure that makes the HHG intensity maximum for a certain laser intensity, as shown in Figures 2d and 3c.This dependence of the gas pressure and laser intensity on the HHG intensity shown in Figures 2d and  3c exhibits similarity to the findings reported in previous studies [48,51].Finally, it is worth noting that the optimal gas pressures in each laser intensity in Figure 3c are much higher than that obtained in Figure 2d, and this can be aAributed to the phase matching.The results in Figure 3c neglected the influence of the phase mismatching of the qth harmonic, which is typically significant at high gas pressure and high laser intensity.However, the results in Figure 2d are obtained by fully solving the harmonic propagation Equation ( 2), which intrinsically contains the phase mismatching of the harmonics.In addition, the optimal laser intensity to maximize the HHG intensity in Figure 3c is 1.2 × 10 15 W/cm 2 , which is different from the optimal intensity (1.5 × 10 15 W/cm 2 ) in Figure 2.This difference in optimal intensity can also be aAributed to phase matching, particularly in the overdriven regime.Therefore, the dynamic phase matching conditions will be investigated in the following.

Phase Matching Analysis of the Harmonic Emissions in the Overdriven Regime
Because the phase matching of the harmonic is dependent on the spatiotemporal evolution of the driving laser in the gas medium, we first examine the electric fields of the driving laser in the reference frame at three selected positions, one at the entrance, one at the middle, and one at the exit of the gas medium, as shown in Figure 4. Figure 4a-d shows the on-axis (r = 0) driving laser electric field of the four laser intensities at the entrance (red line), middle (blue dashed line), and exit (green dash-doAed line), and Figure 4e-h shows the spatiotemporal laser intensity distributions of the four driving lasers at the exit of the gas medium.When the input laser peak intensity is low (7.5 I0), the spatiotemporal evolution of the laser in the gas medium is similar to that in a vacuum, as shown in Figures 4a and 3a.This is because the ionization fraction is low (0.75%, as shown in Figure 3b) in this low-intensity laser case, which leads to a weak plasma defocusing effect and is unable to effectively reshape the spatiotemporal profile of the driving laser.Therefore, in the low-intensity case (7.5 I0), the spatiotemporal distribution of the driving laser during propagation typically demonstrates Gaussian profiles in both the temporal and radial direction, as shown in Figure 4e.But when the input laser peak intensity is high, Finally, it is worth noting that the optimal gas pressures in each laser intensity in Figure 3c are much higher than that obtained in Figure 2d, and this can be attributed to the phase matching.The results in Figure 3c neglected the influence of the phase mismatching of the qth harmonic, which is typically significant at high gas pressure and high laser intensity.However, the results in Figure 2d are obtained by fully solving the harmonic propagation Equation ( 2), which intrinsically contains the phase mismatching of the harmonics.In addition, the optimal laser intensity to maximize the HHG intensity in Figure 3c is 1.2 × 10 15 W/cm 2 , which is different from the optimal intensity (1.5 × 10 15 W/cm 2 ) in Figure 2.This difference in optimal intensity can also be attributed to phase matching, particularly in the overdriven regime.Therefore, the dynamic phase matching conditions will be investigated in the following.

Phase Matching Analysis of the Harmonic Emissions in the Overdriven Regime
Because the phase matching of the harmonic is dependent on the spatiotemporal evolution of the driving laser in the gas medium, we first examine the electric fields of the driving laser in the reference frame at three selected positions, one at the entrance, one at the middle, and one at the exit of the gas medium, as shown in Figure 4. Figure 4a-d shows the on-axis (r = 0) driving laser electric field of the four laser intensities at the entrance (red line), middle (blue dashed line), and exit (green dash-dotted line), and Figure 4e-h shows the spatiotemporal laser intensity distributions of the four driving lasers at the exit of the gas medium.When the input laser peak intensity is low (7.5 I 0 ), the spatiotemporal evolution of the laser in the gas medium is similar to that in a vacuum, as shown in Figures 3a and 4a.This is because the ionization fraction is low (0.75%, as shown in Figure 3b) in this low-intensity laser case, which leads to a weak plasma defocusing effect and is unable to effectively reshape the spatiotemporal profile of the driving laser.Therefore, in the low-intensity case (7.5 I 0 ), the spatiotemporal distribution of the driving laser during propagation typically demonstrates Gaussian profiles in both the temporal and radial direction, as shown in Figure 4e.But when the input laser peak intensity is high, although the laser fields have Gaussian profiles in the temporal and radial direction at the entrance of the medium, the electric fields at the middle and exit of the medium are reduced compared to the fields at the entrance, as shown in Figure 4b-d.Moreover, with an increase in the propagation distance, the spatiotemporal intensity distributions of the driving laser in the cases of high intensity are greatly reshaped, as shown in Figure 4f-h.This spatiotemporal reshaping of the driving laser in high-intensity cases can be attributed to the strong plasma defocusing effect caused by the excessive free electrons at the high ionization level [36,52,53].In addition, one can clearly see the effective wavefront distortions in the high-intensity cases [54,55], as shown in Figure 4f-h.This is because the weaker intensity at a larger radial distance remains the same due to the weak plasma defocusing while the higher intensity at a smaller radial distance is reduced by the strong plasma defocusing, as shown in Figure 4f-h.Moreover, one can also see that with the wavefront distorts during propagation, the laser electric field also undergoes a blue shift in time and mostly in the trailing edge of the beam, resulting in an asymmetric laser pulse distribution in time that the laser intensity in the leading edge is much higher than that in the trailing edge, as shown in Figure 4b-d,f-h).Furthermore, as the input peak laser intensity increases, the effects of spatiotemporal reshaping and frequency blue-shift become more prominent.Therefore, the evolution of the driving laser in the high-intensity cases is different from that in the low-intensity case, which leads to the dynamic phase matching of the harmonics in the overdriven regime different from that in the conventional regime [23] (i.e., weakly ionized gas medium).
although the laser fields have Gaussian profiles in the temporal and radial direction at the entrance of the medium, the electric fields at the middle and exit of the medium are reduced compared to the fields at the entrance, as shown in Figure 4b-d.Moreover, with an increase in the propagation distance, the spatiotemporal intensity distributions of the driving laser in the cases of high intensity are greatly reshaped, as shown in Figure 4f-h.This spatiotemporal reshaping of the driving laser in high-intensity cases can be aAributed to the strong plasma defocusing effect caused by the excessive free electrons at the high ionization level [36,52,53].In addition, one can clearly see the effective wavefront distortions in the high-intensity cases [54,55], as shown in Figure 4f-h.This is because the weaker intensity at a larger radial distance remains the same due to the weak plasma defocusing while the higher intensity at a smaller radial distance is reduced by the strong plasma defocusing, as shown in Figure 4f-h.Moreover, one can also see that with the wavefront distorts during propagation, the laser electric field also undergoes a blue shift in time and mostly in the trailing edge of the beam, resulting in an asymmetric laser pulse distribution in time that the laser intensity in the leading edge is much higher than that in the trailing edge, as shown in Figure 4b-d,f-h).Furthermore, as the input peak laser intensity increases, the effects of spatiotemporal reshaping and frequency blue-shift become more prominent.Therefore, the evolution of the driving laser in the high-intensity cases is different from that in the low-intensity case, which leads to the dynamic phase matching of the harmonics in the overdriven regime different from that in the conventional regime [23] (i.e., weakly ionized gas medium).We next look at the time-frequency analysis of the resulting high harmonics exiting from the gas medium at four laser intensities and the corresponding optimal gas pressures, as shown in Figure 5.One can see that in all the cases of the input peak laser intensities, only short-trajectory emissions contribute to the harmonics.This is because the long-trajectory emissions accumulate more intrinsic phases compared to the short-trajectory emission branches, making phase matching particularly challenging.This is especially true in the overdriven regime, where the spatiotemporal reshaping of the driving laser causes a rapid change in the intrinsic phase.Furthermore, Figure 5 also showed that  (a-d) Simulated on-axis (r = 0) laser electric fields at the entrance (red line), middle (blue dotted line), and exit (green dash-dotted line) of the gas medium for the optimal gas pressure at four laser intensities.(e-h) Spatiotemporal laser intensity distributions at the exit of the gas medium for the optimal gas pressure at four laser intensities.The four laser intensities and their corresponding optimal gas pressures are 7.5 I 0 and 0.3 bar (a,e), 12 I 0 and 0.45 bar (b,f), 15 I 0 and 0.6 bar (c,g), and 20 I 0 and 0.9 bar (d,h).Laser intensities in each panel of (e-h) are normalized to their maximum.
We next look at the time-frequency analysis of the resulting high harmonics exiting from the gas medium at four laser intensities and the corresponding optimal gas pressures, as shown in Figure 5.One can see that in all the cases of the input peak laser intensities, only short-trajectory emissions contribute to the harmonics.This is because the long-trajectory emissions accumulate more intrinsic phases compared to the short-trajectory emission branches, making phase matching particularly challenging.This is especially true in the overdriven regime, where the spatiotemporal reshaping of the driving laser causes a rapid change in the intrinsic phase.Furthermore, Figure 5 also showed that as the input laser intensity increases, harmonic emissions are concentrated gradually towards the leading edge of the laser pulse.This can be explained by two reasons: on the one hand, the wavefront distortion and blue shift in time of the driving laser field in the high-intensity cases lead to the laser pulse intensity in the leading edge being much higher than that in the trailing edge (see Figure 4), which contributed to the harmonic generation.On the other hand, the ionization fraction in the laser pulse trailing edge is much higher than that in the leading edge and limits the phase matching of the harmonics.
as the input laser intensity increases, harmonic emissions are concentrated gradually towards the leading edge of the laser pulse.This can be explained by two reasons: on the one hand, the wavefront distortion and blue shift in time of the driving laser field in the high-intensity cases lead to the laser pulse intensity in the leading edge being much higher than that in the trailing edge (see Figure 4), which contributed to the harmonic generation.On the other hand, the ionization fraction in the laser pulse trailing edge is much higher than that in the leading edge and limits the phase matching of the harmonics.In order to investigate the effect of the spatiotemporal reshaping in the driving laser pulses on the harmonic generation process, we plot the spatial phase matching maps of H80 from the short-trajectory emissions inside the gas medium at t = 0 for the four laser intensities and the corresponding optimal gas pressures, as shown in Figure 6a-d.The phase mismatching in each position was calculated by Equation ( 4), which takes into account the effects of spatiotemporal reshaping and blue-shift of the driving laser in the overdriven regime, as shown in Equation (5).Additionally, the corresponding spatial evolution of harmonic intensity maps of H80 at the four laser intensities inside the gas medium was also investigated, as shown in Figure 6e-f.One can see that with the input peak laser intensity of 7.5 I0, the phase matching area is maximized, which contributed greatly to the harmonic yield buildup.The best region of phase matching is to be found on-axis downstream from the focus and the off-axis positions ranging from 30 µm to 60 µm, as shown in Figure 6a, resulting in a rapid harmonic buildup within the paraxial region in the second half of the medium, as shown in Figure 6e.When the input peak laser intensity increases to 12 I0, the phase matching area is greatly reduced due to the excessive plasma dispersion, and the best region of on-axis phase matching is located near the center of the gas medium, as shown in Figure 6b.This means that the harmonic yield buildup can only occur rapidly near the center of the medium and descends at the rear of the medium, as shown in Figure 6f.However, when the input peak laser intensity increases to 15 I0, despite experiencing an even greater plasma dispersion caused by the higher laser intensity and optimal gas pressure that makes phase matching more challenging, the phase matching area is only slightly reduced compared to the case of an input peak laser intensity of In order to investigate the effect of the spatiotemporal reshaping in the driving laser pulses on the harmonic generation process, we plot the spatial phase matching maps of H80 from the short-trajectory emissions inside the gas medium at t = 0 for the four laser intensities and the corresponding optimal gas pressures, as shown in Figure 6a-d.The phase mismatching in each position was calculated by Equation ( 4), which takes into account the effects of spatiotemporal reshaping and blue-shift of the driving laser in the overdriven regime, as shown in Equation (5).Additionally, the corresponding spatial evolution of harmonic intensity maps of H80 at the four laser intensities inside the gas medium was also investigated, as shown in Figure 6e-f.One can see that with the input peak laser intensity of 7.5 I 0 , the phase matching area is maximized, which contributed greatly to the harmonic yield buildup.The best region of phase matching is to be found on-axis downstream from the focus and the off-axis positions ranging from 30 µm to 60 µm, as shown in Figure 6a, resulting in a rapid harmonic buildup within the paraxial region in the second half of the medium, as shown in Figure 6e.When the input peak laser intensity increases to 12 I 0 , the phase matching area is greatly reduced due to the excessive plasma dispersion, and the best region of on-axis phase matching is located near the center of the gas medium, as shown in Figure 6b.This means that the harmonic yield buildup can only occur rapidly near the center of the medium and descends at the rear of the medium, as shown in Figure 6f.However, when the input peak laser intensity increases to 15 I 0 , despite experiencing an even greater plasma dispersion caused by the higher laser intensity and optimal gas pressure that makes phase matching more challenging, the phase matching area is only slightly reduced compared to the case of an input peak laser intensity of 12 I 0 , as illustrated in Figure 6c.In addition, compared with the case of 12 I 0 , the best region of on-axis phase matching in the case of 15 I 0 is to be found within a 1 mm distance near the exit of the gas medium, which contributed greatly to the rapid harmonic yield buildup at the rear of the medium, as shown in Figure 6g.This difference in phase matching in 12 I 0 and 15 I 0 can be attributed to the effects of spatiotemporal reshaping and blue-shift of the driving laser in the overdriven regime.When the input peak laser intensity further increases to 20 I 0 , a much greater plasma dispersion due to a higher laser intensity and optimal gas pressure makes the phase matching area greatly reduce so that the best region of phase matching is to be found only off-axis, which hinders the harmonic yield buildup along the propagation, as shown in Figure 6h.Additionally, one can also find that as the driving laser intensity increases, the region of maximum harmonic intensity gradually shifts away from the axis, leading to an increase in the divergence angle of the harmonics, as shown in Figure 6e-h.This can be attributed to the wavefront distortion of the driving laser field at high laser intensity due to the spatiotemporal reshaping, as shown in Figure 4c-h.
12 I0, as illustrated in Figure 6c.In addition, compared with the case of 12 I0, the best region of on-axis phase matching in the case of 15 I0 is to be found within a 1 mm distance near the exit of the gas medium, which contributed greatly to the rapid harmonic yield buildup at the rear of the medium, as shown in Figure 6g.This difference in phase matching in 12 I0 and 15 I0 can be aAributed to the effects of spatiotemporal reshaping and blue-shift of the driving laser in the overdriven regime.When the input peak laser intensity further increases to 20 I0, a much greater plasma dispersion due to a higher laser intensity and optimal gas pressure makes the phase matching area greatly reduce so that the best region of phase matching is to be found only off-axis, which hinders the harmonic yield buildup along the propagation, as shown in Figure 6h.Additionally, one can also find that as the driving laser intensity increases, the region of maximum harmonic intensity gradually shifts away from the axis, leading to an increase in the divergence angle of the harmonics, as shown in Figure 6e-h.This can be attributed to the wavefront distortion of the driving laser field at high laser intensity due to the spatiotemporal reshaping, as shown in Figure 4c-h.To further observe the influence of the spatiotemporal reshaping of the driving laser on the four components of the phase matching, we also plot the on-axis (r = 0) phase mismatching of H80 from the short-trajectory emissions at t = 0 for the four laser intensities and the corresponding optimal gas pressures, as shown in Figure 7. Figure 7a-d showed the on-axis (r = 0) phase mismatching of H80 from the total ∆kq (red line), neutral atomic dispersion ∆kn (blue line), plasma dispersion ∆kp (green line), intrinsic dipole phase ∆kd (magenta line), and Gouy phase ∆kg (black line) driven by the four laser intensities, and Figure 7e-h showed the on-axis intensity evolution of H80 driven by the four driving lasers inside the gas medium.One can see that when the input peak laser intensity is 7.5 I0, the ionization fraction is low enough that the phase mismatching caused by the plasma dispersion is nearly compensated by the neutral atom dispersion.Additionally, the positive Gouy phase mismatching ∆kg is also compensated by the negative intrinsic dipole phase mismatching ∆kd downstream from the center of the medium, as shown in Figure 7a.As a result, a favorable phase matching region is formed in the second half of the gas medium, which leads to sustained harmonic buildup along the medium, as shown in Fig- ure 7e.Therefore, the phase matching in this low-intensity laser case is actually a conventional phase matching mechanism where the ∆kg is always positive along the medium and the ∆kd is positive (negative) before (after) the focal position.However, we will see that the phase matching in the high-intensity laser cases is very different from the conventional mechanism, particularly in the overdriven regime.When the input peak laser intensity increases to 12 I0, despite the plasma dispersion ∆kp due to the excessive ionization To further observe the influence of the spatiotemporal reshaping of the driving laser on the four components of the phase matching, we also plot the on-axis (r = 0) phase mismatching of H80 from the short-trajectory emissions at t = 0 for the four laser intensities and the corresponding optimal gas pressures, as shown in Figure 7. Figure 7a-d showed the on-axis (r = 0) phase mismatching of H80 from the total ∆k q (red line), neutral atomic dispersion ∆k n (blue line), plasma dispersion ∆k p (green line), intrinsic dipole phase ∆k d (magenta line), and Gouy phase ∆k g (black line) driven by the four laser intensities, and Figure 7e-h showed the on-axis intensity evolution of H80 driven by the four driving lasers inside the gas medium.One can see that when the input peak laser intensity is 7.5 I 0 , the ionization fraction is low enough that the phase mismatching caused by the plasma dispersion is nearly compensated by the neutral atom dispersion.Additionally, the positive Gouy phase mismatching ∆k g is also compensated by the negative intrinsic dipole phase mismatching ∆k d downstream from the center of the medium, as shown in Figure 7a.As a result, a favorable phase matching region is formed in the second half of the gas medium, which leads to sustained harmonic buildup along the medium, as shown in Figure 7e.Therefore, the phase matching in this low-intensity laser case is actually a conventional phase matching mechanism where the ∆kg is always positive along the medium and the ∆k d is positive (negative) before (after) the focal position.However, we will see that the phase matching in the high-intensity laser cases is very different from the conventional mechanism, particularly in the overdriven regime.When the input peak laser intensity increases to 12 I 0 , despite the plasma dispersion ∆k p due to the excessive ionization that cannot be balanced by the neutral atom dispersion ∆k n , the large positive Gouy phase mismatching ∆k g and the negative intrinsic dipole phase mismatching ∆k d makes the ∆k q decrease gradually and cross the ∆k q = 0 line in the first half of the medium, as shown in Figure 7b.In the second half of the medium, the decrease in ∆k p makes ∆k q decrease gradually again and cross the ∆k q = 0 line for the second time.However, because of the large positive Gouy phase mismatching ∆k g along the medium, the total phase matching ∆k q increases again at the rear of the medium, which hinders the harmonic buildup in the exit plane, as shown in Figure 7f.It can be seen from Figure 7f that the harmonic is accumulated in the first half of the medium due to the phase matching of the first time and then descends gradually.In the second half of the medium, there is a re-accumulation of the harmonic due to the second phase matching and followed by a subsequent decrease again.By comparing the phase mismatching results in Figure 7a,b, one can see that the behavior of the Gouy phase mismatching and the intrinsic dipole phase mismatching in the case of 12 I 0 is very different from that in the case of 7.5 I 0 .This can be attributed to the spatiotemporal reshaping and frequency blue-shift of the driving laser in high laser intensity.Because ∆k g is dependent on the spatial phase of the driving laser and ∆k d is dependent on the intensity gradient and frequency shift of the laser, the spatiotemporal reshaping and frequency blueshift of the driving laser will make these two kinds of phase mismatching very different from that in the ordinary regime (i.e., the low-intensity case) [23].The phase matching mechanism in the overdriven regime is also called the extended phase matching [29,31].
that cannot be balanced by the neutral atom dispersion ∆kn, the large positive Gouy phase mismatching ∆kg and the negative intrinsic dipole phase mismatching ∆kd makes the ∆kq decrease gradually and cross the ∆kq = 0 line in the first half of the medium, as shown in Figure 7b.In the second half of the medium, the decrease in ∆kp makes ∆kq decrease gradually again and cross the ∆kq = 0 line for the second time.However, because of the large positive Gouy phase mismatching ∆kg along the medium, the total phase matching ∆kq increases again at the rear of the medium, which hinders the harmonic buildup in the exit plane, as shown in Figure 7f.It can be seen from Figure 7f that the harmonic is accumulated in the first half of the medium due to the phase matching of the first time and then descends gradually.In the second half of the medium, there is a re-accumulation of the harmonic due to the second phase matching and followed by a subsequent decrease again.By comparing the phase mismatching results in Figure 7a,b, one can see that the behavior of the Gouy phase mismatching and the intrinsic dipole phase mismatching in the case of 12 I0 is very different from that in the case of 7.5 I0.This can be aAributed to the spatiotemporal reshaping and frequency blue-shift of the driving laser in high laser intensity.Because ∆kg is dependent on the spatial phase of the driving laser and ∆kd is dependent on the intensity gradient and frequency shift of the laser, the spatiotemporal reshaping and frequency blue-shift of the driving laser will make these two kinds of phase mismatching very different from that in the ordinary regime (i.e., the low-intensity case) [23].The phase matching mechanism in the overdriven regime is also called the extended phase matching [29,31].With the input peak laser intensity further increasing to 15 I0, although a greater plasma dispersion makes phase matching more challenging, a strong spatiotemporal reshaping of the driving laser results in a large positive Gouy phase mismatching.This large Gouy phase mismatching can almost compensate for the plasma dispersion phase mismatching and the total phase mismatching basically follows the profile of the intrinsic dipole phase mismatching, as shown in Figure 7c.This is because the strong spatiotemporal reshaping and frequency blue-shift of the driving laser results in an enhanced intrinsic dipole phase mismatching that is more dominant than the neutral atomic dispersion phase mismatching, particularly in the fast-dropping stage of the driver pulse intensity.As a result, the balance of ∆kp, ∆kg, and ∆kd makes the total phase mismatching reduce gradually in the second half of the gas medium and realizes a good phase matching near With the input peak laser intensity further increasing to 15 I 0 , although a greater plasma dispersion makes phase matching more challenging, a strong spatiotemporal reshaping of the driving laser results in a large positive Gouy phase mismatching.This large Gouy phase mismatching can almost compensate for the plasma dispersion phase mismatching and the total phase mismatching basically follows the profile of the intrinsic dipole phase mismatching, as shown in Figure 7c.This is because the strong spatiotemporal reshaping and frequency blue-shift of the driving laser results in an enhanced intrinsic dipole phase mismatching that is more dominant than the neutral atomic dispersion phase mismatching, particularly in the fast-dropping stage of the driver pulse intensity.As a result, the balance of ∆k p , ∆k g , and ∆k d makes the total phase mismatching reduce gradually in the second half of the gas medium and realizes a good phase matching near the rear of the medium, which contributed greatly to the harmonic buildup along the medium, as shown in Figure 7g.When the input peak laser intensity further increases to 20 I 0 , even though the significantly greater plasma dispersion can still be counterbalanced by the intensified Gouy phase mismatching caused by the stronger spatiotemporal reshaping of the driving laser, the intrinsic dipole phase mismatching is also greatly enhanced.Consequently, this leads to poor phase matching throughout the gas medium, thereby impeding the harmonic buildup, as depicted in Figure 7d,h.Therefore, we can conclude that the spatiotemporal reshaping of the driving laser in the overdriven regime can significantly influence the phase matching of the harmonics and the harmonic buildup along the gas medium.
It is worth noting that, by comparing the harmonic buildup process in the four laser intensities, as shown in Figure 7e-h, one can find that despite the phase mismatching being much smaller in the case of 7.5 I 0 compared to 12 I 0 and 15 I 0 , the harmonic buildup in the exit plane of the gas medium is still lower in the low-intensity case than that in the high-intensity cases.This is attributed to the high induced atomic dipole moment amplitude of the single-atom response in the high laser intensity and the high atomic number density in the high gas pressure (see Equation ( 5)).One can also see that the harmonic intensity at the exit of the medium in the case of 12 I 0 is lower than that in the case of 15 I 0 , which can be attributed to the phase matching of the harmonics.The harmonic in the case of 15 I 0 has a good phase matching at the rear of the medium, resulting in a continuous harmonic buildup along the medium, as shown in Figure 7c,g.However, the harmonic in the case of 12 I 0 has poor phase matching due to the large Gouy phase mismatching at the rear of the medium, resulting in a decrease in the harmonic at the rear of the medium.This discrepancy in phase matching between the two laser intensities can also be used to explain the difference between the results in Figures 2d and 3c.Therefore, the phase matching of the harmonic can also greatly affect the harmonic yield buildup.We can also conclude that both the induced atomic dipole moment amplitude of the singleatom response and the phase matching of the macroscopic medium should be carefully considered to maximize the harmonic generation.In addition, by comparing the phase matching in the overdriven regime as shown in Figure 7b-d, one can also see that only a suitable spatiotemporal reshaping of the driving laser results in favorable phase matching and harmonic buildup.In the case of 12 I 0 , a weak spatiotemporal reshaping that makes the Gouy phase mismatched just slightly changed compared to the low-intensity case, resulting in poor phase matching at the rear of the medium, as shown in Figure 7b.However, in the case of 20 I 0 , although a strong spatiotemporal reshaping makes the ∆k g change considerably compared to the low-intensity case, the intrinsic dipole phase mismatching is also greatly enhanced, which also results in poor phase matching along the medium.Therefore, it is the suitable spatiotemporal reshaping of the driving laser in the case of 15 I 0 that makes the ∆k g and the ∆k d a relatively large value in the first half of the medium to compensate for the plasma dispersion and a small value at the rear of the medium to realize favorable phase matching.

Conclusions
In this work, we demonstrated an approach to enhance the XUV harmonics in the plateau region by increasing the intensity of the driving laser while keeping the laser energy constant.We theoretically investigated the effects of different driving laser intensities at different gas pressures on the harmonics generation.The simulation results showed that by increasing the laser intensity to the overdriven regime, the average XUV harmonics yield in the plateau region is approximately twice as high as that obtained optimally in the conventional loose focusing geometry scheme by utilizing a relatively low-intensity driving laser with the same laser energy.The analysis of the HHG in the macroscopic medium using a 1D model, along with a further quantitative analysis of the phase matching, has revealed that the observed enhancement in harmonics can be attributed to two key factors.On the one hand, it is the amplification of the induced atomic dipole moment amplitude of the single-atom response in the high-intensity driving laser that results in the improvement of the harmonic generation efficiency in the microscopic, which contributes to the enhancement of the harmonic.On the other hand, the favorable phase matching along the gas medium in the overdriven regime can also contribute to the harmonic buildup in the macroscopic.Furthermore, by investigating the effects of the spatiotemporal reshaping of the driving laser on the harmonic's phase mismatching in the overdriven regime, the results revealed that the spatiotemporal reshaping can greatly change the phase mismatching of the harmonics, particularly the Gouy phase mismatching and the intrinsic dipole phase mismatching.In addition, the results also showed that only a suitable spatiotemporal reshaping of the driving laser can contribute to forming a favorable phase matching and

Figure 1 .
Figure1.Simulated HHG spectra generated for different gas pressures at four laser intensities, where I0 = 1 × 10 14 W/cm 2 .The harmonic yields are integrated over the radial direction in the near field.

Figure 1 .
Figure1.Simulated HHG spectra generated for different gas pressures at four laser intensities, where I 0 = 1 × 10 14 W/cm 2 .The harmonic yields are integrated over the radial direction in the near field.

Figure 2 .Figure 2 .
Figure 2. (a-c) Optimal HHG spectra (a-b) and the statistical harmonic yields of different harmonic orders (c) driven by the four laser intensities at the corresponding optimal gas pressures, where I0 = Figure 2. (a-c) Optimal HHG spectra (a,b) and the statistical harmonic yields of different harmonic orders (c) driven by the four laser intensities at the corresponding optimal gas pressures, where I 0 = 1 × 10 14 W/cm 2 .The harmonic yields are integrated over the radial direction in the near field.The arrows in (b) represent the selected harmonic orders to analyze furtherly in this paper, which are uniformly referred to as H80 for convenience.(d) The harmonic yield of selected H80 versus gas pressure for different laser intensities.

Figure 3 .
Figure 3. (a,b) Simulated on-axis (r = 0) laser intensity evolution (a) and probability evolution (b) along the propagation axis for the optimal gas pressure at four laser intensities: 7.5 I0 and 0.3 bar (black line), 12 I0 and 0.45 bar (blue line), 15 I0 and 0.6 bar (red line), and 20 I0 and 0.9 bar (green line).z0 represents the Rayleigh length for each laser beam waist.(c) The calculated harmonic yield of H80 for different gas pressures at the four laser intensities by numerically integrating Equation (1).

Figure 3 .
Figure 3. (a,b) Simulated on-axis (r = 0) laser intensity evolution (a) and probability evolution (b) along the propagation axis for the optimal gas pressure at four laser intensities: 7.5 I 0 and 0.3 bar (black line), 12 I 0 and 0.45 bar (blue line), 15 I 0 and 0.6 bar (red line), and 20 I 0 and 0.9 bar (green line).z 0 represents the Rayleigh length for each laser beam waist.(c) The calculated harmonic yield of H80 for different gas pressures at the four laser intensities by numerically integrating Equation (1).

Figure 4 .
Figure 4.(a−d) Simulated on-axis (r = 0) laser electric fields at the entrance (red line), middle (blue doAed line), and exit (green dash-doAed line) of the gas medium for the optimal gas pressure at four laser intensities.(e−h) Spatiotemporal laser intensity distributions at the exit of the gas medium for the optimal gas pressure at four laser intensities.The four laser intensities and their corresponding optimal gas pressures are 7.5 I0 and 0.3 bar (a,e), 12 I0 and 0.45 bar (b,f), 15 I0 and 0.6 bar (c,g), and 20 I0 and 0.9 bar (d,h).Laser intensities in each panel of (e−h) are normalized to their maximum.

Figure 4 .
Figure 4.(a-d) Simulated on-axis (r = 0) laser electric fields at the entrance (red line), middle (blue dotted line), and exit (green dash-dotted line) of the gas medium for the optimal gas pressure at four laser intensities.(e-h) Spatiotemporal laser intensity distributions at the exit of the gas medium for the optimal gas pressure at four laser intensities.The four laser intensities and their corresponding optimal gas pressures are 7.5 I 0 and 0.3 bar (a,e), 12 I 0 and 0.45 bar (b,f), 15 I 0 and 0.6 bar (c,g), and 20 I 0 and 0.9 bar (d,h).Laser intensities in each panel of (e-h) are normalized to their maximum.

Figure 5 .
Figure 5. Time-frequency analysis of harmonic emissions (normalized to the maximum in each panel) at the exit plane of the gas medium for the four laser intensities and the corresponding gas pressures.

Figure 5 .
Figure 5. Time-frequency analysis of harmonic emissions (normalized to the maximum in each panel) at the exit plane of the gas medium for the four laser intensities and the corresponding gas pressures.

Figure 6 .
Figure 6.(a-d) Spatial phase matching map of H80 from the short-trajectory emissions inside the gas medium at t = 0 for the four laser intensities and the corresponding optimal gas pressures.(e-h) The corresponding spatial evolution of harmonic intensity maps (normalized to the maximum in each panel) of H80 at the four laser intensities inside the gas medium.z0 represents the Rayleigh length for each kind of laser beam waist.

Figure 6 .
Figure 6.(a-d) Spatial phase matching map of H80 from the short-trajectory emissions inside the gas medium at t = 0 for the four laser intensities and the corresponding optimal gas pressures.(e-h) The corresponding spatial evolution of harmonic intensity maps (normalized to the maximum in each panel) of H80 at the four laser intensities inside the gas medium.z 0 represents the Rayleigh length for each kind of laser beam waist.

Figure 7 .
Figure 7. (a−d) Simulated on-axis (r = 0) phase mismatching of H80 from the total ∆kq (red line), neutral atomic dispersion ∆kn (blue line), plasma dispersion ∆kp (green line), intrinsic dipole phase ∆kd (magenta line), and Gouy phase ∆kg (black line) inside the medium driven by the four laser intensities and their optimal gas pressures.(e−h) On-axis intensity evolution of H80 inside the gas medium driven by the four driving lasers and their optimal pressures.The four laser intensities and their corresponding optimal gas pressures are 7.5 I0 and 0.3 bar (a,e), 12 I0 and 0.45 bar (b,f), 15 I0 and 0.6 bar (c,g), and 20 I0 and 0.9 bar (d,h).

Figure 7 .
Figure 7.(a-d) Simulated on-axis (r = 0) phase mismatching of H80 from the total ∆k q (red line), neutral atomic dispersion ∆k n (blue line), plasma dispersion ∆k p (green line), intrinsic dipole phase ∆k d (magenta line), and Gouy phase ∆k g (black line) inside the medium driven by the four laser intensities and their optimal gas pressures.(e-h) On-axis intensity evolution of H80 inside the gas medium driven by the four driving lasers and their optimal pressures.The four laser intensities and their corresponding optimal gas pressures are 7.5 I 0 and 0.3 bar (a,e), 12 I 0 and 0.45 bar (b,f), 15 I 0 and 0.6 bar (c,g), and 20 I 0 and 0.9 bar (d,h).