## Abstract

In this letter we present a technique for the implementation of Nth-order ultrafast temporal differentiators. This technique is based on two oppositely chirped fiber Bragg gratings in which the grating profile maps the spectral response of the Nth-order differentiator. Examples of 1^{st}, 2^{nd}, and 4^{th} order differentiators are designed and numerically simulated.

© 2007 Optical Society of America

## 1. Introduction

An Nth-order optical temporal differentiator is a device that provides the Nth-time derivative of the complex envelope of an arbitrary input optical signal. This operation is performed on optical devices at operation speeds several orders of magnitude over electronics. These devices may find important applications as basic building blocks in ultrahigh-speed all-optical analog–digital signal processing circuits [1]. Moreover, 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 [2]. Several schemes have been previously proposed based on integrated-optic transversal filter [1], long-period fiber gratings [3], phase-shifted fiber Bragg grating [4], and two-arm interferometer [5].

In this letter we use a technique for temporal differentiation based on the use of two oppositely chirped fiber Bragg gratings (FBG) [6]. As it can be seen in Fig. 1, the system includes two linearly chirped FBGs connected by an optical circulator. The first, FBG_{a}, is the spectral shaper, and provides the spectral response for pulse shaping. The second, FBG_{b}, cancels the dispersion introduced by the first grating. Obviously, the order of the FBGs can be arbitrarily selected.

This scheme has been previously proposed and experimentally demonstrated in [7] and [8]. Specifically, in [7], phase-shifts are introduced in the shaper FBG to generate spectral-phase-encoded bit. In [8], a bandpass Gaussian FBG optical filter in which the bandwidth can be continuously adjusted is presented. Besides the inherent advantages of FBGs (all-fiber approach, low insertion loss, and the potential for low cost), this scheme can provide a direct implementation of an Nth-order differentiator, avoiding the concatenation of N first order differentiators devices. Furthermore, this approach has the possibility of adjusting the bandwidth and tuning the central wavelength [8].

## 2. Theory

The temporal operation of a Nth-order differentiator can be expressed as *f _{out}*(

*t*)=

*d*(

^{N}f_{in}*t*)/

*dt*, where

^{N}*f*(

_{in}*t*) and

*f*(

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

*t*is the time variable. We can also express this in frequency domain as,

*F*(

_{in}*ω*)=(j

*ω*)

*(*

^{N}F_{out}*ω*) where

*F*(

_{in}*ω*) and

*F*(

_{out}*ω*) are the spectral functions of

*f*(

_{in}*t*) and

*f*(

_{out}*t*), respectively (

*ω*is the base-band frequency, i.e.,

*ω*=

*ω*-

_{opt}*ω*, where

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

_{opt}*ω*is the central optical frequency of the signals). Thus, the spectral response of the ideal Nth-order differentiator is:

_{0}$${H}_{N}\left(\omega \right)=\frac{{F}_{\mathit{out}}\left(\omega \right)}{{F}_{\mathit{in}}\left(\omega \right)}={\mathrm{\left(}\mathrm{j}\omega \right)}^{N}$$

Moreover, in a real system we have a finite bandwidth, so we have to window the spectral response function:

where *W*(*ω*) is a window function, which must be selected to meet:

where the operative band is the region where the differentiator operation works with accuracy, and *trans*(*ω*) is a transient function which must have low amplitude values at the edges of the band of interest in order to avoid an abrupt discontinuity. Notice that not any window function verifies this condition on *trans*(*ω*), even in the case of a window function presenting low values at the edges of the band of interest.

The objective is to obtain a spectral response of the whole system (composed by the two FBGs), *H _{syst}*(

*ω*), proportional to the differentiator spectral response:

where *H _{a}*(

*ω*),

*H*(

_{b}*ω*),

*R*(

_{a}*ω*),

*R*(

_{b}*ω*),

*ϕ*(

_{a}*ω*),

*ϕ*(

_{b}*ω*) are the spectral response in reflection, reflectivity and phase of the FBGs. In this approach we assume that FBG

_{b}is a dispersion compensator, so we can consider that

*R*(

_{b}*ω*) presents an ideal flat-top response in the band of interest, so the shape of the reflectivity is influenced by FBG

_{a}solely. Thus, we have:

Regarding the phase, we have two oppositely linearly chirped FBGs, so *$\ddot{\varphi}$ _{a}*(

*ω*) = -

*$\ddot{\varphi}$*(

_{b}*ω*) =

*$\ddot{\varphi}$*, where

_{a}*$\ddot{\varphi}$*(

*ω*) denotes ∂

^{2}

*ϕ*(

*ω*)/∂

*ω*

^{2}, and

*$\ddot{\varphi}$*is a constant value, which is obtained from the FBG

_{a}_{a}design.

At this point, we present the theory to design the spectral shaper, FBG_{a}. The refractive index of FBG_{a} can be written as:

where *n _{av,a}*(

*z*) represents the average refractive index of the propagation mode, Δ

*n*describes the maximum refractive index modulation,

_{max,a}*A*(

_{a}*z*) is the normalized apodization function,

*Λ*is the fundamental period of the grating,

_{0,a}*φ*(

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

*z*∈ [-

*L*/2,

_{a}*L*/2] is the spatial coordinate over the grating, with

_{a}*L*the length of FBG

_{a}_{a}. In the following we consider a constant average refractive index

*n*=

_{av,a}*n*+(Δ

_{eff,a}*n*/2), where

_{max,a}*n*is the effective refractive index of the propagation mode.

_{eff,a}Notice that (1) implies 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 of FBG_{a} at *z*=0 The chirp factor of FBG_{a}, which is defined as *C _{K,a}*=∂

^{2}

*φ*(

_{a}*z*)/∂

*z*

^{2}, and

*L*can be calculated from [9]:

_{a}where *c* is the light vacuum speed, and Δ*ω _{g,a}* is the FBG

_{a}bandwidth. It is well known that when a chirped FBG introduces an enough high dispersion, the spectral response of the grating is a scaled version of its corresponding apodization profile [10]. This high dispersion condition can be expressed as:

where Δ*t _{a}* is the temporal length of the inverse Fourier transform of the FBG

_{a}spectral response without the dispersive term, which can be calculated from the temporal length of ℑ

^{-1}[

*H*(

_{N,w}*ω*)], where ℑ

^{-1}denotes inverse Fourier transform. It is worth noting that the broader (narrower) bandwidth, the shorter (longer) minimum length of the grating required for FBG

_{a}to map properly the spatial profile on the spectral response [6].

If condition (9) is met and Born approximation is applicable, both temporal and spectral envelopes reproduce the shape of the apodization profile function, so we can obtain the apodization profile which corresponds to Ra(*ω*) [11]. In the case of high reflectivity an approximate function [12] must be applied over *R _{a}*(

*ω*). In particular, a logarithmic based function is used in our approach, and we obtain an expression which is valid for both weak and strong gratings:

## 3. Examples and results

Here we give three design examples for 1^{st}, 2^{nd} and 4^{th} order differentiators, which are numerically simulated. For all the examples we assume a carrier frequency (*ω _{0}*/2π) of 193 THz, an effective refractive index

*n*=1.45 for FBG

_{eff,a}_{a}, a band of interest (Δ

*ω*/2π) of 5 THz centred at

*ω*(

_{0}*ω*-Δ

_{0}*ω*/2 ≤

*ω*≤

_{opt}*ω*+Δ

_{0}*ω*/2), a FBG

_{a}bandwidth Δ

*ω*=Δ

_{g,a}*ω*, and a maximum reflectivity for FBG

_{a}of 90 %.

In the first example we design a system which implements a 1^{st}-order differentiator. The corresponding ideal spectral response is *H _{1}*(

*ω*)=j

*ω*, and we choose a function based on a hyperbolic tangent as window,

*W*(

_{th}*ω*)=(1/2)[1+tanh(4-|16

*ω*/

*Δω*|)]:

The spectral shaper (FBG_{a}) must be designed to properly map the desired spectral response. From the temporal length of ℑ^{-1}[*H*
_{1,w}(*ω*)] we obtain Δ*t _{a}*≈2 ps. Using expression (9) we have |

*$\ddot{\varphi}$*|>>1.5915×10

_{a}^{-25}

*s*

^{2}/

*rad*, and choose

*$\ddot{\varphi}$*=-1.6×10

_{a}^{-23}

*s*

^{2}/

*rad*. Moreover, the odd order of 1

^{st}differentiator implies that π-phase shift must be introduced in FBG

_{a}at

*z*=0. The desired reflectivity for FBG

_{a}in the band of interest is obtained from (5):

where *C _{R}*=2.8494 is a normalization constant selected to get a maximum reflectivity max(

*R*(ω))=0.9. Using (10) at the maximum reflectivity and apodization, with max(

_{a}*A*(

_{a}*z*))=1, we obtain Δ

*n*=1.4484 × 10

_{max,a}^{-3},

*n*= 1.45072.

_{av,a}Additionally, using (7) and (8), we obtain *C _{K,a}*=5.8543 × 10

^{6}rad/m

^{2}and

*L*=5.1936 cm. The fundamental period of the grating FBG

_{a}_{a}can be obtained from

*Λ*=π

_{0,a}*c*/(

*n*

_{av}*ω*)=535.36 nm. The period of FBG

_{0}_{a}varies from 542.39 nm to 528.51 nm along the length of the grating. This supposes a relative period variation of 2.591 %, which is within accuracy of currently available fabrication techniques [13].

Finally, using (10) we obtain the apodization profile function:

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

*A*(

_{a}*z*)≤1, and

*N*=1.

Moreover, the dispersion parameter of the dispersion compensator (FBG_{b}) is *$\ddot{\varphi}$ _{b}*(

*ω*) = -

*$\ddot{\varphi}$*= -1.6×10

_{a}^{23}

*s*

^{2}/

*rad*, which must present a flat top spectral response in the band of interest.

As a second example we design a 2^{nd} order differentiator using the same methodology. We obtain again Δ*t _{g,a}*≈2 ps, so we have the same technological parameters as in the first example. The apodization profile which is given by (13), where

*C*= 13.568, and

_{R}*N*=2 (same

*L*and

_{a}*C*as for first example).

_{A}Finally, in a third example, we design a 4^{th} order differentiator. We have again Δ*t _{g,a}*≈2 ps, and the same technological parameters, with an apodization profile described by (13), where

*C*= 243.1, and

_{R}*N*=4 (same

*L*and

_{a}*C*as for first example).

_{A}Figure 2 shows the results from our numerical simulations corresponding to these examples. Figures 2(a), 2(b), and 2(c) show the phase response of the spectral shaper (FBG_{a}), the dispersion compensator (FBG_{b}), and the whole system, for the first, second and third example, respectively. Figures 2(d), 2(e), and 2(f) compare the spectral responses of the spectral shaper (FBG_{a}) and the ideal differentiator, for the first, second and third example, respectively. Figures 2(g), 2(h), and 2(i) show the temporal waveform of the input pulse and the output pulses of the designed system, and the ideal differentiator, for the first, second and third example, respectively. We have applied an input gaussian envelope pulse, described by *f _{in}*(

*t*) ∝ exp(-

*t*

^{2}/(2σ

^{2})), with σ = 500 fs (FWHM= 1.177 ps).

It is worth noting that in our simulations we have supposed ideal cancellation of dispersions of both FBGs. In practice this requires a careful monitoring of the chirp profile of each grating in order to avoid excessive phase ripple. This has been achieved in [7], even with tunable chirp in [8]. Considerations about dispersion and phase ripple tolerance can be found in [8] and [14].

## 4. Conclusion

In this paper, we have presented an Nth-order differentiator based on a pair of oppositely chirped FBGs, and we have analytically designed and numerically simulated three examples, the 1^{st}, 2^{nd}, and 4^{th} order differentiators.

In addition to the inherent advantages of FBGs, we find two main features that could be of practical relevance. Firstly, we can implement a Nth-order differentiator using a single device, which is more energetically efficient than the concatenation of N first order differentiators. Secondly, the proposed scheme allows tuning the central wavelength and adjusting the bandwidth [8] according to the input signal.

## Acknowledgments

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

## References and Links

**1. **N. Q. Ngo, S. F. Yu, S. C. Tjin, and C. H. Kam, “A new theoretical basis of higher-derivative optical differentiators,” Opt. Commun. **230**, 115–129, (2004). [CrossRef]

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

**3. **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]

**4. **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]

**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 “Grating design of oppositely chirped FBGs for pulse shaping,” IEEE Photon. Technol. Lett. **10**, 435–437 (2007). [CrossRef]

**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. **I. Littler, M. Rochette, and B. Eggleton, “Adjustable bandwidth dispersionless bandpass FBG optical filter,” Opt. Express **13**, 3397–3407 (2005). [CrossRef] [PubMed]

**9. **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]

**10. **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]

**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. **B. Bovard, “Derivation of a matrix describing a rugate dielectric thin film,” Appl. Opt. **27**, 1998–2004 (1988). [CrossRef] [PubMed]

**13. **James F. Brennan III and Dwayne L. LaBrake, “Fabrication of chirped fiber bragg gratings of any desired bandwidth using frequency modulation,” US patent 6728444 (April
2004).

**14. **I. C. M. Littler, L. Fu, and B. J. Eggleton, “Effect of group delay ripple on picosecond pulse compression schemes,” Appl. Opt. **44**, 4702–4711 (2005). [CrossRef] [PubMed]