Chirped pulse amplification1 is the main approach for high-power/high-energy ultrashort laser pulse’ generation. It is achieved by stretching ultrashort laser pulses. Those stretched pulses have a much lower value of instantaneous power and, therefore may be amplified by broadband laser amplifiers without damaging the medium. Amplified chirped pulses may then be directed to the compressing element, which allows collection of all spectral components back into ultrashort pulse, but with higher energy. The conventional way for doing this is the use of a pair of diffracting gratings for both stretching and compression (Treacy stretchers and compressors1). An especially elegant idea was suggested and implemented by Galvanauskas et al.2,3 with the use of chirped fiber Bragg gratings, which are gratings with a gradually varied period along the fiber. It is possible to use the same grating for compression, and stretching stage, but to illuminate it from the opposite end. In this case, the influence of smooth heterogeneities of time-delay dispersion (TDD) at the stretching stage, , is compensated for by those of the compression stage, so that . A similar performance was demonstrated with the use of chirped volume Bragg gratings (chirped VBGs or CBGs), which are produced by holographic recording in the bulk of photo-thermo-refractive (PTR) glass and have dramatically higher apertures when compared with those of fiber gratings.4 Recent development of VBGs based on PTR glass5 allowed for operation of stretching-compression schemes at much higher values of power, see for example, Ref. 6.
The present paper is devoted to the development of numerical and analytic tools for the study of stretching and compression by volume CBGs. While this analysis is applicable for any volume CBG, the examples shown in this paper are based on CBGs recorded in PTR glass.
Basic Scheme and System of Equations, Definition of Chirp, and Time-Delay Dispersion
The basic scheme of stretching, amplification of stretched pulse, and subsequent recompression back into short pulse by means of a volume CBG is presented in Fig. 1. Due to gradual variation of CBG period in the direction, different spectral components of an incident pulse are reflected from different parts of the CBG and, therefore, have different delays. After amplification, the stretched pulse is launched to the same CBG from the opposite side and compressed back to its original width. In this work, we consider the amplifier as a linear device that does not affect any parameters of a laser pulse but power. Therefore, we assume that a CBG-reflected stretched pulse is transmitted from the cross-section at the front surface of the CBG () to the opposite end of the same CBG (), without any additional distortions.
We assume that the dielectric permittivity and magnetic permeability of CBG depend on coordinate as
Here, smooth heterogeneity and spatial modulation are small in comparison with and are relatively slow varying functions of ; is a differentiable function, and () has the value close to , see Eq. (13). (It should be emphasized that our notation for the smooth correction to the background refractive index has nothing to do with the nonlinearity of refractive index.) In the standard approach of slow varying envelope approximation (SVEA), one should ignore second derivatives of slowly varying amplitudes and , keep the resonant terms and only, and equalize each of them to zero. As a result, one gets the following system of coupled equations:
It is convenient now to introduce new amplitudes and , defined by
Then for these small amplitudes and , one gets our system in the form
Possible background absorption in the medium with coefficient , for intensity , may be described by purely imaginary , namely .
Solution of the linear system of ordinary differential equations [Eq. (7)] may be represented as the solution of the Cauchy problem, where both boundary values are known at the same point, e.g., at . Linearity of the system allows writing.10). 9), which is present in the numerator of Eq. (11), is unity as the consequence of Eq. (7), even in the presence of absorption. This allows additional checking of the accuracy of numeric integration of those equations.
The amplitude of reflection and transmission for the opposite boundary conditions may be expressed through the elements of the same matrix .
Subsequently, we neglect the term. Then, equating the local value of detuning to zero, one can find the position of the point of exact Bragg resonance for the given vacuum wavelength .
Here is the central wavelength defined by basic parameter of CBG: , with defined by Eq. (4). This condition allows connecting spectral chirp parameter of CBG with the second derivative of the phase correction . Indeed, one has to assume in Eq. (13) and take -derivative of the left- and right-hand sides of that equation. As a result, one gets
Here the expression for is an alternative form of denoting group velocity, and . Constant chirp () corresponds to , where the mid-point of CBG with length has resonant wavelength . Indeed,
In the approximation of unchanged group velocity , the value of time-delay in the stretching process as a function of wavelength is
Analytic Expression of Diffraction Efficiency of CBG
Consider first a weak CBG for which the first order of perturbation theory is valid. In the zeroth approximation and at , one can take for all . Then for the reflection amplitude coefficient, one gets from Eq. (5)
Derivation of this expression, see Refs. 7 and 8, is similar to the calculation of quantum-mechanical transmissivity of parabolic potential barrier, see Ref. 9, and will not be discussed here. To go over the rather heavy derivation of Eq. (19) from Refs. 8 and 7, we suggest using the known structure of Eq. (19) and checking all the coefficients in it via the first-order perturbation result [Eq. (18)]. Actual numerical modeling (see below) confirmed the validity of Eq. (19) with great accuracy, especially for apodized CBG, where spatial refractive index modulation comes smoothly to zero at the ends of CBG in this particular modeling, see Eq. (40).
Approximate Expression for Time-Delay Dispersion
Approximate evaluation of TDD may be done on the basis of assumption of slowly varying behavior of the coefficients and in Eq. (7). We take the SVEA equations in the symmetric form.
In the assumption of , the solution of the form exists. From that we get the result for two eigenvalues of matrix .
We define the positive root as the one for which, in the region , the real square root satisfies the condition ; in other words, the sign of is the same as the sign of in that region where propagation is not forbidden. Eigenvectors for and for are, respectively,
Group velocity of the wave in the absence of coupling, , is taken as , while group velocity of the wave is taken as . The general expression for the group velocity in the presence of both constant values and for a given mode may be intuitively written as
In the approximation of constant (or reasonably slow varying) and , one gets for mode #1 (a predominantly -wave) and for mode #2 (a predominantly -wave)
We take those expressions as approximations for group velocities of incident and reflected waves in the use of CBG for the stretcher-compressor scheme and verify the degree of their validity by comparison with results of direct wave modeling of CBG, Sec. 5.
In this approximate approach (Fig. 2), the incident wave (the one to be stretched by CBG) propagates from the point of the left boundary of CBG to the left point inside CBG, where decisive reflective transformation of into takes place. The points and are defined by the condition . In CBG approximated by constant chirp, , those points are symmetrically positioned around , i.e., around perfect Bragg-resonant point.
After the decisive reflection from the point , incident wave of type gets transformed into reflected wave of type . Total time-delay due to propagation from the entrance to reflection point and back is22) and (24). At the compression stage, the wave enters CBG through the cross-section and propagates in the -form down to , and after the reflection, it propagates back to . The corresponding round-trip delay time of this stage is given by
Equation (15) yields the expression for . Here is the central wavelength of CBG. In that case, , and integrals may be found analytically.
In the case when both and , one can expand in terms of small ratio .
One should take the value of from Eqs. (31) or (32) above and then average the expression over the spectral content of the pulse in question. As an even more crude estimate, one may try to take the value of at , and then the formula for time-delay of the peak of recompressed pulse becomes
It is worth discussing separately two effects that influence the delay time of the cycle stretching–recompression.
The first one is that the round-trip length of the stretching process is shorter than by the thickness of forbidden zone: , see Fig. 2. This leads to a shorter delay time of the stretching–recompression cycle: shorter by approximately , i.e., of the first order in .
The second effect is due to considerable ( and more) decrease of group velocities in the vicinity of reflection points and . The thickness of this vicinity is of the order of . Thus, this second effect results in a longer delay time, also of the first order in .
What is truly remarkable is that these two effects compensate each other in the first order in . The resulting delay time does decrease (in comparison with ), but in the second order in the coupling constant, i.e., proportionally to .
A monochromatic wave stretches in time from to . One can define time delay dispersion for the quasi-monochromatic packet with the wavelength . Many sources advise to one calculate the phase of reflection coefficient and postulate that
Our numeric modeling shows that this is valid for large modulation of only, .
Parameters of Numerical Modeling
We have developed the program for numerical modeling of the scheme presented in Fig. 1 with the use of Mathematica software package. In our approach, the fields of all pulses were decomposed into time Fourier series via discrete Fourier transform (DFT) subroutine. The discrete index in DFT programs takes only non-negative values in the range , where is the total number of points (either in time domain or in frequency domain). For that reason the frequency of an individual component was related to that index by
In that manner, we account for both positive and negative values of frequency detuning. Meanwhile discrete time points were numbered just as
Here is the total time interval under consideration. Average refractive index of PTR glass was taken for each wavelength as either constant value (for central wavelength ) or from Sellmeier formula for that glass.10 The system of ordinary differential equations was integrated numerically for each frequency . Sometimes it was advantageous to divide all the integration length (thickness of CBG) into four separate regions, with accuracy goal achieved for each region independently. Examples below are demonstrated mostly for () and constant chirp parameter from Eq. (15). Value of coupling constant was chosen via dimensionless parameter , with varying from 15 to 150.
Modeling of Stretching–Compressing by CBG with Imperfections
It should be noted that CBGs with large apertures and long stretching times4,6,11 show some spatial variations of phase resulting from optical aberrations in a holographic recording system and optical homogeneity of a recording material (PTR glass). This is why, among other results of our numerical modeling of stretching–compression process, we would like to discuss here the influence of small and very inhomogeneous variations of phase of the grating. Equation (13) allows finding the position , where formal Bragg condition is satisfied for a given wavelength . We took a small oscillatory addition to the phase, , where and is the period of perturbations. Since the term is present in the right-hand side of Eq. (13), the equation for detuning becomes
If is larger than , then there is a possibility that Eq. (37) has several solutions ,, … for the given wavelength . Then application of the simplest formula38) is not important. In particular, Fig. 3 shows z-dependence of Bragg resonant wavelength detuning (in nanometers) for particular perturbation of the phase given by Eq. (39). Actual wave modeling of stretching by perturbed CBG was done for and , coupling coefficient , and . Additional modulation from Eq. (39) was taken with the parameters , . so that , and the condition of multivalued function from Eq. (38) was satisfied. Stretching and compression of incident Gaussian pulse with was modeled and depicted in Fig. 4. It shows spectra of the incident pulse (upper curve) and of the recompressed pulse (lower curve). One can see that for a CBG with spectral width equal to that of a laser pulse at the level of , the spectral width of a recompressed pulse is practically the same. However, there is some modulation of the spectrum in the vicinity of the maximum, which is caused by ripples in the dispersion curve depicted in Fig. 3; cutting off short- and long-wavelength wings of the spectrum and the spectral width of the CBG is due to finite spectral width of CBG.
Figure 5 shows the intensity profiles of the input pulse, stretched pulse (multiplied by factor 80 for illustrative purposes), and recompressed pulse for CBG with the above oscillatory perturbations of the phase of the grating. Input pulse was a transform-limited Gaussian pulse with . The properties of CBG are the same as for Figs. 3 and 4. There is a small precursor in the recompressed pulse (arrow in Fig. 5) containing ~6% of recompressed energy. Nevertheless one observes rather good quality of recompression: diffraction efficiency of full cycle of stretching–compression for energy was 0.92 and ratio of peak intensities was 0.61. It means that steep perturbation of phase with amplitude is not harmful for recompression. Meanwhile, our modeling for perturbations with amplitude had shown considerable hindering recompression.
Study of Time-Delay Dispersion
For the evaluation of our analytic (but approximate) results for time-delay dispersion, TDD(), we had to choose a special approach in numerical modeling. Namely, illumination of CBG by a very short incident pulse ( and therefore with a very broad spectrum) resulted in a stretched pulse with very flat top, for which it was difficult to find TDD at the stretching stage. If any distortions in the temporal profile of stretched pulse were present, the arrival time of the peak of reflected pulse was difficult to determine. On the other hand, if the incident pulse was relatively long (and hence was relatively narrow-band), the reflected pulse was not stretched to full possible duration . Here is the group velocity at in unexposed glass. Still, determining the delay time for stretching stage was difficult due to large duration of the input pulse.
For that reason we calculated from numerical wave modeling the delay time of a recompressed pulse in comparison with the arrival time of an incident pulse. The delivery of the stretched pulse to the back of CBG for recompression was considered to be instantaneous in the modeling. Two types of CBG were considered. One type was with taken as a constant through the whole thickness of CBG; hereafter, we call this type of CBG as uniform or nonapodized. The other type had the profile
This apodization function suppressed not-quite Bragg contributions of the ends of CBG; the latter contribution leads to Fresnel function-like oscillations in the reflection spectrum of nonapodized CBG. For both cases we calculated TDD (at ) for different values of dimensionless parameter (and for constant chirp ). The results for TDD slightly varied as a function of duration of input Gaussian pulse (). The shortest pulse duration , for which this CBG still reflected almost all the spectrum, was . The stretched pulse in both cases had a duration of , i.e., 300 times longer than the input one.
Figure 6 shows the comparison of numerical modeling of time-delay with our analytic approximation, Eqs. (27), (28), and (33). The square points in Fig. 6 show full duration of stretching–recompression process versus dimensionless coupling coefficient for unapodized CBG. Triangular points describe corresponding data for apodized CBG. Dashed and dotted lines are calculated by the analytic model [Eqs. (27) and (28)]. Parabolic curve describing the simple formula [Eq. (33)] (not shown) is almost the same as the dashed curve for unapodized CBG. We see quadratic deviation of those curves from (form horizontal line) versus coupling . All data are for the input pulses with at the central wavelength of our CBG. We see that the approach used in the derivation of the integrals [Eqs. (27), (28), and (33)] yields very reasonable correspondence with the results of numerical modeling.
Here is an important (and originally unexpected) observation. Too large a value of coupling () for that particular CBG and , while yielding 99.4% diffraction efficiency for energy of recompressed pulse, resulted in relatively poor quality of recompression: peak intensity of recompressed pulse constituted 52% of the incident one. Our interpretation of that decrease of recompression quality is the following. Different spectral components of incident pulse exhibit somewhat different values of the sum , from integrals [Eqs. (27) and (28)] with different values of , see Fig. 7 for apodized CBG with . The presence of shorter values of points to the formation of precursor in the recompressed pulse. This effect leads to deteriorated recompression. Meanwhile, the fact that two considerably distant wavelengths have those shorter values of explains interference-type oscillations in the precursor.
We have developed a detailed Mathematica-based numerical tool for modeling the process of stretching–recompression by CBG with arbitrary profiles of the grating’s phase and coupling coefficient. To better understand the results of numerous variants of that modeling, we developed approximate analytical model of time-delay dispersion . An unexpected result of that analytic model is that dependence of TDD on the coupling constant starts with terms proportional to . We have shown that excessively large coupling in CBG leads to deterioration of recompression quality. We show that perturbations of grating phase with small ( top-to-bottom) amplitude do not hinder the recompression quality much, enabling of energy in the recompressed pulse.
The work was partially supported by Navy contracts NN00024-09-C-4134 and N68335-12-C-0239 and HEL JTO contract W911NF-10-1-0441. Authors are grateful to Vadim Smirnov for fruitful discussions.
Sergiy Kaim graduated from Mechnikov Odessa University in 2006, where he received his MS degree in physics. Since 2008, he is pursuing his PhD in physics at the University of Central Florida, where his research interests are in the areas of laser beam quality characterization, efficiency and quality of laser beam combination, and properties of volume Bragg gratings.
Sergiy Mokhov graduated from the department of theoretical physics of Kiev University in Ukraine in 1996. Then he pursued his PhD in nuclear physics at that university. In 2004, he moved to CREOL–the College of Optics at the University of Central Florida. From 2006, he has worked with Professor B. Zeldovich and Professor L. Glebov’s group on the theory of volume Bragg gratings. He earned his PhD in optics from CREOL in 2011.
Boris Y. Zeldovich graduated from Moscow University in 1966 and received his doctor of physical and mathematical sciences degree from the Lebedev Physics Institute in Moscow, Russia. He is currently a professor of optics and physics at the University of Central Florida’s CREOL College of Optics and Photonics. He is the codiscoverer of optical phase conjugation. He is a member of the Russian Academy of Science and won the 1997 Max Born award for physical optics from OSA.
Leonid B. Glebov got his PhD from the State Optical Institute, Leningrad, Russia (1976). Since 1995, he has been at CREOL, University of Central Florida, as a research professor. He is a founder and VP of OptiGrate Corp. He is the coauthor of a book, 300+ scientific articles, and 10 patents. He is a fellow of ACerS, OSA, SPIE, and NAI. He is a recipient of a Gabor award in holography. His main directions of research are glasses, holography, and lasers.