Κ Is Not Satisfied
Total Page:16
File Type:pdf, Size:1020Kb
A Coherent Technical Note August 29, 2018 Propagation, Dispersion and Measurement of sub-10 fs Pulses Table of Contents 1. Theory 2. Pulse propagation through various materials o Calculating the index of refraction . Glass materials . Air . Index of refraction in the 600 -1100 nm range . Example: Broadening of a 7 fs pulse through different materials o Group delay, group velocity and group velocity dispersion . Group velocity . Group velocity dispersion . Group Delay for various materials . Example: Output pulse duration as a function of GDD o Pulse broadening and distortion due to TOD . Example: output pulse duration as a function of TOD . Example: Output pulse shape for SF10 . TOD dispersion curve for various materials o Pulse broadening and distortion due to FOD and higher o Summary of pulse broadening effects in various materials 3. Choosing the proper optics o Some practical tips 4. Appendix A: Pulse duration measurement 5. Appendix B: Refractive index for various materials referenced in this work Any ultrafast laser pulse is fully defined by its intensity Depending on their exact nature, phase distortions may and phase, either in time or frequency domain. broaden the pulse and modify its shape such that the Propagation in any media, inclusive of air, results in pulse is not “transform- limited” anymore. distortions of phase or amplitude. Large distortions in phase can be introduced by propagation of the beam This means that the time-bandwidth relationship even through optical elements with very low absorption ∆t∆ν ≤ κ is not satisfied. Here κ depends on the such as lenses or prisms. This effect is frequently shape of the spectrum (κ =0.441 for a Gaussian pulse, called temporal chirp and is due to chromatic (i.e. κ =0.315 for a hyperbolic secant pulse and κ =0.886 wavelength-dependent) dispersion. for a square pulse). 1 | www.coherent.com | [email protected] | 800-527-3786 | 408-764-4983 A Coherent Technical Note August 29, 2018 1. Theory With χ the dielectric permissibility Here for simplicity we will assume pulses with a By combining the two previous equations and moving temporal Gaussian shape. While this is only a into the frequency domain by Fourier transform we get: convenient approximation it still gives a good idea of what is happening. ( ) ( ) 2 + 2 2 , = 0 The electric field is defined by the expression ϖ � 2 2 2 ϖ � ϖ With ( ) = 1 + ( ) the dielectric constant ( ) ( ) = 2 We assume the susceptibility and dielectric constant − 1+ � � ϖ ϖ − here are real (i.e. there is no absorption). With = , where is the center wavelength (in 2 A general solution for the previous equation is our case 8000 nm). 0 a is defined as the chirp parameter of the pulse and we ( , ) = ( , 0) ( ) will start without chirp so that a=0. − ϖ With defined as = 2ln (2) , being the FWHM ϖ ϖ Where the propagation constant k is determined by the pulse duration. � linear optics dispersion relation The expression for the pulse then becomes ( ) ( ) ( ) = 2 = 2 ( ) 2 = 2 ϖ ϖ 2 �−0−� � � 2 2 ϖ ϖ ϖ n being the refractive index. A wave equation can be derived for the electric field E from the Maxwell equations (in absence of external The main reason to work in the frequency domain is charges and currents, and considering only non- that phases are additive, unlike in the time domain. magnetic permeability and uniform medium). 2. Pulse Propagation Through Various In Cartesian coordinates the equation is Materials 1 ( ) Let’s continue to assume a Gaussian pulse initially 2 + 2 + 2 2 , , , chirp-free (transform limited) with a 800 nm center � 2 2 2 − 2 2� = ( , , , ) wavelength. To calculate the pulse envelope we need to Fourier- transform the electric field into the frequency 0 2 domain, add the frequency-dependent phase, Fourier- With the magnetic permeability of free space transform back in the time domain and calculate the 0 norm of E(t,z). To summarize, the propagated field in µ = + The Polarization P contains two terms, , the frequency domain is given by i.e. the linear and non-linear term respectively. Assuming very weakly focused pulses (with no ( ) substantial changes in the x and y direction) we can ( , ) = ( , 0) neglect . − ϖ ϖ ϖ We then go back in the time domain in order to We then obtain the reduced wave equation calculate the envelope and duration of the output pulse. 1 ( ) ( ) 2 2 , = , � 2 − 2 2� 0 2 From classic electrodynamics we know that ( , ) = ( ) (z, ) 0 ϖ ϖ 2 | www.coherent.com | [email protected] | 800-527-3786 | 408-764-4983 A Coherent Technical Note August 29, 2018 As a first step we need to determine the index of chirp-free 7 fs input pulse. This is representative of the refraction of the different materials we will consider. typical pulse produced by the Coherent Vitara UBB Then we calculate the pulse duration after propagation Ti:Sapphire laser. through 5 mm and 10 mm of material starting with a Vitara UBB <10 fs 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 spectrum intensity 0.2 0.1 0 620 640 660 680 700 720 740 760 780 800 820 840 860 880 900 920 940 960 980 1000 1020 wavelength in nm Typical spectral amplitude of pulse from Vitara UBB (Transform-limited pulse duration = 7 fs) We are going to consider 6 media, including the most Calculating the index of refraction common types of glass used for laser optics, a very low . Glass materials dispersion material (CaF2), a high dispersion material (SF10) and finally air. For all glass materials, we use the well-known - Fused Silica Sellmeier equation: - Sapphire - ( ) CaF2 = 1 + 2 + 2 + 2 1 2 3 - SF10 2 2 2 2 − 1 − 2 − 3 - BK7 with λ in μm. The B and C coefficients are given in the - Air table below. Material 1 B1 2 B2 3 B3 4 C1 5 C2 6 C3 Fused Silica 0.6961663 0.4079426 0.8974794 0.00467914826 0.0135120631 97.9340025 Sapphire 1.43134930 0.650547130 5.34140210 0.00527992610 0.0142382647 325.017834 CaF2 0.5675888 0.4710914 3.8484723 0.00252642999 0.0100783328 1200.555973 SF10 1.62153902 0.256287842 1.64447552 0.0122241457 0.0595736775 147.468793 BK7 1.03961212 0.231792344 1.01046945 0.00600069867 0.0200179144 103.560653 3 | www.coherent.com | [email protected] | 800-527-3786 | 408-764-4983 A Coherent Technical Note August 29, 2018 Air (Edlen Model) An accurate calculation of the refractive index of air The ten constants for air are given in the table below. requires a set of ten constants, in addition to its dependence on temperature, pressure and humidity - Temperature - Pressure - Humidity K1 K2 K3 K4 K5 1.167052145280E+03 -7.242131670320E+03 -1.707384694010E+01 1.202082470250E+04 -3.232555032230E+06 K6 K7 K8 K9 K10 1.491510861350E+01 -4.823265736160E+03 4.051134054210E+05 -2.238555575678E-01 6.501753484448E+02 There are also five coefficients that depend on the temperature (expressed in K) Ωa Ca Aa Ba Xa = + = + + A = + + B = K + + = + 4 9 2 2 2 2 6 7 8 a 1 2 a 3 4 5 Ω 10 Ω Ω Ω Ω Ω Ω − � − − We also need to include: Where RH is the relative humidity and the pressure is - Pressure vapor saturation : expressed in Pascal (101325 Pa = 1 atm). = 10 (2 ) - 6 4 Finally, there are seven additional constants : - Partial vapor humidity : = ( ) 100 Ab Bb Cb Db Eb Fb Gb 8342.54 2406147 15998 96095.43 0.601 0.00972 0.003661 - We can then calculate the coefficient Xb that we can use to calculate the full refractive index of air (1 + 10 ( 273.15) = 1−8+ ( 273.15) � − − � ( ) = 1 + ( ( ) 1) Armed with these parameters, − we calculate an 292 .75 10− 3.7345 expression for ns 0.0401−10 − � � � ( ) = 1 + + + 10 1 1 38.9 130 −8 − 2 � �� � � �� − 2 − 2 4 | www.coherent.com | [email protected] | 800-527-3786 | 408-764-4983 A Coherent Technical Note August 29, 2018 Index of Refraction in the 600 nm to 1100 nm Range By applying the previous formula for air and the different materials, we get the index of refraction curve from 600 nm to 1100 nm. Index of refraction Fused Silica Index of refraction Sapphire 1.464 1.775 1.46 1.771 1.456 1.767 1.452 1.763 1.448 1.759 1.444 index of refraction 1.755 Index of refraction 1.44 1.751 1.436 1.747 600 650 700 750 800 850 900 950 1000 1050 1100 600 650 700 750 800 850 900 950 1000 1050 1100 Wavelength in nm Wavelength in nm 1.728 Index refraction SF10 1.524 Index refraction BK7 1.724 1.52 1.72 1.516 1.716 1.512 1.712 1.508 1.708 1.504 Index of refraction Index of refraction 1.704 1.5 1.7 1.496 600 650 700 750 800 850 900 950 1000 1050 1100 600 650 700 750 800 850 900 950 1000 1050 1100 Wavelength in nm Wavelength in nm Index of refraction CaF2 Index refraction air 1.014 1.444 1.44 1.01 1.436 1.006 1.432 1.002 1.428 0.998 1.424 0.994 index of refraction index of refraction 1.42 0.99 1.416 0.986 600 650 700 750 800 850 900 950 1000 1050 1100 600 650 700 750 800 850 900 950 1000 1050 1100 Wavelength in nm Wavelength in nm 5 | www.coherent.com | [email protected] | 800-527-3786 | 408-764-4983 A Coherent Technical Note August 29, 2018 Index refraction air 1.000268 1.0002675 1.000267 1.0002665 Refractive index of air on a much expanded scale shows its wavelength dependence for given 1.000266 temperature, pressure and humidity index of refraction 1.0002655 1.000265 1.0002645 600 650 700 750 800 850 900 950 1000 1050 1100 Wavelength in nm The refractive indices of all materials are generally The degree of broadening depends both on the change wavelength-dependent which means that each of the refractive index vs.