Doppler effects in the propagation of a few-cycle pulse through a dense medium

This numerical study demonstrates that Doppler redshift exists in the reflected spectrum of a few-cycle pulse, propagating through a dense medium. It manifests itself in two different forms, a sharp low-frequency spike (LFS) located at the red edge of the reflected spectrum and a relatively broader redshift near the carrier frequency. With the variation of the laser and medium parameters, the dominant reflection mechanism changes between bulk generation of backwards propagation waves and nonlinear reflection near the front face. This leads to the manifestation of Doppler effect changing accordingly between the two different forms. This study unifies the physical mechanism behind the LFS and dynamic nonlinear optical skin effect, which enriches the theoretical explanation of the spectral redshift of few-cycle pulse propagation beyond the intrapulse four-wave mixing. © 2015 Optical Society of America OCIS codes: (320.2250) Femtosecond phenomena; (190.0190) Nonlinear optics; (190.7110) Ultrafast nonlinear optics; (320.0320) Ultrafast optics; (190.5530) Pulse propagation and solitons. References and links 1. A. Baltuska, Z. Wei, M. S. Pshenichnikov, and D. A. Wiersma, “Optical pulse compression to 5 fs at a 1-MHz repetition rate,” Opt. Lett. 22(2), 102–104 (1997). 2. M. Nisoli, S. De Silvestri, O. Svelto, R. Szipöcs, K. Ferencz, Ch. Spielmann, S. Sartania, and F. Krausz, “Compression of high-energy laser pulses below 5 fs,” Opt. Lett. 22(8), 522–524 (1997). 3. M. Wegener, Extreme Nonlinear Optics (Springer, 2005). 4. Q. Q. Xu, D. Z. Yao, X. N. Liu, Q. Zhou, and G. G. Xiong, “Solitary propagation effect of a well-defined chirped femtosecond laser pulse in a resonance-absorbing medium,” Phys. Rev. A 86(2), 023853 (2012). 5. Z. D. Wang, T. T. Liu, J. L. Ma, and J. Xiao, “Interference resonant propagation and spectral properties of two-colour femtosecond Gaussian pulses in three-level Λ-type atomic medium,” Opt. Commun. 330, 155–159 (2014). 6. L. Guo, X. T. Xie, and Z. M. Zhan, “Shaping of few-cycle laser pulses via a subwavelength structure,” Chin. Phys. B 22(9), 094212 (2013). 7. J. K. Ranka, R. W. Schirmer, and A. L. Gaeta, “Coherent spectroscopic effects in the propagation of ultrashort pulses through a two-level system,” Phys. Rev. A 57(1), R36–R39 (1998). 8. R. K. Bullough, P. M. Jack, P. W. Kitchenside, and R. Saunders,“Solitons in laser physics,” Phys. Scr. 20(3), 364–381 (1979). 9. J. C. Eilbeck, “Reflection of short pulses in linear optics,” Gen. Phys. 5(9), 1355–1363 (1972). 10. J. C. Eilbeck and R. K. Bullough, “The method of characteristics in the theory of resonant or nonresonant nonlinear optics,” Gen. Phys. 5(6), 820–829 (1972). #240039 Received 5 May 2015; revised 15 Jun 2015; accepted 16 Jun 2015; published 19 Jun 2015 © 2015 OSA 29 Jun 2015 | Vol. 23, No. 13 | DOI:10.1364/OE.23.017046 | OPTICS EXPRESS 17046 11. V. P. Kalosha and J. Herrmann, “Formation of optical subcycle pulses and full Maxwell–Bloch solitary waves by coherent propagation effects,” Phys. Rev. Lett. 83(3), 544–547 (1999). 12. X. T. Xie and M. A. Macovei, “Single-cycle gap soliton in a subwavelength structure,” Phys. Rev. Lett 104(7), 073902 (2010). 13. R. Marskar and U. L. Österberg, “Backpropagation and decay of self-induced-transparency pulses,” Phys. Rev. A 89(2), 023828 (2014). 14. K. S. Yee, “Numerical solution of initial boundary value problem involving Maxwell’s equation in isotropic media,” IEEE Trans. Antennas Propag. 14(3), 302–307(1996). 15. R. W. Ziolkowski, “Ultrafast pulse interactions with two-level atoms,” Phys. Rev. Lett. 52(4), 3082–3094 (1995). 16. R. Marskar and U. L. Österberg, “Multilevel Maxwell–Bloch simulations in inhomogeneously broadened media,” Opt. Express 19(18), 16784–16796 (2011). 17. W. Forysiak, R. G. Flesch, and J. V. Moloney, “Doppler shift of self-reflected optical pulses at an interface: dynamic nonlinear optical skin effect,” Phys. Rev. Lett. 76(20), 3695–3698 (1996). 18. L. R. Franco, “Self-reflected wave inside a very dense saturable absorber,” Phys. Rev. Lett. 55(20), 2149–2151 (1985). 19. L. R. Franco, “Propagation of light in a nonlinear absorber,” J. Opt. Soc. Am. B 4(11), 1878–1884 (1987). 20. L. Allen and J.H. Eberly, Optical Resonance and Two-Level Atoms (Wiley, 1975). 21. J. Xiao, Z. Y. Wang, and Z. Z. Xu, “Spectrum of a few-cycle laser pulse propagation in a two-level atom medium,” Chin. Phys. 10(10), 0941–0945 (2001). 22. J. Xiao, Z. Y. Wang, and Z. Z. Xu, “Area evolution of a few-cycle pulse laser in a two-level-atom medium,” Phys. Rev. A 65(3) R031402 (2002). 23. C. W. Luo, K. Reimann, M. Woerner, T. Elsaesser, R. Hey, and K. H. Ploog, “Phase-resolved nonlinear response of a two-dimensional electron gas under femtosecond intersubband excitation,” Phys. Rev. Lett. 92(4), 047402 (2004). 24. J. Xiao, Z. Y. Wang, and Z. Z. Xu, “Propagation of few-cycle pulse laser in two-level atom medium,” Chin. Phy. Lett. 18(7), 912 (2001). 25. R. K. Bullough, P. J. Caudrey, J. C. Eilbeck, and J. D. Gibbon, “A general theory of self-induced transparency,” Opt. Quantum Electron. 6(1), 121–140 (1974).


Introduction
With the rapid development in ultrafast laser technology [1,2], it has become more attractive to extend the study of the interaction between matter and light into extreme nonlinear regime [3].The fundamental problem in this field is the interaction of few-cycle pulses with resonant twolevel atomic (TLA) ensembles, which lays the foundation for understanding the more complex cases, such as the quantum dot ensembles [4], three-level atoms [5], and specifically designed structures [6].
The transmitted spectrum of few-cycle pulses propagating in vacuum to a TLA material has been studied extensively by the previous works [4][5][6][7].However, few studies have investigated the reflected spectrum.This is because many previous works are based on the envelope form of the Maxwell equations, where forward propagation approximation is implied and the reflected field is ignored [8][9][10].Nevertheless, if the density is high enough to allow the backpropagation to make a difference, the full Maxwell-Bloch (MB) equations [11][12][13] without rotating wave approximation (RWA) and slowly varying envelope approximation (SVEA) have to be integrated.To solve these equations, one usually uses Yee's finite difference timedomain method (FDTD) assisted by the predictor-corrector method [14,15], or a pseudospectral operator-splitting method combined with a hybrid interface [16].
A dense TLA medium prepared in the ground state can be regarded as a saturable absorbing material.An intense pulse incident upon such a medium can build an absorption front within an optical wavelength [17].If the induced spatial inhomogeneity is sufficiently sharp, this absorption front will act as an interface for self-reflection [18,19].Since this reflector moves with the propagating pulse, a doppler redshift appears in the reflection spectrum.This is known as the dynamic nonlinear optical skin (DNOS) effect, which has been used to explain the slight redshift in reflected spectrum for long pulse incident on the interface [17].However, for the few-cycle pulse, theoretical studies about the behavior of the reflected spectrum are far from complete.The existing theories on the origin of the redshift for the few-cycle case are confined to intrapulse four-wave mixing (FWM) [11,12].Whether absorption fronts exist in the medium, whether the Doppler redshift plays a role in the reflected spectrum, and in which way it is manifested still remain unsolved.
In this paper, we study the reflected field of a few-cycle pule propagating through a dense TLA medium by solving the MB equations without the rotating wave and slowly varying envelope approximations.We prove that absorption fronts exist even in the few-cycle cases and that the Doppler redshift induced by the resulting moving front contributes to the formation of the reflected spectrum.We also clarify how the manifestation of Doppler effect in the reflected spectrum depends on the dominant reflection mechanism, which in turn can be modified with laser duration, area, medium density and length.

Theoretical model and results
First, we consider the propagation of a few-cycle pulse along z direction through a resonant TLA.The MB equations take the form [20][21][22] ∂ ∂ The macroscopic coherent polarization is P x (z) = Ndu, where u is the real part of the offdiagonal density matrix element ρ 12 = (u + iv)/2.w = ρ 22 − ρ 11 is the population inversion between the excited state 2 and the ground state 1.The lifetime of the excited state T 1 = 1ps and the dephasing time T 2 = 0.5ps.The Rabi frequency Ω = dE x /h with dipole moment d = 2 × 10 −29 A s m.The carrier frequency is ω p = 2.3fs −1 , equal to the transition frequency ω 0 between ground and excited states.The TLA medium is initially prepared in the ground state, w 0 = −1, u = v = 0.The medium length L = 45μm.The collective frequency parameter ω c = Nd 2 /ε 0 h = 0.05fs −1 , where N = 1.1 × 10 20 cm −3 is the medium density.The incident pulse has a hyperbolic secant envelop, Ω(t = 0, z) = Ω 0 cos[ω p (z−z 0 )/c] sech[1.76(z−z0 )/(cτ p )], where Ω 0 = 0.7fs −1 is the peak Rabi frequency, τ p = 5fs the full width at half maximum (FWHM), and z 0 = 26.25μm is used to avoid the initial disturbance from the input pulse.The corresponding envelope area A = Ω 0 τ p π/1.76 = 2π.
The FDTD method assisted by the predictor-corrector method is used to solve the full MB equations [14,15].The reflected field and the corresponding spectrum obtained with the above parameters are shown in Fig. 1.The reflected field in Fig. 1(a) consists of a leading part produced by front surface reflection, and a tail generated by backpropagation waves.The reflected fields from the back surface and the remaining atomic oscillation are irrelevant to the Doppler redshift in question.Hence, the reflections corresponding to these two parts are omitted.Figure 1(b) shows that, besides the double-peak structure, a low-frequency spike (LFS) centered at 0.215ω 0 with an amplitude of 0.7×10 −4 fs −1 appears at the red edge of the reflected spectrum.To find out the origin of the LFS, we take the respective Fourier transformations to the tail and leading parts.The results show that, this LFS and the double-peak structure correspond to the tail and leading part of the reflected fields, respectively.The spectrum of the reemitted field produced by macroscopic coherent polarization at z = 1.5μm is shown in Fig. 1(b) [23].Apparently, it covers most of the frequency components in reflected spectrum except for the LFS.That is to say, nonlinear polarization is responsible for the spectral broadening of the double-peak structure, but has no contribution to the generation of the LFS.
Nonlinear polarization may include many nonlinear effects, such as intrapulse FWM and self-phase modulation (SPM).According to the previous works on few-cycle pulse propagation, intrapulse FWM, which is ignored in the framework of RWA, is significant for very short pulses.It is widely accepted as the nonlinear effect that determines the spectral transformation due to the ultrashort feature of the few-cycle pulse [11].Moreover, SPM may also contribute to the spectral widening.Refractive index does not change under RWA, but it changes with the carrier in the few-cycle case.Thus, SPM, which is an ignored nonlinear process in the former case, can be very important for the short pulse propagation.However, which nonlinear effects play the dominant role in the spectral broadening is determined by the laser and medium parameters.As for the parameters in Fig. 1, the SPM may be the dominant nonlinear effect that results in a symmetrical broadening, while for the parameters used in Fig. 5, as in [11], intrapulse FWM probably is the dominant nonlinear effect that results in redshift peaks excluding the Doppler redshift.
According to the dynamic nonlinear optical skin (DNOS) effect [17], when a long pulse is incident on a saturable medium, an absorption front is built up near the interface and its motion within the nonlinear skin depth produces a slight redshift that appears in the reflected spectrum.In contrast, for few-cycle pulses, the absorption fronts can travel beyond the nonlinear skin depth with a greater speed, as shown in the following.When the spatial structure of the self-reflected waves are short enough to yield constructive interference in the backwards direction, low-frequency backpropagation waves are generated.Therefore, velocity-dependent low-frequency components are expected in the reflected spectrum.To validate this physical picture, we have determined the absorption front velocity from the numerical simulation of the population inversion.The locations z i at which the population inversion occurred at t i are recorded with a time interval Δt = t i+1 − t i = 12.5fs.The average velocity is then obtained using υ i = z i+1 −z i Δt .Figure 1(c) shows that the velocities of the absorption front of few-cycle pulse are relatively stable and distribute within 1.79 ∼ 2.1×10 8 m/s.The relation between velocity υ and frequency ω of the self-reflected waves fulfills ω = c−υ c+υ ω 0 , with which we get a frequency band [0.178ω 0 , 0.253ω 0 ] corresponding to the velocity range.Surprisingly, the frequencies of backpropagation waves are exactly within the frequency band of the LFS, which justifies the our assumption that a Doppler effect inducing redshift in backwards propagation waves is the source of the LFS, as shown in Fig. 1(c).Thus, besides the FWM, the Doppler effect also plays an important role in the reflected spectra of the few-cycle pulses.In the following, we discuss how the laser and medium parameters influence the self-reflections, and how the Doppler effect manifests itself in the reflected spectrum.To further prove that the physical mechanism behind the LFS is unified with the DNOS effect, we study the reflected spectrum with the parameters in section 2, except for a variable duration τ p .A pulse with duration τ p = 10fs is used to study the influence of pulse duration on the Doppler redshift.The results are shown in the top of Fig. 2. It can be seen from Fig. 2(a) that the LFS produced by backpropagation waves disappears.When a 2π pulse propagates through a TLA medium, the energy transfer from the pulse leading edge to the trailing edge due to Rabi flop, which results in a slowdown of the propagating pulse.Obviously, the longer the incident pulse is, the more slowly it propagates [24].Since the amount of Doppler redshift is proportional to the propagation velocity, the redshift of the LFS for a 10fs pulse decreases compared with that for a 5fs pulse.Moreover, as the pulse area is fixed, increasing duration τ p is equivalent to reducing peak intensity.Since the condition for reflection to play a significant role is Ω ∼ ω c u [25], 10fs pulse is more vulnerable to reflections.The dramatic reflection loss causes deceleration and prohibits pulse propagating through the medium.As shown in Fig. 2(b), both the energy and propagation velocity of the main pulse decrease until settling around 135μm.Compared with the case of 5fs pulse in Fig. 1(c), the absorption front of 10fs pulse travels more slowly and decelerates from 9.09×10 7 m/s to 8.16×10 6 m/s within 625fs, as shown in Fig. 2(c).According to the relation between velocity and frequency, the corresponding frequency band [0.53ω 0 , 0.95ω 0 ] is obtained.This shows that the Doppler redshift no longer manifests itself as an isolated LFS but is merged with the spectrum of the reemitted field.

Influence of pulse duration
As shown in Fig. 2(f), the propagation velocity of a 40fs pulse further decreases from 9.4×10 6 m/s to 1.44×10 6 m/s within 150fs.The corresponding frequency obtained by velocitydependent theory ranges from 0.94ω 0 to 0.99ω 0 .This frequency band is consistent with the LFS obtained by numerical simulation in Fig. 2 (d), in terms of its center frequency and FWHM.For a 40fs pulse, the nonlinear reflection near the front surface becomes so intense that the pulse barely propagates within the medium.This is demonstrated in in Fig. 2 (e), where the penetration depth is less than 2 μm.Therefore, when it comes to the reflection mechanism of a 40fs pulse, the bulk generation of backpropagation waves is trivial, while the nonlinear reflection near the interface becomes dominant, i.e., DNOS effect.
For pulses with a fixed area, the longer the pulse is, the more severely it will suffer from reflections.As the duration increases from 5fs to 40fs, both the nonlinear reflections near front surface and backwards propagation inside the medium increase.Consequently, the propagation velocity of the main pulse decreases, so does the velocity of the moving absorption front following it.According to the velocity-frequency relation, the frequency of the LFS increases until it merges into the redshift directly produced by nonlinear effects.At this point, the LFS disappears and the dominant reflection mechanism changes to be the nonlinear reflection near the front surface instead of the bulk generation of backward propagation waves.Thus, the underlying physics of the LFS is essentially consistent with that of the DNOS effect.That is, the Doppler redshift induced by a moving absorption front manifests itself in two forms: one is the LFS for short pulses, and the other is the DNOS effect for long pulses.

Influence of pulse area
Fig. 3.The reflected (solid line) and reemitted spectra (dashed line) (a)(e).Time evolution of velocities (b)(f) and corresponding frequencies (b)(g).The top and bottom correspond to 3π and 4π pulses, respectively.Subscripts 1 and 2 in (f)(g) stand for the first and second moving fronts respectively.(c) The instantaneous fields of 2π ( dotted line), 3π (dashed line), and 4π(solid line) incident pulses.The insert in (a) is the enlarged view to the LFS for a 3π incident pulse.
Besides pulse duration, pulse area also has a significant influence on the behavior of the Doppler redshift.Here, we give a discussion about the behavior of the Doppler redshift with respect to different pulse areas.For a 3π pulse, the reflected spectrum, as shown in Fig. 3(a), consists of a strong spike centered at ω 0 and a LFS located at 0.054ω 0 .Compared with that of a 2π pulse, the LFS of a 3π pulse has a lower frequency and a larger amplitude of 5.06×10 −4 fs −1 .According to the area theory, when a 3π pulse propagates through a TLA medium, its area decreases to 2π through splitting.Since the 2π pulse obtained in this way is much more intense and shorter than the incident 2π pulse, the former can propagate much faster than the latter, as shown in Fig. 3(c).Accordingly, the absorption front generated by the 3π pulse has a higher speed compared to that of the 2π pulse, which results in a larger redshift of the LFS. Figure 3(b) shows that, the velocity of the moving front ranges from 2.656×10 8 m/s to 2.7×10 8 m/s, and the corresponding frequency is within [0.052ω 0 , 0.06ω 0 ].The obtained frequencies are consistent with the frequency band of the LFS shown in Fig. 3(a).Meanwhile, since the amplitude of the reflection is proportional to the injected energy, the amplitude of the LFS is expected to increase.Moreover, after a complete Rabi flop, the remaining energy produces a small area pulse, which is supposed to propagate forwards following the main 2π pulse [22].However, for the dense medium in question, the small area pulse suffers from dramatic reflections and propagates backwards instead of forwards, which generates the strong spike with center frequency, as shown in Fig. 3(a).
As for a 4π pulse, both the amount of redshift and amplitude of the LFS increase.The LFS shown in Fig. 3(e) has a center frequency of 0.027ω 0 and an amplitude of 5.6×10 −4 fs −1 .This pulse splits into two 2π pulses with different propagation velocities, and two accompanied absorption fronts are induced.It can be seen in Fig. 3(f) that the first split pulse travels with a stable velocity much faster than the second.Hence, the backwards propagation waves produced by these two fronts have distinct frequency bands.From Figs. 3(f) and 3(g) we can see that, for the first pulse, the velocity of the moving front and frequency of the backpropagation waves are stabilized around 2.84×10 8 m/s and 0.027ω 0 , respectively.The velocity of the second front continuously decline from 1.992×10 8 m/s to 1.868×10 8 m/s, generating backpropagation waves with frequency within [0.202ω 0 , 0.233ω 0 ].Therefore, the LFS appeared in Fig. 3(e) is essentially the Doppler redshift induced by the first moving front, while the second moving front is responsible for the redshift near 0.2ω 0 .Since the first 2π pulse split from 4π pulse is more intense and shorter than that split from the 3π pulse, a LFS with a larger redshift is expected in the reflected spectrum of a 4π pulse.Thus, the redshift of the LFS is proportional to the pulse area, if the other parameters fixed.

Influence of medium density
So far, we have investigated the influences of the laser parameters on the manifestation of the Doppler redshift.Next, the influence of medium parameters, such as medium density and length, are discussed in the following sections.First, we take 2π and 4π pulses as examples to analyze the influence of medium density.The changes in frequency and amplitude of the LFS with respect to density for 2π and 4π pulses are shown in Figs.4(a) and 4(b), respectively.With the increase of the density, the changes of amplitude and redshift of the LFS have an upward trend and a downward trend, respectively.As the density increases, the reflections enhance while the propagation velocity of the main pulse reduces.Thus, the LFS has a larger amplitude but a smaller redshift.Figure 4(c) describes the influence of the density on the reflected spectrum of a 2π pulse.As ω c increases from 0.05fs −1 to 0.2fs −1 , the amplitude of the LFS increases as a result of the enhancement in reflections.Meanwhile the amount of redshift decreases due to the slowdown of the moving absorption front.Moreover, the LFS broadens as it moves towards the center frequency ω 0 .This is because, the enhanced reflections prohibit forward propagation and result in deceleration of the moving absorption front.Consequently, Fig. 4. Dependence of amplitude (blue squares) and frequency (red circles) of the LFS on medium density for 2π (a) and 4π (b) pulses, respectively.(c)The variation of the reflected spectra with increasing medium density for a 2π pulse.
the wider velocity range of the absorption front results in a broader frequency band of the reflected field.Therefore, for a relatively dilute medium, the Doppler effect manifests itself as a sharper LFS with a lower frequency and a lower amplitude.However, with the increase of the medium density, the LFS broadens and strengthens as it moves towards ω 0 .When the medium is dense enough to forbid propagating beyond the nonlinear skin depth, as the case of ω c =0.3fs −1 , the LFS is replaced by the redshift induced by the DNOS effect.Note that, the spectrum of ω c =0.3fs −1 in Fig. 4(c) is similar to that shown in Fig. 2(d), implying that the underlying physics is actually the same.That is, the process of increasing density is essentially another case where the nonlinear reflections near the surface gradually take over from the bulk generation of backpropagation waves.

Influence of medium length
Now let us investigate the influence of the medium length.Considering the case where a 4π pulse propagates through a TLA system with a density ω c = 0.1fs −1 , we change the medium length from L = 45μm to L = 150μm.The reflected and transmitted spectra are shown in Figs.5(a) and 5(b), respectively.It can be seen that, when the medium length changes to 150μm, the blueshift in the transmitted spectrum increases by 0.05ω 0 .This z-dependent maximum shift is produced by intrapulse FWM and can be explained with the phase-matching condition [11].However, regarding the redshift shown in the reflected spectrum Fig. 5(a), it barely moves as medium length increases.Thus, it is fair to suggest that the intrapulse FWM, which induces the blueshift in transmitted spectrum, is not responsible for the z-independent redshift.This redshift is actually the Doppler redshfit produced by the moving absorption front.
To prove this, we investigate the variation of the redshift versus medium length, in terms of the maximum redshift and its amplitude.The medium density is set as ω c = 1.0fs −1 .By increasing the medium length from λ 0 /4 to 45μm, the reflected spectrum is obtained, as shown in Fig. 5(c).The maximum redshift in the reflected spectrum increases dramatically within a subwavelength scale then levels off after L = λ 0 .This variation of the maximum redshift against the medium length can be predicted by the velocity-dependent theory.At the beginning of the propagation, the nonlinear reflection near the front surface plays a dominant role.At this stage, the main pulse propagates with a relatively low velocity, so does the induced absorption front near the surface.According to the velocity-dependent theory, a slight Doppler redshift appears in the self-reflected waves.However, once the main pulse crosses the surface, it experiences a dramatic acceleration within a subwavelength scale, which results in the increase of the corresponding redshift.After then, the main pulse propagates with a relatively stable velocity, so does the induced absorption front in the medium.The dominant reflection mechanism changes to be the bulk generation of backpropagation waves.In this case, the stable velocity of the absorption front indicates a constant frequency of the backpropagation wave, forming a LFS in the reflected spectrum.This is why the maximum redshift is insensitive to the change of the medium length after a subwavelength scale.
The amplitude of the maximum redshift is expected to increase as the propagation length increases.However, the fact is that the growth rate of the amplitude is gradually vanishing as the propagation length gets longer.As shown in Fig. 5(c), the growth rate of the amplitude is 178% from z = λ 0 /4 to z = 4λ 0 , while only 22.2% from z = 4λ 0 to z = 54λ 0 (45μm).This is because, within a few wavelength scale, the resonance frequency ω 0 of the main pulse is still available to support the increase of the LFS amplitude, as shown in Fig. 5(d).While for a longer distance, the transmitted pulse is blueshifted and gradually centered on 1.5ω 0 due to the nonlinear propagation effects.Since the large detuning prohibits the energy exchange between atoms and the electric field, the atoms prepared in ground states are rarely excited and the absorption front then disappears.This leads to the decrease of the growth rate of the LFS's amplitude along the propagation path.Thus, only if the medium is long enough to allow the absorption front to reach the stable velocity, the LFS, if it exists for the medium considered, can be observed in the spectrum.Otherwise, the bulk generation of backpropagation is missed due to the limited medium length, leaving a slight redshift produced by the nonlinear reflection near the front surface.Moreover, with the increase of the length, the LFS's amplitude keeps increasing until the center frequency of the main pulse is largely detuned.

Conclusion
In conclusion, we have made a systematic analysis of the Doppler effect of the absorption front with variable laser and medium parameters.The manifestation of the Doppler effect in the reflected spectrum depends sensitively on the laser and medium parameters, and generally takes two forms: a broad redshift near carrier-wave frequency of the incident pulse and a sharp LFS isolated from the main redshift.The former is produced by the nonlinear reflection near the front surface of the medium, while the latter corresponds to the bulk generation of the backwards propagation waves inside the medium.The speed of the laser pulse near the front surface is much lower than that inside the medium, and the amount of the Doppler redshift is dependent on the velocity.Thus, the redshift produced by the nonlinear reflection near the surface is much smaller than that by the backpropagation waves.Moreover, for the case where nonlinear reflections near the front surface is dominant, the absorption front experiences a deceleration and only penetrates a nonlinear skin depth.The corresponding frequency of the self-reflected field is distributed in a relatively broad band smaller than the carrier-wave frequency.The Doppler redshift manifests itself as the DNOS effect.By reducing pulse duration, increasing pulse area, reducing the medium density, or increasing the medium length, the dominant reflection mechanism gradually changes from the nonlinear reflection near the front surface to the bulk generation of backpropagation waves.When the latter plays the dominant role, the velocity of the moving absorption front is relatively stable and close to vacuum light speed c.In this case, the Doppler redshift manifests itself as a sharp LFS located at the red edge of the reflected spectrum.Our study unifies the underlying physics between the DNOS effect and the observed LFS.Moreover, we investigate how the laser and medium parameters affect the performance of the Doppler effect in the reflected spectrum, and how the behavior of Doppler redshift changes against parameters' variation.These results are crucial to a further understanding to the reflected spectrum of a few-cycle pulse propagating through dense media, and provide a basis for assessing the performance of theories based on the forward propagation approximation.

Fig. 1 .
Fig. 1.(a) The reflected field of a few-cycle pulse propagation through a dense TLA medium, which consists of a leading (thick line) and tail parts (thin line).(b) The corresponding reflected (solid line) and reemitted field (dashed line) spectra.(c) The velocity variations of the moving absorption front versus time and the corresponding frequencies obtained with ω/ω 0 = c−υ c+υ .The horizontal dashed line is the central frequency of the LFS, ω/ω 0 = 0.215, and the two horizontal solid linea represent its FWMH τ = 0.063ω 0 .

Fig. 2 .
Fig. 2. (a)(d)The reflected spectrum.(b)(e) The spatial evolution of the absorption front.(c)(f) The time evolution of the absorption front's velocity (circles) and the backpropagation wave's frequency (squares).The top and bottom correspond to τ p = 10fs and τ p = 40fs, respectively.The dashed line in (d) is the spectrum of the incident pulse.

Fig. 5 .
Fig. 5.The left and right are the reflected and transmitted spectra, respectively.The top and bottom are that for ω c = 0.1fs −1 and ω c = 1.0fs −1 , respectively.