## Abstract

An approach to reconstruct a quasi-phase-matching grating by using a discrete layer-peeling algorithm is presented. Experimentally measured output spectra of Šolc-type filters, based on uniform and chirped QPM structures, are used in the discrete layer-peeling algorithm. The reconstructed QPM structures are in agreement with the exact structures used in the experiment and the method is verified to be accurate and efficient in quality inspection on quasi-phase-matching grating.

© 2012 OSA

## 1. Introduction

Quasi-phase-matching (QPM) nonlinear devices, such as periodically poled lithium niobate (PPLN) and lithium tantalite (PPLT), have been widely studied in the field of nonlinear optical interactions because of their capability of achieving the desired phase-matching condition by spatially modulating the quadratic nonlinearities [1]. More important, possible engineering of QPM structures has afforded new degree of freedom in the nonlinear process. Recent attention has been given from uniform structure to Fibonacci optical superlattice [2], linearly chirped [3], apodized [4,5], multi-segmented [6], and aperiodically QPM gratings [7,8], which can provide more spatial reciprocal vectors to compensate for nonlinear phase mismatch in the broadband nonlinear process and multi-wavelength second harmonic generation (SHG). Engineered QPM structures also are used in controlling the cascaded quadratic nonlinearities in the phase mismatched SHG process by varying the local phase mismatch, which were used to achieve the tunable-chirp pulse compression [9], quadratic nonlinearity management [10] and two-color pulse compression based on quadratic nonlinearities [11]. Recently ultralow-power frequency-resolved optical gating (FROG) was demonstrated by using aperiodically poled lithium niobate (APPLN) because of adequate bandwidth for SHG in Ref. [12]. Thus the period qualities of such non-uniform QPM gratings become crucial in the applications. Usually Hydrofluoric (HF) acid etching method is used to characterize the quality of domain-reversed structure but unavoidably destroy the crystal surface.

We show an approach to reconstruct quasi-phase-matching gratings by using a discrete layer-peeling (DLP) algorithm. DLP method was proven to be a simple, accurate and efficient approach to reconstruct the local structure of refractive index of fiber bragg grating [13–16] and other optical devices [17]. However, as far as we know, no work has been reported yet on the synthesis of QPM structure.

Chen et al. first demonstrated narrow-band Šolc-type wavelength filters experimentally based on PPLN crystal, which is sandwiched in the crossed polarizers [18]. By applying the external electric field along the transverse direction or using the misalignment of LN crystal orientation, the optical axes of ferroelectric domains are alternately rotated with the folded azimuth angle *θ* in successive reversal domains, thus PPLN is similar to so-called folded Šolc filter structure with periodically alternating azimuth angles of the crystal axes [19]. Polarization-independent [20] and narrowband multi-wavelength Šolc-type wavelength filters [21–23] in APPLN and bandwidth control by using local temperature gradient technique [24] were demonstrated in Šolc-type wavelength filters. By confining the light in the waveguides, efficient optical wavelength filter were recently reported [25, 26].

The ordinary and extraordinary waves propagating in QPM gratings are phase-matched due to the reciprocal vectors from the modulation of electro-optical (EO) coefficient in Šolc-type filters, which can be considered as a grating-assisted co-directional coupler, such as long period fiber grating based on refractive index modulation in the coupling region. Therefore the polarization coupling in QPM gratings can be described by using co-propagating coupling mode theory [27]. The grating-assisted co-directional coupler is usually treated as a finite-impulse-response (FIR) filter and each point in the impulse response is determined by the entire grating [28, 29]. We successfully use DLP algorithm to reconstruct the poling structure of QPM gratings in this work.

## 2. Numerical method and analysis

The coupling between ordinary and extraordinary waves in QPM gratings is described by the well-known coupled-mode equations as follows [19]:

*A*(

_{o}*z*, Δ

*β*) and

*A*(

_{e}*z*, Δ

*β*) represent the normalized amplitudes of ordinary and extraordinary waves, respectively, which can be considered as slowly varying amplitudes of co-directional propagating waves. $\kappa (z)=-\text{sin}\left(2\theta \right)\frac{\omega}{c}\frac{{n}_{e}^{2}-{n}_{o}^{2}}{4\sqrt{{n}_{o}{n}_{e}}}{e}^{i\left[\frac{\pi}{2}+\varphi (z)\right]}$ is the coupling coefficient along QPM grating [19], while

*ϕ*(

*z*) = [2

*π*/Λ(

*z*) − 2

*π*/Λ

_{0}]

*z*denotes the local phase variation induced by QPM grating chirp in ununiformed QPM structures, and

*n*and

_{o}*n*are refractive indices of ordinary and extraordinary waves. The detuning parameter Δ

_{e}*β*for QPM structure is defined as: Δ

*β*= (

*β*−

_{o}*β*) − 2

_{e}*πm*/Λ

_{0},

*β*and

_{o}*β*are wave vectors of ordinary and extraordinary waves and Λ

_{e}_{0}=

*λ*

_{0}/(

*n*−

_{o}*n*) represents the reference grating period and

_{e}*m*is the ordinary QPM order. The polarization coupling between co-directional propagating waves occurs at the resonance frequency

*ω*

_{0}when satifaying the phase matching condition: Δ

*β*(

*ω*

_{0}) = 0. The local coupling coefficients are determined by the difference of refractive indices Δ

*n*=

*n*−

_{o}*n*and QPM grating period Λ(

_{e}*z*) simultaneously. We assume the dispersion of refractive index negligible and keep Δ

*n*constant in the simulation. Now that we need to figure out the period distribution of QPM gratings from the coupling coefficient

*κ*(

*z*).

The QPM structure is divided into *N* layers of uniform sections separated by a distance Δ*z*. The discretized coupling model is illustrated in Fig. 1. Based on the piecewise uniform coupling model [30], we can separate the distributed coupling between ordinary and extraordinary waves into a stack of discrete instantaneous scattering points.

The matrix *G*_{Δ} describes the pure propagation between the instantaneous scattering points and the other matrix *G _{ρ}* describes the coupling in the

*j*-th layer of QPM structure [31].

*ρ*is the coupling strength of the instantaneous scattering point

_{j}*z*along the structure. The discretized coupling coefficient

_{j}*κ*(

*z*) can be obtained from the coupling strength

_{j}*ρ*by using the following relation [27]:

_{j}The main task now is to obtain the required strength *ρ _{j}* of the instantaneous scattering points from the transmission spectra. If the discretized coupling model is reformulated into the space-time domain, the discrete layer-peeling algorithm (DLP) by using the inverse-scattering method becomes efficient and is briefly introduced as follows (see details in Ref. [13, 16, 27, 31]). We define time delay with respect to the propagation of ordinary wave, i.e., ordinary wave propagates without time delay. The field vectors

*a*

_{o(e),j}(

*τ*) in the time domain are related to the spectral fields by using a discrete-time Fourier transform [13, 16].

*τ*= 0, 1,· · ·,

*N*. The transfer matrix model defined in Eq. (2) can be transformed to time domain by replacing

*A*

_{o(e),j}(

*z*, Δ

*β*) with Eq. (4) [27]. We can obtain the time-domain response in the

*j*-th layer, which reveal the complex coupling strength

*ρ*of QPM grating structure. The time-domain response is given as follows:

_{j}It is known that the time domain field coefficients, *a _{o,j}*(

*τ*) and

*a*(

_{e,j}*τ*) are discrete and correspond to the response of different time delay.

*a*(0) and

_{o,j}*a*(0) denote the two shortest time delay through QPM structure from causality requirements. At the end of QPM structure (layer

_{e,j}*N*), the discrete coupling strength is simply expressed as:

Since the value of *ρ _{N}* is obtained from the output spectra, we can peel this layer off and calculate the impulse responses of layer

*N*– 1 by using transfer matrix method [27] and repeat the procedure until all the layers have been reconstructed. After the coupling strength is known, we can get the coupling coefficient

*κ*(

*z*) of each layer from Eq. (3). The period chirp distribution of QPM grating can be obtained from the derivative of phase of the coupling coefficient by using Eq. (7).

_{j}## 3. Experimental results and discussions

The setup of tunable Šolc-type filter based on QPM structure is shown in Fig. 2(a). A Z-cut MgO doped PPLN crystal is placed in the oven between two crossed polarizers and its temperature is precisely adjusted by using a temperature controller. The temperature instability is less than 0.1°C. The PPLN crystal is fabricated by using electric-poling technique and the dimension is 50 mm (length), 12 mm (width) and 1 mm (thickness). There are one periodic QPM grating (3 mm width) with the nominal period 20.40 *μ*m and two chirped gratings (3 mm width) with period change from 19.90 to 20.40 *μ*m linearly and quadratically, respectively. We use a C + L broadband ASE source with an output wavelength range from approximately 1520 nm to 1610 nm and an optical spectrum analyzer (OSA) to detect the output light. The output spectrum of broadband ASE source used in the experiment is shown in Fig. 2(b). The following experimental transmission spectra of Šolc-type filter are normalized by the spectra of ASE source.

The typical measured transmission spectrum of Šolc-type filter based on uniform QPM grating is shown in Fig. 3(a) as the red dot line. It shows when temperature is 35.5 °C, the center wavelength is measured as 1589.44 nm, which is in good agreement with the theoretical prediction. The central wavelength of Šolc-type filter is tunable due to the temperature-dependent refractive indices [32] and wavelength tuning by changing the temperature is 0.60 nm/°C.

We use the output spectra transmitted from the uniform and chirped QPM gratings to investigate the synthesis of different grating structures by using DLP algorithm. The uniform QPM grating is divided into *N* = 200 layers and the reconstructed period distribution by using DLP method is showed in Fig. 3(b). The result shows a good agreement with the theoretical expectation. The nominal chirped QPM gratings are expressed by Λ(*z*) = Λ_{0} + (Λ* _{N}* − Λ

_{0})(

*z*/

*L*)

*m*, Λ

_{0}, Λ

*is 19.90*

_{N}*μ*m and 20.40

*μ*m, showing the linearly (

*m*= 1) and quadratically (

*m*= 2) chirped structures, respectively. Figures 4(a) and 4(b) show the measured transmission spectra of the linearly and quadratically chirped QPM gratings. We observe band-pass transmission spectra with some fluctuations in the middle wavelength region compared to the theoretical results.

The spectrum fluctuations are attributed to the possible deficiency of chirped QPM structures due to non-uniform electrical field distribution imprinted in QPM gratings during poling process. But the bandwidths of spectra are in good agreement with the theoretical results, which means that the grating periods at the input and output points are same with the nominal grating ones. We reconstruct the grating period distributions by using DLP algorithm and the results are shown in Fig. 5. The reconstructed linearly chirped QPM grating shows better agreement with pre-designed grating distribution than quadratically chirped grating, which indicates that manufacturing accuracy of local QPM period decrease when using the aperiodically grating pattern. The primary results show that it is possible to accurately inspect the period quality of QPM gratings without any etching destroys. We also admit that the output spectra are not fine enough and limit the accuracy of the synthesis on QPM grating.

## 4. Conclusion

We present an approach to retrieve the aperidically QPM grating structures efficiently and accurately by using the discrete layer-peeling algorithm. The measured output spectra of Šolc-type filter based on uniform and chirped QPM structures are used to reconstruct the grating structures, which show that the accuracy of local grating period decrease as the gratings become more complicated. This method has shown a potential application on the QPM grating design or the quality inspection during the fabrication.

## Acknowledgments

This work was supported in part by the National Natural Science Foundation of China ( 60978004, 60937003, and 61132004), Shanghai Municipal Education Commission ( 10YZ14) and Shuguang Program ( 10SG38) and Shanghai Science and Technology Development Funds ( 11XD1402600), SSTCP ( 11511502501-2) and Leading Academic Discipline Project ( S30108), and X.L. Zeng acknowledges the support of Marie Curie IIF( GA-2009-253289).

## References and links

**1. **K. Mizuuchi, K. Yamamoto, M. Kato, and H. Sato, “Broadening of the Phase-Matching Bandwidth in Quasi-Phase-Matched Second-Harmoic Generation,” IEEE J. Quantum Electron. **30**, 1596–1604 (1994). [CrossRef]

**2. **S. N. Zhu, Y. Y. Zhu, and N. B. Ming, “Quasi-phase-matched third-harmonic generation in a quasi-periodic optical superlattice,” Science **278**, 843–846 (1997). [CrossRef]

**3. **M. Charbonneau-Lefort, M. M. Fejer, and B. Afeyan,“Tandem chirped quasi-phase-matching grating optical parametric amplifier design for simultaneous group delay and gain control,” Opt. Lett. **30**, 634–636 (2005). [CrossRef] [PubMed]

**4. **J. Huang, X. P. Xie, C. Langrock, R. V. Roussev, D. S. Hum, and M. M. Fejer, “Amplitude modulation and apodization of quasiphase-matched interactions,” Opt. Lett. **31**, 604–606 (2006). [CrossRef] [PubMed]

**5. **T. Umeki, M. Asobe, T. Yanagawa, O. Tadanaga, Y. Nishida, K. Magari, and H. Suzuki, “Broadband wavelength conversion based on apodized *χ*^{(2)} grating,” J. Opt. Soc. Am. B **26**, 2315–2322 (2009). [CrossRef]

**6. **X. Zeng, S. Ashihara, Z. Wang, T. Wang, Y. Chen, and M. Cha, “Excitation of two-colored temporal solitons in a segmented quasi-phase-matching structure, ” Opt. Express **17**, 16877–16884 (2009). [CrossRef] [PubMed]

**7. **Y. W. Lee, F. C. Fan, Y. C. Huang, B. Y. Gu, B. Z. Dong, and M. H. Chou, “Nonlinear multiwavelength conversion based on an aperiodic optical superlattice in lithium niobate,” Opt. Lett. **27**, 2191–2193 (2002). [CrossRef]

**8. **X. Zeng, X. Chen, F. Wu, Y. Chen, Y. Xia, and Y. Chen, “Second-harmonic generation with broadened flattop bandwidth in aperiodic domain-inverted gratings,” Opt. Commun. **204**, 407–411 (2002). [CrossRef]

**9. **A. M. Schober, G. Imeshev, and M. M. Fejer, “Tunable-chirp pulse compression in quasi-phase-matched second-harmonic generation,” Opt. Lett. **27**, 1129–1131 (2002). [CrossRef]

**10. **K. Beckwitt, F. Ö. Ilday, and F. W. Wise, “Frequency shifting with local nonlinearity management in nonuniformly poled quadratic nonlinear materials,” Opt. Lett. **29**, 763–765 (2004). [CrossRef] [PubMed]

**11. **X. Zeng, S. Ashihara, X. Chen, T. Shimura, and K. Kuroda, “Two-color pulse compression in aperiodically poled lithium niobate,” Opt. Commun. **281**, 4499–4503 (2008). [CrossRef]

**12. **H. Miao, S. Yang, C. Langrock, R. V. Roussev, M. M. Fejer, and A. M. Weiner, “Ultralow-power second-harmonic generation frequency-resolved optical gating using aperiodically poled lithium niobate waveguides,” J. Opt. Soc. Am. B **25**, A41–A53 (2008). [CrossRef]

**13. **J. Shaar, L. G. Wang, and T. Erdogan, “On the synthesis of fiber bragg gratings by layer peeling,” IEEE J. Quantum Electron. **37**, 165–173 (2001). [CrossRef]

**14. **J. Zhang, P. Shum, S. Y. Li, N. Q. Ngo, X. P. Cheng, and J. H. Ng, “Design and fabrication of flat-band long-period grating,” IEEE Photon. Technol. Lett. **15**, 1558–1560 (2003). [CrossRef]

**15. **H. Li, T. Kumagai, and K. Ogusu, “Advanced design of a multichannel fiber Bragg grating based on a layer-peeling method,” J. Opt. Soc. Am. B **21**, 1929–1938 (2004). [CrossRef]

**16. **Y. Choi, J. Chun, and J. Bae, “Numerically extrapolated discrete layer-peeling algorithm for synthesis of nonuniform fiber Bragg gratings,” Opt. Express **19**, 8254–8266 (2011). [CrossRef] [PubMed]

**17. **E. C. Levy and M. Horowitz, “Layer-peeling algorithm for reconstructing the birefringence in optical emulators,” J. Opt. Soc. Am. B **23**, 1531–1539 (2006). [CrossRef]

**18. **X. Chen, J. Shi, Y. Chen, Y. Zhu, Y. Xia, and Y. Chen, “Electro-optic Šolc-type wavelength filter in periodically poled lithium niobate,” Opt. Lett. **28**, 2115–2117 (2003). [CrossRef] [PubMed]

**19. **Y. Q. Lu and Z. L. Wan, “Electro-optic effect of periodically poled optical superlattice LiNbO_{3} and its applications,” Appl. Phys. Lett. **77**, 3719–3721 (2000). [CrossRef]

**20. **Q. Zhang, X. Zeng, F. Pang, X. Chen, and T. Wang, “Tunable polarization-independent Šolc-type wavelength filter based on periodically poled lithium niobate,” Opt. Laser Technol. **44**, 1992–1994 (2012). [CrossRef]

**21. **C. H. Lin, Y. H. Chen, S. W. Lin, C. L. Chang, Y. C. Huang, and J. Y. Chang, “Electro-optic narrowband multi-wavelength filter in aperiodically poled lithium niobate,” Opt. Express **15**, 9859–9866 (2007). [CrossRef] [PubMed]

**22. **X. Gu, X. Chen, Y. Chen, X. Zeng, Y. Xia, and Y. Chen, “Narrowband multiple wavelengths filter in aperiodic optical superlattice,” Opt. Commun. **237**, 53–58 (2004). [CrossRef]

**23. **C. L. Chang, Y. H. Chen, C. H. Lin, and J. Y. Chang, “Monolithically integrated multi-wavelength filter and second harmonic generator in aperiodically poled lithium niobate,” Opt. Express **16**, 18535–18544 (2008). [CrossRef] [PubMed]

**24. **Y. L. Lee, Y. C. Noh, C. S. Kee, N. E. Yu, W. Shin, C. Jung, D. K. Ko, and J. Lee, “Bandwidth control of a Ti:PPLN Šolc filter by a temperature-gradient-control technique,” Opt. Express **16**, 13699–13706 (2008). [CrossRef] [PubMed]

**25. **C. Y. Huang, C. H. Lin, Y. H. Chen, and Y. C. Huang, “Electro-optic Ti:PPLN waveguide as efficient optical wavelength filter and polarization mode converter,” Opt. Express **15**, 2548–2554 (2007). [CrossRef] [PubMed]

**26. **C. S. Kee, Y. L. Lee, and J. Lee, “Electro- and thermo-optic effects on multi-wavelength Šolc filters based on *χ*^{(2)} nonlinear quasi-periodic photonic crystals,” Opt. Express **16**, 6098–6103 (2008). [CrossRef] [PubMed]

**27. **J. K. Brenne and J. Skaar, “Design of grating-assisted codirectional couplers with discrete inverse-scattering algorithms,” J. Lightwave. Technol. **21**, 254–263 (2003). [CrossRef]

**28. **G. Lenz, B. J. Eggleton, and C. R. Giles, “Dispersive properties of optical filters for WDM systems,” IEEE J. Quantum Electron. **34**, 1390–1402 (1998). [CrossRef]

**29. **L. Wang and T. Erdogan, “Layer peeling algorithm for reconstruction of long-period fibre gratings,” Electron. Lett. **37**, 154–156 (2001). [CrossRef]

**30. **T. Erdogan, “Fiber grating spectra,” J. Lightwave Technol. **15**, 1277–1294 (1997). [CrossRef]

**31. **R. Feced and M. N. Zervas, “Efficient inverse scattering algorithm for the design of grating-assisted co-directional mode couplers,” J. Opt. Soc. Am. A **17**, 1573–1582 (2000). [CrossRef]

**32. **O. Gayer, Z. Sacks, E. Galun, and A. Arie, “Temperature and wavelength dependent refractive index equations for MgO-doped congruent and stoichiometric LiNbO_{3},” Appl. Phys. B **91**, 343–348 (2008). [CrossRef]