Abstract

In integrated optics the radiation modes represent a negative aspect regarding the propagation of guided modes. They characterize the losses of the substrate region but can contribute to enhance the guided modes by considering the coupling through properly designed gratings arranged at the core/substrate interface. By tailored gratings, the radiation modes become propagating modes and increase the guided power inside the waveguide guiding region. This enhancement is useful especially in low intensity processes such as second harmonic χ(2) conversion process. For this purpose, we analyze accurately the radiation modes contribution in a χ(2) GaAs/AlGaAs nonlinear waveguide where second harmonic signal is characterized by a low power intensity. This analysis considers a new design approach of multiple grating which enhances a fundamental guided mode at λFU =1.55 μm and a codirectional second harmonic guided mode at λSH =0.775 μm. In particular we analyze the second harmonic conversion efficiency by studying the coupling effect of three gratings. The combined effects of the gratings provide an efficient second harmonic field conversion. Design considerations, based on the coupled mode equations analysis, are theoretically discussed. A good agreement between analytical and numerical results is observed.

© 2009 Optical Society of America

1. Introduction

Studies of χ(2) and the quasi-phase matching (QPM) technique promote research works on nonlinear optic (NLO) devices for ultrafast signal processing based on second-order nonlinearity. After the development of high-efficiency QPM- second harmonic generation (SHG) devices using ferroelectric crystal waveguides, the research was directed to implementation of a compact and efficient coherent light source by combining the SHG device with a semiconductor laser. Moreover a strong motivation of research on NLO devices has been the need for the development of all-optical wavelength converters for dense wavelength division multiplexing (DWDM) optical communication systems. In particular the second harmonic generation in χ(2) nonlinear process represents a good solution for DWDM applications in telecommunication systems. But this process requires high intensity input power and large interaction lengths for a good second harmonic (SH) conversion. For this reasons the introduction of tailored gratings which enhance the SH field is important.

The novel design approach proposed in this work takes into account the multiple grating effects which enhance the guided SH signal along the symmetrical waveguide shown in Fig. 1. The new accurate analysis concerns the coupling of the radiation modes with the propagating fundamental and SH modes: properly designed gratings, arranged at the core/substrate interface, generate scattered waves which increase the guided modal power inside the waveguide guiding region. By tailored gratings, the radiation modes supply energy to the low intensity SH mode by providing an efficient SH conversion in short structures. In the presented structure three grating are considered: the first grating (grating 1 of Fig. 1), considered at one interface between core and substrate, is designed by quasi phase matching QPM condition. The QPM technique is a practical method of substantially increasing the second-harmonic power by effectively reducing the phase mismatch between the fundamental and the SH fields: the technique relies on introducing the grating 1 in the waveguide to compensate for the difference of the propagation refractive indexes in the waveguide.

A second grating (grating 2 of Fig. 1), placed at second the interface between the core and the substrate region, increases the fundamental power through the coupling between the substrate modes and the fundamental guided mode. On the same second interface, after the grating 2, we consider a third grating (grating 3 of Fig. 1) which enhances the power transferred to the SH mode through the coupling between the substrate modes and the SH guided mode by optimizing the conversion efficiency.

The substrate modes are excited through the diffraction effect of the grating 2 and 3 of Fig. 1 which supply the energy to both the propagating fields through the coupled energy of the radiation modes. The coupling effect is analyzed by combining the coupled mode equations systems related to the three grating. The novelty of this work is in the combined effect of the three gratings applied to a symmetrical GaAs/AlGaAs waveguide typically used in NLO devices [1]–[3]. In particular we analyze the enhancement of a SH signal (λSH= 0.775 μm) generated by a codirectional fundamental mode (λFU= 1.55 μm). The goal of the presented work is to find the optimum geometrical configuration which overcomes the problems associated to the total waveguide length such as reflections and losses, by providing a good SH field conversion. With the introduction of the multiple grating configuration it is possible to decrease the total waveguide length by obtaining a good conversion efficiency according to real input powers.

We resume the analysis of this work in the following steps: i) we analyze the modes involved in the designed structure including substrate modes and diffraction theory; ii) we study the coupled mode equations systems by combining the effect of each grating; iii) we provide the design criteria steps of the periodic structure; iv) finally we model the tailored waveguide through the time domain Hertzian potential formulation (HPF) used to solve nonlinear χ(2) processes [4], and by the finite difference time domain (FDTD) modelling. A good agreement between numerical and analytical results is checked.

 

Fig. 1. Multiple gratings design in nonlinear optical symmetrical waveguide.

Download Full Size | PPT Slide | PDF

2. Modes of the discontinuous periodic waveguide.

As shown in Fig. 2, the symmetrical waveguide is characterized by a discrete set of orthogonal guided modes and by radiation modes. The radiation field (related to the radiation modes) thus qualifies in all aspects as a mode, except that is not confined to the waveguide core but reaches undiminished to infinite distance in x direction normal to the core. The modes of this type are called radiation modes. Their propagation constant β are not constrained to a discrete set of values, since they are related to the angle of the incident can be chosen arbitrarily. The values of the propagation constant thus form a continuum, so that we also speak of the radiation modes as modes of the continuum. The guided and radiation modes form a complete orthogonal set of modes. Any field distribution can be expressed by series expansion into these modes. All the modes are classified as transverse electric (TE) and as transverse magnetic (TM) modes. According to the assumption of TE polarized source, we limit the discussion to TE waves.

For the transverse electric (TE) field with no variation along the y-direction, only three field components exists, and they are related as follows

Hx(x,z)=1jωμzEy(x,z)
Hz(x,z)=1jωμxEy(x,z)

where Ey(x,z) may be expressed as a modal expansion [5]

Ey(x,z)={kakψk(x)+0b(kx)φ(x;kx)dkx}ejβz

and ak and b(kx) are the amplitudes of the guided and continuum modes respectively. In the case of two TE guided modes (fundamental and SH modes) the guided field are represented by [5]

Eyω=aωψω(x)ejβωz
Ey2ω=a2ωψ2ω(x)ejβ2ωz

with

ψω,2ω(x)={cos(σω,2ωx)corecos(σω,2ωd/2)exp[Δω,2ω(xd/2)]substrate

and

aω,2ω=d1+1Δω,2ω

where the parameters σω,2ω and Δω,2ω are defined by the following equations

Δω,2ω=(n22(kω,2ω)2(βω,2ω)2)1/2
σω,2ω=(n12(kω,2ω)2(βω,2ω)2)1/2
 

Fig. 2. Modes in symmetrical planar waveguide.

Download Full Size | PPT Slide | PDF

The radiation modes are necessary to describe radiation phenomena in region around the core. Imperfections, such as step discontinuities and sinusoidal core interface deformations, cause some of the guided mode power to radiate into the space outside and inside the core. The diffraction theory [6] of plane wave scattering from a sinusoidally deformed dielectric interface explains the radiation modes excitation mechanism. We assume that only two TE modes (fundamental and SH mode) are propagated inside the guiding core. The diffraction theory is applied at interface related to the grating 2 and 3 of Fig. 1. The theory is based on approximate plane wave analysis [7] and on the sinusoidal approximation of one core-substrate interface [6],[7]. The function f(z) which represents the sinusoidal interface deformation is given as

f(z)=tsin(2πΛ2,3z)

As reported in Fig. 3 the sinusoidally deformed dielectric interface scatters the incident radiation. The arrows shown in the Fig. 3 indicate the propagation vector of the incident, reflected, and transmitted plane waves labeled i, r, and t. The shorter arrows labeled 1 and -1 indicate scattered plane waves. Two scattered waves emerge on either side of the interface. One set of scattered waves accompanies the reflected wave in the core region and another set accompanies the transmitted wave in the substrate region. The angle ϕ of the direction of each wave indicated in Fig. 3 is obtained from the equation

cos(ϕ)=β/nikω,2ω

where the propagation constant β indicates either, βω,2ω, β (+), or β (-), ni stands for n1 and n2, and kω,2ω= 2π/λFU,SH. We observe that if the period Λ2,3 satisfies the substrate mode condition

2πΛ2,3±βω,2ω<nikω,2ω

for i=1, and i=2, propagating scattered waves will exist. Far from the surface only the zero-order reflected and transmitted waves can be observed. Only in the immediate vicinity of the distorted surface is there any field distortion. Though the grating, this distorted field exchanges power with the zero-order field by increasing the guided power. The field of zero-order are the incident i, reflected r, and transmitted t plane waves corresponding to the fundamental and SH fields. In particular in the TE case the zero-order field is given by

eyω,2ω=Aiexp[j(σω,2ωx+βω,2ωz)]+Arexp[j(σω,2ωxβω,2ωz)]core region
eyω,2ω=Atexp[j(Δω,2ωx+βω,2ωz)]substrateregion

The total scattered field includes in each region the contributions of the scattered field a + γi and a - γi of Fig. 3, and is represented by [7]

EyS=t2Aiσω,2ω(Δω,2ωσω,2ω)[exp[j(σ()xβ()z)]σ()+Δ()exp[j(σ(+)xβ(+)z)]σ(+)+Δ(+)]
coreregion
EyS=t2Aiσω,2ω(Δω,2ωσω,2ω)[exp[j(Δ()x+β()z)]σ()+Δ()exp[j(Δ(+)x+β(+)z)]σ(+)+Δ(+)]
substrateregion

with the definitions [7]

β(±)=βω,2ω±2πΛ2,3
σ(±)=(n12(kω,2ω)2(β(±))2)1/2
Δ(±)=(n22(kω,2ω)2(β(±))2)1/2

where (15) indicates the phase matching between the zero order field and the scattered waves.

 

Fig. 3. The sinusoidally distorted core-substrate interface acts like a phase grating. The incident i, reflected r, and diffracted rays are shown. The sinusoidal approximation is applied to the grating 2 and 3 of Fig. 1.

Download Full Size | PPT Slide | PDF

3. Coupled mode theory.

In our analysis three coupled equations system are considered. The first set of coupled equations is related to grating 1 of Fig. 1. The function of this grating is to couple the fundamental field with the SH field. In particular by indicating with Aω the amplitude of the pump fundamental field and with A the amplitude of the SH wave we obtain [8],[9]

{dAω(z)dz+j(2kLωcosKz)Aω(z)=j[kNLexp(j2δz)]*[Aω(z)]*A2ω(z)dA2ω(z)dz+j(2kL2ωcosKz)Aω(z)=jkNLexp(j2δz)[Aω(z)]2

with

kLω=ωε04[Eyω(x)]*ΔεEyω(x)dx
kL2ω=2ωε04[Ey2ω(x)]*ΔεEy2ω(x)dx
kNL=2ωε04[Ey2ω(x)]*χ˜NL[Eyω(x)]2dx

where

2δ=β2ω(2βω+K),K=2π/Λ1

is the phase mismatch between the fundamental and SH field,

Δε=(ε1ε2)sin(aπ)πexp[jKz]

is the first term of the permittivity Fourier series, and finally

χ˜NL=(χ1(2)χ2(2))sin(aπ)πexp[jKz]

is the first term of the nonlinear coefficient Fourier series, and a denotes the duty ratio of the period.

The grating 2 of Fig. 1 enhances the low pump fundamental field coupled into the waveguide. In this case coupling between the fundamental field and the scattered fields (generated by grating 2) is considered. In particular the equations to solve are the following [10]

{dAω(z)dz=jiaγi±(z)kγi,Aω*exp[j(β±(βω2π/Λ2)z]dβ±daγi±(z)dz=jAω(z)kγi,Aω*exp[j(β±(βω2π/Λ2)z]

with the coupling coefficient defined as

kγi,Aω=kSω=Eyω(x)ΔεEyS(x)dx

Finally the grating 3 optimize the SH conversion by coupling the SH field with the scattered field. In particular we take into account the system

{dA2ω(z)dz=jiaγi±(z)kγi,A2ω*exp[j(β±(βω2π/Λ3)z]dβ±daγi±(z)dz=jA2ω(z)kγi,A2ω*exp[j(β±(βω2π/Λ3)z]

where

kγi,A2ω=kS2ω=Ey2ω(x)ΔεEyS(x)dx

4. Design and results

We assume that the fundamental pump source (λFU=1.55 μm) is a polarized TE field. The GaAs core (n1FU=1.55 μm)=3.374 and n1SH=0.755 μm)=3.691) is characterized by a thickness d=0.22 μm and supports only a TE fundamental mode at λFU=1.55 μm and only a SH TE mode at λFU=0.775 μm. In this way all the TE power is matched with the two propagating modes and the modal dispersion is low. In order to minimize the reflections along the z- propagating the direction and to conserve the single mode condition, the parameters h and t of Fig. 1 are fixed to the low value of 0.025 μm. The substrate material is Al0.4Ga0.6As (n2FU=1.55 μm)=3.2 and n2SH=0.755 μm)=3.4). The GaAs core, combined with the Al0.4Ga0.6As, characterizes the periodically switched nonlinearity (PSN) [1]. For the PSN effect to occur, a medium is required in which susceptibility coefficient χ(2) is modulated periodically along the direction of light propagation. We observe that in a generic case χ(2) is a tensor, however most semiconductor, which crystallize in zinc-blende structures, have a symmetry and their second-order susceptibility has a single nonzero independent component (for GaAs at the transparency region below the optical energy gap χx,y,z (2)= 200 pm/V). The grating 1 couples the fundamental mode with the SH one. The period Λ1=2.31*10-6 μm satisfies the QPM condition (phase mismatch δ=0):

2βω+2πΛ1=β2ωorΛ1=(λFU)/(N2ωNω)

where N ω,2ω are the effective refractive indices of both modes. The duty ratio a is assumed equal to ½. In Fig. 4 we show the normalized coupling coefficient kNL versus z obtained by applying the QPM condition. Usually the input power is of the order of 1–100 mW [1], but the SH conversion requires high intensity fundamental power. For this purpose the grating 2 is designed in order to enhance the fundamental power inside the waveguide through the coupling with the scattered field generated by the fundamental mode Aω(z). In this case the set of equations to solve are (18) and (25) with the initial conditions Aω(0)=(P0)1/2 and A2m(0)= 0 where P0 is the input pump power. Moreover the grating 3 is added in order to optimize the SH conversion efficiency by increasing the SH coupled power. Concerning the grating 3 the set of equations to solve are (18) and (27) with the initial conditions Aω(0)= Aω(L1) and A2m(0)= A2m(L1). The period Λ2 and Λ3 satisfy the propagation condition (10) and provide a constructive interference given by the phase matching condition (15) between the incident and the scattered field. In particular we choose a strong coupling condition by analyzing the coupling coefficient versus the period and versus the z-direction. Figure 5 shows that regular coupling coefficients kSω,2ω (strong coupling) are obtained with Λ2=0.95*10-6 μm and Λ3=1*10-6 μm. The plots of Fig. 5(a) and Fig. 5(b) are obtained by evaluating the scattered fields EyS (see Fig. 6(a) and Fig. 6(b)) generated by the fundamental and SH modes. Figures 7 and 8 illustrates the normalized substrate coupling coefficients kSω and kS versus z-propagating direction. The SH conversion efficiency is given by η(z)=|A(z)|2/|Aω(0)2| [8]. Figure 9 shows the analytical solution of the conversion efficiency η by considering the grating 1, the grating 1 with the effect of the grating 2, and finally the effect of the three grating together by fixing a reference length of L1=2mm. The input power used during the calculus is P0 =100 mW. A conversion efficiency of about 70 % is obtained after a length of z=10 mm. It is clear from the Fig. 9 that the grating 3 increases the efficiency of about 20 % (optimization) respect to the case of grating 1 combined with the grating 2. In this way is possible to reduce the total length of the device by overcoming the problems of losses and reflections along the propagating direction. Concerning the total length of the structure, we observe that the theoretical effective interaction length [8] Leff≅2.4(Λ1/2a)1/2 =3.6 mm not considers the variation of the fundamental power along the longitudinal direction, and not takes into account the grating parameters h and t which characterize the power reflected on each step discontinuity, thus, in our practical case, a high interaction efficiency is obtained with grating lengths larger than Leff. Moreover the choice of the parameter h=t=0.025 μm is justified by the following considerations. As reported in Fig. 10, low reflection coefficients [5] Γω,2ω=(βd ω,2ω - βD ω,2ω)/(βd ω,2ω + βD ω,2ω) at each step discontinuity are obtained in the single SH mode region by decreasing the values of h. But a strong coupling is obtained in the single mode SH region by increasing the h values (see Fig. 11 where are reported the coupling coefficients versus h). Figure 12 justifies the final choice of h=t=0.025 μm: the figure analyzes the numerical and analytical efficiency η versus z for different h values by showing a good convergence between numerical HPF, numerical two dimensional (2D) FDTD, and analytical results. A good matching between the HPF and 2D FDTD spectra is observed in Fig. 13 where we calculate the discrete Fourier transform (DFT). The convergence of the solutions confirms the accuracy of the presented theoretical model. The numerical simulations are performed by parallel calculus with time steps Δt=1.66*10-6 sec., dimension of spatial cell Δz=0.05 μm, and sinusoidal source at λ0= 1.55 μm given by

E0fundamental=sin(ωfundamental·t·Δt)

The used numerical boundary conditions are the absorbing boundary conditions (ABCs) [13].

 

Fig. 4. Grating 1: normalized nonlinear coupling coefficient versus z.

Download Full Size | PPT Slide | PDF

 

Fig. 5. Substrate coupling coefficient: (a) coupling coefficient versus Λ2; (b) coupling coefficient versus Λ3. Strong couplings are obtained with Λ2 = 0.95 μm and Λ3 = 1 μm, where the coupling coefficients are regular respect to the z- direction.

Download Full Size | PPT Slide | PDF

 

Fig. 6. (a) Scattered field generated by the fundamental mode. (b) Scattered field generated by the SH field. The coordinate system is the same of Fig. 1.

Download Full Size | PPT Slide | PDF

 

Fig. 7. Grating 2: substrate coupling coefficient versus L1.

Download Full Size | PPT Slide | PDF

 

Fig. 8. Grating 3: substrate coupling coefficient versus L2.

Download Full Size | PPT Slide | PDF

 

Fig. 9. Analytical and HPF numerical results of the efficiency SH conversion versus z.

Download Full Size | PPT Slide | PDF

 

Fig. 10. Reflection coefficients Γω,2ω=(βd ω,2ω - βD ω,2ω)/(βd ω,2ω + βD ω,2ω) of the propagating modes at each step discontinuity. Inset: schematic diagram of a step discontinuity. The propagation constants βω,2ω d and βω,2ω D are related to region with core thickness d and D, respectively.

Download Full Size | PPT Slide | PDF

 

Fig. 11. Coupling coefficients versus h calculated by assuming z=0.

Download Full Size | PPT Slide | PDF

 

Fig. 12. Structure with grating 1, 2 and 3: analytical, 2D-FDTD and HPF results of the efficiency SH conversion versus z for different h values. The analytical results are averaged values.

Download Full Size | PPT Slide | PDF

 

Fig. 13. Comparison between HPF and 2D-FDTD normalized spectra of the fundamental and second harmonic modes. Inset: time evolution of the fundamental and second harmonic normalized fields generated inside the waveguide of Fig. 1.

Download Full Size | PPT Slide | PDF

5. Conclusion

We presented in this work a new theoretical model for low intensity linear and nonlinear processes. The novel accurate model analyzes the effect of the radiation modes coupling on modal conversion efficiency. In particular we applied the theory in a nonlinear GaAs/AlGaAs waveguide with discontinuous core-substrate interfaces by providing design criteria for second harmonic enhancement through the use of multiple gratings. A good agreement between analytical and numerical results proves the second harmonic enhancement obtained through the radiation modes interaction. The design approach can be applied also to complex 3D optical waveguides such as ridge and channel waveguides which behave as multiplexers or polarizers. Technological and experimental works are now under study.

References and links

1. E. U. Rafailov, P. L. Alvarez, C. T. A. Brown, W. Sibbett, R. M. De la Rue, P. Millar, D. A. Yanson, J. S. Roberts, and P. A. Houston, “Second-harmonic generation from a first-order quasi-phase-matched GaAs/AlGaAs waveguide crystal,” Opt. Lett. 26, 1984–1986 (2001).. [CrossRef]  

2. X. Yu, L. Scaccabarozzi, J. S. Harris, P. S. Kuo, and M. M. Fejer, “Efficient continuous wave second harmonic generation pumped at 1.55 μm in quasi-phase-matched AlGaAs waveguides,” Opt. Express 13, 10742–10748 (2005). [CrossRef]   [PubMed]  

3. X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, “Growth of GaAs with orientation-patterned structures for nonlinear optics,” J. Cryst. Growth 301, 163–167 (2007). [CrossRef]  

4. A. Massaro, V. Tasco, M. T. Todaro, R. Cingolani, M. De Vittorio, and A. Passaseo, “Scalar time domain modeling and coupling of second harmonic generation process in GaAs discontinuous optical waveguide,” Opt. Express 16, 14496–14511 (2008). [CrossRef]   [PubMed]  

5. T. Rozzi and M. Mongiardo, Open Electromagnetic Waveguides, (IEE Electromagnetic Waves Series 43, London1997). [CrossRef]  

6. D. Marcuse, Theory of Dielectric Optical Waveguides (Academic Press, New York1974).

7. D. Marcuse, “Hollow dielectric waveguides for distributed feedback lasers,” IEEE J. Quantum Electron. 26, 1265–1276 (1972).

8. T. Suhara and M. Fujimura, Waveguide Nonlinear-Optic Devices (Berlin: Springer, 2003).

9. T. Suhara and H. Nishihara, “Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped grating,” IEEE J. Quantum Electron. 8, 661–669 (1972).

10. S. Ura, S. Murata, Y. Awtsuji, and K. Kintaka, “Design of resonance grating coupler,” Opt. Express 16, 12207–12213 (2008). [CrossRef]   [PubMed]  

11. A. Massaro and T. Rozzi, “Rigorous time-domain analysis of dielectric optical waveguides using Hertzian potentials formulation,” Opt. Express 14, 2027–2036 (2006). [CrossRef]   [PubMed]  

12. A. Taflove and S. C. Hagness, Computational Electrodynamic: the Finite-difference Time-domain Method, (Arthec House Publishers, sec. ed., London2000), ch. 2,3,4,7.

13. G. Mur, “Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations,” IEEE Trans. Electromagn. Compat. 23, 377–382 (1981). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. E. U. Rafailov, P. L. Alvarez, C. T. A. Brown, W. Sibbett, R. M. De la Rue, P. Millar, D. A. Yanson, J. S. Roberts, and P. A. Houston, "Second-harmonic generation from a first-order quasi-phase-matched GaAs/AlGaAs waveguide crystal," Opt. Lett. 26, 1984-1986 (2001)..
    [CrossRef]
  2. X. Yu, L. Scaccabarozzi, J. S. Harris, P. S. Kuo, and M. M. Fejer, "Efficient continuous wave second harmonic generation pumped at 1.55 ?m in quasi-phase-matched AlGaAs waveguides," Opt. Express 13, 10742-10748 (2005).
    [CrossRef] [PubMed]
  3. X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
    [CrossRef]
  4. A. Massaro, V. Tasco, M. T. Todaro, R. Cingolani, M. De Vittorio, and A. Passaseo, "Scalar time domain modeling and coupling of second harmonic generation process in GaAs discontinuous optical waveguide," Opt. Express 16, 14496-14511 (2008).
    [CrossRef] [PubMed]
  5. T. Rozzi and M. Mongiardo, Open Electromagnetic Waveguides, (IEE Electromagnetic Waves Series 43, London 1997).
    [CrossRef]
  6. D. Marcuse, Theory of Dielectric Optical Waveguides (Academic Press, New York 1974).
  7. D. Marcuse, "Hollow dielectric waveguides for distributed feedback lasers," IEEE J. Quantum Electron. 26, 1265-1276 (1972).
  8. T. Suhara, and M. Fujimura, Waveguide Nonlinear-Optic Devices (Berlin: Springer, 2003).
  9. T. Suhara, and H. Nishihara, "Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped grating," IEEE J. Quantum Electron. 8, 661-669 (1972).
  10. S. Ura, S. Murata, Y. Awtsuji, and K. Kintaka, "Design of resonance grating coupler," Opt. Express 16, 12207-12213 (2008).
    [CrossRef] [PubMed]
  11. A. Massaro, and T. Rozzi, "Rigorous time-domain analysis of dielectric optical waveguides using Hertzian potentials formulation," Opt. Express 14, 2027-2036 (2006).
    [CrossRef] [PubMed]
  12. A. Taflove, S. C. Hagness, Computational Electrodynamic: the Finite-difference Time-domain Method, (Arthec House Publishers, sec. ed., London 2000), ch. 2,3,4,7.
  13. G. Mur, "Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations," IEEE Trans. Electromagn. Compat. 23, 377-382 (1981).
    [CrossRef]

2008

2007

X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
[CrossRef]

2006

2005

2001

1981

G. Mur, "Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations," IEEE Trans. Electromagn. Compat. 23, 377-382 (1981).
[CrossRef]

1972

D. Marcuse, "Hollow dielectric waveguides for distributed feedback lasers," IEEE J. Quantum Electron. 26, 1265-1276 (1972).

T. Suhara, and H. Nishihara, "Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped grating," IEEE J. Quantum Electron. 8, 661-669 (1972).

Alvarez, P. L.

Awtsuji, Y.

Brown, C. T. A.

Cingolani, R.

De la Rue, R. M.

De Vittorio, M.

Fejer, M. M.

X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
[CrossRef]

X. Yu, L. Scaccabarozzi, J. S. Harris, P. S. Kuo, and M. M. Fejer, "Efficient continuous wave second harmonic generation pumped at 1.55 ?m in quasi-phase-matched AlGaAs waveguides," Opt. Express 13, 10742-10748 (2005).
[CrossRef] [PubMed]

Harris, J. S.

X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
[CrossRef]

X. Yu, L. Scaccabarozzi, J. S. Harris, P. S. Kuo, and M. M. Fejer, "Efficient continuous wave second harmonic generation pumped at 1.55 ?m in quasi-phase-matched AlGaAs waveguides," Opt. Express 13, 10742-10748 (2005).
[CrossRef] [PubMed]

Houston, P. A.

Kintaka, K.

Kuo, P. S.

Lin, A. C.

X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
[CrossRef]

Marcuse, D.

D. Marcuse, "Hollow dielectric waveguides for distributed feedback lasers," IEEE J. Quantum Electron. 26, 1265-1276 (1972).

Massaro, A.

Millar, P.

Mur, G.

G. Mur, "Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations," IEEE Trans. Electromagn. Compat. 23, 377-382 (1981).
[CrossRef]

Murata, S.

Nishihara, H.

T. Suhara, and H. Nishihara, "Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped grating," IEEE J. Quantum Electron. 8, 661-669 (1972).

Passaseo, A.

Rafailov, E. U.

Roberts, J. S.

Rozzi, T.

Scaccabarozzi, L.

X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
[CrossRef]

X. Yu, L. Scaccabarozzi, J. S. Harris, P. S. Kuo, and M. M. Fejer, "Efficient continuous wave second harmonic generation pumped at 1.55 ?m in quasi-phase-matched AlGaAs waveguides," Opt. Express 13, 10742-10748 (2005).
[CrossRef] [PubMed]

Sibbett, W.

Suhara, T.

T. Suhara, and H. Nishihara, "Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped grating," IEEE J. Quantum Electron. 8, 661-669 (1972).

Tasco, V.

Todaro, M. T.

Ura, S.

Yanson, D. A.

Yu, X.

X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
[CrossRef]

X. Yu, L. Scaccabarozzi, J. S. Harris, P. S. Kuo, and M. M. Fejer, "Efficient continuous wave second harmonic generation pumped at 1.55 ?m in quasi-phase-matched AlGaAs waveguides," Opt. Express 13, 10742-10748 (2005).
[CrossRef] [PubMed]

IEEE J. Quantum Electron.

D. Marcuse, "Hollow dielectric waveguides for distributed feedback lasers," IEEE J. Quantum Electron. 26, 1265-1276 (1972).

T. Suhara, and H. Nishihara, "Theoretical analysis of waveguide second-harmonic generation phase matched with uniform and chirped grating," IEEE J. Quantum Electron. 8, 661-669 (1972).

IEEE Trans. Electromagn. Compat.

G. Mur, "Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations," IEEE Trans. Electromagn. Compat. 23, 377-382 (1981).
[CrossRef]

J. Cryst. Growth

X. Yu, L. Scaccabarozzi, A. C. Lin, M. M. Fejer, and J. S. Harris, "Growth of GaAs with orientation-patterned structures for nonlinear optics," J. Cryst. Growth 301, 163-167 (2007).
[CrossRef]

Opt. Express

Opt. Lett.

Other

T. Suhara, and M. Fujimura, Waveguide Nonlinear-Optic Devices (Berlin: Springer, 2003).

T. Rozzi and M. Mongiardo, Open Electromagnetic Waveguides, (IEE Electromagnetic Waves Series 43, London 1997).
[CrossRef]

D. Marcuse, Theory of Dielectric Optical Waveguides (Academic Press, New York 1974).

A. Taflove, S. C. Hagness, Computational Electrodynamic: the Finite-difference Time-domain Method, (Arthec House Publishers, sec. ed., London 2000), ch. 2,3,4,7.

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (13)

Fig. 1.
Fig. 1.

Multiple gratings design in nonlinear optical symmetrical waveguide.

Fig. 2.
Fig. 2.

Modes in symmetrical planar waveguide.

Fig. 3.
Fig. 3.

The sinusoidally distorted core-substrate interface acts like a phase grating. The incident i, reflected r, and diffracted rays are shown. The sinusoidal approximation is applied to the grating 2 and 3 of Fig. 1.

Fig. 4.
Fig. 4.

Grating 1: normalized nonlinear coupling coefficient versus z.

Fig. 5.
Fig. 5.

Substrate coupling coefficient: (a) coupling coefficient versus Λ2; (b) coupling coefficient versus Λ3. Strong couplings are obtained with Λ2 = 0.95 μm and Λ3 = 1 μm, where the coupling coefficients are regular respect to the z- direction.

Fig. 6.
Fig. 6.

(a) Scattered field generated by the fundamental mode. (b) Scattered field generated by the SH field. The coordinate system is the same of Fig. 1.

Fig. 7.
Fig. 7.

Grating 2: substrate coupling coefficient versus L1.

Fig. 8.
Fig. 8.

Grating 3: substrate coupling coefficient versus L2.

Fig. 9.
Fig. 9.

Analytical and HPF numerical results of the efficiency SH conversion versus z.

Fig. 10.
Fig. 10.

Reflection coefficients Γω,2ω=(βd ω,2ω - βD ω,2ω)/(βd ω,2ω + βD ω,2ω) of the propagating modes at each step discontinuity. Inset: schematic diagram of a step discontinuity. The propagation constants βω,2ω d and βω,2ω D are related to region with core thickness d and D, respectively.

Fig. 11.
Fig. 11.

Coupling coefficients versus h calculated by assuming z=0.

Fig. 12.
Fig. 12.

Structure with grating 1, 2 and 3: analytical, 2D-FDTD and HPF results of the efficiency SH conversion versus z for different h values. The analytical results are averaged values.

Fig. 13.
Fig. 13.

Comparison between HPF and 2D-FDTD normalized spectra of the fundamental and second harmonic modes. Inset: time evolution of the fundamental and second harmonic normalized fields generated inside the waveguide of Fig. 1.

Equations (34)

Equations on this page are rendered with MathJax. Learn more.

H x ( x , z ) = 1 jωμ z E y ( x , z )
H z ( x , z ) = 1 jωμ x E y ( x , z )
E y ( x , z ) = { k a k ψ k ( x ) + 0 b ( k x ) φ ( x ; k x ) d k x } e jβz
E y ω = a ω ψ ω ( x ) e j β ω z
E y 2 ω = a 2 ω ψ 2 ω ( x ) e j β 2 ω z
ψ ω , 2 ω ( x ) = { cos ( σ ω , 2 ω x ) core cos ( σ ω , 2 ω d / 2 ) exp [ Δ ω , 2 ω ( x d / 2 ) ] substrate
a ω , 2 ω = d 1 + 1 Δ ω , 2 ω
Δ ω , 2 ω = ( n 2 2 ( k ω , 2 ω ) 2 ( β ω , 2 ω ) 2 ) 1 / 2
σ ω , 2 ω = ( n 1 2 ( k ω , 2 ω ) 2 ( β ω , 2 ω ) 2 ) 1 / 2
f ( z ) = t sin ( 2 π Λ 2,3 z )
cos ( ϕ ) = β / n i k ω , 2 ω
2 π Λ 2,3 ± β ω , 2 ω < n i k ω , 2 ω
e y ω , 2 ω = A i exp [ j ( σ ω , 2 ω x + β ω , 2 ω z ) ] + A r exp [ j ( σ ω , 2 ω x β ω , 2 ω z ) ] core region
e y ω , 2 ω = A t exp [ j ( Δ ω , 2 ω x + β ω , 2 ω z ) ] substrate region
E y S = t 2 A i σ ω , 2 ω ( Δ ω , 2 ω σ ω , 2 ω ) [ exp [ j ( σ ( ) x β ( ) z ) ] σ ( ) + Δ ( ) exp [ j ( σ ( + ) x β ( + ) z ) ] σ ( + ) + Δ ( + ) ]
core region
E y S = t 2 A i σ ω , 2 ω ( Δ ω , 2 ω σ ω , 2 ω ) [ exp [ j ( Δ ( ) x + β ( ) z ) ] σ ( ) + Δ ( ) exp [ j ( Δ ( + ) x + β ( + ) z ) ] σ ( + ) + Δ ( + ) ]
substrate region
β ( ± ) = β ω , 2 ω ± 2 π Λ 2,3
σ ( ± ) = ( n 1 2 ( k ω , 2 ω ) 2 ( β ( ± ) ) 2 ) 1 / 2
Δ ( ± ) = ( n 2 2 ( k ω , 2 ω ) 2 ( β ( ± ) ) 2 ) 1 / 2
{ d A ω ( z ) dz + j ( 2 k L ω cos Kz ) A ω ( z ) = j [ k NL exp ( j 2 δz ) ] * [ A ω ( z ) ] * A 2 ω ( z ) d A 2 ω ( z ) dz + j ( 2 k L 2 ω cos Kz ) A ω ( z ) = j k NL exp ( j 2 δz ) [ A ω ( z ) ] 2
k L ω = ω ε 0 4 [ E y ω ( x ) ] * Δ ε E y ω ( x ) dx
k L 2 ω = 2 ω ε 0 4 [ E y 2 ω ( x ) ] * Δ ε E y 2 ω ( x ) dx
k NL = 2 ω ε 0 4 [ E y 2 ω ( x ) ] * χ ˜ NL [ E y ω ( x ) ] 2 dx
2 δ = β 2 ω ( 2 β ω + K ) , K = 2 π / Λ 1
Δ ε = ( ε 1 ε 2 ) sin ( a π ) π exp [ jKz ]
χ ˜ NL = ( χ 1 ( 2 ) χ 2 ( 2 ) ) sin ( a π ) π exp [ jKz ]
{ d A ω ( z ) dz = j i a γi ± ( z ) k γi , A ω * exp [ j ( β ± ( β ω 2 π / Λ 2 ) z ] d β ± d a γi ± ( z ) dz = j A ω ( z ) k γi , A ω * exp [ j ( β ± ( β ω 2 π / Λ 2 ) z ]
k γi , A ω = k S ω = E y ω ( x ) Δ ε E y S ( x ) dx
{ d A 2 ω ( z ) dz = j i a γi ± ( z ) k γi , A 2 ω * exp [ j ( β ± ( β ω 2 π / Λ 3 ) z ] d β ± d a γi ± ( z ) dz = j A 2 ω ( z ) k γi , A 2 ω * exp [ j ( β ± ( β ω 2 π / Λ 3 ) z ]
k γi , A 2 ω = k S 2 ω = E y 2 ω ( x ) Δε E y S ( x ) dx
2 β ω + 2 π Λ 1 = β 2 ω or Λ 1 = ( λ FU ) / ( N 2 ω N ω )
E 0 fundamental = sin ( ω fundamental · t · Δ t )

Metrics