Evolution of the frequency chirp of Gaussian pulses and beams when passing through a pulse compressor

Evolution of the frequency chirp of Gaussian pulses and beams when passing through a pulse compressor Derong Li, 3, Xiaohua Lv *, Pamela Bowlan, Rui Du, Shaoqun Zeng, Qingming Luo Britton Chance Center for Biomedical Photonics, Wuhan National Laboratory for Optoelectronics, Huazhong University of Science & Technology, Wuhan, 430074, China Georgia Institute of Technology, School of Physics 837 State Street NW, Atlanta, Georgia 3033 USA 3 Key Lab for Biomedical Informatics and Health Engineering, Institute of Biomedical and Health Engineering, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, 58055, China *xhlv@mail.hust.edu.cn Abstract: The evolution of the frequency chirp of a laser pulse inside a classical pulse compressor is very different for plane waves and Gaussian beams, although after propagating through the last (4th) dispersive element, the two models give the same results. In this paper, we have analyzed the evolution of the frequency chirp of Gaussian pulses and beams using a method which directly obtains the spectral phase acquired by the compressor. We found the spatiotemporal couplings in the phase to be the fundamental reason for the difference in the frequency chirp acquired by a Gaussian beam and a plane wave. When the Gaussian beam propagates, an additional frequency chirp will be introduced if any spatiotemporal couplings (i.e. angular dispersion, spatial chirp or pulse front tilt) are present. However, if there are no couplings present, the chirp of the Gaussian beam is the same as that of a plane wave. When the Gaussian beam is well collimated, the introduced frequency chirp predicted by the plane wave and Gaussian beam models are in closer agreement. This work improves our understanding of pulse compressors and should be helpful for optimizing dispersion compensation schemes in many applications of femtosecond laser pulses. 009 Optical Society of America OCIS codes: (30.030) Ultrafast optics; (30.550) Pulse compression; (60.030) Dispersion References and links. E. B. Treacy, Optical pulse compression with diffraction gratings, IEEE J. Quantum Electron. 5(9), 454 458 (969).. R. L. Fork, O. E. Martinez, and J. P. Gordon, Negative dispersion using pairs of prisms, Opt. Lett. 9(5), 50 5 (984). 3. J. Squier, F. Salin, G. Mourou, and D. Harter, 00-fs pulse generation and amplification in Ti:AI O 3, Opt. Lett. 6(5), 34 36 (99). 4. C. L. Blanc, G. Grillon, J. P. Chambaret, A. Migus, and A. Antonetti, Compact and efficient multipass Ti:sapphire system for femtosecond chirped-pulse amplification at the terawatt level, Opt. Lett. 8(), 40 4 (993). 5. O. E. Martinez, 3000 times grating compressor with positive group velocity dispersion: application to fiber compensation in.3-.6 um region, IEEE J. Quantum Electron. 3(), 59 64 (987). 6. S. Zeng, D. Li, X. Lv, J. Liu, and Q. Luo, Pulse broadening of the femtosecond pulses in a Gaussian beam passing an angular disperser, Opt. Lett. 3(9), 80 8 (007). 7. D. Li, X. Li, S. Zeng, and Q. Luo, A generalized analysis of femtosecond laser pulse broadening after angular dispersion, Opt. Express 6(), 37 47 (008). 8. D. Li, S. Zeng, Q. Luo, P. Bowlan, V. Chauahan, and R. Trebino, Propagation dependence of chirp in Gaussian pulses and beams due to angular dispersion, Opt. Lett. 34(7), 96 964 (009). 9. S. Szatmari, G. Kuhnle, and P. Simon, Pulse compression and traveling wave excitation scheme using a single dispersive element, Appl. Opt. 9(36), 537 5379 (990). 0. S. Szatmári, P. Simon, and M. Feuerhake, Group velocity dispersion compensated propagation of short pulses in dispersive media, Opt. Lett. (5), 56 58 (996). #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7070

. S. Akturk, X. Gu, M. Kimmel, and R. Trebino, Extremely simple single-prism ultrashort- pulse compressor, Opt. Express 4(), 00 008 (006).. M. Nakazawa, T. Nakashima, and H. Kubota, Optical pulse compression using a TeO acousto-optical light deflector, Opt. Lett. 3(), 0 (988). 3. S. Zeng, X. Lv, C. Zhan, W. R. Chen, W. Xiong, S. L. Jacques, and Q. Luo, Simultaneous compensation for spatial and temporal dispersion of acousto-optical deflectors for two-dimensional scanning with a single prism, Opt. Lett. 3(8), 09 093 (006). 4. Y. Kremer, J. F. Léger, R. Lapole, N. Honnorat, Y. Candela, S. Dieudonné, and L. Bourdieu, A spatiotemporally compensated acousto-optic scanner for two-photon microscopy providing large field of view, Opt. Express 6(4), 0066 0076 (008). 5. O. E. Martinez, Grating and prism compressors in the case of finite beam size, J. Opt. Soc. Am. B 3(7), 99 934 (986). 6. O. E. Martinez, Pulse distortions in tilted pulse schemes for ultrashort pulses, Opt. Commun. 59(3), 9 3 (986). 7. Z. L. Horváth, Z. Benkö, A. P. Kovács, H. A. Hazim, and Z. Bor, Propagation of femtosecond pulses through lenses, gratings, and slits, Opt. Eng. 3(0), 49 500 (993). 8. K. Varjú, A. P. Kovács, K. Osvay, and G. Kurdi, Angular dispersion of femtosecond pulses in a Gaussian beam, Opt. Lett. 7(), 034 036 (00). 9. J. C. Diels, and W. Rudolph, Ultrashort Laser Pulse Phenomena (Academic, San Diego, Calif., 996). 0. C. Fiorini, C. Sauteret, C. Rouyer, N. Blanchot, S. Seznec, and A. Migus, Temporal aberrations due to misalignments of a stretcher-compressor system and compensation, IEEE J. Quantum Electron. 30(7), 66 670 (994).. K. Osvay, A. P. Kovács, Z. Heiner, G. Kurdi, J. Klebniczki, and M. Csatári, Angular dispersion and temporal change of femtosecond pulses from misaligned pulse compressors, IEEE J. Sel. Top. Quantum Electron. 0(), 3 0 (004).. K. Osvay, A. P. Kovacs, G. Kurdi, Z. Heiner, M. Divall, J. Klebniczki, and I. E. Ferincz, Measurement of noncompensated angular dispersion and the subsequent temporal lengthening of femtosecond pulses in a CPA laser, Opt. Commun. 48(-3), 0 09 (005). 3. A. E. Siegman, Lasers, (University Science, Mill Valley, CA, 986). 4. X. Gu, S. Akturk, and R. Trebino, Spatial chirp in ultrafast optics, Opt. Commun. 4(4-6), 599 604 (004). 5. S. Akturk, M. Kimmel, P. O Shea, and R. Trebino, Measuring spatial chirp in ultrashort pulses using single-shot Frequency-Resolved Optical Gating, Opt. Express (), 68 78 (003). 6. S. Akturk, X. Gu, E. Zeek, and R. Trebino, Pulse-front tilt caused by spatial and temporal chirp, Opt. Express (9), 4399 440 (004). 7. S. Akturk, X. Gu, P. Gabolde, and R. Trebino, The general theory of first-order spatio-temporal distortions of Gaussian pulses and beams, Opt. Express 3(), 864 866 (005). 8. P. Gabolde, D. Lee, S. Akturk, and R. Trebino, Describing first-order spatio-temporal distortions in ultrashort pulses using normalized parameters, Opt. Express 5(), 4 5 (007).. Introduction A classical pulse compressor (commonly including four gratings or prisms, or a pair of gratings and prisms with a double-pass configuration) [,] can stretch or compress a femtosecond laser pulse by introducing variable amounts of positive or negative frequency chirp (also referred to as just the chirp ). Pulse compressors are ubiquitous in ultrafast optics because of their ability to tailor the pulse s temporal duration, which is essential for making and maintaining intense, short pulses. Important applications include chirped-pulse amplification (CPA) and material-dispersion compensation, the latter of which is necessary for generating ultrashort pulses [3 5]. Since the duration of the pulse after the compressor depends on its group-delay dispersion, it is vital to able to accurately calculate this quantity [6 8]. In some practical applications, the classical 4-dispersive element compressor is not always suitable, and instead, a single angular dispersion element [9 ] or a pair of angular dispersion elements (also just referred to as an element ) in a single-pass configuration [ 4] are used to for dispersion control. In these cases, spatiotemporal couplings (meaning that there are x-ω, or equivalently x-t cross terms in the field) such as angular dispersion, spatial chirp, and pulse front tilt are present in the output pulse, making it even more difficult to calculate the chirp that is added to the pulse by the compressor. To model compressors, either a plane wave [,,9 ], or a Gaussian beam [5 8,5 8]) model is usually used. Martinez first showed that these two models sometimes give very different results. Namely, the propagation dependence of the chirp of a pulse while propagating inside a classical pulse compressor is very different for the two types of beams. #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 707

Interestingly, however, after propagating through the 4th dispersive element (i.e., once all of the angular dispersion is removed), both models predict the same results [5]. Similarly it has been shown several times, that, when a pulse passes through a single angular disperser, the chirp of a Gaussian beam increases nonlinearly with propagation distance away from the disperser, while that of a plane wave increases linearly [8]. Though compressors are very commonly used, it seems that there is still more to learn about how they affect ultrashort pulses. Usually, the diffraction integral is used to investigate the propagation dependence of an ultrashort pulse, but this approach is quite complex making it difficult to use to understand how compressors work and why they effect plane waves differently than they do Gaussian beams [5]. Actually, the evolution of the chirp is entirely caused by a spatio-spectral phase φ(x,z,ω) that is added to the pulse by the compressor (since the chirp is defined as the second order derivative of the pulse s spectral phase with respect to angular frequency) [9]. Therefore, it is sufficient to determine the chirp s evolution by calculating this phase change that is acquired by the pulse due to propagation through the compressor. In this paper, we have derived expressions for the propagation dependence of the frequency chirp of a femtosecond Gaussian laser pulse after passing through each element of a classical four-element compressor. These expressions were derived by calculating the phase change acquired by the Gaussian beam after propagating through each element of the compressor. These calculations reveal the physical mechanism by which chirp is acquired by an ultrashort pulse in a pulse compressor. As the effects of misalignment of the dispersive elements are very complex [0 ], we consider only perfectly aligned compressors in this paper.. Phase acquired by a Gaussian pulse in a compressor A steadily propagating electromagnetic field in free space is governed by the scalar approximation of the Helmholtz equation. The plane wave is the simplest solution to this equation, and the Gaussian beam is also a special solution which is obtained if the slowly varying amplitude (SVA) approximation is made. The Gaussian beam is a very good model for realistic laser beams, so it is frequently adopted for modeling optics experiments. The phase of a Gaussian beam propagating in free space is given by [3]: φ kr z = + () R( z) zr ( r, z) kz tan ( ). where k is the wave-number, z is the propagation distance away from the beam waist on the axis, r = (x + y ) / is the distance from the z axis, R(z) = z + z R /z is the radius of curvature of the wave front, z R = kw o / is the Rayleigh range of the Gaussian beam, and w 0 is the beam waist size. Equation () describes the phase shift of a Gaussian beam at the point (r, z) relative to the original point (0, 0). In this equation, the first term is the geometrical phase shift which is the same as the phase of a plane wave. The second term represents a phase shift relative to the radial position. This is due to the finite beam size of the Gaussian beam. No such term exists in the expression for plane waves because they are infinite in space. The third term is the Gouy phase shift which is relative to the geometrical phase shift, and is also unique to a Gaussian beam [3]. #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 707

Fig.. (a) A classical pulse compressor, which consists of four identical prisms (or other angular dispersers, such as gratings). (b) Optical path diagram of a femtosecond laser pulse passing through a classical pulse compressor. For two adjacent elements, the spacings are L, L, L 3, respectively. Due to angular dispersion, in the system, different spectral components have different paths. Taking the entrance vertices (O ~O 4) of each element as the original point, four reference distances can be established: they are x -z through x 4-z 4. The distances for any spectral component are defined as x ω-z ω through x 4 ω-z 4 ω. In this figure, the dashed line between the second and third element is the wave front of the pulse. When the alignment is perfect, θ = θ 3 and L = L 3. BW is the beam waist, and d is the distance between the beam waist and the first dispersive element. As shown in Fig., in order to study the chirp evolution of a femtosecond Gaussian laser pulse propagating through a pulse compressor, we first need to define all of the system parameters. For a single angular dispersion element, the deflection angle of a femtosecond laser pulse after passing through the element is described by [5,6]: θ θ θ ( γ, ω) = α γ + β ω = γ + ω. γ ω where α is the angular magnification, β is the angular dispersion, γ is the angle of incidence, and ω is the frequency. For the four elements in a perfectly aligned pulse compressor, the parameters must obey the following equations [5]:, β, α α () β α = = (3) α = α, β = β, (4) 3 3 β α = = (5) 3 4, β4. α3 α3 where the subscripts on each parameter correspond to the number of the element. We will use α and β to denote these parameters for the first element. It is assumed that the angular dispersion only occurs in the x-z plane and the Gaussian beam still follows the propagation rules of free space in the y direction. In this case, we omit the information for the y axis in order to simplify our calculation. When an arbitrary spectral component ω of a femtosecond Gaussian laser pulse passes through the first element and arrives at the point (x ω, z ω ), the #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7073

waist position is at a distance d/α + z ω relative to the apex of the element for the opposite propagation direction for which the equivalent Rayleigh range is z R = z R /α, where d is the distance between the beam waist and the first dispersive element [8]. Thus the phase expression for the Gaussian beam can be written as: kx d / z φ (,, ω) ( / α ) tan ( ). ω α + ω x z = k d + z + ω ω ω R( d / α + z ) z ω R In the same way, when an arbitrary spectral component ω passes through the second, third and fourth element and arrives at the locations (x ω, z ω ), (x 3ω, z 3ω ) and (x 4ω, z 4ω ), the corresponding waist position is equivalent to d + α z ω + z ω, d/α + z ω + z ω /α + z 3ω, and d + α z ω + z ω + α z 3ω + z 4ω with respect to the apex of the first element for the opposite propagation direction, and the corresponding Rayleigh ranges are of z R = z R,z R3 = z R /α, and z R4 = z R. Then the corresponding phase expressions for the Gaussian beam are: (6) kx d + α z + z φ (,, ω) ( α ) tan ( ), ω ω ω x z = k d+ z + z + ω ω ω ω R( d + α z + z ) z ω ω R (7) 3ω 3( x3 ω, z3 ω, ) = k( d / + z ω + zω / + z3 ω ) + R d α + z ω + zω α + z3 ω φ ω α α d / α + z ω + zω / α + z3 ω tan ( ), z R3 kx ( / / ) (8) kx4ω φ4 ( x4 ω, z4ω, ω) = k( d+ α z ω + zω + α z3 ω + z4 ω ) + R( d+ α z ω + zω + α z3 ω + z4ω ) (9) d+ α z ω + zω + α z3 ω + z4 ω tan ( ). z R4 where the subscripts of the phase functions represent the number of elements. From the above analysis, we obtain the spatio-spectral phase of the Gaussian beam at any position when passing through a pulse compressor. The corresponding frequency chirp can be obtained by calculating the second order derivative of the phase with respect to the frequency ω. 3. Chirp evolution of a Gaussian pulse in a compressor As shown in Fig., after passing through the first angular dispersion element, each spectral component of the femtosecond laser pulse is separated in space, and propagates a different distance. The reference spectral component propagates a distance z, while an arbitrary spectral component ω propagates a distance z ω, which can be related to z using the angle θ. The relevant formulae are [7]: z ω = z cosθ+ x sin θ, (0) x ω = z sinθ+ x cos θ. () And the corresponding first and second order derivatives are: ' z ω = 0, () '' dθ ω z, z = dω (3) x dθ = z (4) d ω ' ω, #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7074

d θ = z (5) d ω ''. ω x where β = dθ /dω is the angular dispersion introduced by the element. The second order derivative of the expression forφ with respect to spectral frequency ω (which describes the introduced chirp of the Gaussian beam passing through the first element) can then be obtained as follows: kβ z α z ( d+ α z ) φ ( ) β β [ ]. '' G z = k z + = k z R( d / α + z ) ( d+ α z) + zr Comparing the above expression with the phase function φ (Eq. (6)), we can see that the first term is the chirp introduced by the geometrical phase shift, which is negative. The second term is the chirp introduced by the radially dependent phase term which has an x-ω coupling, and this is positive. The chirp introduced by the Guoy phase shift is usually far less than the first two terms, so for simplicity we have left it out of Eq. (6) and will neglect it for the rest of our discussion. After the pulse has passed through the second element (which is anti-parallel with the first element and separated by a distance L ), the angular dispersion is totally eliminated but different spectral component of the pulse are still separated transversely in space, or spatial chirp, and also the x- ω coupling term in the phase (wave front tilt dispersion) are present [4 8]. In this interval, the relevant distances are: (6) z = z + ω const, (7) x = x + ξω+ const (8) ω. where ξ = αβl is the spatial chirp present in the pulse before entering the second element. Then the corresponding first and second derivatives of z ω and x ω with respect to ω can be obtained as follows: x ' z ω = 0, (9) '' z ω = 0, (0) = αβ () ' L, ω '' x ω = () 0. And thus the chirp of the pulse after passing through the second element is: '' α L ( d + α L + z ) φg ( z ) = kβ L + kα β L = kβ L [ ]. (3) R( d + α L + z ) ( d + α L + z ) + z R Comparing Eq. (3) with the phase function φ (Eq. (7)), we can see that the first term is the chirp introduced by the geometrical phase shift, which is negative. As the second element removes angular dispersion, no extra chirp will be introduced by the geometrical phase shift after passing through the second element. The second term represents the frequency chirp introduced by the radially dependent phase when the pulse travels the distance z, which is positive. Here we can see that though the angular dispersion has been eliminated after passing through the second element, x-ω couplings, namely, spatial chirp and wave-front-curvature dispersion are still present, and these combined with a changing beam spot size introduce additional frequency chirp as the pulse propagates through this region of the compressor. After the pulse has passed through the third element, angular dispersion is present again. As shown in Fig., the distance for any spectral component between the third and the fourth elements is L 3 cosθ 3 [], and thus the relevant formulae are: #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7075

z3 ω = z3 cosθ3 x3 sin θ3, (4) x3 ω = ( L3 z3)sinθ3+ x3 cos θ3. (5) Then the first and second derivatives of z 3ω and x 3ω with respect to frequency ω can be obtained as follows: ' z3 ω = 0, (6) z '' dθ3 3ω z3, = dω (7) ' dθ3 x3 ω = ( L3 z3), (8) dω '' d θ3 x ω = ( L3 z3). (9) dω where the angular dispersion introduced by the third element is dθ 3 /dω = -β. Also, the chirp of the femtosecond laser pulse after passing through the third element is: '' ( ) = + ( ) 3 3 3 3 φ z k G β L k β z k β L z R( d / a + L + L / α + z ) ( ) 3 α ( d + α L + L + α z 3) ( d + α L + L + α z 3) + z 3 3 3 R = kβ L kβ z + kβ L z Comparing the above expression with the phase function φ 3 in Eq. (8), we see that the first term is the chirp from the geometrical phase term acquired from propagating between the first and the second elements, which is negative. The second term is the chirp from the geometrical phase shift due to propagation between the third and the fourth elements, which is also negative. Finally, the third term is the chirp introduced by the radially dependent phase shift due to propagation through the same distance, which is positive. If the propagation distance z 3 = L 3, i.e. when the pulse arrives at the entrance of the fourth element, the chirp from the last term becomes zero. Provided that z 3 = L 3, and the angular dispersion from the first and third elements are equal and opposite, after the pulse has passed through the fourth element, all of the spectral components overlap transversely in space (i.e. no spatiotemporal couplings are present), and both the angular dispersion and all other spatiotemporal couplings have been removed, giving:. (30) z = z, (3) 4ω 4 x = x. (3) 4ω 4 So, after passing through the fourth element, neither the geometrical phase shift nor the radially dependent phase shifts introduce any chirp and then we get [5]: φ '' ( z G 4 ) = kβ L kβ L 3 = kβ L. (33) In this formula we can see that the final frequency chirp of the Gaussian beam is the same as that of a plane wave after passing through the pulse compressor. This final expression depends only on the geometrical phase shift and has no dependence on the radially dependent phase shift. From the above analysis of the chirp of the Gaussian beam, we can easily obtain the corresponding chirp of the plane wave by considering only the geometrical phase term. As the term describing the geometrical phase shift of the Gaussian beam is equal to the phase term of the plane wave, the chirp introduced by this term is the corresponding chirp of the plane wave: #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7076

φ ( z P ) = kβ z, (34) '' φ ( z P ) = kβ L, (35) '' φ ( z P ) = kβ L kβ z, (36) '' 3 3 φ '' ( z P 4 ) = kβ L. (37) From the above analysis, we can see that, for the Gaussian beam, the chirp introduced by the geometrical phase term is negative while that introduced by the radially dependent phase shift is positive at any position after passing through the element; it becomes zero when no ω- x coupling terms are present in the phase. Also, the radially dependent phase shift in the phase function of the Gaussian beam is the fundamental reason for the difference between the frequency chirp of a Gaussian beam and plane wave when passing through the classical pulse compressor. The frequency chirp evolution of the Gaussian beam and plane wave over the whole propagation process are shown in Fig., using the parameters α =, β = 0. rad/µm, d = z R = L = L = L 3 = m. Fig.. Comparison of the chirp evolution of a Gaussian beam (green) and plane wave (black) when passing through a 4-dispersive element pulse compressor. The red arrows in the figure show the difference in the two models. As shown in Fig., though the frequency chirp evolution of Gaussian beams is different from that of plane waves, the final result is the same for both models. In the following section, we provide a detailed analysis of the physical mechanism of this phenomenon. 4. Comparison of the Gaussian beam and plane wave models for a pulse compressor Here we compare the propagation dependent frequency chirps predicted by the two models in order to better understand their differences. For the plane wave model, in each interval of propagation, the following expressions are obtained: '' φ ( ) = β, (38) P z k z φ '' z P ( ) = 0, (39) φ P ( z ) = kβ z, (40) '' 3 3 #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7077

( ) = 0. (4) φ '' z P 4 where the z s are the corresponding propagation intervals as shown in Fig.. We can see that the frequency chirp of the plane waves only changes in the interval where the angular dispersion is nonzero (e.g. after the pulse has passed through the first and third elements, as shown in Fig. ). The corresponding chirp evolution of the Gaussian beam in each interval is: '' α z ( d + α z ) φ ( z ) = k z [ ], G β ( d + α z ) + z R (4) φ '' ( d + α L + z ) ( d + α L ) G α β ( z ) = k L [ ] ( d + α L + z ) + z ( d + α L ) + z R = [ ], ( ) ( ) kα β L R d + α L + z R d + α L R (43) '' α ( d+ α L + L + α z ) α L ( d+ α L + L ) φg ( z ) = k z + k ( L z ) k L ( ) ( ) 3 β β β 3 3 3 3 d+ α L + L + α z + z d+ α L + L + z 3 R R L z = kβ z kα β L [ ( ) ], ( ) ( ) 3 3 3 R d+ α L + L L R d+ α L + L + α z 3 3 (44) φ '' G ( z ) = 0. (45) In contrast to the plane wave model, for Gaussian beams, the chirp changes as the pulse propagates even when the angular dispersion is zero. When the angular dispersion is zero but the x-ω coupling in the phase is nonzero, the frequency chirp of the Gaussian beam still changes along the propagation distance (as shown in Fig., after the pulse has passed through the second element). After the pulse has passed through the first and third elements (where the angular dispersion is not zero), the changes in the frequency chirp consist of two parts; one is introduced by the geometrical phase shift (corresponding to the first term in Eqs. (4) and (44)) and the other is due to the radially dependent phase shift (corresponding to the second term in Eqs. (4) and (44)). Calculating the difference in the chirp predicted by the two models over corresponding intervals (Eq. (4) is substracted from Eq. (38), and, etc.), we get the following results: 4 φ φ = '' PG ( z ) k α β z R d + α z ( ), '' ( z ) k L φ = α β, PG R ( d L z α ) R ( d α L + + + ) ( ), R( d + α L + L ) L R( d + α L + L + α z ) 3 3 '' L z 3 3 PG z = kα β L 3 (46) (47) (48) φ '' PG ( z ) = 0. (49) 4 From the above equations we can clearly see that over any propagation interval, the fundamental reason for the difference between frequency chirp of the Gaussian beam and plane wave is due to the radially dependent, or the x-ω coupling term in the phase of the Gaussian beam. As shown in Fig., (A) In the interval 0-L, the chirp of the Gaussian beam increases less rapidly than that of the plane wave because the radially dependent phase shift introduces a positive change in the chirp (see Eq. (46)). (B) In the interval L to L, there is no further #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7078

change of the chirp of the plane wave while the chirp of Gaussian beam continues to change, due to the spatiotemporal couplings in the radially dependent term. However, these changes can be either positive or negative, as determined by the parameters of the beam itself (such as the Rayleigh length z R ) and the propagation distances (such as d, L, L ) (see Eq. (47)). (C) In the interval L to L 3, the chirp of the Gaussian beam increases faster than that of the plane wave, because in this interval, different spectral components of the pulse show a tendency to converge, that is to say, the spatial chirp is becoming smaller, which is the opposite of what happens between 0 to L (where the beam is diverging and the spatial chirp is increasing), so the chirp introduced by radially dependent phase shift becomes negative (see Eq. (48)). When the propagation distance z 3 = L 3 (at the entrance of the fourth element), the changes in the difference of the chirp between the Gaussian beam and the plane wave in this interval offsets their difference in the intervals 0 to L and L to L exactly, leading to the final chirp of Gaussian beam and plane wave being the same. Actually, Eqs. (46) ~(48) also reveal this rule: '' '' '' φpg ( L ) + φ ( ) ( ) 0. PG L + φ PG L = (50) 3 Thus, for both the Gaussian beam and the plane wave, the final chirps are the same. The fundamental reason for this difference is that the Gaussian beam has an x-ω coupling (or a spatiotemporal coupling) in its phase which introduces some chirp, and this term is absent in the phase of plane wave. If the spot size of the Gaussian beam increases, or if the Rayleigh range is longer and the beam is better collimated, the propagation dependence of the chirp predicted by the Gaussian beam is closer to that of the plane wave, as shown in Fig. 3. Fig. 3. The chirp evolution of a femtosecond laser pulse when passing through the pulse compressor. The Gaussian beam (green line) compared with a plane wave (black line). Each of the subfigures has different Rayleigh range: from (a) to (d), they are m, 3m, 5m, and 0m, respectively. The other parameters are the same as those of Fig.. When the Rayleigh range increases, the chirp evolution of the Gaussian beam becomes closer to that of a plane wave. This phenomenon is also shown in Eqs. (46) ~(49) where you can see that if the Gaussian beam is well collimated, that is to say, d, L, L, L 3 << z R,(which is called the approximate '' condition for good collimation [5,6]), = 0. Also, we can see that the chirp changes of φ PG #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7079

the Gaussian beam can be either positive (see subfigures b and c) or negative (see figure a) in the interval L to L, and depends on the parameters of the Gaussian beam and the propagation distance. 5. Discussion Using the Kirchhoff-Fresnel diffraction integral, we can obtain the complete electric field of the pulse while propagating through compressor which will tell us not only the frequency chirp, but also the spatial chirp, pulse front tilt, spot size, and etc [6 8,5,6]. However, for studying the evolution of the chirp, an analysis of the whole electric field, is not necessary, as the chirp is directly determined by the phase shift of the laser pulse. In this paper, we used a much simpler and more straightforward method, by analyzing the phase acquired by a Gaussian pulse due to propagation through a compressor, and directly calculated the second order derivative of the phase with respect to the spectral frequency ω. This gives us the chirp evolution of the pulse when passing through a pulse compressor. We found the x-ω coupling term in the phase of the Gaussian beam to be the fundamental reason for its chirp evolution inside a pulse compressor being different from that of a plane wave. When the spot size of the Gaussian beam increases, or when it is well collimated and has a longer Rayleigh range, the chirp acquired by the Gaussian beam is much closer to that of a plane wave, as shown in Fig. 3. The reason for this is that the radially dependent phase term of the Gaussian beam becomes less significant (i.e. the curvature of the wave front increases). If the Rayleigh range is long enough, no radially dependent phase term exists, leading to equivalent chirps for the Gaussian beam and the plane wave models. Although the spot size of the Gaussian beam has no influence on the final chirp after propagation through the 4th dispersive element (assuming a well aligned compressor), it could significantly influence the value of the chirp when the pulse is passing through a single angular dispersion element, a pair of elements in a single-pass structure, or when the pulse compressor is not perfectly aligned. This of course affects the duration of the output laser pulse. Another phenomenon worth noting is that for Gaussian beams, as long as spatiotemporal couplings exist (i.e. x-ω cross terms), even though there is no angular dispersion, an additional frequency chirp will be introduced as the pulse propagates. Martinez implied such a phenomenon but provided no explanation [5]. In essence, propagation of an ultrashort pulse in the presence of spatiotemporal couplings is a three-dimensional (x, ω and z) effect, and it causes the spatial terms of the Gaussian to mix with the frequency terms, and phase terms are transferred to the intensity (and vice versa). Namely, the first dispersive element introduces angular dispersion which propagation changes into spatial chirp, wave-front-curvature dispersion (or pulse front tilt, if viewed in the time domain), frequency chirp, and others. And once the angular dispersion vanishes after the second prism, there are still coupling terms remaining the spatial chirp and the wave-front-tilt dispersion so propagation again transfers these terms into frequency chirp. Specificially it is the radially dependent phase term in the Gaussian model that allows for this mixing, and if this term vanishes, then the frequency chirp of Gaussian beam will be the same as that of plane wave (such as after the pulse has passed through the fourth element, which also makes this term vanish). The spatial chirp and frequency chirp are commonly regarded as independent parameters; the spatial chirp describes the transverse separation (perpendicular to the propagation direction) of different spectral components while the frequency chirp describes the longitudinal delay (along the propagation direction). Here we again illustrate that this is not the case [6 8], and we have shown that these two quantities can be coupled by propagation. 6. Conclusion In this paper, we studied the chirp evolution of a Gaussian beam when passing through a classical pulse compressor by directly calculating the acquired spatio-spectral phase. Compared with the chirp evolution predicted by the plane wave model, we found that a spatiotemporal coupling or an x-ω dependent term in the phase of the Gaussian beam is the fundamental reason for the difference in these two models predictions for pulse compressors. #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 7080

For a Gaussian beam, the existence of spatiotemporal couplings also introduces an additional frequency chirp even when no angular dispersion exists. If the Gaussian beam is well collimated, the frequency chirp evolution is closer to that of a plane wave after passing through the angular dispersion elements. This work provides a deeper understanding of the physical mechanism of the frequency chirp evolution when a laser pulse passes through a classical pulse compressor. Our analysis will also be helpful for optimization of dispersion compensation schemes in many applications of femtosecond laser pulses. Acknowledgments This work was supported by the National Natural Science Foundation (NSFC) (3090033, 309700), and Program for Changjiang Scholars and Innovative Research Team in University. #0 - $5.00 USD Received Jun 009; revised 4 Sep 009; accepted 5 Sep 009; published 0 Sep 009 (C) 009 OSA 4 September 009 / Vol. 7, No. 9 / OPTICS EXPRESS 708