Theory of Photonic Crystal Slabs by the Guided-Mode Expansion Method Dario Gerace* Lucio Claudio Andreani Dipartimento di Fisica Alessandro Volta, Università degli Studi di Pavia, Pavia (Italy) *Quantum Photonics Group, Institute of Quantum Electronics, ETHZ, gerace@phys.ethz.ch Workshop on Numerical Methods for Optical Nanostructures, ETHZ, July 2007
Where is Pavia??? Europe Italy Pavia is 25 Km far from Milano A nice small town with an ancient (~XIII century) university tradition
Collaboration: Acknowledgements Mario Agio (now at Physical Chemistry Laboratory, ETH-Zurich) Support: MIUR (Italian Ministry for Education, University and Research) under Cofin program Silicon-based photonic crystals and FIRB project Miniaturized electronic and photonic systems INFM-CNR (National Institute for the Physics of Matter) under PRA (Advanced Research Project) PHOTONIC
2D photonic crystals embedded in planar waveguides (photonic crystal slabs) GaAs membrane (air bridge) Chow et al., Nature 407, 983 (2000) GaAs/AlGaAs PhC slab with internal source H. Benisty et al, IEEE J. Lightwave Technol. 17, 2063 (1999) W1 waveguide in Si membrane S.J. McNab et al, Opt. Express 11, 2923 (2003)
Novel effects with 3D control of light Guiding light through a 90 o bent waveguide: Loncar et al., APL (2000) 3D confinement of light ultra-high Q nanocavities (Q~10 6 recently measured) Song et al., Nature Mat. (2005)
Novel effects with 3D control of light Strong light-matter coupling with single QD emitters: Hennessy et al., Nature (2007) Electrically-driven PhC-based laser: Park et al., Science (2004)
Outline The Guided Mode Expansion (GME) method Theory, convergence tests, comparisons with other methods Overview of applications Disorder-induced losses in PhC slab waveguides PhC slab cavities Q-factors Triangular holes in PhC slabs Quantum theory of radiation-matter interaction
Reminder: eigenmodes in planar dielectric slabs 1D Helmholtz equation: d=260 nm Exact solution!
Photonic crystal slabs: dispersion, losses,, and the light-line line issue Periodically textured planar waveguide Modified light dispersion Truly guided mode
Basic ideas of guided-mode expansion (GME) Motivation: fast and accurate frequency-domain solver Photonic crystal slabs combine the features of 2D photonic crystals: control of light propagation in the xy plane by Bragg diffraction y Slab waveguides: control of light propagation in the vertical (z) direction by total internal reflection z x Electromagnetic field is a combination of 2D plane waves in xy and guided modes along z
Guided-Mode Expansion method Equation for the magnetic field: Expand magnetic field in a finite basis set, linear eigenvalue problem Hamiltonian 1 ω H( r) = ε ( r) c 2 H ν 2 μνcν = 2 2 H( r), Hr ( ) = μ c μ Hμ ( r) ω c c 1 H μν = ε ( )( * r H μ ( r )) ( H ν ( r )) dr μ H( r) = 0 Basis states H μ (r): guided modes of effective waveguide, where each layer is taken to have an average dielectric constant approximation 1 Off-diagonal components folding and splitting of PhC bands ε GG, ' Modes above the light line are coupled to radiative waveguide modes: Fermi s golden rule radiative losses
Effective waveguide Consider a PhC slab with semi-infinite claddings. For the basis states H μ (r), we choose the guided modes of an effective homogeneous waveguide with an average dielectric constant ε jin each layer j=1,2,3. z d/2 d/2 3 2 1 ε 3 ε 2 ε 1 We usually take ε j to be the spatial average of ε j (xy) ε j (ρ) in each layer: 1 ε j = cellε j( ρ) dρ, A=unit cell area A The average dielectric constants must fulfill ε 2 > ε1, ε3.
Diffraction losses by perturbation theory Modes above the light line are coupled to radiative PhC slab modes, yielding an imaginary part of mode frequencies: ω Im 2 k = π 2 c G', λ, j H k,rad 2 ρ j k + 2 ω G'; 2 c The matrix element between a quasi-guided and a radiative mode is 1 * rad H k, rad = ( Hk ( r)) ( H r r r k G λ j ( )) d ε ( ) + ',, The sum is over G reciprocal lattice vectors (out-of-plane diffraction channels) λ=te, TM polarization of radiative PhC slab modes j=1,3 radiative modes that are outgoing in medium 1,3 Approximation: radiative PhC slab modes are replaced with those of the effective waveguide..
Quality factor and propagation losses The quality factor of a photonic mode is defined as Spatial attenuation is described by an imaginary part of the wavevector: where Q = Im( k) v g = ω. 2Im( ω) Im( ω) =, dω dk is the mode group velocity. Propagation losses in db are obtained as v g α loss = 4.34 2 Im( k). To treat line and point defects (linear waveguides and nanocavities): use an in-plane supercell
Discussion of GME method Approximations (besides numerical ones): - For dispersion: the basis set of guided modes of the effective waveguide is not complete, since leaky modes are not included. - For losses: radiative PhC slab modes are replaced with those of the effective homogeneous waveguide. Advantages: - Guided and quasi-guided photonic modes are obtained, without any artificial layers (PML ) in the vertical direction. - For thin slabs a few guided modes are sufficient numerical effort is comparable to that of a 2D plane-wave calculation. - Diffraction losses by perturbation theory very efficient procedure, little additional numerical effort beyond dispersion. - Calculated values for losses are more accurate when they are small good for line defects and for high-q nanocavities.
2D periodicity: : the triangular lattice d/a=0.5, r/a=0.3 Symmetry plane: horizontal mid-plane (x,y) Even modes w.r.t. reflection through (x,y) have a band gap Defect modes with parity σ xy =+1 can be engineered!
Frequency ωa/(2πc) 0.5 0.4 0.3 0.2 10-2 Convergence tests Triangular lattice on membrane, ε=12.11, d/a=0.5, r/a=0.3 Γ M K (a) Frequency ωa/(2πc) 0.5 0.4 0.3 0.2 10-2 Γ M K (a) Im(ω)a/(2πc) 10-3 10-4 10-5 Im(ω)a/(2πc) 10-3 10-4 10-5 10-6 (b) 1 2 3 4 5 6 7 8 10-6 4 5 6 7 8 9 (b) 10 11 12 Number of guided modes N α Core effective ε Few guided modes are sufficient for convergence Use of average ε is justified
Comparison with MIT PBG code Triangular lattice on membrane, ε=12, d/a=0.6, r/a=0.45 0.7 0.6 Frequency ωa/(2πc) 0.5 0.4 0.3 0.2 0.1 odd, σ xy = 1 even, σ xy =+1 0.0 Γ M K Γ S.G. Johnson et al., PRB 60, 5751 (1999) Guided-mode expansion
Comparison with FDTD* Triangular lattice on membrane, ε=11.56, d/a=0.65 r/a=0.25 FDTD 0.5 GME 0.025 ωa/2πc 0.4 0.3 0.2 0.1 σ kz =+1 σ kz =-1 Q -1 (2 Im(ω)/ω) 0.020 0.015 0.010 0.005 4 even, σ kz = +1 odd, σ kz = 1 3 1 2 2 1 3 4 0.0 Γ K 0.000 M Γ K *T. Ochiai and K. Sakoda, PRB 63, 125107 (2001)
Losses in W1 waveguide: comparisons 1.1 Membrane, ε=11.56, h=0.6a, r=0.3016a Energy (ev) 1.0 0.9 0.8 0.7 σ kz = 1 σ kz =+1 0.6 0.0 0.2 0.4 0.6 0.8 1.0 Wavevector ka/π Loss (db/mm) 100 80 60 40 20 0 FDTD & Fourier: M. Cryan et al., IEEE-PTL 17, 58 (2005) GME -20 1300 1350 1400 1450 1500 1550 Wavelength (nm)
Applications of GME method Photonic band dispersion and gap maps: 1D and 2D lattices Linear waveguides: propagation losses Nanocavities: Q-factors optimization Geometry optimization, design, comparison with expts Extrinsic losses: effects of disorder Radiation-matter interaction: exciton-polaritons, weak/strong coupling Main reference on the method (see also refs. therein): L.C. Andreani and D. Gerace, Phys. Rev. B 73, 235114 (2006) The code is freely available on the web at http://fisicavolta.unipv.it/dipartimento/ricerca/fotonici/
Disorder-induced induced losses in W1 waveguides ωa/2πc 0.32 0.30 0.28 d/a=0.5 r/a=0.28 Δr/a=0.005 Δr/a=0.01 Δr/a=0.02 Δr/a=0.04 Δr/a=0.08 0.26 0.0 0.5 1.0 k x a/π 0.0 0.1 0.2 0.3 v g /c 10-6 10-5 10-4 10-3 10-4 10-3 10-2 10-1 10 0 Im(ωa/2πc) α loss a (db) Δr r.m.s. deviation from average hole radius Δr/a<<1 The losses show a quadratic behavior as a function of the disorder parameter typical of Rayleigh scattering Gerace & Andreani, Optics Letters (2004)
Comparison with experimental results* Experiment Transmission (db) 0-10 -20-30 Present theory 0.180 mm 0.300 mm 0.029 mm -40 2 mm 0.700 mm 0,320 0,295 0,300 0,305 0,310 Frequency ωa/2πc membrane W1 waveguide, d/a=0.494, r/a=0.37, a=445 nm, Δr=5 nm Frequency ωa/2πc 0,315 0,310 0,305 0,300 0,295 0,290 0,2 0,4 0,6 0,8 1,0 Wavevector ka/π *S.J. McNab, N. Moll, Yu. Vlasov, Optics Express 11, 2923 (2003) 2,7 db/mm 1 10 100 1000 Diffraction losses (db/mm)
Photonic crystal nanocavities Point defects in PC slabs behave as 0D cavities with full photonic confinement in very small volumes. The Q-factor can be increased by shifting the position of nearby holes. T. Akahane, T. Asane, B.S. Song, and S. Noda, Nature 425, 944 (2003). V=6.5 10 10-17 cm cm 3 52 26 0 Q (10 3 )
Q-factor of photonic crystal cavities 10 5 L3 defect L2 defect L1 defect expt.(*) L3 10 4 Q-factor L2 10 3 L1 10 2 0.0 0.1 0.2 0.3 0.4 Displacement Δx/a Supercell in two directions + Golden Rule Im(ω) and Q=ω/2Im(ω) *Expt: T. Akahane et al., Nature 425, 944 (2003). Andreani, Gerace, Agio, Photonics and Nanostructures 2, 103 (2004)
Effect of size disorder on cavity Q-factors 10 5 σ/a=0 σ/a=0.005 σ/a=0.005 σ/a=0 10 5 Q-factor 10 4 σ/a=0.010 σ/a=0.010 10 4 σ/a=0.015 σ/a=0.020 10 3 0.0 0.1 0.2 0.3 Displacement Δx/a σ/a=0.015 σ/a=0.020 0.0 0.1 0.2 Shrinkage Δr/a 10 3 Gerace & Andreani, Photon. Nanostructures 3, 120 (2005)
Triangular lattice with triangular holes 0.5 a L Frequency ωa/(2πc) 0.4 0.3 0.2 0.1 0.34 0.68 σ 0.0 xy =+1 σ xy = -1 Γ M K Γ M K Γ 0.66 Frequency ωa/(2πc) 0.32 0.30 0.28 even (σ xy =+1) gap 0.26 0.24 Γ M K Γ 0.64 0.62 0.60 0.58 0.56 odd (σ xy =-1) gap 0.54 0.52 0.50 0.48 Γ M K Γ d/a=0.68, L/a=0.85: no complete photonic gap (closed by 2 nd -order mode) d/a=0.5, L/a=0.8: even gap at ω, odd gap at 2ω Takayama et al., APL 87, 061107 (2006); Andreani&Gerace, PRB 73, 235114 (2006)
PhC slab with embedded quantum wells Quantum-well exciton resonance in interaction with photonic modes possibility of strong-coupling regime or photonic crystal polaritons (analogous to polaritons in bulk crystals or in planar microcavities)
Weak and strong coupling regimes: radiative PhC slab polaritons Energy (ev) 1.60 1.55 1.50 1.45 1.40 (a) 5 4 3 1.490 1.488 1.486 1.484 (b) 4 5 4 1.35 1.30 M Photonic dispersion QW Exciton Light lines Γ X 1.482 1.480 M Γ Photon Exciton Polaritons X Strong coupling regime with vacuum Rabi splitting occurs when photonic mode linewidth is smaller than exciton-photon coupling ~ a few mev Gerace & Andreani, PRB 75, 235325 (2007)
Experimental configuration and parametric processes: : polariton stimulation? 4 3 E-E ex (mev) 2 1 0-1 Photonic mode LP UP QW exciton -2-3 -20-10 0 10 20 Angle, θ (deg) Gerace & Andreani, PRB 75, 235325 (2007)
More on GME method
Photonic crystal slabs: structures Air bridge GaAs/AlGaAs Silicon-on-Insulator (SOI)
Expansion in the basis of guided modes H k ( r) guided, = G, α c( k + G, α ) H k+ G α ( r) k G α guided H k+ G, α ( r) = Bloch vector, chosen to be in the first Brillouin zone = reciprocal lattice vectors of the 2D Bravais lattice = index of guided mode = guided modes of effective waveguide at k+g Dimension of linear eigenvalue problem is (N PW N α ) (N PW N α ). Since the guided modes are simple trigonometric functions, the matrix elements can be calculated analytically. Both TE and TM guided modes appear in the expansion of a PhC slab mode. The index α includes mode order and polarization.
Matrix elements: the dielectric tensor The matrix elements between guided modes depend on the inverse dielectric matrices in the three layers, 1 η j ( G, G') = ε j A 1 ( ρ) e i( G' G) ρ Like in 2D plane-wave expansion, they are evaluated from numerical inversion of the direct dielectric matrices: η j ( G, G') = [ ε ( G, G')] 1 ε j ( G, G') = ε j ( ρ) e A 1 i( G' G) ρ Convergence properties, possible optimization beyond the Ho-Chan- Soukoulis procedure are similar to 2D plane-wave expansion. The off-diagonal components of η j (G,G ) are responsible for the folding of photonic bands in the first BZ, gap formation, splittings etc. j dρ dρ
Discussion of GME method (dispersion) Main drawback: the basis set of guided modes of the effective waveguide is not complete, since leaky modes are not included the GME method is an APPROXIMATE one Main advantage: folded photonic modes in the first Brillouin zone may fall above the light line guided and quasi-guided photonic modes are obtained, without the need of introducing any artificial layers (PML ) in the vertical direction. When a few guided modes are sufficient, as it often happens, the numerical effort is comparable to that of a 2D plane-wave calculation. Very suited for parameter optimization, design The dispersion of photonic bands in a PhC slab can be compared with that of the ideal 2D system and with the dispersion of free waveguide modes. Multimode slabs are naturally treated.
Symmetry properties, TE/TM mixing When the PhC slab is symmetric under reflection in the xy plane: separation of even (σ xy =+1) and odd (σ xy =+1) modes k = kˆ x Low-lying photonic modes are dominated by lowest-order modes of the effective waveguide: σ xy =+1 modes are dominated by TE waveguide modes quasi-te σ xy = 1 modes are dominated by TM waveguide modes quasi-tm However, in general, any PhC mode contains both TE and TM waveguide modes only the symmetry labels σ xy =±1 have a rigorous meaning.
Coupling matrix element Using the guided-mode expansion for the magnetic field of a PhC slab mode, the matrix element with a radiative mode is where * k, rad = k G, α ) G, α H c( + H 1 * rad guided, rad H k, rad = ( Hk ( r)) ( H r r r k G λ j ( )) d ε ( ) + ',, are matrix elements between guided and radiation modes of the effective waveguide, which are calculated analytically. Both guided and radiation modes are normalized according to H * μ ( r) Hν ( r) dr = δμν (within a large box for radiation modes).
Micro-roughness model Each hole is subdivided into N sections of length L c =2πr/N, where L c is a correlation length of size fluctuations. Micro-roughness is characterized by two parameters: σ and L c, in analogy to the Marcuse-Payne-Lacey model for strip waveguides. σ r L c x z y y W Si strip x σ d SiO 2 cladding W L c