31 December 2013 Stretching and compressing of short laser pulses by chirped volume Bragg gratings: analytic and numerical modeling
Author Affiliations +
Optical Engineering, 53(5), 051509 (2013). doi:10.1117/1.OE.53.5.051509
Abstract
The theory of stretching and compression of short light pulses by chirped volume Bragg gratings (CBGs) is reviewed based on spectral decomposition of short pulses and on wavelength-dependent coupled wave equations. Analytic theory of diffraction efficiency of CBG with constant chirp and approximate theory time-delay dispersion is presented. Comparison of approximate analytic results with exact numeric coupled-wave modeling shows excellent agreement for smooth heterogeneities of CBGs.
Kaim, Mokhov, Zeldovich, and Glebov: Stretching and compressing of short laser pulses by chirped volume Bragg gratings: analytic and numerical modeling

1.

Introduction

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, Tstr(λ), is compensated for by those of the compression stage, so that Tstr(λ)+Tcompr(λ)const. 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.

2.

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 z 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 (z=0) to the opposite end of the same CBG (z=L), without any additional distortions.

Fig. 1

Basic scheme of stretching–amplification–compression via Chirped Bragg grating (CBG).

OE_53_5_051509_f001.png

We assume that the dielectric permittivity and magnetic permeability of CBG depend on coordinate z as

(1)

ε(z)=(n02+2n0δn(z)+O(δn)2)·εvac,μμvac,
where δn(z) is the variation of refractive index. Then the equation for the complex amplitudes A(z) and B(z) of monochromatic component of electric field

(2)

Ereal(z,t)=12{exp(iωt)[A(t)exp(ikz)+B(t)exp(ikz)]+compl.conj.}
with k=ωn/c,n0n(ω) is

(3)

d2dz2[A(z)eikz+B(z)eikz]+(ωc)2[n02+2n0δn(z)]·[A(z)eikz+B(z)eikz]=0,
where the terms of the order of O(δn)2 were ignored. We also assume that the spatial variation of the refractive index δn(z) has the form

(4)

δn(z)=n2(z)+n1(z)cos[Qz+ψ(z)].

Here, smooth heterogeneity n2(z) and spatial modulation n1(z) are small in comparison with n0 and are relatively slow varying functions of z; ψ(z) is a differentiable function, and Q (rad/m) has the value close to 2k, see Eq. (13). (It should be emphasized that our notation n2(z) 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 A(z) and B(z), keep the resonant terms exp(ikz) and exp(ikz) only, and equalize each of them to zero. As a result, one gets the following system of coupled equations:

(5)

dAdz=i(ω/c)n2(z)·A(z)+i(ω/2c)n1(z)exp[iQz2ikz+iψ(z)]·B(z),dBdz=i(ω/2c)n1(z)exp[iQz+2ikziψ(z)]·A(z)i(ω/c)n2(z)·B(z).

It is convenient now to introduce new amplitudes a(z) and b(z), defined by

(6)

a(z)=A(z)exp{i2[(2kQ)zψ(z)]},b(z)=B(z)exp{i2[(2kQ)zψ(z)]}.

Then for these small amplitudes a(z) and b(z), one gets our system in the form

(7)

ddz[a(z)b(z)]=V^(z)[a(z)b(z)],V^(z)=[iν(z)iκ(z)iκ(z)iν(z)],
Here ν(z) is the local value of detuning, measured in (rad/m), and κ(z) is the local strength of coupling, of dimensions (1/m)

(8)

ν(z)=(ωc)n2(z)+2kQdψ(z)/dz2,κ(z)=(ωc)n1(z)2.

Possible background absorption in the medium with coefficient α(1/m), for intensity I(z)exp(αz), may be described by purely imaginary n2(z), namely n2(z)=iλα/4π.

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 z=0. Linearity of the system allows writing.

(9)

[a(z)b(z)]=[Maa(z)Mab(z)Mba(z)Mbb(z)]·[a(0)b(0)],dM^(z)dz=V^(z)M^(z),M^(0)=[1001].
Actual boundary conditions for the reflection problem are usually prescribed at the opposite ends of CBG; for example

(10)

a(z=0)=1,b(z=L)=0,
where we assumed unit amplitude of incident wave a at the beginning z=0 of CBG and no incident wave b at the end z=L of CBG. Then the reflection and transmission amplitudes are found by imposing the conditions in Eq. (10).

(11)

r(ba)b(0)a(0)=Mba(L)Mbb(L),t(aa)a(L)a(0)=Mbb(L)Maa(L)Mba(L)Mab(L)Mbb(L).
The determinant of the matrix M^(z) from Eq. (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 M^(L).

(12)

r(ba)a(L)b(L)=Mab(L)Mbb(L),t(bb)b(0)b(L)=1Mbb(L),
which is especially convenient for the case when one uses the same CBG both for stretching and for recompression.

Subsequently, we neglect the n2(z) term. Then, equating the local value of detuning ν(z) to zero, one can find the position zB(λ) of the point of exact Bragg resonance for the given vacuum wavelength λ.

(13)

ν(zB)=02πn(λ)λ=2πn(λ0)λ0+12dψdz.

Here λ0 is the central wavelength defined by basic parameter of CBG: Q=4πn(λ0)/λ0, with Q defined by Eq. (4). This condition allows connecting spectral chirp parameter dλ/dzB of CBG with the second derivative of the phase correction ψ(z). Indeed, one has to assume λ=λ(z) in Eq. (13) and take z-derivative of the left- and right-hand sides of that equation. As a result, one gets

(14)

dλdzB=λ24πngr(λ)d2ψdz2,ngr(λ)=n(λ)λdndλ.

Here the expression for ngr is an alternative form of denoting group velocity, vgr=dω/dk=c/ngr and ngr=n(ω)+ω(dn/dω)n(λ)λ(dn/dλ). Constant chirp (dλ/dzB=const) corresponds to ψ(z)=β·(z0.5L)2, where the mid-point z=0.5L of CBG with length L has resonant wavelength λ0. Indeed,

(15)

ψ(z)=β·(z0.5L)2dλdzB=λ22πngrβ.

In the approximation of unchanged group velocity v0=c/ngr, the value of time-delay in the stretching process as a function of wavelength λ is

(16)

T1(λ)TDD(λ)2[zB(λ)z1]ngr/c.

We will see below (Fig. 7 in Sec. 5.3) that this formula yields a reasonable result for CBG with low efficiency only.

3.

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 n20, one can take A(z)1 for all z. Then for the reflection amplitude coefficient, one gets from Eq. (5)

(17)

r(AB)=B()=iκexp[iψ(z)]dz,
which, for ψ(z)=β·(z0.5L)2, yields explicit value of

(18)

|r|2π|κ|2/|β|.

An analytic expression has been derived in Ref. 8, which is valid for any strenth of coupling. It may be simplified for constant chirp profile of the grating, with the use of Eq. (15) to:

(19)

|r|2=1exp(π|κ|2/|β|),|t|2=1|r|2.

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).

4.

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 κ=k(z) and ν=ν(z) in Eq. (7). We take the SVEA equations in the symmetric form.

(20)

ddz(ab)=V^(ab),V^=(iνiκ*iκiν).

In the assumption of ν=const,κ=const, the solution of the form (a,b)=exp(pz)(a0,b0) exists. From that we get the result for two eigenvalues of matrix V^.

(21)

p1,2=±iμ,μ=+ν2|κ|2.

We define the positive root as the one for which, in the region |κ|<|ν|, the real square root μ=+ν2|κ|2 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 p1=+iμ and for p2=iμ are, respectively,

(22)

(a1b1)=const1(1κ/(μ+ν)),(a2b2)=const2(κ*/(μ+ν)1).

Group velocity of the wave a in the absence of coupling, |κ|=0, is taken as va=v0, while group velocity of the wave b is taken as vb=v0. The general expression for the group velocity in the presence of both constant values ν and κ for a given mode may be intuitively written as

(23)

vgr=P[Watt/m2]/U[Joule/m3],
where P is the Poynting vector and U is the energy density. Assuming P=|a|2|b|2, U=const·(|a|2+|b|2) and checking the case of |κ|=0, we get const=1/v0. Thus, in general case, group velocity in the presence of grating becomes

(24)

vgroup=v0|a|2|b|2|a|2+|b|2.

In the approximation of constant (or reasonably slow varying) ν and κ, one gets for mode #1 (a predominantly a-wave) and for mode #2 (a predominantly b-wave)

(25)

v1=(va)=v0|μν|,v2=(vb)=v0|μν|.

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 a (the one to be stretched by CBG) propagates from the point z=z1 of the left boundary of CBG to the left point z=z2 inside CBG, where decisive reflective transformation of a into b takes place. The points z=z2 and z=z3 are defined by the condition ν(z2,3,λ)=±|κ|. In CBG approximated by constant chirp, ν(z,λ)=β·[zzB(λ)], those points are symmetrically positioned around zB(λ), i.e., around perfect Bragg-resonant point.

(26)

z2=zB(λ)z0,z3=zB(λ)+z0,z0=|κ|/|β|.

Fig. 2

Transmitting and reflecting wave propagation in CBG. Here L is the total thickness of CBG; z=z2 is the point where the input wave a of the given frequency ω=2πc/λ. hits the region forbidden for propagation (marked by curved dashed line).

OE_53_5_051509_f002.png

After the decisive reflection from the point z=z2, incident wave of type a gets transformed into reflected wave of type b. Total time-delay due to propagation from the entrance z=z1 to reflection point z=z2 and back is

(27)

T1(λ)T(z1z2z1)[second]=2z1z2[v(z)]1dz,
where v(z)=v0·|μ(z)/ν(z)| as a consequence of Eqs. (22) and (24). At the compression stage, the wave enters CBG through the cross-section z=z4 and propagates in the b-form down to z=z3, and after the reflection, it propagates back to z=z4. The corresponding round-trip delay time of this stage is given by

(28)

T2(λ)T(z4z3z4)[second]=2z3z4[v(z)]1dz,
where again v(z)=v0·|μ(z)/ν(z)|. For a general z-dependence of ν(z) and κ(z), one should numerically find zB(λ),z2(λ),z3(λ) and numerically calculate the corresponding integrals. Calculation can be done analytically for constant chirp and constant κ.

(29)

ν(z)=β·[zzB(λ)],z1=0,z4=L,zB(λ)=L2+(dλdz)1·(λλ0).

Equation (15) yields the expression for dλ/dz. Here λ0 is the central wavelength of CBG. In that case, z2=zB(λ)z0,z3=zB(λ)+z0,z0=|κ|/β, and integrals may be found analytically.

(30)

T(z1z2z1)=2v0[zB(λ)]2z02,T(z4z3z4)=2v0[LzB(λ)]2z02,

(31)

TcycleT(z1z2z1)+T(z4z3z4)=2v0{[zB(λ)]2z02+[LzB(λ)]2z02}.

In the case when both z0zB(λ) and z0LzB(λ), one can expand Tcycle(λ) in terms of small ratio z0/L.

(32)

Tcycle(λ)1v0{2L4z02LL2L24[zB(λ)L/2]2}.

One should take the value of Tcycle 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 Tcycle at |zBL/2|z, and then the formula for time-delay of the peak of recompressed pulse becomes

(33)

Tcycle1v0(2L4z02L).

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 2[zB(λ)z1] by the thickness of forbidden zone: 2z0=2|κ/β|, see Fig. 2. This leads to a shorter delay time of the stretching–recompression cycle: shorter by approximately (4z0/v0), i.e., of the first order in z0.

The second effect is due to considerable (50% and more) decrease of group velocities in the vicinity of reflection points z=z2 and z=z3. The thickness of this vicinity is of the order of z0. Thus, this second effect results in a longer delay time, also of the first order in z0.

What is truly remarkable is that these two effects compensate each other in the first order in z0. The resulting delay time does decrease (in comparison with 2L/v0), but in the second order in the coupling constant, i.e., proportionally to z02.

A monochromatic wave stretches in time from t= to t=+. One can define time delay dispersion TDD(λ) for the quasi-monochromatic packet with the wavelength λ. Many sources advise to one calculate the phase of reflection coefficient φ(λ)=arg(r(λ)) and postulate that

(34)

dωdλ=2πcλ2,TDD(λ)=dφdωdφdλdλdωλ22πcdφ(λ)dλ.

Our numeric modeling shows that this is valid for large modulation of φ(λ) only, |φ|1.

5.

Numerical Modeling

5.1.

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 j in DFT programs takes only non-negative values in the range 0<j<N1, where N 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 j by

(35)

ωj=ωcentral+(2π/Ttotal)×phys(j),phys(j)={j,if0<j<N/2;jN,otherwise.

In that manner, we account for both positive and negative values of frequency detuning. Meanwhile discrete time points were numbered just as

(36)

tj=tinitial+j·Δ,Δ=Ttotal/N.

Here Ttotal is the total time interval under consideration. Average refractive index of PTR glass was taken for each wavelength λj=2πc/ωj as either constant value n=1.4891 (for central wavelength λcentral=1.06μm) or from Sellmeier formula for that glass.10 The system of ordinary differential equations was integrated numerically for each frequency ωj. Sometimes it was advantageous to divide all the integration length L (thickness of CBG) into four separate regions, with accuracy goal achieved for each region independently. Examples below are demonstrated mostly for L=0.1m (10cm) and constant chirp parameter β=8.104rad/m2 from Eq. (15). Value of coupling constant |κ| was chosen via dimensionless parameter S=|κ|L, with S varying from 15 to 150.

5.2.

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 zB(λ), where formal Bragg condition is satisfied for a given wavelength λ. We took a small oscillatory addition to the phase, δψ=acosγz, where γ[1/meter]=2π/Λ1 and Λ1 is the period of perturbations. Since the term dψ/dz is present in the right-hand side of Eq. (13), the equation for detuning zB(λ) becomes

(37)

2πn(λ)λ=2πn(λ0)λ0β(zBb)+aγsinγzB.

If aγ2 is larger than 2|β|, then there is a possibility that Eq. (37) has several solutions (zB)1,(zB)2, … for the given wavelength λ. Then application of the simplest formula

(38)

TDD(λ)2(zB(λ)z1)/v0,
yields multivalued TDD(λ) and gives a suspicion of low-quality recompression, to say nothing about peculiar oscillations, which are even more severe in the stretched pulse. However, the results of our numerical modeling show for phase modulation

(39)

δψ=acos(γz)·exp((z0.5L)2/d2),d=L/4
that if the amplitude 2a (top to bottom) is moderately small, e.g., 2a0.4rad, the influence of formally multivalued feature of TDD(λ) from Eq. (38) 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 L=0.1m and β=8·104rad/m2, coupling coefficient |κ|=S/L, and S=30. Additional modulation δψ from Eq. (39) was taken with the parameters 2a=0.4rad, γ=1405rad/m. so that aγ2/2β=2.47, and the condition of multivalued function TDD(λ) from Eq. (38) was satisfied. Stretching and compression of incident Gaussian pulse with τ0(HWe2IM)=3.0ps 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 Δλ(FW)0.1nm equal to that of a laser pulse at the level of 1/e2, 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.

Fig. 3

Dependence of Bragg resonant wavelength detuning on position inside CBG. Particular perturbation of the phase is given by Eq. (39) with top-to-bottom modulation 2a=0.4rad. Given wavelength corresponds to one or three resonant points.

OE_53_5_051509_f003.png

Fig. 4

Spectra of Fourier transform–limited incident pulse (upper curve) and of recompressed pulse (lower curve); time duration of input pulse Δt (HWe2IM)=3ps.

OE_53_5_051509_f004.png

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 δψ(z) of the grating. Input pulse was a transform-limited Gaussian pulse with Δt(HWe2IM)=3ps. 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 2a0.4rad is not harmful for recompression. Meanwhile, our modeling for perturbations with amplitude 2a0.4rad had shown considerable hindering recompression.

Fig. 5

Intensity profiles of input pulse (left), stretched pulse (multiplied by factor 80 for illustrative purposes, center), and recompressed pulse (right) for CBG with oscillatory perturbations of phase δψ(z) of the grating. Arrow points to a precursor in the recompressed pulse containing 6% of its energy.

OE_53_5_051509_f005.png

5.3.

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 (<1ps 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 2L/vgr. Here vgr is the group velocity at λ0 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 |κ|=S/L taken as a constant through the whole thickness L of CBG; hereafter, we call this type of CBG as uniform or nonapodized. The other type had the profile

(40)

|κ(z)|=SL·apod(z);apod(z)=1(2zL1)10.

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 λ0=1.06μm) for different values of dimensionless parameter S (and for constant chirp β=8·104rad/m2). The results for TDD slightly varied as a function of duration of input Gaussian pulse τ0 (HWe2IM). The shortest pulse duration τ0, for which this CBG still reflected almost all the spectrum, was τ0=3.0ps. The stretched pulse in both cases had a duration of 1ns, i.e., 300 times longer than the input one.

Figure 6 shows the comparison of numerical modeling of time-delay Tcycle 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 S=|κ|L 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 2L/v0 (form horizontal line) versus coupling |κ|=S/L. All data are for the input pulses with τ0=12.0ps at the central wavelength λ0=1.06μm 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.

Fig. 6

Dependence of total time-delay between incident short pulse [Gaussian, with Δt [(HWe2IM)=12ps] and recompressed pulse on the dimensionless coupling coefficient S=|κ|L. Square and triangular points correspond to numerical modeling of unapodized and apodized CBG, respectively. Dashed and dotted lines for corresponding CBG yield the results calculated by analytic model Eqs. (27) and (28). Parabolic curve describing simple formula Eq. (33) (not shown) is almost the same as dashed curve for unapodized CBG. We see quadratic deviation of those curves from 2L/v0 (from horizontal line) versus coupling |κ|=S/L.

OE_53_5_051509_f006.png

Here is an important (and originally unexpected) observation. Too large a value of coupling (S70) for that particular CBG and τ0=3.0ps, 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 Tcycle(λ)=T1(λ)+T2(λ), from integrals [Eqs. (27) and (28)] with different values of zB(λ), see Fig. 7 for apodized CBG with S=90. The presence of shorter values of Tcycle(λ) 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 Tcycle(λ) explains interference-type oscillations in the precursor.

Fig. 7

Dependence of stretching time T1(λ) calculated by Eq. (27) for apodized CBG with S=90, compression time T2(λ), Eq. (28), time of cycle stretching–compression Tcycle(λ)=T1(λ)+T2(λ), time of double pass through blank glass 2L/v0.

OE_53_5_051509_f007.png

6.

Conclusions

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 TDD(λ). An unexpected result of that analytic model is that dependence of TDD on the coupling constant |κ|=S/L starts with terms proportional to |κ|2. We have shown that excessively large coupling in CBG leads to deterioration of recompression quality. We show that perturbations of grating phase with small (<0.4rad top-to-bottom) amplitude do not hinder the recompression quality much, enabling >90% of energy in the recompressed pulse.

Acknowledgments

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.

References

1. 

E. B. Treacy, “Optical pulse compression with diffraction gratings,” IEEE J. Quantum Electron. 5(9), 454–458 (1969).IEJQA70018-9197http://dx.doi.org/10.1109/JQE.1969.1076303Google Scholar

2. 

A. Galvanauskaset al., “All-fiber femtosecond pulse amplification circuit using chirped Bragg gratings,” Appl. Phys. Lett. 66(9), 1053–1055 (1995).APPLAB0003-6951http://dx.doi.org/10.1063/1.113571Google Scholar

3. 

W. Changet al., “Femtosecond pulse spectral synthesis in coherently-spectrally combined multi-channel fiber chirped pulse amplifiers,” Opt. Express 21(3), 3897–3910 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.003897Google Scholar

4. 

K.-H. Liaoet al., “Large-aperture chirped volume Bragg grating based fiber CPA system,” Opt. Express 15(8), 4876–4882 (2007).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.15.004876Google Scholar

5. 

L. B. Glebov, “Volume holographic elements in a photo-thermo-refractive glass,” J. Hologr. Speckle 5(1), 77–84 (2009).1546-900Xhttp://dx.doi.org/10.1166/jhs.2009.011Google Scholar

6. 

G. Changet al., “Femtosecond Yb-fiber chirped-pulse-amplification system based on chirped-volume Bragg gratings,” Opt. Lett. 34(19), 2952–2954 (2009).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.34.002952Google Scholar

7. 

O. V. BelaiE. V. PodivilovD. A. Shapiro, “Group delay in Bragg grating with linear chirp,” Opt. Commun. 266(2), 512–520 (2006).OPCOB80030-4018http://dx.doi.org/10.1016/j.optcom.2006.05.032Google Scholar

8. 

L. Poladian, “Graphical and WKB analysis of nonuniform Bragg gratings” Phys. Rev. E 48(6), 4758–4767 (1993).PLEEE81063-651Xhttp://dx.doi.org/10.1103/PhysRevE.48.4758Google Scholar

9. 

L. D. LandauE. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, §50, problem 4, Pergamon Press, Oxford (1977).Google Scholar

10. 

L. Glebov, “Fluorinated silicate glass for conventional and holographic optical elements,” Proc. SPIE 6545, 654507 (2007).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.720928Google Scholar

11. 

L. Glebovet al., “Volume chirped Bragg gratings—monolithic components for stretching and compression of ultrashort laser pulses,” Opt. Eng.OPEGAR0091-3286Google Scholar

Biography

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.

Sergiy Kaim, Sergiy Mokhov, Boris Y. Zeldovich, Leonid B. Glebov, "Stretching and compressing of short laser pulses by chirped volume Bragg gratings: analytic and numerical modeling," Optical Engineering 53(5), 051509 (31 December 2013). http://dx.doi.org/10.1117/1.OE.53.5.051509
Submission: Received ; Accepted
JOURNAL ARTICLE
7 PAGES


SHARE
KEYWORDS
Fiber Bragg gratings

Numerical modeling

Glasses

Wave propagation

Picosecond phenomena

Modulation

Optical engineering

Back to Top