Twin-Capture Rydberg State Excitation Enhanced with Few-Cycle Laser Pulses
-
Abstract
Quantum excitation is usually regarded as a transient process occurring instantaneously, leaving the underlying physics shrouded in mystery. Recent research shows that Rydberg-state excitation with ultrashort laser pulses can be investigated and manipulated with state-of-the-art few-cycle pulses. We theoretically find that the efficiency of Rydberg state excitation can be enhanced with a short laser pulse and modulated by varying the laser intensities. We also uncover new facets of the excitation dynamics, including the launching of an electron wave packet through strong-field ionization, the re-entry of the electron into the atomic potential and the crucial step where the electron makes a U-turn, resulting in twin captures into Rydberg orbitals. By tuning the laser intensity, we show that the excitation of the Rydberg state can be coherently controlled on a sub-optical-cycle timescale. Our work paves the way toward ultrafast control and coherent manipulation of Rydberg states, thus benefiting Rydberg-state-based quantum technology. -
Due to their significant dipole–dipole interactions and enduring stability, Rydberg atoms have emerged as a promising avenue in the realm of detecting microwave electric fields, high-speed terahertz (THz) imaging and quantum information applications.[1–4] Conventionally, population of Rydberg states involves a sequential single-photon resonant transition through a predetermined pathway using a series of narrow-band laser fields.[5] In contrast, coherent Rydberg wave packets, which consist of a superposition of multiple excited states, can be efficiently generated and persist under intense laser pulses.[6] Achieving coherent excitation and manipulation of Rydberg states depends on a thorough understanding of the mechanisms involved in strong-field excitation, facilitated by advanced ultrafast laser techniques.
For atoms exposed to strong laser fields, the excitation of Rydberg wave packets has long been predicted.[7] However, in the tunneling regime, valence electrons were previously believed to undergo ionization and have little opportunity to settle on the excited states. Tunneling ionization can be well-described theoretically using Keldysh theory[8] and can be easily investigated experimentally through charged particles or coherent emissions, such as terahertz waves[9] and high-order harmonics,[10] resulting from the rescattering of liberated electrons by the ionic core. It was later discovered that Rydberg states played roles in both the phase-matching[11] and emitting processes of high-order harmonic generation (HHG).[12] Additionally, Rydberg states were found to lead to two kinds of independent emission of HHG: extreme-ultraviolet free induced decay (XFID) and excited state HHG.[13] Notably, bright XFID emission was observed[14,15] in excited atoms and the XFID intensity strongly depended on the ellipticity and carrier-envelope phase of the laser pulses.[16]
With increasing experimental observations, more insight has been gained into the buildup dynamics of Rydberg states under strong laser fields. A stabilization mechanism through destructive interference transitions or the formation of Kramers–Henneberger states[17,18] was proposed to explain the survival of Rydberg states under strong laser fields. Channel closing induced by the ponderomotive shift of the ionization potential was identified to be responsible for the oscillations in the intensity dependence of excitation probabilities.[19–22] As a complement to the rescattering model, frustrated tunneling ionization[14,23] provided a more intuitive picture, where the liberated electron can be recaptured by the parent ion driven by the laser field. Regardless of which mechanism dominates, the excitation of Rydberg states can be considered as interference between wave packets created in every half laser cycle.[24,25] Therefore, the phases, including the transition phase and propagation phase, play crucial roles. The electron capture occurs every half cycle and its coherent superposition determines the final occupation of Rydberg states, depending on constructive or destructive interferences. However, how and when the Rydberg state is populated determines the phase of amplitude and, consequently, the final population of Rydberg states after the laser field, which remains unexplored.
In this work, we conducted numerical investigations on the coherent development of excited states in neon atoms driven by intense femtosecond pulses with varying pulse durations and intensities. We found that the excitation probability of Rydberg states is more efficient with shorter laser pulses and can be modulated by varying the laser intensities. Within the framework of strong-field approximation, we propose an analytical model to confirm that the freed electron after tunnel ionization surfs over the atomic potential and the Rydberg state is dominantly populated twice when the electron makes a U-turn within one optical cycle. By varying the peak intensity of the pulse, we observed modulations in the excitation probabilities of Rydberg states, which can be attributed to the intensity-dependent interference of the two capture events during the U-turn.
Numerical simulations were performed with the time-dependent Schrödinger equation (TDSE) describing the interaction of a one-dimensional model neon atom with an intense laser pulse (atomic units are used unless stated otherwise),(1) with α = 0.65 to reproduce the ionization potential of the real neon atom. The electric field of the laser pulse isE (t ) =E 0f (t )cos(ωt +φ ) with a Gaussian envelopef (t ) and carrier-envelope phase (CEP)φ , andω denotes the frequency of 800 nm. The calculations were performed with Gaussian laser pulses with different peak intensities and pulse durations. The populations of excited states (n = 2–20) as a function of driving laser intensity are illustrated in Fig.1(a) . The excitation probability exhibits pronounced oscillations with increasing laser intensity, which was observed due to the fact that the production of Rydberg states is associated with the occurrence of channel closing.[19 ,20 ,22 ,26 ] The higher excited states are easier to be coherently populated, in agreement with a previous study.[27 ] The energy of the high excited states is shifted approximately by the intensity-dependent ponderomotive energy of the electronU p =I /4ω 2. For the states with specific parity, if the energy levels of the excited states are shifted by the energy of two photons, the channel is open again and provides the condition for the channel closing as ΔU p = 2ω . For laser wavelength centered at 800 nm, the oscillation period of the excitation probability for a single excited state is about 0.5 × 1014 W/cm2, shown clearly with the populations of Rydberg states withn = 8 in Fig.1(c) andn = 9 in Fig.1(d) . The state-resolved population measurement is difficult in the experiment. The oscillation period for the total excitation yields of 0.25 × 1014 W/cm2 has been observed experimentally by measuring the intensity dependence of the excited atoms.[28 ] Moreover, we present the derivative of the Rydberg state phase after the laser pulse with respect to the laser intensity, as shown by the green lines in Figs.1(c) and1(d) . As the resonant excitation channel is closed, the phase of the Rydberg state changes dramatically. When the resonant excitation channel is open, i.e., the energy difference between the ground state and the laser-dressed Rydberg state equals multiple photon energies, the phase of the Rydberg state varies smoothly, where phase matching is easier to achieve and it leads to efficient excitation.Fig. 1. TDSE calculations for the coherent excitation probabilities of Rydberg states with 5 fs laser pulse and CEP of φ = 0. Results of excited states (n = 2–20) are shown in (a), and n = 8, 9 are shown in (b), with the laser intensity from 2 × 1014 W/cm2 to 2.5 × 1015 W/cm2. Phase variation dependence on the laser intensity for the Rydberg states of (c) n = 8 and (d) n = 9.Interestingly, except for the fast oscillations, it can be seen from Fig. 1(b) that a much slower oscillation appears as a function of the laser intensity, which depends on the quantum number. To highlight the slow oscillation appearing in the intensity dependence, populations of the Rydberg states are averaged over the laser intensities. The populations of Rydberg state n = 8 using laser pulses with different pulse duration are shown in Fig. 2. For few-cycle pulses, the population reached a maximum above 1015 W/cm2 and was followed by an oscillation dependence on the laser intensity. As the pulse duration increases, the total population decreased and the peak position shifted to a smaller intensity. Different Rydberg state populations show similar oscillation structures and pulse duration dependences, and the slight difference originates from their orbital properties.
As will be discussed in the following, the excitation behavior of the Rydberg states can be explained by the interference of multiple recapture events for each ionization time within the recollision picture. Considering a neon atom exposed to the strong field, an electron freed at timet i will have the probability to revisit the Rydberg state by acceleration of the laser field. Electrons revisiting at different times can be coherently recaptured into the Rydberg states. Interference between the two adjacent recapture times within a half laser cycle contributes to the low oscillation structure dependence on the laser intensity observed in the experiments. The recapture probability for each ionization timet i can be calculated based on the strong-field approximation,(2) W (t i) is the ionization rate given by the Ammosov–Delone–Krainov theory[29 ] and the ground-state depletion is taken into account using. The phase accumulated by the electron from the ionization time t i to them -th recapture timetm is, with the energy of the Rydberg state En and the vector potential of the laser fieldA (t ). The phase describes the Rydberg electron evolving with instantaneous energy with a laser-induced ac-Stark shift. The electron motion after ionization is calculated classically with zero initial velocity and position. Once the electron wave packet passes through the spatial range of the Rydberg state, the recapture occurs with the probabilitya r(t ). The total amplitude for the Rydberg excitation is obtained by the coherent sum over all ionization events,(3) acquired by the Rydberg electron from the capture time t ″ to the end of the laser pulse at timet f. For a slowly varying pulse envelope, we can write the final amplitude of occupying the Rydberg state as a sum over contributions from every half cycle (assume that the pulse duration isN /2 times laser period),(4) denotes the accumulated phase during each half cycle. Contributions from subsequent half cycles lead to the π phase difference because the electrons are recaptured in opposite directions of the same Rydberg state. Finally, the occupation probability can be written as(5) From the first term in Eq. (5), it can be seen that, when multiphoton resonance occurs, that is, ϕ0 = mπ, the ac-Stark shifted energy of the Rydberg state equals multiple photon energy En + Up = mω. For a state with specified parity, the occupation probability peaks are separated by ΔUp = 2ω, corresponding to the channel-closing effects observed with a fine adjustment of the laser intensity discussed in previous work. The second term arising from the intra-cycle phase accumulation Δϕ is attributed to the slow modulation observed in Fig. 2.
By analyzing the electron trajectories within each half laser cycle, we find that electrons born around the peaks of the electric field would have a typical ‘surfing’-like trajectory, where the electron turns around in the spatial range of the Rydberg state such that it can be recaptured coherently into the Rydberg state twice within a half laser cycle as indicated by t2 and t3 in Fig. 3. Here, we consider only the first surfing recapture event due to the quantum dispersion of the electron wave packet. The interference effect in the surfing trajectory contributes to the oscillatory structure in the laser intensity dependence of the Rydberg state population. Let us consider an electron born at the peak of the electric field during the rising edge of the laser pulse. The electron trajectory is calculated classically by neglecting the Coulomb force of the ion. The time interval between these two capture times, Δt = t3 – t2, in the surfing trajectory can be obtained by cos(ωΔt/2) = R0ωα/E0, with R0 being the radius of the Rydberg orbital. Here, E0 is the field amplitude of the laser pulse and α is the slope of the electric field envelope at the ionization time, which is inversely proportional to the pulse duration. The phase difference of the twin-capture event is approximated to be
. Constructive interference appears when Δϕ = 2nπ, leading to an effective Rydberg state population. Destructive interference occurs when Δϕ = (2n+1)π, decreasing the coherent excitation of the Rydberg state. Similar intra-cycle interference has been observed in above-threshold ionization spectra[30] and in the photoelectron holography[31,32] arising from a pair of electron trajectories released within one laser cycle. Fig. 3. The twin-capture model of Rydberg state excitation. Illustration of electron wave packet launch, reentry and twin events of capture into Rydberg orbitals at times t2 and t3. Electron wave packets are partially captured into the Rydberg orbital during each intersection. However, only those wave packets captured at times t2 and t3 can be in phase due to the smaller phase accumulation during the quick U-turn.To verify our proposed model, the coherent excitation probability calculated with Eq. (4) is shown in Fig. 4, considering only the contribution from ionization within one half cycle before the peak intensity of the pulse envelope. It can be seen that both the oscillation period and the optimal intensity decrease as the pulse duration increases, which agrees quite well with the TDSE simulations. The interference structure is determined by the phase difference that accumulates during the time interval Δt between the twin capture events, which is proportional to the laser intensity and the pulse duration. Therefore, the oscillation period is inversely proportional to the pulse duration. As the capture dynamics repeats at every half cycle, the coherent superposition of the Rydberg state also leads to the intra-cycle interference effect. However, for a few-cycle pulse with high intensity, the phase difference between adjacent half cycles can be as large as 38π due to ∼ 60 eV ponderomotive energy at 1 × 1015 W/cm2 pulse intensity, which makes the amplitudes out of phase. Therefore, the populating efficiency decreases as the pulse duration increases, confirmed by the TDSE simulations in Fig. 2.
In conclusion, the reported data demonstrate that the Rydberg population can be enhanced by one order of magnitude when the pulse duration is reduced to nearly a single-cycle laser pulse. Furthermore, the Rydberg population can be significantly modulated by changing the intensity of the driving pulses. The underlying physics is attributed to the crucial step where the returned electron makes a U-turn in the vicinity of Rydberg orbital, leading to constructively interfered twin-capture events. Our work advances new means of coherent Rydberg states excitation by tailoring laser waveforms, which benefits future Rydberg-state-based precision measurements and quantum detection.
Acknowledgments: This work was supported by the National Key Research and Development Program of China (Grant No. 2019YFA0307703), the National Natural Science Foundation of China (Grant Nos. 12234020, 11874066, 12274461, and 11974426), and the Science and Technology Innovation Program of Hunan Province (Grant No. 2022RC1193). -
References
[1] Jing M Y, Hu Y, Ma J, Zhang H, Zhang L J, Xiao L T, Jia S T 2020 Nat. Phys. 16 911 doi: 10.1038/s41567-020-0918-5[2] Facon A, Dietsche E K, Grosso D, Haroche S, Raimond J M, Brune M, Gleyzes S 2016 Nature 535 262 doi: 10.1038/nature18327[3] Luo L, Men L, Liu Z, Mudryk Y, Zhao X, Yao Y, Park J M, Shinar R, Shinar J, Ho K M, Perakis I E, Vela J, Wang J 2017 Nat. Commun. 8 15565 doi: 10.1038/ncomms15565[4] Ebadi S, Wang T T, Levine H, Keesling A, Semeghini G, Omran A, Bluvstein D, Samajdar R, Pichler H, Ho W W, Choi S, Sachdev S, Greiner M, Vuletić V, Lukin M D 2021 Nature 595 227 doi: 10.1038/s41586-021-03582-4[5] Vassen W, Cohen-Tannoudji C, Leduc M, Boiron D, Westbrook C I, Truscott A, Baldwin K, Birkl G, Cancio P, Trippenbach M 2012 Rev. Mod. Phys. 84 175 doi: 10.1103/RevModPhys.84.175[6] Zimmermann H, Buller J, Eilzer S, Eichmann U 2015 Phys. Rev. Lett. 114 123003 doi: 10.1103/PhysRevLett.114.123003[7] Grochmalicki J, Lewenstein M 1988 J. Phys. B: At. Mol. Opt. Phys. 21 3285 doi: 10.1088/0953-4075/21/19/017[8] Keldysh L V 1965 Sov. Phys. Usp. 8 496 doi: 10.1070/PU1965v008n03ABEH003058[9] Zhang D W, Lü Z H, Meng C, Du X Y, Zhou Z Y, Zhao Z X, Yuan J M 2012 Phys. Rev. Lett. 109 243002 doi: 10.1103/PhysRevLett.109.243002[10] Zhao J, Liu J L, Wang X W, Yuan J M, Zhao Z X 2022 Chin. Rev. Lett. 39 123201 doi: 10.1088/0256-307X/39/12/123201[11] Chini M, Wang X, Cheng Y, Wang H, Wu Y, Cunningham E, Li P C, Heslar J, Telnov D A, Chu S I, Chang Z 2014 Nat. Photon. 8 437 doi: 10.1038/nphoton.2014.83[12] Xiong W H, Geng J W, Tang J Y, Peng L Y, Gong Q 2014 Phys. Rev. Lett. 112 233001 doi: 10.1103/PhysRevLett.112.233001[13] Beaulieu S, Camp S, Descamps D, Comby A, Wanie V, Petit S, Legare F, Schafer K J, Gaarde M B, Catoire F, Mairesse Y 2016 Phys. Rev. Lett. 117 203001 doi: 10.1103/PhysRevLett.117.203001[14] Yun H, Mun J H, Hwang S I, Park S B, Ivanov I A, Nam C H, Kim K T 2018 Nat. Photon. 12 620 doi: 10.1038/s41566-018-0255-8[15] Tao W K, Wang L, Song P, Xiao F, Wang J C, Zheng Z G, Zhao J, Wang X W, Zhao Z X 2023 Chin. Phys. Lett. 40 063201 doi: 10.1088/0256-307X/40/6/063201[16] Chetty D, Glover R D, Tong X M, deHarak B A, Xu H, Haram N, Bartschat K, Palmer A J, Luiten A N, Light P S, Litvinyuk I V, Sang R T 2022 Phys. Rev. Lett. 128 173201 doi: 10.1103/PhysRevLett.128.173201[17] Matthews M, Morales F, Patas A, Lindinger A, Gateau J, Berti N, Hermelin S, Kasparian J, Richter M, Bredtmann T, Smirnova O, Wolf J P, Ivanov M 2018 Nat. Phys. 14 695 doi: 10.1038/s41567-018-0105-0[18] Ivanov I A, Kheifets A S, Kim K T 2022 Sci. Rep. 12 17048 doi: 10.1038/s41598-022-21549-x[19] Li Q G, Tong X M, Morishita T, Wei H, Lin C D 2014 Phys. Rev. A 89 023421 doi: 10.1103/PhysRevA.89.023421[20] Li Q G, Tong X M, Morishita T, Jin C, Wei H, Lin C D 2014 J. Phys. B 47 204019 doi: 10.1088/0953-4075/47/20/204019[21] Zimmermann H, Patchkovskii S, Ivanov M, Eichmann U 2017 Phys. Rev. Lett. 118 013003 doi: 10.1103/PhysRevLett.118.013003[22] Chetty D, Glover R D, deHarak B A, Tong X M, Xu H, Pauly T, Smith N, Hamilton K R, Bartschat K, Ziegel J P, Douguet N, Luiten A N, Light P S, Litvinyuk I V, Sang R T 2020 Phys. Rev. A 101 053402 doi: 10.1103/PhysRevA.101.053402[23] Nubbemeyer T, Gorling K, Saenz A, Eichmann U, Sandner W 2008 Phys. Rev. Lett. 101 233001 doi: 10.1103/PhysRevLett.101.233001[24] Liu M Q, Xu S P, Hu S L, Becker W, Quan W, Liu X J, Chen J 2021 Optica 8 765 doi: 10.1364/optica.418636[25] Hu S L, Hao X L, Lv H, Liu M Q, Yang T X, Xu H F, Jin M X, Ding D J, Li Q G, Li W D, Becker W, Chen J 2019 Opt. Express 27 31629 doi: 10.1364/OE.27.031629[26] Piraux B, Mota-Furtado F, O’Mahony P F, Galstyan A, Popov Y V 2017 Phys. Rev. A 96 043403 doi: 10.1103/PhysRevA.96.043403[27] Zhang B, Chen W, Zhao Z 2014 Phys. Rev. A 90 023409 doi: 10.1103/PhysRevA.90.023409[28] Xu S P, Liu M Q, Hu S L, Shu Z, Quan W, Xiao Z L, Zhou Y, Wei M Z, Zhao M, Sun R P, Wang Y, Hua L, Gong C, Lai X, Chen J, Liu X 2020 Phys. Rev. A 102 043104 doi: 10.1103/PhysRevA.102.043104[29] Ammosov M V, Delone N B, Krainov V P 1986 Sov. Phys. JETP 64 1191 doi: 10.1117/12.938695[30] Arbó D G, Ishikawa K L, Schiessl K, Persson E, Burgdörfer J 2010 Phys. Rev. A 82 043426 doi: 10.1103/PhysRevA.82.043426[31] Meckel M, Staudte A, Patchkovskii S, Villeneuve D M, Corkum P B, Dorner R, Spanner M 2014 Nat. Phys. 10 594 doi: 10.1038/nphys3010[32] Bian X B, Bandrauk A D 2012 Phys. Rev. Lett. 108 263003 doi: 10.1103/PhysRevLett.108.263003 -
Related Articles
[1] Man Xing, Jun Wang, Xi Zhao, Shushan Zhou. The Role of Multi-Electron and Multi-Orbital Effects in High-Order Harmonic Generation of Benzonitrile Molecules [J]. Chin. Phys. Lett., 2025, 42(4): 043201. doi: 10.1088/0256-307X/42/4/043201 [2] LI Na-Na, ZHAI Zhen, LIU Xue-Shen. High-Order Harmonic Generation from a Model of Ar+ Ionized Clusters [J]. Chin. Phys. Lett., 2008, 25(7): 2508-2510. [3] GE Yu-Cheng. Laser-Duration Dependence of Emission Properties of High-Order Harmonic Generation [J]. Chin. Phys. Lett., 2008, 25(4): 1305-1308. [4] GE Yu-Cheng. Laser Phase Relations of High-Order Harmonic Generation [J]. Chin. Phys. Lett., 2006, 23(9): 2461-2464. [5] PI Liang-Wen, SHI Ting-Yun, QIAO Hao-Xue. Enhancement of Bichromatic High-Order Harmonic Generation by a Strong Laser Field and Its Third Harmonic [J]. Chin. Phys. Lett., 2006, 23(6): 1490-1493. [6] WANG Bing-Bing, CHEN Jing, LIU Jie, LI Xiao-Feng, FU Pan-Ming. Carrier Envelope Phase Controlled High-Order Harmonic Generation in Ultrashort Laser Pulse [J]. Chin. Phys. Lett., 2005, 22(9): 2237-2240. [7] WANG Run-Hai, JIANG Hong-Bing, YANG Hong, WU Cheng-Yin, GONG Qi-Huang. High-Order Harmonic Generation by Two Non-Collinear Femtosecond Laser Pulses in CO [J]. Chin. Phys. Lett., 2005, 22(8): 1913-1915. [8] CHEN Jing, CHEN Shi-Gang, LIU Jie. High-Order Harmonic Generation in the Ionization Process [J]. Chin. Phys. Lett., 2000, 17(10): 723-725. [9] GAO Liang-hui, LI Xiao-feng, GUO Dong-sheng, FU Pan-ming. Formal Scattering Approach to High-Order Harmonic Generation [J]. Chin. Phys. Lett., 1999, 16(7): 502-504. [10] WANG Bing-bing, LI Xiao-feng, FU Pan-ming. High-Order Harmonic Generation in the Presence of Static Electric Field [J]. Chin. Phys. Lett., 1998, 15(3): 195-197.