## 1. Introduction

Dispersion compensation plays a key role in generation, amplification and propagation of femetosecond pulses. In the dispersive medium, the pulse can be broaden or compressed depending on the sign of chirp and dispersion. To obtain the ultrashort pulses, the pulse group delay (GD) should have about frequency independence after the dispersion compensation. Especially to compress a pulse to near the transform limit one should not only compensate the GD but also eliminate the high order dispersion term. In this chapter, we briefly review the calculation formulae and design rules of the various dispersive devices including grating and prism pairs and chirped mirrors. An integrative optimization of chirped mirror and its function in mode locked Ti:sapphire laser is described in details. Finally, a useful tool ─ dispersion map is introduced.

## 2. Grating pairs with negative and positive dispersion

A parallel grating pair as Treacy demonstrated (Fig. 1) can offer a negative group delay dispersion (GDD) (Treacy, 1969). The commonly used phase shift formula through the parallel grating pair is

where *b* and *G* are the slant and the perpendicular distance path length between the two gratings respectively, *d* the grating constant, *c* the speed of light in vacuum, and

For a grating pair with the positive dispersion, Zhang et al (Zhang et al. 1997) proposed a grating- mirror model in contrast to the grating-lens model of Martinez (Martinez, 1987), to demonstrate that the phase shift of a grating-mirror system is a conjugation of that of parallel grating pair. In the model, a mirror with the radius of curvature of *R* is placed so that the incident point on the grating is right at the circular center of the mirror, as shown in Fig. 2. The simple ray-tracing would result a phase shift as

The subsequent differentiation with respect to frequency results in GD, GDD and the third order dispersion (TOD). Obviously, the only difference between the phase shifts of those two configurations is a constant (4R, depends on the starting point of calculation) and a opposite sign before the term of *b*(1+co*s*). The differentiation would vanish the constant term and leaves an opposite sign. This proves that the grating-mirror configuration is a phase conjugation of the parallel grating pair and can be used as a pulse “stretcher”.

However, this grating-mirror stretcher cannot be realized because, in practice, the incident beam must have a finite size, while the single-mirror system gives no beam collimation. Therefore a telescope is required for the beam collimation.

To derive the similar formula to Equation for the telescope stretcher, we would prefer to start with the folded-grating-telescope stretcher and to adapt it in the coordinate shown in Fig.3. For explicitness, only the incident and the output rays are displayed, the folding mirror and some optical paths have been omitted.

where

All variables are defined in reference (Zhang et al. 1997).It is noticed that, although the formula is complicated, it takes the form of:

indicating that it is indeed the phase conjugation of the pulse compressor, except for the mirror induced stigmatism. Similarly, we can have the formulism for Öffner type stretcher, which is (Jiang et al. 2002)

The only difference with Martinez stretcher is the subscript of the angles. All variables are defined in reference (Jiang et al. 2002).

The above formulas, although they are two dimensional, are very useful in optimization of chirped pulse amplification systems.

## 3. Prism pairs

Prism pairs have also been widely used for dispersion control inside laser oscillators since they can be very low loss in contrast to grating pairs. In 1984 Fork et al. derived formulism of GDD for pair of Brewster-angled prisms (Fork et al. 1984). Nakanuma and Mori noticed that the commonly used formulas developed by Fork were not convenient to use in the case where the adjustable prism insertion is required. They proposed to separate the prism separation and prism insertion (Nakanuma & Mori, 1991) so that those parameters can be easily measured. The simplified formulism is found in reference (Zhang & Yagi, 1993). The GDD and TOD can be expressed to functions of prism separation and prism insertion:

where the coefficients are approximately

where n is the refractive index of the prism material. An important advantage of using the above formulism is that they can be presented in a so called dispersion map, as described in section 5.

## 4. Chirped mirrors

Although grating and prism pairs are widely used as the compressor schemes, they suffer from high order dispersion. By contrast, chirped mirrors (CMs) which are dispersive multilayer structures designed by optimizing the initial design are usually used for pulses compression by providing dispersion without much material in the beam path. They can provide the precisely controlled GD, GDD, even higher order dispersion such as third order dispersion (TOD), as well as the high reflectivity in a broadband or ultra-broadband wavelength range. They are key components for extremely short pulse compression and especially have played a critical role in the development of octave-spanning lasers and the generation of mono-cycle pulses. In 2007, 1.5-octave chirped mirrors have been proposed and demonstrated, for extra-cavity compression to 3 fs pulses (Pervak et al. 2007). Furthermore, they are important for compression of attosecond soft-x-ray pulses (Wonisch et al. 2006).

There are always ripples in GDD curves of broadband chirped mirrors, due to the interference between the reflections from different layers of a multilayer structure. Although the ripples can be minimized by the optimization procedure such as with commercially available software, it is difficult to reduce the ripples to acceptable level when the design of a single chirped mirror has to operate over the wavelength ranges up to an octave spanning or more. Therefore, several other approaches to reduce GDD ripples have been proposed, such as double-chirped mirrors (Kärtner et al. 1997; Matuschek et al. 1998; Kärtner et al. 2001), complementary CM pairs (Pervak et al. 2005; Pervak et al. 2007). By pairing the CMs, the dispersion ripples of CMs are opposite to each other. Thus, the residual ripple of the CM pairs could be very small.

The design rule of CMs consists of two main issues. The first issue is the initial multilayer design. The design procedure of chirped mirrors is usually to start with a favorable initial structure, then to perform optimization. In this procedure, the initial design is very crucial. If the initial design is not close enough to the target, the optimization procedure is impossible to reach a satisfactory result. The commonly used initial designs are double chirp (DC) and modulated layer thickness (MLT) (Tempea et al. 1998).

The second issue is the optimization process including the choices of targets and optimization methods. Although a CM is characterized by a certain value of GDD that is the second derivative of the phase shift on reflection with respect to the angular frequency, GD was mostly chosen as the target, because the second derivative of the phase has two problems: the accuracy and the time consumption in the optimization. However, we found that just taking GD as the target cannot ensure a smooth GDD. Besides the targets setting, we should also consider the optimization methods. Although paring two CMs is a good idea to eliminate the ripples but in some cases such as generating octave-spanning in four-mirror Ti:Sapphire ring cavity it is not enough. Then we introduce the intracavity optimization.

### 4.1. Initial design of chirped mirrors

Here, we summarize two different designs as the initial layer structure for the optimization and compare the differences. First was the DC, and the second was the MLT.

The conception of DC mirrors has been proposed by Käertner (Kärtner et al. 1997; Matuschek et al. 1998) where in addition to the chirp of the Bragg wave number, the thickness of the high refractive index layers is also chirped.

Assuming the double chirped mirrors composed of 2N index steps with high refractive index medium such as TiO_{2} and low refractive index medium such as SiO_{2}. The thicknesses of low index layer are linearly chirped while the thicknesses of high index layer are given by

Where η=1, 2 for linear and quadratic chirp. The clear reduction of oscillations in the GD by the double chirped technique is visible (Matuschek et al. 1998). Moreover, the index-matching can be extended by a properly high quality AR-coating.

MLT method (Tempea et al. 1998) is based on the recognition that in the two component structure, the effects of modulating the layer thickness is similar to that of a corresponding modulation of refractive index. Composing the mirror with narrow bandwidth filters with gradually shifted central wavelength may avoid some resonant interference. The general formula described the variation of the layer thickness is given by

where x is the distance of the respective layer from the substrate and the number of layers within a modulation period Λ must be kept constant,

where the parameter

For an early constant negative GDD in the high reflectance range we need

where d is the total optical thickness of the coating. One problem in CM design is to obtain negative GDD at the shortest wavelength of the high-reflectance range. This difficulty is related to that the shortest wavelength components have a limited penetration depth which tends to result in a constant GD. This problem can be alleviated by increasing the amplitude of the modulation depth by

where A_{1}, A_{2} (>A_{1}) are limits between which the amplitude is allowed to vary.

The above equations consist of a complete recipe of the design of chirped mirror by using MLT method.

We simulated two chirped mirror pairs to compare these two above initial designs—DC (pair 1) and MLT (pair 2) (Chen et al. 2007). The target was CM pairs providing negative GDD of around -50 fs^{2} and TOD of around -35 fs^{3} at 800 nm with the reasonable residual GDD ripples and high reflectivity in wavelength range of 600-1200 nm. In the DC model, we chose the first order. In the MLT model, the number of layer within a modulation period was 5, the minimum modulation amplitude was 12.5 nm and maximum modulation amplitude was 62.5 nm, same as in the reference (Tempea et al. 1998). The total layer number of each mirror in a pair for both models was 80.

After optimization, from Fig. 6 and 7, we can see that both of the CM pairs meet the target at 800 nm with the residual GDD ripples <15 fs^{2} and the reasonable high reflectivity. Meanwhile, we find that the time consuming on optimization for pair 1 is little longer than pair 2, the GDD curve is closer to what we required and the GD is matched better in the spectral range. Particularly, the reflectivity is relatively higher.

It is noticed that, from Fig. 8 we can see that the reflectivity of pair 1 is above 99.4% in most of spectral range except the 1170~1200 nm. However, reflectivity of pair 2 is unexpectedly low (< 99%) in the range of 810~910 nm. If these CM pairs are going to be used inside laser cavity, the loss of mirrors is one of the most important factors we should consider. In this regard, pair 1 is better than pair 2.

The above simulation is just one example to design the CM pairs. We can see each initial design has its own advantages and disadvantages; thus, the choice mainly depends on the specific necessities.

### 4.2. Optimization of chirped mirrors

#### 4.2.1. Optimization targets

It is noticed that in some cases, the two CMs can be better paired with each other with GD as the target than with GDD as the targets (Pervak et al. 2007). For complementary pair mirrors, the merit function is modified as below when we take GD as the target:

where the subscripts j = 1, 2 represents each individual mirror. R_{t} (i) and GD_{t}(i) are the desired values at the wavelength_{j}(i) and GD_{j}(i) are the calculated instantaneous values. R and GD are the corresponding tolerances and n is the number of wavelength points. The weighting parameters α, β and , have to be adjusted according to the error in the process to handle the optimization.

We should also notice that CMs are finally characterized by the GDD. Therefore, in the final optimization step, GDD can be used to be the target. The merit function thus takes the following form in this case:

where the variables are the same as in Eq. (1), except GDD in the place of GD.

Considering these two aspects, we find that the optimization process can be divided into three steps. In the first and second optimization, the GD and GDD are set to be the optimization target respectively, as described in Eq. (1) and Eq. (2), where the derivative of the phase is calculated by the Birge and Kaertner method (Birge et al. 2006), since this method reduces the calculation time. During the third optimization, the exact differential process is applied for fine adjustment for the target of GDD.

#### 4.2.2. Optimization methods

As we said before, to minimize the ripples in GDD curves of broadband chirped mirrors, especially ones operate over the wavelength ranges up to an octave spanning or more, people developed complementary CM pairs. However, in the octave-spanning lasers, particularly in the high repetition rate (~1 GHz) ring cavity lasers (Fortier et al. 2006; Nogueira et al. 2006), there are three CMs and one output coupler. Two of the CMs are usually paired, but the other one has no mirror to pair with. Therefore, the total dispersion is usually heavily oscillated, and the compensation is incomplete in such a laser (Fortier et al. 2006; Nogueira et al. 2006). The same difficulty also happens to the linear cavity, where the end mirrors offer one bounce reflection while other CMs experience twice which can be designed in pair. Silver mirror has been used as the end mirror to avoid this problem (Matos et al. 2004). However, the high reflection loss and the weak durability of the silver mirror make the laser inefficient and short-lived. Though octave spanning spectrum has been obtained by some groups with such a cavity, it requires quite a critical alignment. Therefore, more precise overall intracavity dispersion compensation scheme has to be developed.

In 2008, a new cavity for octave-spanning laser is proposed by the MIT group (Crespo et al. 2008), where two pairs of chirped mirrors (four chirped mirrors) were designed and combined to obtain quite a low ripple of the intracavity GDD. The output coupler, however, was performed by the coating on a fused-silica wedge. Sub-two-cycle pulses were obtained.

In 2009, another method for precise intracavity dispersion compensation: the integrative optimization of all cavity mirrors, provided that the output coupler has a flat GDD (Chen et al. 2009) is proposed by us. Both GDD and TOD were considered. This technique allows minimizing the ripples of the whole GDD for ultra-broadband femtosecond lasers.

Take an example as in four-mirror Ti:Sapphire ring cavity for this intracavity optimization to obtain the octave spanning. The four-mirror ring cavity for optimization is shown in Fig. 9. In an octave-spanning Ti:Sapphire laser, one or two mirrors have to have a low reflectivity at 532 nm for the pumping beam to go through while to maintain a high reflectivity from 600 nm to 1200 nm. This constrict makes the oscillation of the GDD larger than the other mirrors. Therefore, we first designed this pump mirror, and then matched the dispersion of this mirror with the other two mirrors. After that, we optimized these three mirrors together with all requirements (reflection and GDD). In Fig. 9, M1 and M2 are concave mirrors. M1 is the pump mirror. M3 may be a convex mirror (Fortier et al. 2006; Nogueira et al. 2006). The optimization procedure for two mirrors was similar to that in reference (Crespo et al. 2008), but with more steps as described below.

The reflectivity of the pump mirror M1 was designed to be less than 20% in the wavelength range of 500-540 nm, particularly less than 5% at 532 nm. The reflectivity of the other two mirrors was set without the restriction on the transmission at 532 nm. The reflectivity of all the mirrors was maintained to be more than 99.8% in the desired spectral range of 600-1200 nm. The target GDD and TOD of each CM were set to be -50 fs^{2} and -35 fs^{3} at the central wavelength of 800 nm so as to compensate the dispersion employed by the 2.3 mm Ti:Sapphire crystal (The GDD employed by the output coupler is assumed to be zero). For all other wavelengths the GD and GDD were calculated respectively, as for the optimization target of 351 points in the wavelength range of 500-1200 nm. The initial design of all mirrors was adapted from the double chirped mirror. The optimization process was adapted as we said in 4.2.1, which was divided into 3 stages and merit functions were modified as below:

The final optimization results for all CMs are shown in Fig. 10 and Fig. 11. Figure 10 shows that the reflectivity of M1 has a high transmission at 532 nm, and all three mirrors have the reflectivity more than 99% for the wavelength range of 600-1200 nm. It also shows that the combined GD of the three mirrors agrees well with the desired GD curve. This in turn results in a best fit with the GDD curve (Fig. 11). The GDD ripples of the combined three mirrors are less than 25 fs^{2} for most wavelengths from 600 to 1200 nm. This is in comparison with the unbalanced mirror dispersion in reference (Fortier et al. 2006), where ripples of net GDD were as high as more than 100 fs^{2}.

Although the combined GDD is smooth and flat, the large ripple of GDD in each individual mirror is still a concern that may destroy the pulse evolution. To testify the pulse profile changes on those mirrors, we introduce a pulse on the three mirrors in sequence and look at the pulse evolution. The pulse is assumed to be secant hyperbolic and chirp-free with spectrum of 696–904 nm, corresponding to a transform-limited pulse of 8 fs (FWHM). For each CM, a constant GDD is added so that the GDD at 800 nm is based to be 0 fs^{2}. In this way, the effect of the GDD ripples on the pulse could be clearly seen. When this pulse goes to those three mirrors in the order of M1, M2, and M3, the pulse duration becomes 8.21 fs, 8.80 fs, 9.31 fs and the corresponding energy in the main pulse is 97.14%, 92.89%, and 92.24% of the initial value respectively (Fig.12). When three mirrors are combined together, the

pulse width becomes 9.31 fs after one bounce and 10.53 fs after two bounces. The energy in the main pulse is 85.86% of the initial value after two bounces. It is noticed that slight satellite pulses appear after two bounces but the pulse still maintains its shape (Fig. 12), indicating that the residual GDD ripples of these three CMs do not destroy the pulse.

In a laser cavity, the pulse is circulating and the accumulated GDD ripples may destroy the pulse formation. However, the amplitude modulation caused by Kerr lens mode locking plays also an important role in pulse shaping. We simulated the pulse evolution in the Ti: Sapphire laser cavity with the master equation (Haus et al. 1992). The Kerr lens mode locking was modeled as a saturable absorber described by

where a is the electrical field, q_{0} the modulation depth, and P_{A} the saturation fluence of the absorber. By solving the master equation, we obtained a stable solution for the pulse profile. The intracavity pulse spectrum expands from 600 nm to 1200 nm at about 20 dB, shown in Fig. 13. The simulation demonstrated that the ultrashort pulses in the Ti:Sapphire laser would ultimately build up in such a cavity without pulse splitting or collapse and the intracavity GDD would be compensated well with these three CMs.

The inconvenience of this scheme in practice is that one additional mirror has to be coated separately, in addition to the chirped mirror pair. This will increase the manufacturing cost as well as the risk of phase errors resulted from imperfect coating. Nevertheless, the new scheme should be able to significantly reduce uncontrolled ripples and hopefully to achieve octave-spanning spectrum easier than using unbalanced mirrors.

In summary, we noticed from the simulation that this multi-mirror simultaneous optimization technique can reduce the total cavity dispersion ripples in a ring cavity. This technique allows minimizing the ripples of the whole GDD for ultra-broadband femtosecond lasers. This technique is able to easily extend to the design of chirped mirrors in linear cavities, where the number of bounces is difficult to match for paired mirrors. Another logical extension of this technique is the design of multi-chirped-mirror set for over one octave-spanning spectrum for extra-cavity dispersion compensation.

## 5. Mapping dispersion for compensation

In order to show how the dispersion is compensated in a laser system, or in any other extracavity dispersion compensation scheme, Nakanuma et al introduced a so-called dispersion “map” where the second order and the third order dispersion are expressed in a orthogonal coordinates (Nakanuma et al. 1991). The second order dispersion is the horizontal while the third order dispersion is in the vertical coordinate respectively. This map is very useful in selecting dispersion compensation schemes.

Let us take an example. In a Ti:Sapphire laser, if we only consider the dispersion introduced by the sapphire crystal, its second and third order dispersion are fs^{2} and fs^{3} respectively. In the map, the dispersion can be expressed in a vector with the horizontal and vertical coordinates of (Nakanuma et al. 1991). The purpose of using this map is to basically bring the vector back to zero, by introducing prism pairs. The prism pair dispersion depends on the prism separation and insertion, which are not independent. However, by the geometry introduced in section 3, for the crystal length of 20 mm, the dispersions could be

that is, in the map, the prism pair dispersion can be expressed by two vectors: one represent the prism separation dependent and the other is the prism insertion dependent. One can chose the prism separation and the prism insertion to make the total dispersion zero. As indicated in reference (Zhang et al. 1997), this can be made by solving equations

To obtain finite and positive solutions of l and d, D_{2t,} D_{3t,} D_{2d} and D_{3d} should satisfy either condition 1:

or condition 2:

where

In Ti:sapphire crystal, at the wavelength of 800 nm, the slope of the crystal is 0.7 fs; if we chose SF10 glass prism pair (A_{1}-B_{1}-C_{1}), the slope R_{l} and R_{d} are 1.71 fs and 0.69 fs and respectively. They cannot satisfy either condition, therefore the vector is open. One can choose prism material with R_{l} and R_{d} as small as possible, such as fused silica (A_{1}-B_{2}-C_{2}). However, because the unit dispersion is smaller too, the required prism separation would be enormously long that is not acceptable in a laser cavity. The only choice is to shorten the crystal length so that the residual third order dispersion is minimized (A_{2}-B_{3}-C_{3}) for the shortest pulse generation.

However, as the laser operation wavelength changes, the corresponding material dispersion changes, as well as the dispersion slope does. In this case, the map can help choosing the prism materials, if a prism pair is preferred. For example, a Cr:forsterite laser, which the operation wavelength is 1300 nm, the slope of the dispersion is 3.8 fs, much higher than at

800nm. Figure x shows a typical map for forsterite at 1290 nm. The gain material dispersion D_{2g} and D_{3g} are drawn as a vector OA, which shows a significantly steep slope in the map. By investigation, most of the prism materials available do not fulfill condition 1: R_{l} > R_{g}, but fit to R_{l} < R_{g}. Therefore, we have to search for the materials which satisfy the second part of condition 2: R_{d} < R_{l}. Some of the materials do satisfy condition 2 such as SFS01, SF58, and SF59. In the map, the dispersion of a SFS01 prism pair is drawn as another vector AB and BC which bring the gain dispersion vector back to origin so that the net D_{2} and D_{3} are simultaneously zero, to complete a triangle loop. The prism separation l and the insertion d can be derived from section 3. The dispersion vector of the SF14 prism pair is also plotted in dashed lines in Fig. 15 as a comparison. Apparently it is not possible to complete the triangle using the SF14 prism pair at this wavelength and a net positive TOD will remain in the cavity. The vectors for SF6 prisms (not plotted) will also end up with a net positive TOD because of R_{d} ≈ R_{l}. Checking the wavelength dependence of the dispersion material, we found that SF14 at a wavelength shorter than 1230 nm satisfies the condition 2: R_{d} < R_{l}. Therefore, it is reasonable to choose SF14 prisms as the dispersion compensator for the laser operated at a shorter wavelength than 1230 nm. If either condition 1 or condition 2 can be satisfied, there should be no severe restriction on the gain crystal length. This is particularly good for Cr:forsterite because its low thermal conductivity does not allow to use a thin crystal with high doping density. It is reasonable to use a relatively long crystal to ease the thermal problem (Ivanov et al 1995). This is in contrast to the case of Ti:sapphire lasers where neither of the conditions can be fulfilled by available prism materials so that thin crystal is the only choice.

In the wavelength of 1550 nm, the dispersion slope of the Cr:YAG crystal increases to >11 fs. There are no more glass materials making the prism pair have such a high slope. One can draw the map as in Fig. 15 and will find that any glass material will leave a positive and high TOD. Rather, simply a piece of fused silica itself would do the compensation of second order only (Ishida & Naganuma, 1994), because its bulk material has a negative second order and positive third order, while the fused silica prism pair does not offer a negative third order by increasing the prism separation. Therefore, either prism pair or bulk material will ultimately make the third order positive.

It should be noted that, in an inhomogeneously broadened laser medium, such as Nd:glass, this map does not work well, because it requires much more dispersion to make it mode locked (Lu et al. 2001).

The similar method can be applied into the chirped pulse amplifier (CPA). In a grating stretcher (O-A_{1}) compressor (A_{2}-A_{3}) system with identical incident angle and slant distance,

the material dispersion (A_{1}-A_{2}) plays a key role in breaking the balance (i.e., bring the dispersion to the origin) of dispersion introduced by grating pairs. The necessity of compensating the third order dispersion (A_{2}-O) can be very obvious in the map.

## 6. Conclusions

We have reviewed three important dispersion compensation devices including grating pair, prism pair and chirped mirrors. We also summarized the detailed calculation formulism for readers’ convenience. The emphasis was given to the novel design of chirped mirrors. Finally, a dispersion map was introduced for conveniently mapping intracavity and extracavity dispersions of a femtosecond laser or an amplification system.