Theory of Slow Light Enhanced Four-Wave Mixing in Photonic Crystal Waveguides

doi:10.1364/OE.18.021024

Theory of Slow Light Enhanced Four-Wave Mixing in Photonic Crystal Waveguides

Marco Santagiustina, Carlo G. Someda, Giovanni Vadala', Sylvain Combrie', and Alfredo De Rossi »View Author Affiliations

Optics Express, Vol. 18, Issue 20, pp. 21024-21029 (2010)
http://dx.doi.org/10.1364/OE.18.021024

View Full Text Article

Acrobat PDF (687 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Derivation of the nonlinear propagation equations
3. Numerical results
4. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (25)
Cited By
Figures (2)
Metrics

Abstract

The equations for Four-Wave-Mixing in a photonic crystal waveguide are derived accurately in the hypotesis of negligible nonlinear absorption. The dispersive nature of slow-light enhancement, the impact of Bloch mode reshaping in the nonlinear overlap integrals and the tensor nature of the third order polarization are therefore taken into account. Numerical calculations reveal substantial differences from simpler models, which increase with decreasing group velocity. We predict that the gain for a 1.3 mm long, un-optimized GaInP waveguide will exceed 10 dB if the pump power exceeds 1 W.

1. Introduction

Slow light (SL) can enable interesting applications in photonics [1

Slow Light: Science and Applications , J. B. Khurgin and R. S. Tucker, (CRC Press, Boca Raton, 2009).

] and microwave-photonics [2

M. Santagiustina, “Governing the speed of light: recent advances and future perspectives of slow and fast light in microwave-photonics”, in Proc. 2009 Intern. Top. Meet. on Microwave Photonics , (Valencia, Spain, 2009) Th3.1.

] and it is expected to enhance nonlinear phenomena. In particular, photonic crystal waveguides (PhCWs) present SL propagation [3

T. Baba, “Slow light in photonic crystals”, Nat. Phot. 2, 465–473 (2008). [CrossRef]

] which is predicted to enhance self-phase modulation (SPM) [4

N. A. R. Bhat and J. E. Sipe, “Optical pulse propagation in nonlinear photonic crystals”, Phys. Rev. E 64, 056604 (2001). [CrossRef]

, 5

M Soljacic, M Ibanescu, SG Johnson, Y Fink, and JD Joannopoulos, “Optimal bistable switching in nonlinear photonic crystals”, Phys. Rev. E 66, 055601(R) (2002). [CrossRef]

, 6

T Kamalakis and T Sphicopoulos, “A new formulation of coupled propagation equations in periodic nanophotonic waveguides for the treatment of Kerr-induced nonlinearities”, IEEE J. Quantum Electron. 43, 923–933 (2007). [CrossRef]

]. One intuitive but powerful picture of this effect represents pulses subject to a spatial compression which locally increases the power density [7

T. F. Krauss, “Slow light in photonic crystal waveguides”, J. Phys. D: Appl. Phys. 40, 2666–2670 (2007). [CrossRef]

], pretty much as cars in a highway get closer to each other as their speed is decreased. The dependence of nonlinearity on group velocity has been recently observed in PhCW for third harmonic generation [8

B. Corcoran, C. Monat, C. Grillet, D. Moss, B. J. Eggleton, T. White, L. O'Faolain, and T. Krauss, “Green light emission in silicon through slow-light enhanced third-harmonic generation in photonic-crystal waveguides”, Nat. Phot. 3, 206–210 (2009). [CrossRef]

, 9

B. Corcoran, C. Monat, M. Pelusi, C. Grillet, T. P. White, L. O'Faolain, T. F. Krauss, B. J. Eggleton, and D. J. Moss, “Optical signal processing on a silicon chip at 640Gb/s using slow-light”, Opt. Express 18, 7770–7781 (2010). [CrossRef] [PubMed]

], for SPM and for three photon absorption [10

S. Combrié, Q. Vy Tran, C. Husko, P. Colman, and A. De Rossi, “High quality GaInP nonlinear photonic crystals with minimized nonlinear absorption”, Appl. Phys. Lett. 95, 221108 (2009). [CrossRef]

, 11

C. Husko, S. Combrié, Q. Tran, F. Raineri, C. Wong, and A. De Rossi, “Non-trivial scaling of self-phase modulation and three-photon absorptionin III-V photonic crystal waveguides,” Opt. Express 17, 22442–22451 (2009). [CrossRef]

Efficient four-wave mixing (FWM) was reported recently in 1.3 mm long, GaInP PhCW [12

V. Eckhouse, I. Cestier, G. Eisenstein, S. Combrié, P. Colman, A. De Rossi, M. Santagiustina, C. G. Someda, and G. Vadalà, “Highly efficient four wave mixing in GaInP photonic crystal waveguides”, Opt. Lett. 35, 1440–1142 (2010). [CrossRef] [PubMed]

], with a conversion efficiency comparable to that of about 1 m-long highly nonlinear fiber [13

T. Hasegawa, T. Nagashima, and N. Sugimoto, “Determination of nonlinear coefficient and group-velocity dispersion of bismuth-based high nonlinear optical fiber by four-wave mixing”, Opt. Commun. 281, 782–787 (2008). [CrossRef]

] and about 2 cm-long chalcogenide fiber photonic wire [14

M. D. Pelusi, F. Luan, E. Magi, M. R. E. Lamont, D. J. Moss, B. J. Eggleton, J. S. Sanghera, L. B. Shaw, and I. D. Aggarwal, “High bit rate all-optical signal processing in a fiber photonic wire”, Opt. Express 16, 11506–11512 (2008). [CrossRef] [PubMed]

]. That result confirms the theoretical prediction [15

M. Ebnali-Heidari, C. Monat, C. Grillet, and M. K. Moravvej-Farshi, “A proposal for enhancing four-wave mixing in slow light engineered photonic crystal waveguides and its application to optical regeneration”, Opt. Express 17, 18340–18353 (2009). [CrossRef] [PubMed]

] of enhanced (FWM) on the basis of the square group index scaling factor. However, it must be pointed out that a simple square group index scaling, to model the FWM enhancement, does not actually take into account several fundamental features of the phenomenon.

The first feature is that the group index in the waveguide is a function of frequency and, therefore, it is not the same for the various waves involved in the FWM. Observe in fact in Fig. 1a the numerically calculated [16

http://ab-initio.mit.edu/photons/

] group index (red dots) of a GaInP membrane PhCW. The calculation is carried out with pumps placed at N = 6 different wavelengths approaching the band edge: λ _1M = [1570+10 (M − 1)] nm (M = 1, …, 6). The second pump wavelength is λ _2M = λ _1M + 2 nm; the signal and idler frequency are calculated according to ω _3M = ω _1M − Δω_M, ω _4M = ω _2M + Δω_M , where Δω_M = ω _1M − ω _2M, thus satisfying the FWM frequency condition ω _1M + ω _2M = ω _3M + ω _4M. The waveguide parameters are: a = 480 nm (crystal period), d = 0.38a (hole diameter) and h = 170 nm (PhCW slab height).

Fig. 1 Comparison with model in ref. [15

]. (a) Square of the group index

n_{g}^{2}

and SL enhancement factor S, as a function of the wavelength. In the insets the intensity of electric field |e _i |² of the Bloch mode within an elementary cell is shown at three different wavelengths. (b) Wavelength dependence of the SL scaling factor and of the effective FWM nonlinear coefficient, both normalized to their respective value at λ ₁₁ = 1570 nm: γ_norm = |γ_F ₄(λ _1M)/γ _F4(λ ₁₁)|,S_norm =S(λ _1M)/S(λ ₁₁).

Moreover, the modal superposition of the interacting fields must be considered. In PhCWs the overlap integrals are functions of frequency; in particular, as the frequency approaches the bandgap, the mode spreads into the hole region as shown in Fig. 1a, where the intensity distribution of the electric field within a cell of the PhCW is shown at three different wavelengths. Finally, the PhCWs mode is not constant but periodic along the propagation direction and so, differently from other photonic waveguides (slabs, fibers etc.), the careful determination of the nonlinear effective coefficients is more complicated [4

N. A. R. Bhat and J. E. Sipe, “Optical pulse propagation in nonlinear photonic crystals”, Phys. Rev. E 64, 056604 (2001). [CrossRef]

, 6

, 17

N. C. Panoiu, J. F. McMillan, and C. W. Wong, “Theoretical analysis of pulse dynamics in silicon photonic crystal wire waveguides”, IEEE J. Sel. T. Quantum Electron. 16, 257–266 (2010). [CrossRef]

In this paper, an accurate calculation of all nonlinear effective coefficients necessary to evaluate the FWM interaction in the SL regime of a PhCW is carried out. Similarly to [18

D. Michaelis, U. Peschel, C. Wächter, and A. Braüer, “Reciprocity theorem and perturbation theory for photonic crystal waveguides”, Phys. Rev. E 68, 065601(R) (2003). [CrossRef]

], the derivation is performed through a perturbation approach directly from Maxwell's equations. The resulting SPM coefficient corresponds to that obtained by previous derivations [4

N. A. R. Bhat and J. E. Sipe, “Optical pulse propagation in nonlinear photonic crystals”, Phys. Rev. E 64, 056604 (2001). [CrossRef]

, 6

]; the XPM coefficient is consistent with the one obtained for multimode propagation [17

]. The FWM effective coefficients are determined for the first time to the best of our knowledge. Finally we consider the specific case of GaInP around 1.5 μm, where two and three photon absorption can be neglected [10

S. Combrié, Q. Vy Tran, C. Husko, P. Colman, and A. De Rossi, “High quality GaInP nonlinear photonic crystals with minimized nonlinear absorption”, Appl. Phys. Lett. 95, 221108 (2009). [CrossRef]

], and we evaluate numerically the coefficients.

2. Derivation of the nonlinear propagation equations

The starting point of our analysis are Maxwell's equations in the frequency domain where the linear permittivity ε(r) is a spatial function describing the PhCWs structure and P _NL (r,t) accounts for the nonlinear response (ι is the imaginary unit):

\nabla \times E (r, ω) = ι ω μ H (r, ω), \nabla \times H (r, ω) = - ι ω ε (r) E (r, ω) - ι ω P_{NL} (r, ω) .

(1)

It is assumed that four signals are propagating in the fundamental TE mode of the PhCW, at frequencies that satisfy the FWM condition ω ₁ + ω ₂ = ω ₃ + ω ₄. The electric and magnetic fields are then expanded as:

E (r, ω) = \frac{1}{2} \sum_{i = - 4, \neq 0}^{4} A_{i} e (r, ω_{i}) exp [ι (β_{i} z)] δ_{i}, H (r, ω) = \frac{1}{2} \sum_{i = - 4, \neq 0}^{4} A_{i} h (r, ω_{i}) exp [ι (β_{i} z)] δ_{i},

(2)

where A_i are the complex amplitudes, the pairs: e(r, ω_i ) exp(ιβ_iz), h(r, ω_i ) exp(ιβ_iz) are the Bloch modes at frequencies ω_i (ω_−i = −ω_i ) with β_i = β(ω_i ) the propagation constant and δ_i = δ(ω − ω_i ). Bloch modes satisfy linear Maxwell's equations individually:

\nabla \times [e_{i} exp (ι β_{i} z)] = ι μ ω_{i} h_{i} exp (ι β_{i} z), \nabla \times [h_{i} exp (ι β_{i} z)] = - ι ε ω_{i} e_{i} exp (ι β_{i} z) .

(3)

and also obey the relations:

e_{- i} = e_{i}^{*}, h_{- i} = - h_{i}^{*}

. For the sake of simplicity, we have omitted the explicit dependence on space and replaced the frequency ω_i with the subscript i.

We assume that the nonlinearity is small enough so that the complex amplitude, A_i = A_i (z), in the direction ẑ is slowly varying in comparison to exp(ιβ_iz) and to the Bloch mode within the cell {e _i , h _i }. Furthermore, the Bloch modes are normalized:

\int_{V} (e_{i} \times h_{i}^{*} + e_{i}^{*} \times h_{i}) \cdot \hat{z} dV = 4 a

. Let us stress that, with this choice, |A_i |² = P_i is the active power propagating in the z direction at frequency ω_i [6

]. In the following, we will show that this is the natural choice for normalizing the Bloch modes when group velocity is substantially different from the phase velocity. Similarly to ref. [18

D. Michaelis, U. Peschel, C. Wächter, and A. Braüer, “Reciprocity theorem and perturbation theory for photonic crystal waveguides”, Phys. Rev. E 68, 065601(R) (2003). [CrossRef]

] we consider Eqs. (1) calculated at frequency ω_i ; the second equation is scalarly multiplied by

e_{i}^{*} exp (- ι β_{i} z)

and then subtracted from the first, multiplied by

h_{i}^{*} exp (- ι β_{i} z)

. The result is integrated over the volume of the PhCW unit cell to obtain:

\frac{\partial A_{i}}{\partial z} \int_{V} \hat{z} \cdot [e_{i} \times h_{i}^{*} + e_{i}^{*} \times h_{i}] dV = j ω_{i} \int_{V} e_{i}^{*} exp (- ι β_{i} z) \cdot P_{NL} (r, ω_{i}) dV

(4)

Here, we also used the hypothesis that A_i , ∂A_i/dz are slowly varying functions of z and therefore can be taken constant over one unit cell.

We now introduce the explicit form of the third-order nonlinear polarization P _NL (r, ω_i ). Using the notation of [19

R. Boyd, Nonlinear Optics , Chapt. 4 (Academic Press, San Diego, 2003).

] it reads:

\begin{array}{r} P_{NL} (r, ω_{i}) = ε_{0} χ^{(3)} (r; ω_{i}; ω_{i}, - ω_{i} ω_{i}) ⋮ E_{i} E_{i}^{*} E_{i} + \\ + ε_{0} \sum_{j = 1, \neq i}^{4} [χ^{(3)} (r; ω_{i}; ω_{j}, - ω_{j}, ω_{i}) ⋮ E_{j} E_{j}^{*} E_{i}] + ε_{0} χ^{(3)} (r; ω_{i}; ω_{j}, - ω_{l}, ω_{k}) ⋮ E_{j} E_{l}^{*} E_{k}, \end{array}

(5)

where i,j,k,l = {1, 2, 3, 4}, with the constraint that all indexes are different in the last term. The susceptibility tensor is real because multi-photon absorption can be neglected [10

S. Combrié, Q. Vy Tran, C. Husko, P. Colman, and A. De Rossi, “High quality GaInP nonlinear photonic crystals with minimized nonlinear absorption”, Appl. Phys. Lett. 95, 221108 (2009). [CrossRef]

]. On the right hand side (RHS) of Eq. (5), the first term represents SPM, the summation term XPM, and the last term the non-degenerate FWM. Here, third and other harmonic generations are neglected by assuming that they will not be phase matched. For the sake of brevity the tensor explicit dependence on position and frequencies will also be omitted till the end of the derivation. Inserting Eq. (5) in Eq.(4) we obtain:

\begin{array}{r} 4 a \frac{\partial A_{i}}{\partial z} = \frac{ι ω_{i} ε_{0}}{4} [{| A_{i} |}^{2} A_{i} \int_{V} e_{i}^{*} \cdot χ^{(3)} ⋮ e_{i} e_{i}^{*} e_{i} dV + \\ + \sum_{j = 1, \neq i}^{4} {| A_{j} |}^{2} A_{i} \int_{V} e_{i}^{*} \cdot χ^{(3)} ⋮ e_{j} e_{j}^{*} e_{i} d V + A_{j} A_{l}^{*} A_{k} exp (- ι σ_{i} Δ β z) \int_{V} e_{i}^{*} \cdot χ^{(3)} ⋮ e_{j} e_{l}^{*} e_{k} dV] . \end{array}

(6)

Here Δβ = β ₃ + β ₄ − β ₁ − β ₂ is the linear phase mismatch and σ_i = ±1, where the plus (minus) sign applies for i = 3, 4 (i = 1, 2). In Eq. (6), the SL enhancement of the nonlinear response is hidden in the integrals. In order to make this dependence explicit, we use the identity between the electromagnetic energy velocity v _e and the group velocity that holds in lossless homogeneous media, in periodic ones [20

P. Yeh, “Electromagnetic propagation in birefringent layered media” J. Opt. Soc. Am. 69, 742 (1979). [CrossRef]

] and in PhCWs [21

K. Sakoda, Optical Properties of Photonic Crystals , Chapt. 2 (Springer, Berlin, 2005).

]. By projecting the energy velocity along the axis unit vector ẑ and by using the property that the space-time average magnetic and electric energies are equal for Bloch modes [22

B. Lombardet, L. A. Dunbar, R. Ferrini, and R. Houdre, “Bloch wave propagation in two-dimensional photonic crystals: Influence of the polarization”, Opt. Q. Electr. 37, 293–307 (2005). [CrossRef]

μ_{0} / 4 \int_{V} h_{i} \cdot h_{i}^{*} = 1 / 4 \int_{V} e_{i} \cdot d_{i}^{*}

, the following is obtained:

v_{ei} \cdot \hat{z} = \frac{1 / 4 \int_{V} (e_{i} \times h_{i}^{*} + e_{i}^{*} \times h_{i}) \cdot \hat{z} dV}{1 / 4 \int_{V} (ε_{0} ε_{r} (r) | e_{i} |^{2} + μ_{0} | h_{i} |^{2}) dV} = \frac{4 a}{2 \int_{V} (ε_{0} ε_{r} (r) | e_{i} |^{2} dV} = \frac{2 a}{ε_{0} W_{i}} = v_{gi} .

(7)

Note that energies appearing above are normalized consistently with the choice |A_i |² = P_i . We can now normalize the terms on the RHS of Eq. (6) multiplying them by the factors

η_{i}^{4}

(SPM),

η_{i}^{2} η_{j}^{2}

(XPM) and η_iη_jη_kη_l (FWM), with

η_{i} = \sqrt{2 a / (ε_{0} W_{i} v_{g i})} = 1, \forall_{i}

. The nonlinear coefficients of Eq. (6) are now cast in their canonical form so the equations governing the FWM in PhCW are obtained in a form similar to nonlinear fiber optics [23

G. P. Agrawal, Nonlinear fiber optics , Chapt. 10 (Academic Press, San Diego, 2001).

\frac{d A_{i}}{d z} = ι γ_{i} {| A_{i} |}^{2} A_{i} + 2 ι \sum_{j = 1, \neq i}^{4} γ_{i j} {| A_{j} |}^{2} A_{i} + 2 ι γ_{Fi} A_{l}^{*} A_{j} A_{k} e^{- ι σ_{i} Δ β z}, i = 1, 2, 3, 4,

(8)

with the effective nonlinear coefficients and the relative effective volumes taking the form:

γ_{i} = \frac{n_{2} ω_{i} a}{c V_{i}}; \frac{1}{V_{i}} = \frac{n_{g i}^{2}}{W_{i}^{2}} \int_{V} \frac{ε_{r}}{3 χ_{xxxx}^{(3)}} e_{i}^{*} \cdot χ^{(3)} (r; ω_{i}; ω_{i}, - ω_{i}, ω_{i}) ⋮ e_{i} e_{i}^{*} e_{i} dV;

(9)

γ_{ij} = \frac{n_{2} ω_{i} a}{c V_{ij}}; \frac{1}{V_{ij}} = \frac{n_{gi} n_{gj}}{W_{i} W_{j}} \int_{V} \frac{ε_{r}}{6 χ_{xxxx}^{(3)}} e_{i}^{*} \cdot χ^{(3)} (r; ω_{i}; ω_{j}, - ω_{j}, ω_{i}) ⋮ e_{j} e_{j}^{*} e_{i} dV;

(10)

γ_{Fi} = \frac{n_{2} ω_{i} a}{c V_{Fi}}; \frac{1}{V_{F i}} = {\prod_{n = 1}^{4} (\frac{n_{gn}}{W_{n}})}^{1 / 2} \int_{V} \frac{ε_{r}}{6 χ_{xxxx}^{(3)}} e_{i}^{*} \cdot χ^{(3)} (r; ω_{i}; ω_{j}, - ω_{l}, ω_{k}) ⋮ e_{j} e_{l}^{*} e_{k} dV;

(11)

and where

n_{2} = 3 χ_{xxxx}^{(3)} / (4 ε_{r} ε_{0} c)

is the bulk, nonlinear refractive index coefficient for a linear state of polarization [19

R. Boyd, Nonlinear Optics , Chapt. 4 (Academic Press, San Diego, 2003).

So, for each nonlinear effect (SPM, XPM and FWM) we determined: 1) the correct enhancement factor due to SL; 2) the correct overlap integral. The obtained SPM and XPM coefficients are consistent to those previously found [4

N. A. R. Bhat and J. E. Sipe, “Optical pulse propagation in nonlinear photonic crystals”, Phys. Rev. E 64, 056604 (2001). [CrossRef]

, 6

, 17

]; the FWM coefficient is derived for the first time to the best of our knowledge. By observing Eqs. (9,10,11) a general rule can be remarked: the enhancement factor due to the SL is always given by the geometric mean of the group indexes of the waves interacting through the tensor χ⁽³⁾.

3. Numerical results

The theoretical findings of the previous section are applied to the PhCW of Fig. 1a. To the aim of explicitly calculating the tensor products in Eqs. (9,10,11) the theory of ref. [24

D. C. Hutchings and B. S. Wherrett, “Polarisation dichroism of nonlinear refraction in zinc-blende semiconductors”, Opt. Commun. 111, 507–512 (1994). [CrossRef]

], that can be generally applied to zincblend semiconductors (group symmetry 43m), is exploited and then

χ_{xxxx}^{(3)} = 2 χ_{xyxy}^{(3)} = 2 χ_{xyyx}^{(3)} = χ_{xxyy}^{(3)}

. For signal wavelengths 1.52 μm < λ < 1.62 μm and given that E_g ≃ 1.9 eV for GaInP, the frequencies ω_i,j,k,l at which the tensor elements are to be calculated satisfy the condition 0.4 < h̄ω_i,j,k,l/E_g < 0.43. Then, although Kleinmann symmetry [19

R. Boyd, Nonlinear Optics , Chapt. 4 (Academic Press, San Diego, 2003).

] is not satisfied, it is found that the dichroism parameter, defined by

χ_{xxyy}^{(3)} / χ_{xxxx}^{(3)}

, can be approximated to 0.28 [24

D. C. Hutchings and B. S. Wherrett, “Polarisation dichroism of nonlinear refraction in zinc-blende semiconductors”, Opt. Commun. 111, 507–512 (1994). [CrossRef]

]. Thus all tensors can be determined from the above relations, from the knowledge of the nonlinear refractive index in GaInP, n ₂ = 10⁻¹⁷ m²/W and through [19

R. Boyd, Nonlinear Optics , Chapt. 4 (Academic Press, San Diego, 2003).

e_{i}^{*} \cdot χ^{(3)} (r) ⋮ e_{j} e_{l}^{*} e_{k} = \sum_{m} [e_{im}^{*} D \sum_{nop} χ_{mnop}^{(3)} (r; ω_{i}; ω_{j}, - ω_{l}, ω_{k}) e_{jn} e_{lo}^{*} e_{kp}]

(12)

where the summations over the indexes m,n,o,p are made on all possible values of the coordinate axes {x,y,z} and D is the frequency degeneracy factor [19

R. Boyd, Nonlinear Optics , Chapt. 4 (Academic Press, San Diego, 2003).

] which represents the number of distinct permutations of the three frequencies {ω_j, − ω_l, ω_k } (D = 3 for SPM, D = 6 for XPM and FWM).

We numerically determined, for all previously defined pump and signal wavelengths the Bloch mode electrical field and the dispersion relation [16

http://ab-initio.mit.edu/photons/

], then calculating the FWM SL scaling factor

S = \prod_{k = 1}^{4} n_{gk}^{1 / 2}

and the effective nonlinear coefficient γ _F4 according to Eq. (11). In Fig. 1a, S can be compared to the square of the group index at the mean frequency, which is the approximation used in refs. [12

, 15

]; a slight discrepancy appears at the band edge (where the SL effect is strong). Note that no particular dispersion engineering of the PhCW has been realized on purpose; Δβ could still be reduced by design [25

J. Li, T. P. White, L. O'Faolain, A. Gomez-Iglesias, and T. F. Krauss, “Systematic design of flat band slow light in photonic crystal waveguides”, Opt. Express 16, 6227–6232 (2008). [CrossRef] [PubMed]

], to increase FWM efficiency. To evaluate the effective nonlinear enhancement of PhCW, in Fig. 1b the wavelength dependence of S is compared to that of γ, revealing that the enhancement of FWM coefficient does not follow the pure SL scaling S the large difference deriving from the decrease in the modal overlap integrals.

By following [23

G. P. Agrawal, Nonlinear fiber optics , Chapt. 10 (Academic Press, San Diego, 2001).

] we can finally determine the nonlinear phase matching and gain coefficient for optical parametric amplification of the wave at frequency ω ₃:

Δ κ = Δ β - \sum_{i, j = 1}^{2} [γ_{i} + 2 | i - j | γ_{ji} - 2 γ_{3 i} - 2 γ_{4 i}] P_{i} = Δ β + γ_{pm} P_{0}, g = {[γ_{gm}^{2} P_{0}^{2} - \frac{Δ κ^{2}}{4}]}^{1 / 2},

(13)

where the last terms of each of Eqs. (13) are obtained for P ₁ = P ₂ = P ₀/2 and

γ_{gm}^{2} = γ_{F 3} γ_{F 4}^{*}

. The effective coefficients γ_pm and γ_gm describe the strength of the nonlinearity contribution, respectively, to the phase matching and to the maximum gain. It is remarkable that, differently from fiber optics [23

G. P. Agrawal, Nonlinear fiber optics , Chapt. 10 (Academic Press, San Diego, 2001).

], γ_pm ≠ γ_gm . Figure (2a) compares the maximum achievable gain coefficient γ_gmP ₀ to the actual one g, which is limited by the phase mismatch. As the pump power increases from 0.9 W to 1.4 W the phase mismatch is almost completely canceled by the nonlinear phase terms and the maximum gain is approached. This fact leads to a dramatic increase of the FWM gain

G_{3} = P_{3} (L) / P_{3} (0) = 1 + γ_{gm}^{2} P_{0}^{2} / g^{2} {sinh}^{2} (gL)

which is shown in fig. (2b) for L = 1.3 mm.

Fig. 2 (a) Comparison of the maximum achievable parametric gain coefficient, γ_gmP ₀ (solid curves), and the actually achieved coefficient, g (dashed curves), for two different powers (P ₀ = 0.9 W circles, P ₀ = 1.4 W triangles) as a function of SL scaling factor S. (b) Parametric gain as a function of λ₁ for a waveguide L = 1.3 mm long, for three different pump powers.

4. Conclusions

We have derived the nonlinear equations which describe four-wave mixing in photonic crystal waveguides directly from Maxwell's equations. These equations are exact in the limit in which we can neglect the changes that the nonlinearity induces in the Bloch modes describing the field in the photonic crystal waveguide, a situation by far verified in practice. We demonstrated rigorously the explicit dependence of the nonlinear enhancement on the group index and that in four-wave mixing (where the fields involved have different group indexes) the dispersive nature is rigorously accounted for by the geometric mean of the group indexes of the modes involved. Moreover, we demonstrated a substantial correction arising from Bloch mode reshaping in the nonlinear field overlap. Finally we accounted for the tensor nature of the nonlinear polarization. As an example, we calculated the gain for a 1.3 mm long GaInP waveguide operated at moderately small group velocity (ν_g > c/20). Particularly, when n_g ≃ 12 and the coupled pump power > 1 W, the expected gain exceeds 10 dB even if the waveguide is un-optimized.

Acknowledgments

The research leading to these results has received funding from EC FET-Open programme (FP7/2007-2011) under grant agreement n. 219299 GOSPEL.

References and links

1.	Slow Light: Science and Applications , J. B. Khurgin and R. S. Tucker, (CRC Press, Boca Raton, 2009).
2.	M. Santagiustina, “Governing the speed of light: recent advances and future perspectives of slow and fast light in microwave-photonics”, in Proc. 2009 Intern. Top. Meet. on Microwave Photonics , (Valencia, Spain, 2009) Th3.1.
3.	T. Baba, “Slow light in photonic crystals”, Nat. Phot. 2, 465–473 (2008). [CrossRef]
4.	N. A. R. Bhat and J. E. Sipe, “Optical pulse propagation in nonlinear photonic crystals”, Phys. Rev. E 64, 056604 (2001). [CrossRef]
5.	M Soljacic, M Ibanescu, SG Johnson, Y Fink, and JD Joannopoulos, “Optimal bistable switching in nonlinear photonic crystals”, Phys. Rev. E 66, 055601(R) (2002). [CrossRef]
6.	T Kamalakis and T Sphicopoulos, “A new formulation of coupled propagation equations in periodic nanophotonic waveguides for the treatment of Kerr-induced nonlinearities”, IEEE J. Quantum Electron. 43, 923–933 (2007). [CrossRef]
7.	T. F. Krauss, “Slow light in photonic crystal waveguides”, J. Phys. D: Appl. Phys. 40, 2666–2670 (2007). [CrossRef]
8.	B. Corcoran, C. Monat, C. Grillet, D. Moss, B. J. Eggleton, T. White, L. O'Faolain, and T. Krauss, “Green light emission in silicon through slow-light enhanced third-harmonic generation in photonic-crystal waveguides”, Nat. Phot. 3, 206–210 (2009). [CrossRef]
9.	B. Corcoran, C. Monat, M. Pelusi, C. Grillet, T. P. White, L. O'Faolain, T. F. Krauss, B. J. Eggleton, and D. J. Moss, “Optical signal processing on a silicon chip at 640Gb/s using slow-light”, Opt. Express 18, 7770–7781 (2010). [CrossRef] [PubMed]
10.	S. Combrié, Q. Vy Tran, C. Husko, P. Colman, and A. De Rossi, “High quality GaInP nonlinear photonic crystals with minimized nonlinear absorption”, Appl. Phys. Lett. 95, 221108 (2009). [CrossRef]
11.	C. Husko, S. Combrié, Q. Tran, F. Raineri, C. Wong, and A. De Rossi, “Non-trivial scaling of self-phase modulation and three-photon absorptionin III-V photonic crystal waveguides,” Opt. Express 17, 22442–22451 (2009). [CrossRef]
12.	V. Eckhouse, I. Cestier, G. Eisenstein, S. Combrié, P. Colman, A. De Rossi, M. Santagiustina, C. G. Someda, and G. Vadalà, “Highly efficient four wave mixing in GaInP photonic crystal waveguides”, Opt. Lett. 35, 1440–1142 (2010). [CrossRef] [PubMed]
13.	T. Hasegawa, T. Nagashima, and N. Sugimoto, “Determination of nonlinear coefficient and group-velocity dispersion of bismuth-based high nonlinear optical fiber by four-wave mixing”, Opt. Commun. 281, 782–787 (2008). [CrossRef]
14.	M. D. Pelusi, F. Luan, E. Magi, M. R. E. Lamont, D. J. Moss, B. J. Eggleton, J. S. Sanghera, L. B. Shaw, and I. D. Aggarwal, “High bit rate all-optical signal processing in a fiber photonic wire”, Opt. Express 16, 11506–11512 (2008). [CrossRef] [PubMed]
15.	M. Ebnali-Heidari, C. Monat, C. Grillet, and M. K. Moravvej-Farshi, “A proposal for enhancing four-wave mixing in slow light engineered photonic crystal waveguides and its application to optical regeneration”, Opt. Express 17, 18340–18353 (2009). [CrossRef] [PubMed]
16.	http://ab-initio.mit.edu/photons/
17.	N. C. Panoiu, J. F. McMillan, and C. W. Wong, “Theoretical analysis of pulse dynamics in silicon photonic crystal wire waveguides”, IEEE J. Sel. T. Quantum Electron. 16, 257–266 (2010). [CrossRef]
18.	D. Michaelis, U. Peschel, C. Wächter, and A. Braüer, “Reciprocity theorem and perturbation theory for photonic crystal waveguides”, Phys. Rev. E 68, 065601(R) (2003). [CrossRef]
19.	R. Boyd, Nonlinear Optics , Chapt. 4 (Academic Press, San Diego, 2003).
20.	P. Yeh, “Electromagnetic propagation in birefringent layered media” J. Opt. Soc. Am. 69, 742 (1979). [CrossRef]
21.	K. Sakoda, Optical Properties of Photonic Crystals , Chapt. 2 (Springer, Berlin, 2005).
22.	B. Lombardet, L. A. Dunbar, R. Ferrini, and R. Houdre, “Bloch wave propagation in two-dimensional photonic crystals: Influence of the polarization”, Opt. Q. Electr. 37, 293–307 (2005). [CrossRef]
23.	G. P. Agrawal, Nonlinear fiber optics , Chapt. 10 (Academic Press, San Diego, 2001).
24.	D. C. Hutchings and B. S. Wherrett, “Polarisation dichroism of nonlinear refraction in zinc-blende semiconductors”, Opt. Commun. 111, 507–512 (1994). [CrossRef]
25.	J. Li, T. P. White, L. O'Faolain, A. Gomez-Iglesias, and T. F. Krauss, “Systematic design of flat band slow light in photonic crystal waveguides”, Opt. Express 16, 6227–6232 (2008). [CrossRef] [PubMed]

OCIS Codes
(130.5990) Integrated optics : Semiconductors
(190.4380) Nonlinear optics : Nonlinear optics, four-wave mixing
(130.5296) Integrated optics : Photonic crystal waveguides

ToC Category:
Nonlinear Optics

History
Original Manuscript: June 9, 2010
Revised Manuscript: July 22, 2010
Manuscript Accepted: July 27, 2010
Published: September 20, 2010

Citation
Marco Santagiustina, Carlo G. Someda, Giovanni Vadala', Sylvain Combrie', and Alfredo De Rossi, "Theory of Slow Light Enhanced Four-Wave Mixing in Photonic Crystal Waveguides," Opt. Express 18, 21024-21029 (2010)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-18-20-21024

Sort: Author | Year | Journal | Reset

References

. Slow Light: Science and Applications, J. B. Khurgin and R. S. Tucker, Eds., (CRC Press, Boca Raton, 2009).
. M. Santagiustina, “Governing the speed of light: recent advances and future perspectives of slow and fast light in microwave-photonics”, in Proc. 2009 Intern. Top. Meet. on Microwave Photonics, (Valencia, Spain, 2009) Th3.1.
. T. Baba, “Slow light in photonic crystals.” Nat. Phot. 2, 465–473 (2008). [CrossRef]
. N. A. R. Bhat, and J. E. Sipe, “Optical pulse propagation in nonlinear photonic crystals”. Phys. Rev. E 64, 056604 (2001). [CrossRef]
. M. Soljacic, M. Ibanescu, S. G. Johnson, Y. Fink, and J. D. Joannopoulos, “Optimal bistable switching in nonlinear photonic crystals,” Phys. Rev. E 66, 055601(R) (2002). [CrossRef]
. T. Kamalakis, and T. Sphicopoulos, “A new formulation of coupled propagation equations in periodic nanophotonic waveguides for the treatment of Kerr-induced nonlinearities,” IEEE J. Quantum Electron. 43, 923-933 (2007). [CrossRef]
. T. F. Krauss, “Slow light in photonic crystal waveguides,” J. Phys. D: Appl. Phys. 40, 2666-2670 (2007). [CrossRef]
. B. Corcoran, C. Monat, C. Grillet, D. Moss, B. J. Eggleton, T. White, L. O’Faolain, and T. Krauss, “Green light emission in silicon through slow-light enhanced third-harmonic generation in photonic-crystal waveguides,” Nat. Phot. 3, 206-210 (2009). [CrossRef]
. B. Corcoran, C. Monat, M. Pelusi, C. Grillet, T. P. White, L. O’Faolain, T. F. Krauss, B. J. Eggleton, and D. J. Moss, “Optical signal processing on a silicon chip at 640Gb/s using slow-light,” Opt. Express 18, 7770-7781 (2010). [CrossRef] [PubMed]
. S. Combrié, Q. Vy Tran, C. Husko, P. Colman, and A. De Rossi, “High quality GaInP nonlinear photonic crystals with minimized nonlinear absorption,” Appl. Phys. Lett. 95, 221108 (2009). [CrossRef]
. C. Husko, S. Combrié, Q. Tran, F. Raineri, C. Wong, and A. De Rossi, “Non-trivial scaling of self-phase modulation and three-photon absorptionin III-V photonic crystal waveguides,” Opt. Express 17, 22442-22451 (2009). [CrossRef]
. V. Eckhouse, I. Cestier, G. Eisenstein, S. Combrié, P. Colman, A. De Rossi, M. Santagiustina, C. G. Someda, and G. Vadalà, “Highly efficient four wave mixing in GaInP photonic crystal waveguides,” Opt. Lett. 35, 1440-1442 (2010). [CrossRef] [PubMed]
. T. Hasegawa, T. Nagashima, and N. Sugimoto, “Determination of nonlinear coefficient and group-velocity dispersion of bismuth-based high nonlinear optical fiber by four-wave mixing,” Opt. Commun. 281, 782-787 (2008). [CrossRef]
. M. D. Pelusi, F. Luan, E. Magi, M. R. E. Lamont, D. J. Moss, B. J. Eggleton, J. S. Sanghera, L. B. Shaw, and I. D. Aggarwal, “High bit rate all-optical signal processing in a fiber photonic wire,” Opt. Express 16, 11506-11512 (2008). [CrossRef] [PubMed]
. M. Ebnali-Heidari, C. Monat, C. Grillet, and M. K. Moravvej-Farshi, “A proposal for enhancing four-wave mixing in slow light engineered photonic crystal waveguides and its application to optical regeneration,” Opt. Express 17, 18340-18353 (2009). [CrossRef] [PubMed]
. http://ab-initio.mit.edu/photons/
. N. C. Panoiu, J. F. McMillan, and C. W. Wong, “Theoretical analysis of pulse dynamics in silicon photonic crystal wire waveguides,” IEEE J. Sel. T. Quantum Electron. 16, 257-266 (2010). [CrossRef]
. D. Michaelis, U. Peschel, C. Wächter, and A. Braüer, “Reciprocity theorem and perturbation theory for photonic crystal waveguides,” Phys. Rev. E 68, 065601(R) (2003). [CrossRef]
. R. Boyd, Nonlinear Optics Chapt. 4 (Academic Press, San Diego, 2003).
. P. Yeh, “Electromagnetic propagation in birefringent layered media,” J. Opt. Soc. Am. 69, 742 (1979). [CrossRef]
. K. Sakoda, Optical Properties of Photonic Crystals, Chapt. 2 (Springer, Berlin, 2005).
. B. Lombardet, L. A. Dunbar, R. Ferrini, and R. Houdre, “Bloch wave propagation in two-dimensional photonic crystals: Influence of the polarization,” Opt. Q. Electr. 37, 293-307 (2005). [CrossRef]
. G. P. Agrawal, Nonlinear fiber optics, Chapt. 10 (Academic Press, San Diego, 2001).
. D. C. Hutchings, and B. S. Wherrett, “Polarisation dichroism of nonlinear refraction in zinc-blende semiconductors,” Opt. Commun. 111, 507–512 (1994). [CrossRef]
. J. Li, T. P. White, L. O’Faolain, A. Gomez-Iglesias, and T. F. Krauss, “Systematic design of flat band slow light in photonic crystal waveguides,” Opt. Express 16, 6227-6232 (2008). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1	Fig. 2

« Previous Article | Next Article »

OSA is a member of CrossRef.

OpticsInfoBase

Optics Express

Theory of Slow Light Enhanced Four-Wave Mixing in Photonic Crystal Waveguides

Article Tools

Article Navigation

Abstract

1. Introduction

2. Derivation of the nonlinear propagation equations

3. Numerical results

4. Conclusions

Acknowledgments

References and links

References

Appl. Phys. Lett.

IEEE J. Quantum Electron.

IEEE J. Sel. T. Quantum Electron.

J. Opt. Soc. Am.

J. Phys. D: Appl. Phys.

Nat. Phot.

Opt. Commun.

Opt. Express

Opt. Lett.

Opt. Q. Electr.

Phys. Rev. E

Other

Cited By

Figures

Select a Conference
Session or Paper #
Year