RESEARCH ARTICLE

Mode overlap analyses of propagated waves in direct bonded PPMgLN ridge waveguide

  • Yujie ZHOU ,
  • Liqun FENG ,
  • Qian HU ,
  • Junqiang SUN
Expand
  • Wuhan National Laboratory for Optoelectronics, College of Optoelectronic Science and Engineering, Huazhong University of Science and Technology, Wuhan 430074, China

Received date: 25 Apr 2011

Accepted date: 17 May 2011

Published date: 05 Sep 2011

Copyright

2014 Higher Education Press and Springer-Verlag Berlin Heidelberg

Abstract

Direct bonded periodically poled MgO doped lithium niobate (PPMgLN) ridge waveguide is a new wavelength converter with high conversion efficiency. The optical field distribution of the ridge waveguide is simulated by employing finite-difference method (FDM), the mode overlap of propagated waves in the ridge waveguide is calculated and the relationship between the overlap coefficient and the waveguide structure sizes is also investigated. The overlap coefficient to difference frequency generation (DFG) process conversion efficiency calculation is firstly introduced.

Cite this article

Yujie ZHOU , Liqun FENG , Qian HU , Junqiang SUN . Mode overlap analyses of propagated waves in direct bonded PPMgLN ridge waveguide[J]. Frontiers of Optoelectronics, 2011 , 4(3) : 343 -347 . DOI: 10.1007/s12200-011-0133-1

Introduction

Quasi-phase-matched (QPM) waveguide devices, employed as wavelength converter based on periodically poled lithium niobate (PPLN) with high efficiency have attracted much attention. Ti diffusion or annealed proton exchange (APE) has been mainly engaged in waveguide fabrication [1,2]. Ti-diffused waveguide exhibits low conversion efficiency due to weak light confinement and photorefractive damage, although low propagation losses are available. On the other hand, waveguides fabricated with APE method can achieve stronger confinement and higher resistance to photorefractive damage than those of Ti-diffused waveguide. However, it exhibits several drawbacks. First, it is polarization-dependent because proton exchange only increases its refractive index for extraordinary wave, propagating with transverse magnetic (TM) mode in a z-cut LiNbO3 substrate. A polarization diversity approach, either with two waveguides or with counter-propagating waves in a single waveguide is needed for APE waveguides to achieve polarization independent operation. Second, excessive proton exchange results in the formation of a “dead layer” whose nonlinearity is degraded. Third, these devices produce asymmetric mode profiles at corresponding wavelengths in wavelength conversion, reducing mode overlap due to the mismatch of mode peak positions.
Direct-bonded ridge-waveguide devices have been studied to overcome these limitations described as above [3-6]. The ridge waveguide neither shows the degradation of the nonlinear coefficient nor photorefractive damage because no ion-exchange process employed in the process. It exhibits strong confinement owing to step-index profile. Therefore, high conversion efficiency can be obtained by using ridge waveguide structures. In addition, this structure confines both transverse electric (TE) and transverse magnetic (TM) modes. It is easy to realize polarization-insensitive wavelength converter using a single waveguide. Therefore, a little of attention has been paid to the optical field distribution of the ridge waveguide. And the mode overlap of the propagated waves in the ridge waveguide is an essential factor related to the conversion efficiency. In this paper we simulated the optical field distribution of the ridge waveguide by employing finite-difference method (FDM). We also calculated the mode overlap of the propagated waves in the ridge waveguide and detailed the overlap coefficient variations along with the structure sizes. We firstly introduced the overlap coefficient to difference frequency generation (DFG) process conversion efficiency calculation.

Ridge waveguide structure

The major structures of ridge waveguide produced by direct bonding consist of three dielectric layers: a ridge region of high refractive index, n1, surrounded by substrate and cladding layers with lower refractive indices n2 and n3 respectively. The cross-section of dimensions W×H of direct-bonded ridge waveguide is illustrated in Fig. 1, where W defines the ridge width and H defines the ridge height. The method of effective dielectric constants originally developed by Knox and Toulios [7] is applied to basic design factors used for direct-bonded ridge waveguide device, from the choice of material layers with suitable refractive indices to the optimization of the geometry and dimensions of the waveguide core. Because ZnO or MgO-doped lithium niobate exhibits strong resistance to photorefractive damage, particular attentions are paid to the examples of ridge dielectric waveguide structures based on periodically poled ZnO doped lithium niobate and periodically poled MgO doped lithium niobate (PPMgLN) as guide layers, MgO doped lithium niobate and LiTaO3 (LT) as substrate layers respectively. The cladding layer is air. Here discussion will be limited to the ridge waveguide based on PPMgLN as guide layer and LT as substrate layer.
Fig.1 Ridge waveguide cross-section

Full size|PPT slide

Mode characteristics of ridge waveguide

Optical field distribution

The modes supported by this waveguide are of two types: quasi-TE mode that polarized parallel to crystal surface and quasi-TM mode that polarized perpendicular to crystal surface. Here discussion will be limited to quasi-TM modes, but a similar methodology can also be applied to find the quasi-TE modes. To minimize the transmission and coupling loss, the fabricated waveguide should allow single-mode propagation of 1550 nm band light.
The waveguide parameters of the waveguide are n1=2.1299@1.55 μm, 2.172@0.775 μm [8], n2 =2.126@1.55 μm, 2.16@0.775 μm [9], n3 =1, and the wavelength of fundamental light is 1.55 μm, second harmonic light is 0.775 μm. According to effective dielectric constants method, it is known that only one mode of fundamental light can propagate in the waveguide when the height ranges 8.4 μm and 13.2 μm with core aspect ratio of width and height (WH) 1∶1. Here we choose the waveguide height H being 9 μm. And in that case, the ridge waveguide in this structure is a multi-mode waveguide for 0.775 μm light. But the periodically poled period is fit for the first-order mode to satisfy the quasi-phase matching condition in PPMgLN. So the first-order mode is the main element of the light transmission in the waveguide. Here we calculate the fundamental mode for each wavelength in the ridge waveguide.
The scalar wave equation evolved from the Maxwell’s equations [10] can be written as
2E(x,y)x2+2E(x,y)y2+(k02n12-β2)E(x,y)=0,
where E represents a wave function and designates a scalar field, n1 is the refractive index of the waveguide core, and β is the propagation constant.
We use FDM to solve scalar wave equations and derive electric field matrix E and the propagation constant [11]. Optical intensity distribution of TM mode for the above examples of lithium-niobate-based waveguide structures is presented in Figs. 2(a) and 2(b), showing the near-field distributions for fundamental wave 1550 nm and second harmonic (SH) wave 775 nm respectively. We obtain the mode sizes at 1/e2 intensities of 8.65 μm (vertical) and 7.05 μm (horizontal) at 1550 nm, and also for values of 7.35 μm (vertical) and 6.95 μm (horizontal) at 775 nm. Due to the step-index profile of the waveguide, the mode size of the fundamental light is almost as the same as that of the second harmonic light. And the effective index of TM mode is 2.12712 at 1550 nm and 2.16887 at 775 nm.
Fig.2 Optical intensity distribution of TM mode for PPMgLN ridge waveguide. (a) Fundamental wave; (b) second harmonic wave

Full size|PPT slide

Mode overlap

The overlap coefficient between the fundamental light and the second harmonic light is given by
Γ=|E1*(x,y)E2(x,y)dxdy|2|E1(x,y)E1*(x,y)dxdy||E2(x,y)E2*(x,y)dxdy|,
where E1(x,y) and E2(x,y) denote the electric field distributions of different propagated waves in the waveguide
Figures 3(a) and 3(b) show the normalized optical intensity overlap along the width and height direction respectively. It is obvious that the intensity overlap along the width direction is greater than that along the height direction due to symmetric geometry. It is proven that an asymmetric geometry can lead to a poor overlap of the optical intensity profiles of a propagating mode at different wavelengths.
Fig.3 Normalized optical intensity along width (a) and height (b) direction

Full size|PPT slide

We calculate a value of 92.75% for the power overlap integral between the 775 and 1550 nm light. Figure 4 shows the power overlap coefficient versus the core width W between the fundamental light and the second harmonic light when the core height is 9 μm. The power overlap coefficient as a function of the core height when the core width retains 9 μm is presented in Fig. 5. We find that the overlap coefficient increases with the core width and height. The growth is faster with the core height and smaller with the width. It is due to the asymmetric geometry in height direction which induces poor overlap of the power intensity. Figure 6 demonstrates the power overlap coefficient between the fundamental light and the second harmonic light when the core height and width of the waveguide are varied. The maximal overlap coefficient is 0.96 at the waveguide single-mode condition.
Fig.4 Power overlap coefficient curve along with core width

Full size|PPT slide

Fig.5 Power overlap coefficient curve along with core height

Full size|PPT slide

Fig.6 Power overlap coefficient plots with ridge width and ridge height

Full size|PPT slide

DFG process conversion efficiency

It was recognized that multi domain ferroelectrics could enhance the efficiency of nonlinear interactions, due to the sign change in the nonlinear susceptibility accompanying domain reversal. Under the assumption of validity of Kleinman’s symmetry, the approximation of slowly varying envelope, the Maxwell equations for electric fields at the three frequencies can be reduced to the following coupled-mode equations with the Fourier components [12]:
Epz=iωpΓsinpcdeffEsEiexp(-i(Δβ-2π/Λ)z)-0.5αpEp,
Esz=iωsΓpinscdeffEpEi*exp(i(Δβ-2π/Λ)z)-0.5αsEs,
Eiz=iωiΓpsnicdeffEpEs*exp(i(Δβ-2π/Λ)z)-0.5αiEi,
where c is the speed of light in the vacuum; nj and Ej are the index of refraction, and the electric field under light-frequencies ωj (j = p,s,i; and ωp, ωs, ωi denote the pump, signal and idler waves), respectively. Δββ = βp-βs-βi) represents the propagation constant mismatching, Λ is the grating period. Equations (3a)–(3c) are the expressions described the QPM DFG process. Гsi is the mode overlap of the signal light and the idler light, Гpi is the mode overlap of the pump light and the idler light, Гps is the mode overlap of the pump light and the signal light. αp, αs, αi represent the waveguide loss of pump, signal and idler lights, respectively.
In order to access the largest nonlinear coefficient in LiNbO3, d33, the periodically poled orientation should follow its optic axis [001] direction (z-direction) and the polarizations of the optical waves in the waveguide should parallel to the domain inversions. Here, the waveguide layer and the substrate layer are both z-cut, and the waveguide periodical pole structure is along the z-direction, as shown in Fig. 7. The optical waves propagate along the [100] direction (x-direction). And the input waves are z-polarized, the modes existing in the waveguide are TM mode.
Fig.7 PPMgLN ridge waveguide

Full size|PPT slide

Assuming the input wavelengths of the pump light and the signal light are 775 and 1545 nm, the generated idler wavelength is 1555 nm. The intrinsic waveguide losses were estimated to be 0.35 dB/cm for 1550 nm band light and 0.7 dB/cm for 775 nm band light. In case of a 9-μm-wide, 9-μm-high ridge, the effective indexes are 2.16887, 2.12732, and 2.12701 for the pump light, the signal light and the idler light respectively, and the poling period is 18.62 μm. The mode overlap of the propagated waves Гsi, Гpi, Гps are 0.9996, 0.9166, 0.9269, respectively.
We calculated the conversion efficiency of the waveguide by measuring the DFG characteristics. The power of the input pump and signal waves are 100 and 1 mW, the length of the waveguide is 50 mm. Figure 8 shows a typical DFG tuning curve for the PPMgLN waveguide considering the mode overlap of the propagated waves. We successfully obtained the DFG conversion efficiency of 260% at a pump wavelength of 775 nm.
Fig.8 DFG conversion efficiency for a 50mm-long waveguide

Full size|PPT slide

Conclusion

We firstly calculate the mode overlap of the propagated waves in the ridge waveguide and discussed the overlap coefficient relationship with the structure sizes. We also introduced the overlap coefficient to the DFG process conversion efficiency calculation, which can help us simulate the propagated waves interaction in the direct bonded waveguide better.

Acknowledgements

This work was supported by the National High Technology Research and Development Program of China (No. 2009AA03Z410)
1
LeeY L, SucheH, MinY H, LeeJ H, GrundkotterW, QuiringV, SohlerW. Wavelength- and time-selective all-optical, channel dropping in periodically poled Ti:LiNbO3 channel waveguides.IEEE Photonics Technology Letters, 2003, 15(7): 978–980

DOI

2
BortzM L, FejerM M. Annealed proton-exchanged LiNbO3 waveguides.Optics Letters, 1991, 16(23): 1844–1846

DOI PMID

3
GawithC B E, ShepherdD P, AbernethyJ A, HannaD C, RossG W, SmithP G R. Second-harmonic generation in a direct-bonded periodically poled LiNbO3 buried waveguide.Optics Letters, 1999, 24(7): 481–483

DOI PMID

4
UmekiT, AsobeM, NishidaY, TadanagaO, MagariK, YanagawaT, SuzukiH. Highly efficient+5-dB parametric gain conversion using direct-bonded PPZnLN ridge waveguide.IEEE Photonics Technology Letters, 2008, 20(1): 15–17

DOI

5
UmekiT, TadanagaO, AsobeM. Highly efficient wavelength converter using direct-bonded PPZnLN ridge waveguide.IEEE Journal of Quantum Electronics, 2010, 46(8): 1206–1213

DOI

6
AsobeM, MiyazawaH, TadanagaO, NishidaY, SuzukiH. A highly damage-resistant Zn:LiNbO3 ridge waveguide and its application to a polarization-independent wavelength converter.IEEE Journal of Quantum Electronics, 2003, 39(10): 1327–1333

DOI

7
KnoxR M, TouliosP P. Integrated circuits for the millimeter through optical frequency range. In: Proceedings of MRI Symposium on Submillimeter Waves. Polytechnic Press,1970, 497

8
ZelmonD E, SmallD L, JundtD. Infrared corrected sellmeier coefficients for congruently grown lithium niobate and 5 mol.% magnesium oxide-doped lithium niobate.Optical Society of America B,1997, 14(12): 3319–3322

9
AbedinK S, ItoH. Temperature-dependent dispersioon relation of ferroelectric lithium tantalate. Journal of Applied Physics, 1996, 80(11): 6561–6563

10
KatsunariO. Fundamentals of Optical Waveguides. 2nd<Eds/>ed. Japan: Elsevier, 2005

11
KawanoK, KitohT. Introduction to Optical Waveguide Analysis: Solving Maxwell’s Equations and the Schrödinger Equation. Japan: John Wiley and Sons Inc., 2001

12
SuharaT, NishiharaH. Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped gratings. IEEE Journal of Quantum Electronics, 1990, 26(7): 1265–1276

DOI

Outlines

/