## Abstract

The letter presents a technique for Nth-order differentiation of periodic pulse train, which can simultaneously multiply the input repetition rate. This approach uses a single linearly chirped apodized fiber Bragg grating, which grating profile is designed to map the spectral response of the Nth-order differentiator, and the chirp introduces a dispersion that, besides space-to-frequency mapping, it also causes a temporal Talbot effect.

© 2007 Optical Society of America

## 1. Introduction

The letter presents a technique for Nth-order differentiation of periodic pulse train. In addition to the interest in optical computing and information systems, Nth-order differentiators are of immediate interest for generation of Nth-order Hermite-Gaussian (HG) temporal waveform from an input Gaussian pulse, which can be used to synthesize any temporal shape by superposition [1]. Several all-fiber schemes have been previously proposed based on long-period fiber gratings [2], phase-shifted fiber Bragg gratings (FBGs) [3,4], two-arm interferometer [5], and two oppositely chirped FBGs [6].

A schematic of the proposed general architecture is shown in Fig. 1. This approach exploits the well-known property of linearly-chirped FBGs, which apodization profile maps its spectral response [6–11]. The dispersion introduced by the FBG must meet two conditions: first, it must be high enough so the grating profile of the FBG maps the spectral response of the differentiator [8,9]. Second, it must meet the temporal Talbot condition [12], so dispersion does not affect to the waveform of the output pulses.

Besides the inherent advantages of FBGs (all-fiber approach, low insertion loss, and the potential for low cost), this scheme avoids the concatenation of N first order differentiator devices, which reduce energetic efficiency and increase the implementation complexity. Two different FBG based approaches have been previously demonstrated for all optical N-order time differentiation, namely multiple-phase-shifted FBG [4] and two oppositely chirped FBGs [6]. Concerning the first solution [4], it provides optical operation bandwidths in the tens-of- GHz, while our approach is specially suited for differentiating ultra-broadband optical waveforms (e.g. picosecond and sub-picosecond optical pulses). Regarding the second one [6], it requires two oppositely chirped FBGs, while this approach only requires one FBG. Furthermore, this approach can simultaneously multiply the input repetition rate. As a drawback, the system depends on the repetition rate of the input pulse train.

## 2. Theory

The analytic expression of an Nth-order differentiator in temporal domain is *f _{out}(t*)=

*d*, where

^{N}_{fin}(t)/dtN*f*and

_{in}(t)*f*are the complex envelopes of the input and output of the system respectively, and

_{out}(t)*t*is the time variable. In frequency domain,

*F*=(

_{out}(*ω*)*jω*)

^{N}*F*(

_{in}*ω*), where

*F*(

_{in}*ω*) and

*F*are the spectral functions of

_{out}(ω)*f*and

_{in}(t)*f*, respectively (ω is the base-band frequency, i.e.,

_{out}(t)*ω*=

*ω*, where

_{opt}-ω_{0}*ω*is the optical frequency, and

_{opt}*ω*is the central optical frequency of the signals). Thus, an Nth-order differentiator is essentially a linear filtering device providing a spectral transfer function of the form

_{0}*H*(

_{N}*ω*)=

*F*(

_{out}*ω*)/

*F*(

_{in}*ω*)=(

*jω*)

*. We are interested in obtaining an analytic expression of a feasible spectral response, so the ideal spectral response function must be windowed,*

^{N}*H*(

_{N,w}*ω*)=

*H*(

_{N}*ω*)

*W*(ω)=(

*jω)*, where

^{N}W(ω)*W(ω)*is a window function.

The objective is to obtain a spectral response proportional to the differentiator spectral response. The chirped FBG introduces a dispersion term and we have *H _{r}*(

*ω*)=(

*R*(

*ω*))

^{1/2}exp(

*jϕ*(ω))∝|

_{r}*H*(ω)|exp(

_{Nw}*jϕ*(

_{N}*ω*)+

*jω*

^{2}$\ddot{\varphi}$

*/2) where*

_{r}*H*(

_{r}*ω*),

*R(ω)*, and

*ϕ*are the spectral response in reflection, reflectivity and phase of the FBG,

_{r}(ω)*ϕ*=phase(

_{N}(ω)*H*(

_{N}*ω*))=phase(

*H*(

_{N,w}*ω*)) and $\ddot{\varphi}$ =∂2

*ϕ*(

_{r}*ω*)/

*∂ω*

^{2}is the first order dispersion coefficient of the FBG, which is a constant value for linearly chirped FBGs. Regarding the reflectivity, we have:

The refractive index of the FBG can be written as:

where *n _{av}(z)* represents the average refractive index of the propagation mode,

*Δn*describes the maximum refractive index modulation,

_{max}*A(z)*is the normalized apodization function,

*Λ*is the fundamental period of the grating,

_{0}*φ(z)*describes the additional phase variation (chirp), and

*z*∈[-

*L/2,L/2*] is the spatial coordinate over the grating, with

*L*the length of FBG. In the following we consider a constant average refractive index

*n*=

_{av}*n*+(Δ

_{eff}*/2), where*

_{nmax}*n*is the effective refractive index of the propagation mode.

_{eff}Notice that when *N* is odd, the differentiator spectral response presents a π-phase shift at *ω*=0. In our approach, this condition is attained by introducing a *π*-phase shift in the grating at z=0.

The chirp factor of the FBG, which is defined as *C _{K}=∂^{2}φ(z)/∂z^{2}*, and the length of the grating

*L*, can be calculated from [13]:

where c is the light vacuum speed, and *Δω _{g}* is the grating bandwidth.

The phase filtering of the FBG must be designed to cause a Talbot effect. In general, the pulses waveform of the periodic pulse train is affected by dispersion, but under Talbot condition [12] the pulses are reflected without undergoing distortion (self-image effect). This condition can be expressed as:

where *s/m* must be an irreducible rational fraction. Repetition rate multiplication can be achieved for *m*>1. As a result the reflected signal has a repetition rate *m* times that of the input signal.

Regarding the amplitude filtering, we apply an apodization profile to the grating that is accurately mapped on the spectral response under high dispersion condition [9], which can be expressed as:

where *Δt _{g}* can be calculated from the temporal length of ℑ

^{-1}[

*H*], and ℑ

_{N,w}(ω)^{-1}denotes inverse Fourier transform.

It is worth noting that in most cases *T*
^{2}≫(Δ*t _{g}*)

^{2}, and it is probable that (5) not only satisfies (6), but greatly exceeds it, so from (4) we can deduce that it is necessary to use a longer FBG than strictly required for space-to-frequency mapping.

The apodization profile can be obtained from the expression [8]:

3. Examples and results

In this section, examples of 1st and 4th order differentiators are designed and numerically simulated. We assume a carrier frequency (*ω _{0}*/2π) of 193 THz, an effective refractive index

*n*=1.45, a band of interest (Δ

_{eff}*ω*/2π) of 5 THz centred at

*ω*

_{0}(

*ω*

_{0}-Δ

*ω*/2≤

*ω*

_{opt}≤

*ω*

_{0}+Δ

*ω*/2), a grating bandwidth Δ

*ω*=Δ

_{g}*ω*, a maximum reflectivity of 50 %, and a pulse train period

*T*=40 ps.

For the first example (1st order differentiator) the corresponding ideal spectral response is *H _{1}(ω)=jω*. We choose a function based on a hyperbolic tangent as window,

*W(ω)=W*=(1/2)[1+tanh(4-|16

_{th}(ω)*ω*/Δ

*ω*|)], and we have

_{g}*H*=

_{1,w}(ω)*H*. The desired reflectivity is obtained from (1):

_{1}(ω)W_{th}(ω)where *C _{R}*=2.1238 is a normalization constant to get a maximum reflectivity of 50 %.

From the temporal length of ℑ^{-1}[*H _{1,w}(ω)*] we obtain Δ

*t*≈2 ps. Using expressions (5)and (6) we have $\mid {\stackrel{\u0308}{\varphi}}_{r}\mid =\frac{s}{m}2.5464\times {10}^{-22}{s}^{2}\u2044\mathrm{rad}$ and |$\ddot{\varphi}$

_{g}_{r}|≫1.5915×10

^{-25}

*s*

^{2}/rad.

We choose $\ddot{\varphi}$
_{r}=-1.2732×10^{-22}
*s*
^{2}/*rad*, where *s*=1 and *m*=2 have been selected. This implies that the input repetition rate is multiplied by two, so we have an output period of repetition *T _{out}*=

*T*/2=20 ps. The odd order of 1st differentiator implies that π-phase shift must be introduced in the grating at

_{in}*z*=0.

Using (7), we obtain Δ*n _{max}*=2.8160×10

^{-4},

*n*=1.45014. Additionally, using (3) and (4), we obtain

_{av}*C*=-7.3508×10

_{K}^{5}rad/m

^{2}and

*L*=41.346 cm. The fundamental period of the grating can be obtained from

*Λ*

_{0}=π

*c*/(

*n*)=535.574 nm. Finally, using (7) we obtain the apodization profile function:

_{av}ω_{0}where *C _{A}*=1.2011 is a normalization constant selected to get a normalized apodization profile function 0≤

*A(z)*≤1, and

*N*=1.

As a second example we design a 4th-order differentiator using the same methodology. We obtain again Δ*t _{g}*≈2 ps, and we design the same technological parameters as in the first example (so repetition rate is also doubled). The apodization profile is given by (9), where

*C*=181.02, and

_{R}*N*=4 (same

*L*and

*C*as for first example).

_{A}Figures 2(a), 2(b), and 2(c) show the results from our numerical simulations corresponding to the first example, and Fig. 2(d), 2(e), and 2(f) show the corresponding to second example. The spectral responses of the designed FBG and the ideal differentiator are showed in Fig. 2(a) and 2(d) for first and second example, respectively. The output pulse corresponding to an input gaussian pulse described by *f _{in,1}(t)*∝exp(-

*t*

^{2}/(2

*σ*

^{2})) with σ=800 fs, are showed in Fig. 2(b), and 2(e) for first and second example, respectively. The output pulse corresponding to an antisymmetric HG pulse described by

*f*(

_{in,2}*t*)∝∂

*f*∝

_{in,1}(t)/∂t*t*·exp(-

*t*

^{2}/(2

*σ*

^{2})) are showed in Fig. 2(c), and 2(f) for first and second example, respectively.

From (4) and (6) we obtain that the length of the grating requires *L*≫0.051 cm for space-to- frequency, which can be satisfied from *L*>5 cm. The length of the grating designed (41.346 cm) is much longer, but is within the accuracy of currently available fabrication techniques, as it can be seen in [11], where Talbot effect was achieved in a 96 cm linearly chirped FBG.

Considerations about the effect of group delay ripple can be found in [14], where it has been associated to the amplitude jitter of the rate-multiplied pulse train.

## 4. Conclusion

In this letter, an apodized linearly chirped FBG is designed and numerically simulated to simultaneously perform amplitude filtering to obtain the spectral response of the Nth-order differentiator, and phase filtering to cause temporal Talbot effect. The amplitude filtering exploits the fact that apodization profile of linearly chirped FBG maps its spectral response. This idea was used in a previous article published recently [6]. The main difference is that second chirped FBG for compensation of dispersion introducing by first FBG proposed in [6] is not required. This method reduces the number of required FBG to one and in addition the use of the Talbot effect allows repetition rate multiplication. These advantages must be weighed against the disadvantage that the system depends on the input repetition rate, and that the FBG length required for temporal Talbot effect can be much longer than the strictly required for space-to-frequency mapping.

It is worth noting that, unlike other pulse shaping techniques based on chirped FBGs [7–11], where the objective is shaping a specific output pulse from a known input pulse waveform, this system performs an operation (Nth-order temporal differentiation) than can be applied over different input waveforms to get the respective output waveforms.

## Acknowledgements

This work was supported by the Spanish Ministerio de Educacion y Ciencia under Project “Plan Nacional de I+D+I TEC2004-04754-C03-02” and “Plan Nacional de I+D+I TEC2007- 68065-C03-02”.

## References and Links

**1. **H. J. A. da Silva and J. J. O’Reilly, “Optical pulse modeling with Hermite - Gaussian functions,” Opt. Lett. **14**, 526- (1989). [CrossRef] [PubMed]

**2. **R. Slavík, Y. Park, M. Kulishov, R. Morandotti, and J. Azaña, “Ultrafast all-optical differentiators,” Opt. Express **14**, 10699–10707 (2006). [CrossRef] [PubMed]

**3. **N. K. Berger, B. Levit, B. Fischer, M. Kulishov, D. V. Plant, and J. Azaña, “Temporal differentiation of optical signals using a phase-shifted fiber Bragg grating,” Opt. Express **15**, 371–381 (2007). [CrossRef] [PubMed]

**4. **M. Kulishov and J. Azaña, “Design of high-order all-optical temporal differentiators based on multiple-phase-shifted fiber Bragg gratings,” Opt. Express **15**, 6152–6166 (2007). [CrossRef] [PubMed]

**5. **Y. Park, R. Slavik, and J. Azaña “Ultrafast all-optical first and higher-order differentiators based on interferometers” Opt. Lett. **32**, 710–712 (2007). [CrossRef] [PubMed]

**6. **M. A. Preciado, V. García-Muñoz, and M. A. Muriel “Ultrafast all-optical Nth-order differentiator based on chirped fiber Bragg gratings,” Opt. Express **15**, 7196–7201 (2007). [CrossRef] [PubMed]

**7. **A. G. Jepsen, A. E. Johnson, E. S. Maniloff, T. W. Mossberg, M. J. Munroe, and J. N. Sweetser, “Fibre Bragg grating based spectral encoder/decoder for lightwave CDMA,” Electron. Lett. **35**, 1096–1097 (1999). [CrossRef]

**8. **M. A. Preciado, V. García-Muñoz, and M. A. Muriel , “Grating design of oppositely chirped FBGs for pulse shaping,” IEEE Photon. Technol. Lett. **19**, 435–437 (2007). [CrossRef]

**9. **J. Azaña and L. R. Chen, “Synthesis of temporal optical waveforms by fiber Bragg gratings: a new approach based on space-to-frequency-to-time mapping,” J. Opt. Soc. Am. B **19**, 2758–2769 (2002). [CrossRef]

**10. **S. Longhi, M. Marano, P. Laporta, and V. Pruneri, “Multiplication and reshaping of high-repetition-rate optical pulse trains using highly dispersive fiber Bragg gratings,” IEEE Photon. Technol. Lett. **12**, 1498–1500 (2000). [CrossRef]

**11. **S. Longhi, M. Marano, P. Laporta, and O. Svelto, “Propagation, manipulation, and control of picosecond optical pulses at 1.5 µm in fiber Bragg gratings, J. Opt. Soc. Am. B **19**, 2742–2757 (2002). [CrossRef]

**12. **J. Azaña and M. A. Muriel, “Temporal Talbot effect in fiber gratings and its applications,” Appl. Opt. **38**, 6700–6704 (1999). [CrossRef]

**13. **J. Azaña and M. A. Muriel, “Real-time optical spectrum analysis based on the time-space duality in chirped fiber gratings,” IEEE J. Quantum Electron. **36**, 517–527 (2000). [CrossRef]

**14. **J. T. Mok and B. J. Eggleton, “Impact of group delay ripple on repetition-rate multiplication through Talbot self-imaging effect,” Opt. Commun. **232**, 167–178, (2004). [CrossRef]