## Abstract

We propose a high-repetition-rate soliton-train source based on adiabatic compression of a dual-frequency optical signal in nonuniform fiber Bragg gratings. As the signal propagates through the grating, it is reshaped into a train of Bragg solitons whose repetition rate is predetermined by the frequency of initial sinusoidal modulation. We develop an approximate analytical model to predict the width of compressed soliton-like pulses and to provide conditions for adiabatic compression. We demonstrate numerically the formation of a 40-GHz train of 2.6-ps pulses and find that the numerical results are in good agreement with the predictions of our analytical model. The scheme relies on the dispersion provided by the grating, which can be up to six orders of magnitude larger than of fiber and makes it possible to reduce the fiber length significantly.

©1998 Optical Society of America

## 1. Introduction

High-repetition-rate optical pulse sources are key components in designing high-speed fiberoptic communication systems. Recently, several soliton-pulse sources based on lithium niobate modulators have been proposed, and nearly transform-limited pulses at repetition rates of up to 15 GHz have been generated [1]. However, it is difficult to achieve repetition rates > 20 GHz by this electronic technique. An all-optical method, capable of generating soliton-like pulse trains at high repetition rates, makes use of modulational instability in optical fibers [2]. Pulse trains at repetition rates of up to 300 GHz have been generated by using such a method [3]. However, individual pulses contain a significant pedestal, leading to nonlinear interactions between neighboring solitons.

A novel all-optical technique which was proposed a few years ago makes use of adiabatic compression of a dual-frequency signal inside an optical fiber and has been used to generate a stable train of pedestal-free, non-interacting solitons [4–11]. It has been shown that a dispersion-decreasing fiber (DDF) can be used for adiabatic pulse compression and, in fact, is mathematically equivalent to a uniform fiber with distributed gain [5,6]. If the group-velocity dispersion (GVD) inside a DDF decreases adiabatically, the soliton self-adjusts to preserve the balance between dispersion and nonlinearity, resulting in pulse compression and enhancement of peak power. Soliton pulse trains with a repetition rate of up to 200 GHz have been generated by using a DDF [6]. However, because of the relatively small GVD of optical fibers, this technique requires long fiber lengths (> 1 km). Moreover, fabrication of fibers with complex dispersion profiles usually involves splicing of several different fibers or drawing the fiber with an axially varying core diameter.

A more attractive solution consists of adiabatic compression of pulses in a highly dispersive nonlinear medium such as a fiber Bragg grating. Grating dispersion just outside the stop band is up to six orders of magnitude larger than that of silica fiber [12,13] and can be tailored simply by changing the grating profile. Indeed, a pulse-compression scheme based on a fiber grating with slowly decreasing dispersion was recently proposed [14]. Moreover, almost any grating profile can be manufactured using the state-of-the-art grating-writing techniques [15,16]. To change the dispersion along the grating length one can vary the average refractive index, the period of the grating (chirped grating), or the index modulation depth. Experimental results on adiabatic soliton evolution in chirped fiber Bragg grating are reported in Ref. [17] elsewhere in this issue. In this paper we consider only changes in the index modulation depth but the theory is general.

We propose a high-repetition-rate soliton source based on adiabatic compression of a dual-frequency signal in a nonuniform fiber Bragg grating operating in transmission. As the signal propagates through the grating whose index modulation depth decreases along its length, it is reshaped into a train of solitons through adiabatic compression, as shown in Fig. 1. In Section II we develop a simple analytical model based on the nonlinear Schrödinger equation (NLSE) for optimization of the grating parameters required for generating high-quality, soliton-like, pulse trains. In Section III we demonstrate numerically that a soliton train at a 40-GHz repetition rate can be generated by this technique and find that the numerical results are in good agreement with the predictions of our analytical model. Finally, in Section IV we summarize our results.

## 2. Analytical model

High-intensity pulse propagation through a fiber Bragg grating is governed by the nonlinear coupled-mode equations of the form [18]

Here *E*
_{+} and *E*
_{-} are the slowly varying amplitudes of forward and backward propagating waves, respectively, *V* = *c*/*n* is the velocity of light in fiber, *n* is the average refractive index, *κ* = *πηΔn*/*λ*
_{B} is the coupling coefficient, *Δn* is the index modulation depth, *η* is the fraction of the energy in the fiber core, *λ*
_{B} is the Bragg wavelength. Γ_{s} and Γ_{×} are the nonlinear parameters responsible for self- and cross-phase modulation such that Γ_{s} = Γ_{×} ≡ Γ_{0} = 2 *π*
*n*
_{2}/*λ*, *n*
_{2} is the nonlinear refractive index, and *λ* is the signal wavelength. The nonlinear coupled mode equations have solitary-wave solutions, called Bragg solitons [19], that have been observed in recent experiments [20–22].

If the grating parameters are nearly constant over the spectral bandwidth of the pulse and the central frequency of the incident pulse is close to but outside the grating stop band, the nonlinear coupled-mode equations can be reduced to a simple NLSE having the form [23–25]

Here *A* is the envelope of the Bloch wave associated with the grating, β_{2} represents the GVD of the grating, Γ is the effective nonlinear coefficient, and *t* is the reduced time. The parameters β_{2} and Γ are frequency dependent inside the grating and are related to the grating parameters through the relations [23–25]

where *δ* = *n*/*c*(*ω* - *ω*
_{B}) represents detuning of the incident signal from the Bragg frequency, *v* = [1-(*κ*/*δ*)^{2}]^{1/2}, and *γ* = (1-*v*
^{2})^{-1/2}. We note that if the grating parameters *κ* or *δ* vary with *z*, both β_{2} and *Γ* also vary along the grating length.

The fundamental soliton, by definition, obeys the following relation [26]

where *L*_{D}
= *T*
^{2}/|β_{2}| is the dispersion length, *L*_{NL}
= (Γ*I*_{in}
)^{-1} is the nonlinear length, *I*_{in}
=
|*A*|^{2} is the peak intensity of the pulse inside the grating, and *T* is a measure of the pulse width. For a pulse of the form sech(*t*/*T*), *T* = *T*
_{FWHM}/1.763, where *T*
_{FWHM} is the intensity full width at half
maximum. The soliton energy *E*
_{S} = 2*I*_{in}
*Tσ*
_{eff}, where *σ*_{eff}
is the effective core area. Note that the intensity inside the uniform fiber grating is enhanced with respect to that outside the grating by a factor 1/*v*, i.e. *I*_{in}
= *I*_{out}
/*v* [25].

We now consider soliton propagation in an axially nonuniform grating for which both β_{2} and Γ are *z* dependent. Combining equations (4) and (5) we obtain the following condition for maintaining *N*=1 soliton in a nonuniform grating:

where ${\mathrm{\beta}}_{2}^{\mathit{\text{eff}}}$ is an effective GVD parameter defined to include the frequency dependence of both β_{2} and Γ (see Eq. (4)) and is given by

From Eq. (6) we find that if ${{\beta}_{2}}^{\mathit{eff}}$ decreases adiabatically with *z* such that the soliton energy *E*_{s}
remains constant, *T* must follow changes in ${{\beta}_{2}}^{\mathit{eff}}$ to maintain the condition (5), while the peak intensity must increase. Note also that as ${{\beta}_{2}}^{\mathit{eff}}$ decreases, the soliton accelerates, i.e. *v* increases with *z*. The soliton width at the end of the grating is given by

We next determine the constrains on the parameters of the grating, *κ* and *L*, for adiabatic compression of a short pulse inside the grating. Following Refs. [5,6], we reduce the NLSE Eq. (3) with z-dependent parameters β_{2} and Γ to a standard NLSE with an effective gain. With the transformation

Eq. (3) becomes

where *τ* = *t*/*T*(0) and assume β_{2} is negative,

The effective gain coefficient *g*(*ξ*) takes into account axial variations of both β_{2} and Γ.

Following Ref. [6], we define the total effective gain as

Using Eq. (11) it can be reduced to

With the help of Eq. (8) and (12) the pulse width at the grating output can be written as

where *G*_{eff}
represents the factor by which the input pulse is compressed. Equations (12) and (13) show that the compression factor is determined by the ratio of the effective dispersion at the end points of the grating.

To ensure that the compression remains adiabatic, we impose the condition that the gain coefficient g is small enough that little amplification occurs over one dispersion length, i.e.

If *g*(ξ) varies slowly with ξ, it can be averaged over the grating length. From Eq. (11), *G*_{eff}
(*L*) = exp(2*gL*) for *z* = *L*. The condition (14) can thus be written as ln(*G*_{eff}
) ≪ *L*/*L*_{D}
or,
using Eq.(12), as

We note that the simple theory developed in this Section is based on the assumption that the second-order dispersion of the grating is dominant, i.e. the impact of higher-order dispersion terms is negligible. This assumption is valid only for input signals whose wavelength is tuned close to a stop-band edge but remains far enough from the edges so that the third- and higher-order dispersion terms are small [13].

3. Numerical simulations

We apply the results of the previous section to design a high-repetition-rate soliton source based on adiabatic compression of a dual-frequency signal in a nonuniform fiber Bragg grating whose index modulation depth Δ*n* decreases along its length (see Fig. 1). Although the theory developed in Section 2 assumes that the input pulses are solitons, it allows to predict the width of compressed soliton-like pulses even if the input signal is sinusoidal and to provide conditions for adiabatic compression. As Δ*n* decreases, grating stop band narrows, as shown in Fig. 1. The dual-frequency optical signal can be obtained from a laser operating simultaneously in two longitudinal modes or by using two distributed feedback semiconductor lasers, each operating in a single longitudinal mode. Coherent beating between the two modes generates a sinusoidally modulated optical signal such that

where *T*_{S}
= (Δ*v*)^{-1} is the modulation period for a mode spacing of Δ_{ν}. We consider a fiber grating whose coupling coefficient is nonuniform with a functional form

where *κ*
_{0} is maximum coupling coefficient, and C is a constant. The effective dispersion parameter corresponding to this profile is given by

Figure 2 shows the axial variations of the coupling coefficient and the effective GVD as functions of *z*. Note that a linear variation of *κ* in Eq. (17) was chosen only for simplicity. In principle, *κ* can have any functional form as long as the conditions (6) and (15) are satisfied.

To test the validity of our simple analytic model, numerical simulations were performed using the nonlinear coupled-mode equations. As an example, we consider propagation of a 40-GHz sinusoidal signal with *T*_{S}
= 25 ps and *I*_{out}
= 7.04 GW/cm^{2} through a 70-cm-long apodized fiber grating. The *κ* profile is given by Eq. (15) with *κ*
_{0} = 80 cm^{-1}, δ = 160 cm^{-1}, and *C* is chosen such that *G*_{eff}
= 5. Since *gL*_{D}
≈ 0.1 for this choice of parameters, the adiabatic-compression condition (14) is reasonably well satisfied. Figures 3 and 4 show pulse shapes and spectra respectively at the input (upper plot) and output (lower plot) ends of the grating. These results show that the 12.5-ps sinusoidal pulses evolve and compress by a factor of about 4.8 to become 2.6 ps soliton-like pulses. This behavior is in good agreement with the prediction of our simple analytic theory. We note that each pulse in the generated train contains a very small pedestal which is not visible on a linear scale. The origin of the pedestal may be related to the fact that the compression is not perfectly adiabatic, resulting in an asymmetric spectrum, as seen in Fig. 4. The time-bandwidth product for each pulse in the train is found to be about 0.34 (FWHM), while it is 0.315 for a soliton, which means that the pulses are very close to transform-limited. Finally, we note that in this example the ratio between the soliton separation and soliton pulse width is 9.6. Since pulses are widely separated, neighboring solitons practically do not interact with each other [27].

## 4. Conclusions

In this paper we have proposed a high-repetition-rate soliton-train source based on a nonuniform fiber Bragg grating. In our scheme, a dual-frequency signal is reshaped into a train of soliton-like pulses as it is compressed adiabatically during its propagation in a fiber grating with a nonuniform coupling coefficient. We develop an simple analytical model to predict the width of compressed soliton-like pulses and to provide conditions for adiabatic compression. We demonstrate numerically that a 40-GHz train of 2.6 ps solitons can be generated using a fiber grating with linearly decreasing coupling coefficient. Higher repetition rates can be achieved using stronger and longer gratings or other types of photonic crystals. Since the proposed scheme relies on the dispersion provided by the grating, which is many orders of magnitude larger than that of silica fibers, the device length can be reduced significantly, compared to fiber-based devices. However, note that the intensities required when silica-fiber gratings are used are very large (~10 GW/cm^{2}). The required intensity can be scaled down using a medium with high nonlinearity. For example, the nonlinear coefficient of As_{2}S_{3} chalcogenide fibers is 100 times larger than that of standard silica fibers [28] and, therefore, the peak intensity of 7 GW/cm^{2} used in our numerical example in Section III will scale down to 70 MW/cm^{2} (peak power ~ 10 W) if such fibers are used.

## Acknowledgments

N.M. Litchinitser thanks the Aileen S. Andrew Foundation for postdoctoral fellowship grant.

## References and links

**1. **J. J. Veselka and S. K. Korotky, “Pulse generation for soliton systems using lithium niobate modulators,” IEEE J. Sel. Top. Quantum Electron. **2**, 300–310 (1996). [CrossRef]

**2. **A. Hasegawa, “Generation of a train of soliton pulses by induced modulational instability in optical fibers,” Opt. Lett. **9**, 288–290 (1984). [CrossRef] [PubMed]

**3. **K. Tai, A. Tomita, J. L. Jewell, and A. Hasegawa, “Generation of subpicosecond solitonlike optical pulses at 0.3 THz repetition rate by induced modulational instability,” Appl. Phys. Lett. **49**, 236–238 (1986). [CrossRef]

**4. **E. M. Dianov, P. V. Mamyshev, A. M. Prokhorov, and S. V. Chernikov, “Generation of a train of fundamental solitons at a high repetition rate in optical fibers,” Opt. Lett. **14**, 1008–1010 (1989). [CrossRef] [PubMed]

**5. **V. A. Bogatyrev, M. M. Bubnov, E. M. Dianov, A. S. Kurkov, P. V. Mamyshev, A. M. Prokhorov, S. D. Rumyantsev, V. A. Semenov, S. L. Semenov, A. A. Sysoliatin, S. V. Chernikov, A. N. Gur’yanov, G. G. Devyatykh, and S. I. Miroshnichenko, “A single-mode fiber with chromatic dispersion varying along the length,” J. Lightwave Technol. **9**, 561–566 (1991). [CrossRef]

**6. **P. V. Mamyshev, S. V. Chernikov, and E. M. Dianov, “Generation of fundamental soliton trains for high-bit-rate optical fiber communication lines,” IEEE J. Quantum Electron. **27**, 2347–2355 (1991). [CrossRef]

**7. **S.V. Chernikov, D.J. Richardson, R.I. Laming, E.M. Dianov, and D.N. Payne, “70 Gbit/s fiber based source of fundamental solitons at 1550nm,” Electron. Lett. **28**, 1210–1212 (1992). [CrossRef]

**8. **M. Romagnoli, S. Trillo, and S. Wabnitz, “Soliton switching in nonlinear couplers,” Optical and Quantum Electon. **24**, S1237–S1267 (1992). [CrossRef]

**9. **S. V. Chernikov, J. R. Taylor, and R. Kashyap, “Comblike dispersion-profiled fiber for soliton pulse train generation,” Opt. Lett. **19**, 539–541 (1994). [CrossRef] [PubMed]

**10. **E. A. Swanson and S. R. Chinn, “23-GHz and 123-GHz soliton pulse generation using two cw lasers and standard single-mode fiber,” IEEE Photon. Technol. Lett. **6**, 796–799 (1994). [CrossRef]

**11. **N. Akhmediev and A. Ankiewicz, “Generation of a train of solitons with arbitrary phase difference between neighboring solitons,” Opt. Lett. **19**, 545–547 (1994). [CrossRef] [PubMed]

**12. **B. J. Eggleton, T. Stephens, P. A. Krug, G. Dhozi, Z. Brodzeli, and F. Ouellette, “Dispersion compensation over 100 km at 10 Gbit/s using a fiber grating in transmission,” Electron. Lett. **32**, 1610–1611 (1996). [CrossRef]

**13. **N. M. Litchinitser, B. J. Eggleton, and D. B. Patterson, “Fiber Bragg gratings for dispersion compensation in transmission: Theoretical model and design criterion for nearly ideal pulse compression,” J. Lightwave Technol. **15**, 1303–1313 (1997). [CrossRef]

**14. **G. Lenz and B. J. Eggleton, “Adiabatic Bragg soliton compression in nonuniform grating structures,” J. Opt. Soc. Am. B, in press.

**15. **A. Asseh, H. Storay, B. E. Sahlgren, S. Sandgren, and R. A. H. Stubbe, “A writing technique for long fiber Bragg gratings with complex reflectivity profiles,” J. Lightwave Technol. **15**, 1419–1423 (1997). [CrossRef]

**16. **L. Dong, M. J. Cole, M. Durkin, M. Ibsen, and R. I. Laming, “40Gbit/s 1.55um transmission over 109km of non-dispersion-shifted fiber with long continuously chirped fiber gratings,” in 1996Opt. Fiber Commun. Conf. (OFC'96), postdeadline paper PD6.

**17. **R. E. Slusher, B. J. Eggleton, C. M. de Sterke, and T. A. Strasser, “Nonlinear pulse reflections from chirped fiber gratings,” Opt. Express **3**, (1998). http://epubs.osa.org/oearchive/source/6996.htm [CrossRef] [PubMed]

**18. **C. M. de Sterke and J. E. Sipe, in *Progress in Optics* XXXIII, ed. by E. Wolf (North-Holand, Amsterdam, 1994), pp. 203-260.

**19. **A. B. Aceves and S. Wabnitz, “Self-induced transparency solitons in nonlinear refractive periodic media,” Phys. Lett. A **141**, 37–42 (1989). [CrossRef]

**20. **B. J. Eggleton, R. E. Slusher, C. M. de Sterke, P. A. Krug, and J. E. Sipe, “Bragg grating solitons,” Phys. Rev. Lett. **76**, 1627–1630 (1996). [CrossRef] [PubMed]

**21. **D. Taverner, N.G.R. Broderick, D.T. Richardson, R.I. Laming, and M. Ibsen, “Nonlinear self-switching and multiple gap-soliton formation in a fiber Bragg grating,” Opt. Lett. **23**, 328–330 (1998). [CrossRef]

**22. **B. J. Eggleton, C. M. de Sterke, and R. E. Slusher, “Bragg solitons in the nonlinear Schrödinger limit: theory and experiment,” submitted to J. Opt. Soc. Am. B (1998).

**23. **C. M. de Sterke and J. E. Sipe, “Coupled modes and the nonlinear Schrödinger equation,” Phys. Rev. A **42**, 550–555 (1990). [CrossRef]

**24. **C. M. de Sterke and B. J. Eggleton, “Bragg solitons and the nonlinear Schrödinger equation,” Phys. Rev. E. , in press.

**25. **B. J. Eggleton, C. M. de Sterke, A. Aceves, J. E. Sipe, T. A. Strasser, and R. E. Slusher, “Modulational instabilities and tunable multiple soliton generation in apodized fiber gratings,” Opt. Commun. **149**, 267–271 (1998). [CrossRef]

**26. **G. P. Agrawal, *Nonlinear Fiber Optics*, (Academic, New York, 1995).

**27. **J. P. Gordon, “Interaction forces among solitons in optical fibers,” Opt. Lett. **8**, 596–598 (1983). [CrossRef] [PubMed]

**28. **M. Asobe, “Nonlinear optical properties of chalcogenide glass fibers and their application to all-optical switching,” Opt. Fiber Technol. **3**, 142–148 (1997). [CrossRef]