Mid-infrared laser-driven broadband water-window supercontinuum generation from pre-excited medium

We theoretically investigate the broadband water-window supercontinuum generation from pre-excited medium with a mid-infrared pulse. We find that the wavelength scaling of the harmonic yield from near-visible (0.8 μm) to mid-infrared (1.8 μm) in single-atom level is λ−2.7. Using an intense phase-stabilized few-cycle 1.6 μm laser pulse, a broadband water window supercontinuum with bandwidth of approximately 140 eV is obtained. We also investigate the macroscopic effects and find that large initial population of the excited state leads to the high-density of free electrons, which shift the carrier-envelop phase of the driving pulse and further diminish the water-window supercontinuum generation. The highly-ionized medium also results in poor temporal and spatial properties of the attosecond pulse. Instead, small initial population of the excited state can produce well phase-matched xuv supercontinuum in water-window region and an 100-as pulse with central wavelength of 2.8 nm and pulse energy of 0.15 nJ can be filtered out. © 2011 Optical Society of America OCIS codes: (190.7110) Ultrafast nonlinear optics; (190.4160) Multiharmonic generation; (300.6560) Spectroscope, x-ray. References and links 1. F. Krausz, “Attosecond physics.” Rev. Mod. Phys. 81, 163-234 (2009). 2. E. Goulielmakis, M. Schultze, M. Hofstetter, V. S. Yakovlev, J. Gagnon, M. Uiberacker, A. L. Aquila, E. M. Gullikson, D. T. Attwood, R. Kienberger, F. Krausz, U. Kleineberg, “Single-Cycle Nonlinear Optics,” Science, 320, 1614-1617 (2008). 3. E. takahashi, P. Lan, O. D. Mucke, Y.Nabekawa, and K. Midorikawa, “Infrared two-color multicycle laser field synthesisi for generating an intense attosecond pulse,” Phys. Rev. Lett. 104, 233901 (2010). 4. Y. Zheng, Z. Zeng, P. Zou, L. Zhang, X. Li, P. Liu, R. Li, and Z. Xu, “Dynamic chirp control and pulse compression for attosecond high-order harmonic emission,” Phys. Rev. Lett. 103, 043904 (2009). 5. P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. 71, 1994-1997 (1993). 6. M. Lewenstein, Ph. Balcou, M. Yu. Ivanov, Anne L’Huillier, and P. B. Corkum, “Theory of high harmonic generation by low-frequency laser fields,” Phys. Rev. A 49, 2117 (1994). 7. Y. Zheng, Z. Zeng, X. Li, X. Chen, P. Liu, H. Xiong, H. Lu, S. Zhao, P. Wei, L. Zhang, X. Wang, J. Liu, Y. Cheng, R. Li, and Z. Xu, “Enhancement and broadening of extreme-ultraviolet supercontinuum in a relative phase controlled two-color laser field,” Opt. Lett. 33, 234-236 (2008). #156105 $15.00 USD Received 6 Oct 2011; revised 27 Oct 2011; accepted 3 Nov 2011; published 14 Nov 2011 (C) 2011 OSA 21 November 2011 / Vol. 19, No. 24 / OPTICS EXPRESS 24376 8. W. Becker, A. Lohr, M. Kleber, and M. Lewenstein, “A unified theory of high-harmonic generation: Application to polarization properties of the harmonics,” Phys. Rev. A 56, 645 (1997). 9. Z. Chang, “Single attosecond pulse and xuv supercontinuum in the high-order harmonic plateau,” Phys. Rev. A 70, 043802 (2004). 10. H. Mashiko, S. Gilbertson, C. Li, S. D. Khan, M. M. Shakya, E. Moon, and Z. Chang, “Double optical gating of high-order harmonic generation with carrier-envelope phase stabilized laser,” Phys. Rev. Lett. 100, 103906 (2008). 11. E. Takahashi, T. Kanai, K. Ishikawa, Y. Nabekawa, and K. Midorikawa, “Coherent Water Window X Ray by Phase-Matched High-Order Harmonic Generation in Neutral Media,” Phys. Rev. Lett. 101, 253901 (2008). 12. W. Hong, P. Lu, Q. Li, and Q. Zhang, “Broadband water window supercontinuum generation with a tailored mid-IR pulse in neutral media,” Opt. Lett. 34, 2102-2104 (2009). 13. Q. Zhang, P. Lu, P. Lan, W. Hong, and Z. Yang, “Multi-cycle laser-driven broadband supercontinuum with a modulated polarization gating, ” Opt. Express 16, 9795-9803 (2008). 14. W. Hong, P. Wei, Q. Zhang, S. Wang, and P. Lu, “Mid-infrared modulated polarization gating for ultra-broadband supercontinuum generation,” Opt. Express 18, 11308-11315 (2010). 15. J. Tate, T. Auguste, H. Muller, P. Salières, P. Agostini, and L. F. DiMauro, “Scaling of Wave-Packet Dynamics in an Intense Midinfrared Field,” Phys. Rev. Lett. 98, 013901 (2007). 16. A. Sanpera, J. B. Watson, M. Lewenstein, and K. Burnett, “Harmonic-generation control,” Phys. Rev. A 54, 4320 (1996). 17. E. Priori, G. Cerullo, M. Nisoli, S. Stagira, S. De Silvestri, P. Villoresi, L. Poletto, P. Ceccherini, C. Altucci, R. Bruzzese, and C. de Lisio, “Nonadiabatic three-dimentional model of high-order harmonic generation in the few-optical cycle regime,” Phys. Rev. A 61, 063801 (2000). 18. M. V. Frolov, N. L. Manakov, and Anthony F. Starace, “Wavelength Scaling of High-Harmonic Yield: Threshold Phenomena and Bound State Symmetry Dependence,” Phys. Rev. Lett. 100, 173001 (2008). 19. K. Schiessl, K. L. Ishikawa, E. Persson, and J. Burgdorfer, “Quantum Path Interference in theWavelength Dependence of High-Harmonic Generation,” Phys. Rev. Lett. 99, 253903 (2007). 20. P. Lan, P. Lu, P. Lan, F. Li, Y. Li, and Z. Yang, “Carrier-envelope phase measurement from half-cycle high harmonics,” Opt. Express 16, 5868-5873 (2008). 21. W. Hong, Q. Zhang, X. Zhu, and P. Lu, “Intense isolated attosecond pulse generation in pre-excited medium,” Opt. Express 19, 4728-4739 (2011).


Introduction
High-order harmonic generation (HHG) is one of the major topics in the interaction of intense laser fields with atoms.It has provided an attractive way towards attosecond extreme ultraviolet (XUV) light source, making a breakthrough in attosecond science [1] and nonlinear optics in the extreme ultraviolet region [2].One of the major applications of HHG is to produce isolated attosecond pulses.In the past few years, how to broaden the bandwidth of the supercontinuum and shorten the pulse duration of the isolated attosecond pulse have attracted great interests, because the generation of shorter attosecond pulses enables a higher time resolution for ultrafast measurements [3,4].Recently, Cavalieri and co-workers [2] have broke through the 100-asbarriers.In their experiment, a sub-4-fs near one optical cycle driving pulse has been employed to generated a 40 eV supercontinuum and an 80-as pulse has been filtered out, with the pulse energy of 0.5 pJ.
The HHG process can be well explained by the quasi-classical three-step model [5,6].First, the electron is ionized by tunnelling through the potential barrier, then it is driven by the laser field treated as a free electron, finally, it may return to the atomic core and emit a harmonic photon in the transition back to the ground state.The maximum kinetic energy acquired by a free electron from the laser field when it recombines with the atomic core is 3.17U p (where U p = I/4ω 2 is the ponderomotive energy of the laser field with frequency ω and intensity I), which coincides with the cutoff energy of the emitted harmonic photons.Manipulating the three steps in HHG process (ionization, acceleration and recombination) offers a possibility to control the electron trajectories for the generation of a broadband supercontinuum.For instance, two-color schemes have been proposed to control electron dynamics to generate broadband attosecond pulses [7].An alternative way to broaden the bandwidth of the supercontinuum is applying polarization gating technique (PG) [8,9] and double-optical-gating technique (DOG) [10].It has been demonstrated that the electron only returns efficiently in a short time period when PG is adopted and then an isolated attosecond pulse can be produced.Recently, a very broadband XUV continuum, which supports a 16-as isolated pulse generation, has been produced with double-optical-gating technique (DOG).
The well-known I p + 3.17U p law for the position of the HHG cutoff (where I p is the atomic binding energy) shows that using mid-infrared lasers allows the generation of higher order harmonics [11] and offers the possibility for the generation of broader supercontinuum in waterwindow [the spectra range between the K-absorption edges of carbon (284 eV) and oxygen (543 eV)] region.Hong et al. have found that the two-color schemes in the mid-IR regime are much more effective in electron trajectory control than those in the near-IR regime [12].The DOG and PG techniques have also been applied in the mid-IR region to produce broadband supercontinuum [13,14].However, a theoretical investigation about the wavelength scaling from near-visible (0.8 μm) to mid-infrared (1.8 μm) of high order harmonics shows that the integration of the harmonic yield follows a λ − (5−6) scaling at constant intensity [15].The spreading of the electron wave packet leads to a decrease of the HHG yield in mid-infrared regime [15], which leads to low harmonic efficiency thus low energy of attosecond pulse.
An effective way to increase the harmonic efficiency of harmonics driven by mid-infrared pulses is to populate part of the electron wave packet at the excited state.It has been demonstrated that the ionization rate can be significantly enhanced when atomic systems are prepared in a coherent superposition of ground and excited states [16], since part of the electron wave packet on the excited state can be much more easily projected to the continuum due to its relatively small ionization potential.Thereby the harmonic efficiency can be significantly increased compared with HHG generated from the ground state alone.In this paper, we theoretical investigate the generation of a broadband water-window supercontinuum in pre-excited medium.A better wavelength scaling of the harmonic yield is demonstrated in infrared region.Using a few-cycle phase-stabilized 1.6 μm driving pulse when the target is prepared in a coherent superposition of two electron states, a supercontinuum in water-window region with bandwidth of approximately 140 eV is straightforwardly obtained.The macroscopic effects are also investigated associated with the initial population of the excited state.The results show that in large initial population case the harmonic spectra is deeply modulated and the temporal and spatial properties of the generated pulses are poor.Instead, small initial population of the excited state can produce well phase-matched water-window xuv supercontinuum and a 100 as attosecond pulse can be filtered out with pulse energy of 0.15 nJ and central wavelength of 2.8 nm.

Single-atom response
At single-atom level, we perform a quantum simulation by solving the time-dependent Schrödinger equation.In the single-active-electron (SAE) and dipole approximations, the Hamiltonian is (in atomic units), where V (x) is the atomic Coulomb potential and E(t) is the electric field.We choose the soft-core model within SAE approximation where a and b are parameters describing different ions.For a He + ion, we choose a = 2 and b = 0.5 to retrieve the ionization energy of the ground state (1s), i.e. 54.4 eV, and the first excited state (2s) ionization energy is 13.6 eV.In our simulation, a 8-fs 1600-nm few-cycle laser pulse with intensity of 6 × 10 14 W/cm 2 is used.The electric field of the driving pulse can be written as where E 0 , φ and ω are the amplitude, carrier envelope phase (CEP) and laser frequency of the driving pulse, respectively.The initial state is prepared as a coherent superposition of the 1s and 2s states where α and β are the amplitudes of the ground and first excited states and ϕ is the initial relative phase between two states.ϕ is set to be zero because it has no influence on the results.The initial population of the excited state is described as p = β 2 .
The instantaneous wave function can be expressed as where ψ g (r,t) and ψ e (r,t) are the instantaneous wave functions originating from the ground state and the excited state, respectively.In order to obtain the harmonic spectrum, we should calculate the time-dependent dipole acceleration The harmonic spectrum is then obtained by Fourier transforming of the time-dependent dipole acceleration In order to describe the features of harmonic spectra, we split up different contributions to the dipole acceleration for the coherent superposition as [16] If we assume that the ground and excited state are not coupled to any other bound states during the pulse, we can rewrite the time-dependent wave functions as ψ g (r,t) = αe −iω g t |g >, ( where |c > is the continuum state.γ e and γ c , are the time-dependent amplitudes of the excited and continuum states.Here we assume the ground state is not depleted, so the term < ψ g (r,t)| x|ψ g (r,t) > doesn't contribute to the harmonic generation.Then we can rewrite the relevant contributions to the acceleration as The first term is the dipole acceleration one would obtain starting from the excited state.The second term can be thought of as an interference term, i.e. the dipole transition between the continuum and the ground state, where the excited state is responsible for the ionization.Since the excited state is almost depleted at the intensity of 6 × 10 14 W/cm 2 (as shown in Fig. 1), the ionization rate can be highly increased and so are the conversion efficiency.As a result, this part of dipole acceleration would play an irreplaceable role on the harmonic spectra.

Macroscopic propagation response
The macroscopic response can be described by the co-propagation of the laser and harmonics beams which is simulated by numerically solving Maxwell wave equations for the fields of the laser pulse E f and the harmonics E h in cylindrical coordinates separately.
where ω p is the plasma frequency and is given by The nonlinear polarization of gas is , where n 0 is the gas density and n e is the free-electron density.This model takes into account both temporal plasma-induced phase modulation and spatial plasma-lensing effects on the laser pulse.The procedures of solution of above equations have been clearly described in [17].Firstly , we investigate the population variations of ground and excited states.Figure 1 shows the normalized populations of the excited state (red dashed line) and the ground state (blue dashed line) as functions of time when the initial state is a coherent superposition of the ground and excited states with equally weighted initial populations.The electric field (green line) of the driving pulse is also shown in this figure.The excited state ionizes completely within the first 1.5 optical cycles while the ground state is hardly ionized.Since harmonic generation originates from the coherent dipole transition between the continuum and the bound states, only those states that remain populated during the pulse will contribute to the harmonic generation.The interference term < ψ g (r,t)| x|ψ e (r,t) >, where the excited state is responsible for the ionization and the dipole transition is between the continuum and the ground state, plays a more important role than other two terms in HHG process.Therefore, conversion efficiency can be highly increased in the coherent superposition case.Fig. 2. High-order harmonic spectra generated from He + at the intensity of 6×10 14 W/cm 2 for cases of ground state (blue dashed line) and superposition state (red thick line).In superposition state case, the initial populations of ground and excited state are equal.The inset is the harmonic spectra for both cases when the wavelength of the driving pulse is 0.8 μm.

Result and discussion
To conform the harmonic efficiency increasing, we calculate the harmonic spectra at singleatom level.Figure 2 presents the harmonic spectra of the He + ion in coherent superposition state case (red thick line) and ground state case (blue dotted line), respectively.The CEP of the driving pulse is set to be 0.6 π.The harmonic efficiency of the superposition state case is about 6 orders of magnitude higher than that of the ground case with the same cutoff harmonic order of approximately 620th.A water-window supercontinuum with the bandwidth of approximately 140 eV is straightforwardly obtained.The inset is the harmonic spectra when the wavelength of the driving pulse is 0.8 μm (near-visible) which is shown for comparison.We can see that the bandwidth of the supercontinuum is further broadened and more energetic harmonic photons in water-window region are generated using a 1.6 μm driving pulse than that of a 0.8 μm driving pulse with equal intensity.
A deeper insight is obtained by investigating the temporal and spectral structures of HHG using wavelet time-frequency analysis method.Figure 3 presents the time profiles of the harmonics when the initial state is superposition state.The maximum harmonic order of the highest peak is about 620th, agreeing well with harmonic spectra cutoff harmonic order.We can also see that multi-recombinations of the electron wave packet are remarkable.This is due to the high ionization rate of the excited state.
To further demonstrate the advantage of using coherent superposition state, we investigate the wavelength scaling of HHG yield from near-visible (0.8 μm) to mid-infrared (1.8 μm).The wavelength scaling is evaluated by integrating the HHG yield (defined as radiated energy per unit time) from 50 to 80 eV, shown in Fig. 4. As pointed out by J.Tate and co-workers [15], the harmonic yield follows a λ − (5−6) scaling at constant intensity when the initial state is only ground state.However, in our simulation, when using coherent superposition state, the harmonic yield scales as λ −2.7 , which falls more slowly as the increasing of the driving laser wavelength.Although the physics underlying the scaling law is not yet fully understood: the scaling law of the HHG yield (λ − (5−6) ) at constant intensity has discrepancy on the theory predicted λ −3 and high-order returns of the electron wave packet have not been paid enough attention to in previous analysis [18,19].It is expected that the superposition state case may have potential to help us understand this owing to its noteworthy multi-recombinations of the electron wave packet.
As for few-cycle driving pulses, the CEP plays an important role [20].Figure 5(a) shows the population variations of the excited state as functions of CEP.Since the ionization rate of excited state has a significant influence on HHG process in superposition state case, we can predict that the different variations of the population of the excited state will change the struc- ture of the harmonic spectra.As shown in Fig. 5(b), the harmonic efficiency and bandwidth of the generated xuv supercontinuum sensitively depend on the CEP of driving pulse.This is quite different from the ground state case, of which the CEP only influences the modulation of the supercontinuum and slightly changes its bandwidth and the harmonic cutoff.A qualitative explanation to the results is that the ionization time of the electron trajectory coincides with the time when the population of the excited state decreases fastest at some CEPs, which forms a broadband supercontinuum with high conversion efficiency While at other CEPs this condition can't be satisfied and the bandwidth are relatively narrow and the harmonic efficiencies are low.Next, we investigate the macroscopic effects.As described above, high ionization rate of the excited state is responsible for the high harmonic efficiency in superposition state case.Since the population of the excited state is almost totally depleted during the laser pulse while the ground state is hardly ionized, more populations of the excited state mean higher conversion efficiency of the harmonics at single-atom level.However, considering macroscopic effects, more populations of the excited state bring high density of free electrons into the medium, which may deteriorate the phase-matching condition and change the spatiotemporal propagation dynamics of the laser field including the spatial defocusing of the beam and temporal distortion of the electric field [21].To demonstrate this issue, we first perform three-dimensional propagation simulations for driving and harmonic fields in the gas medium.The parameters are as below: the beam waist of the focused laser beam is 50 μm and the gas target is 500 μm long with the density of 2.6 × 10 18 /cm 2 .The gas jet is placed 1 mm after the laser focus.Figure 6 shows the macroscopic harmonic spectra generated from He + in coherent superposition state with different initial excited state populations of p = 0.05 (red line), p = 0.1 (blue line) and p = 0.3 (green line).We can see that the modulation in the water-widow supercontinuum region is weaken in the case of p = 0.05 after propagation, which implies the phase-matching conditions of one electron trajectory is satisfied.While in the cases of p = 0.1 and p = 0.3, the supercontinua are more deeply modulated than that at single-atom level.Moreover, in these two cases, the harmonic cutoffs are changed in comparison with the single-atom response, which may be due to the distortion of the electric field after propagation.Figure 7 shows the electric field of the driving pulses after propagation with the excited state populations of p = 0.05 (red line), p = 0.1 (blue line) and p = 0.3 (green line), respectively.The dashed line shown in this figure is the initial driving field at the entrance of the gas jet.It is found that the electric field distortion is unconspicuous for the case of p = 0.05.On the contrary, it is considerable for the cases of p = 0.1 and p = 0.3.This distortion changes the temporal profiles of the electric field, which can be regarded as CEP shifts.As shown in Fig. 5, the efficiencies and cutoffs of the supercontinua strongly depend on the CEP, so the CEP shifts caused by the distortion of the electric field after propagation change the macroscopic xuv supercontinua.
We further investigate the temporal profiles of the supercontinuum after propagation by means of applying a square window with width of 50 harmonics to the supercontinuum from 550th to 600th harmonic orders for each case.The results are shown in Figs.8(a)-(c).The bandwidth is 39 eV, approximately 100 as to the Fourier-transform limit.For the case of p = 0.05, a pure attosecond pulse is filtered out from the double-peak structure, which indicates that the harmonics from short quantum path are well phase-matched while that from the other electron trajectories are suppressed.On the contrary, the satellite pulses become noticeable for the case of p = 0.1 and multi-peaks attosecond pulses with equal weighted intensities are appeared for the case of p = 0.3 because of the deep modulation of the supercontinua.The spatial distribution is also an important characteristic of the generated attosecond pulse.Figures 8(d)-(f) show the spatiotemporal profiles of the attosecond pulse generation in the cases of p = 0.05, p = 0.1 and p = 0.3.This characteristic is related to the high-density of free electrons.In the spatial domain, the refractive index introduced by the free electrons causes the laser phase front to advance faster on axis,which leads to defousing of the beam.The higher density of free electrons in the cases of p = 0.1 and p = 0.3 make the spatiotemporal distribution much more complicated than that for the case of p = 0.05.We further calculate the emitted power of the attosecond pulses and find it is estimated to be approximately 0.15 nJ in the case of p = 0.05.

Conclusion
In summary, we theoretically investigate the generation of a broadband water-window supercontinuum with high efficiency.When the medium (He + ) is prepared in a coherent superposition of ground and excited states, a better wavelength scaling of the harmonic yield is demonstrated: the harmonic yield scales as λ −2.7 at constant intensity, which falls more slowly as the increase of the driving laser wavelength yet the physical mechanism of this scaling law should be further discussed.It should be emphasize that the λ −2.7 scaling keeps unchanged when the initial population of the excited is higher than 2%.Using a phase-stabilized few-cycle 1.6 μm laser pulse to irradiate the prepared excited medium, the enhancement of the harmonic efficiency is achieved due to the high ionization rate of the excited state and a water-window supercontinuum of about 140 eV is obtained straightforwardly.We also investigate the macroscopic effects.The results show that large initial population of the excited state causes the medium to be high-ionized and introduces high density of the free electrons, which lead to the distortion of the electric field and make the supercontinuum deeply modulated.The highlyionized medium also results in poor temporal and spatial properties of the attosecond pulse.On the other hand, small initial population of the excited state can produce well phase-matched xuv supercontinuum and an isolated 100-as attosecond pulse with the pulse energy of 0.15 nJ and central wavelength of 2.8 nm can be filtered out.We also calculated different initial populations of the excited state and find that 2% − 7% initial population of the excited state can produce an isolated attosecond pulse.When it increases to 8% − 15%, the supercontinuum is modulated and a satellite pulse is emerged.When the initial population of the excited state becomes more, for instance, 30%, the supercontinuum becomes deeply modulated and multi-peaks attosecond pulses with equal weighted intensities are appeared.

Fig. 1 .
Fig.1.Normalized populations of excited state (red dashed line) and ground state (blue dashed line) as functions of time when the initial state is a coherent superposition of ground and excited states with equally weighted populations.The electric field (green line) of the driving pulse with intensity of 6 × 10 14 W/cm 2 and CEP φ = 0.6π.

Fig. 3 .Fig. 4 .
Fig. 3. Time-frequency distributions of the harmonics when the initial state is coherent superposition state.

Fig. 5 .
Fig. 5. (a) The population variations of the excited state as functions of CEP and (b) the harmonic spectra as functions of CEP.

Fig. 6 .Fig. 7 .
Fig.6.The macroscopic harmonic spectra generated from He + in coherent superposition state with different initial excited state population of p = 0.05 (red line), p = 0.1 (blue line) and p = 0.3 (green line).

Fig. 8 .
Fig. 8. Normalized temporal profiles of the attosecond pulses with the excited-state populations of (a) p = 0.05, (b) p = 0.1 and (c) p = 0.3.Spatiotemporal profiles of the attosecond pulses at the end of the medium with the excited-state populations of (d) p = 0.05, (e)p = 0.1 and (f) p = 0.3.