Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Design and analysis of two-dimensional high-index-contrast grating surface-emitting lasers

Open Access Open Access

Abstract

The use of a two-dimensional (2D) high-index-contrast grating (HCG) with square periodic lattice is proposed to realize surface-emitting lasers. This is possible because the use of 2D HCG, in which multiple resonant leaky modes are excited by the 2 orthogonal directions of the grating, causes the high reflective zone to be split into two regions. Hence, a dip of the reflectivity is formed to support the excitation of a resonant cavity-mode inside the 2D HCG. With suitable design on the dimensions of the 2D HCGs, Q factor as high as 1032 can be achieved.

©2009 Optical Society of America

1. Introduction

Recently, extensive studies have been concentrated on the use of one-dimensional (1D) high-index-contrast sub-wavelength gratings (HCGs) as the broadband reflectors [1, 2]. This is because the HCGs can provide a high reflectivity (i.e., > 99%) over a broadband region (i.e., > 15% of the operating wavelength). Furthermore, the realization of HCGs, which consist of a single-layer 1D-grating (i.e., periodic stripes with high-index-contrast) sandwiched in-between two low-index cladding layers, can be relatively simple. In contrast to HCGs, the fabrication of highly reflective distributed Bragg reflectors (DBRs) is difficult due to the epitaxial growth constraints of semiconductor dielectric layers with low-index-contrast. Hence, the use of HCG was proposed to replace the top DBR mirror of vertical-cavity surface-emitting lasers (VCSELs) [3]. In fact, the design and fabrication of VCSELs can be further simplified if the use of DBRs can be completely avoided.

High reflectivity of 1D-HCGs can be explained by the mechanism of guided-mode resonance (GMR). The periodic modulation of refractive index generates counter-propagating leaky modes inside the grating structure. Under phase matching conditions, these counter-propagating leaky modes establish standing waves and give rise to GMR. [4–6]. As a result, the leaky modes re-radiate reflectively through the resonant interaction with the grating. It is noted that the high reflective zone of the 1D HCGs is arisen from the blend of the adjacent leaky modes [7]. Hence, if the number and spectral location of the leaky modes is controlled, it is possible to introduce a dip to the high reflective zone (i.e., a cavity mode of the HCGs). In this paper, we propose the use of a two-dimensional (2D) HCG, in which a single-layer grating is constructed by a 2D high-index-contrast of square periodic lattice, to realize a dip in the high reflective zone of the 2D HCG. Hence, it can be shown that VCSEL without using DBRs can be formed by a 2D HCG with the introduction of an active layer.

2. Design and analysis

Figure 1 shows the schematic of a̱ proposed 2D HCG considered in our studies. It is assumed that the square periodic lattice is constructed by a layer of high-index periodic lattice with refractive index of nH. The square periodic lattice is supported by a layer of low-index buffer layer, which has refractive index and thickness of n L and t buf respectively, laid on a substrate. The top surface of the square periodic lattice is surrounded by air. In order to understand how laser cavity can be formed within the 2D HCGs, the corresponding reflection and transmission characteristics are analyzed and compared with that of the 1D HCGs. It is assumed that the structure of 1D HCGs is identical to that of the 2D HCGs except that the 2D grating (i.e., square periodic lattice) is replaced by infinite long periodic stripes with the same period.

 figure: Fig. 1.

Fig. 1. Schematic of a 2D HCG. Thickness of grating layer, t g, was set to 460 nm. The refractive indices of grating layer, buffer layers and substrate were set to n H = 3.48, n L = 1.45 and n sub = 3.48 respectively.

Download Full Size | PDF

Figure 2 shows the reflection and transmission spectra of the 1D and 2D HCGs. The period, ∧, and height, tg, of the 1D- and 2D-gratings were set to 700 and 460 nm respectively and the corresponding thickness of buffer layer, t buf, was set to infinite in the calculation. In addition, Poly-silicon (n H = 3.48) and SiO2 (n L = 1.45) were assumed to be high-index layer and buffer layer respectively with negligible absorption loss. A 3-D finite-difference time-domain (FDTD) method was used to investigate the reflection and transmission characteristics of the HCGs [8]. In addition, periodic boundary condition was used in the FDTD calculation to approximate the periodicity of the 1D- and 2D- gratings.

 figure: Fig. 2.

Fig. 2. Plot of (a) reflection and (b) transmission spectra of 1D and 2D HCGs. The HCGs have same period (∧ = 700 nm) and thickness of gratings (t g = 460 nm) and were supported by an infinite thick buffer layer (i.e., t buf = ∞).

Download Full Size | PDF

1D HCG with fill factor of 75% shows broadband reflection for the illumination of TM polarized light. Two closely spaced leaky modes (at 1.49 and 1.63 μm) can be seen from the transmission spectrum. These modes interact with the waveguide grating leading to efficient reflection [9]. Each mode corresponds to a high reflective mode (i.e., transmittance approaching zero) due to GMR. For the illumination of TE polarized light, the reflection spectrum is a narrow band. This is evident from the corresponding transmission spectrum which shows wide spacing of the weak leaky modes. For 1D HCG with a fill factor of 65%, one extra isolated leaky mode appears at a shorter wavelength (1.229 μm). As a result, the wide spectral spacing of the leaky modes contributes to a shallow and wide dip in the reflection spectrum. On the other hand, for the case of 2D HCG with a fill factor of 65%, four leaky modes appeared in the transmission spectrum. It is noted that the leaky modes #1 (1.239 μm), #2 (1.306 μm) and #3 (1.48 μm), #4 (1.612 μm) contribute to the high reflective regions centered at ~1.3 and ~1.6 μm respectively. Mode #2 is an extra leaky mode appearing in the transmission spectra of the 2D HCG when compared to that of the 1D HCG with the same fill factor. Furthermore, the presence of a mode #2 leads to a significant drop of the reflectivity at ~1.343 μm (i.e., the reflection dip). In other words, the presence of multiple modes, which cause the separation of GMR location, is due to 2D nature of the grating. As a result, a dip in the reflection spectrum is obtained. This dip of the reflectivity represents the formation of a resonant cavity mode of the 2D HCG.

In order to verify that the formation of cavity mode is due to the 2D nature of the grating instead of Fabry-Perot resonance, the values of t buf is allowed to vary between 775 and 925 nm (see Fig. 3). I̱t is observed that when t buf increases from 775 to 925 nm, the reflectivity at the dip wavelength of ~1.343 μm increases from 7% to ~ 25%. However, the minimum value of the reflection dip (~3%) occurs only when t buf = ∞. This suggests that the dip in the reflection spectrum is due to the 2D nature of the grating and the cavity mode is not formed inside the buffer layer. However, the buffer layer has influence on the fineness of the reflection dip.

 figure: Fig. 3.

Fig. 3. Reflection spectra of 2D HCGs with varying thickness of buffer layer (i.e., t buf varies between 775 and 925 nm).

Download Full Size | PDF

It can be shown that there is a proportional relationship between the wavelength of the reflection dip and the dimensions of the grating. For a 1D HCG, the resonant conditions of resonant leaky modes can be determined by eigen-equations [10]:

tan(κitg)=κi(γi+δi)/(κi2γiδi)(TEmode)
tan(κitg)=ng2κi(nL2γi+δi)/(nL2κi2ng4γiδi)(TMmode)

where i is an integral denoting the diffraction order of grating, βi (= k[n gsinθ - /∧]) is effective propagation constant of the waveguide grating, k (= 2π/λ) is the wavevector in free space, θ is the incidence angle to the normal of the grating, λ is the resonant wavelength, and n g is the effective reflective index of the grating layer. The parameters κi(=ng2k2βi2),γi(=βi2k2), and δi(=βi2nL2k2) are defined as the wave vectors within the grating layer, air and buffer layer respectively. If the resonant wavelength λ increases to by a factor of Y, it can also be shown that the values of β i, κ i, γ i and δ i can be reduced to β i/Y, κ i/Y γ i/Y and δ i/Y respectively if θ = 0 and ∧ increases by a factor of Y. In order to satisfy Eq. (1) for this new resonant wavelength, t g has to be increased to Yt g. Hence, there is a proportional relationship between resonant wavelength and the dimensions of the grating. For our 2D HCGs with the dimensions of grating identical in both x and y directions (see Fig. 1), the corresponding eigen-equations and phase matching conditions are similar to that of the 1D grating. Therefore, the proportional relationship should also be held for the 2D HCGs. In order to demonstrate this scalable property of the 2D HCGs numerically, a design with a reflection dip located at a wavelength of 1.343 μm was used for the calculation. The corresponding values of, ∧, fill factor, t g and t buf, were set to 700 nm, 65%, 460 nm and 825 nm respectively. Reflection dip at wavelength of 1.55 and 2 μm can be achieved by scaling up ∧, t g and t buf, altogether by a factor 1.157 and 1.489 respectively as shown in Fig. 4. It can be seen from Fig. 4 that the wavelength of reflection dips and reflection patterns are similar for all the three cases under investigation.

 figure: Fig. 4.

Fig. 4. Plot of the reflection spectra of the HCGs versus normalized wavelength (i.e., normalized to the wavelength of the reflection dip). For the reflection dip at 1.343 μm, the dimensions of 2D HCG were set to ∧ = 700 nm, fill factor = 65%, t g = 460 nm and t buf = 825 nm. For a reflection dip at 1.55 and 2 μm, the corresponding dimensions of 2D HCG were scaled by a factor of 1.157 and 1.489 respectively

Download Full Size | PDF

3. Surface emitting laser with 2D HCG

Conventional 1D 2nd order grating structure using low-index contrast has been used to fabricate surface emitting lasers [11]. The proposed 2D grating can also be applied to realize surface-emitting lasers by inserting a thin active layer between the grating and buffer layers. However, the 2D grating structure, which need not to be a 2nd order grating, is necessary to have high-index-contrast. This is because high-index-contrast between the grating and buffer layers is essential to achieve high reflective zones and a sharp reflection dip. Reduced index contrast below 2 (i.e., n H/n L < 2) may lead to loss in reflectivity and poor Q factor of the device. Hence, refractive index and peak gain wavelength of the active layer should be close to the buffer layer and wavelength of the reflection dip respectively. If the buffer layer was assumed to be made by SiO2, which has a refractive index of 1.45 at a wavelength of 1.343 μm, Nd:YAG can be used as an active layer. Hence, the refractive index, n act, and the peak gain wavelength of the active layer can be set to be 1.82 and 1.343 μm respectively [12] in the calculation. In addition, the thickness of the active layer, t act, was assumed to be 100 nm. In our analysis, a Lorentz model [8] was used to implement the gain material into FDTD algorithm. In this case, the corresponding permittivity, ε, can be expressed as

ε(f)=εo+εlorωo2/(ωo2+22πf(2πf)2)

where εo = 3.3124, ωo = 1.4×103 THz (~1.343 μm), and δ = 31.4 THz (~ 50 nm). The emission spectra observed from the surface of the 2D HCGs can be calculated by a Fourier transform of the electric fields.

Two configurations of the 2D-HCG laser, i.e., with 1) t buf = 725 nm and 2) t buf = ∞, are considered in the studies. Fig. 5 plots the emission spectra of the 2D-HCG laser at threshold. Sharp resonance is observed at ~1.34 μm for the case t buf = ∞. However, the resonance wavelength shifted to ~1.342 μm for the case t buf = 725 nm. This is because the effective refractive index of the buffer layer is slightly increased due to the influence of substrate (which is assumed to have refractive index larger than the buffer layer) so that the resonant frequency of the standing wave inside the 2D grating is slightly shifted. For the case of t buf = ∞, Q factor is found to be ~925. This large value of Q is expected from the reflection spectrum of the 2D HCG with t buf = ∞ as the corresponding reflectivity is the lowest at the dip wavelength (see Fig. 3). On the other hand, although the presence of buffer layer will increase the reflectivity of the reflection dip, the reduction of 10% linewidth of the reflection dip (see Fig. 3) can compensate for the increase of reflectivity so that the value of Q factor can be improved. The inset of Fig. 5 shows the Q factor of the 2D HCGs with different values of t buf. The corresponding reflectivity and 10% linewidth of the reflection dip are also plotted in the figure. It is observed that the magnitude of reflection increases with t buf but the linewidth of the reflection dip reaches a minimum for t buf between 725 and 775 nm. As a result, a maximum value of Q factor equal to 1032 can be obtained at t buf = 725 nm. However, the Q factor does not improve for the value of t buf other than 725 nm.

 figure: Fig. 5.

Fig. 5. Surface emission spectra of 2D HCG lasers with t buf of 725 nm and infinite thick. The inset shows the reflectivity and 10% linewidth (defined as the width of the dip spectrum at a level 10% above the reflection dip, see Fig. 3) of the reflection dip, as well as the Q factor of the 2D HCGs versus t buf.

Download Full Size | PDF

4. Conclusions

In conclusion, the reflection characteristics of 2D HCGs were investigated. With suitable selection of ∧, fill factor and t buf of the 2D HCGs, cavity resonant (i.e., a dip of the reflectivity) can be observed from the corresponding reflection spectrum. The excitation of cavity mode (i.e., the dip) in 2D HCGs can be explained by the excitation of extra resonant leaky mode (#2 mode in Fig. 2(b)) causing the formation of two high reflective zones in the reflection spectrum. It is noted that there is a very narrow range of ∧, fill factor and t g to generate the cavity mode. By inserting a thin active layer between the grating and buffer layers, a 2D HCG surface-emitting laser with Q factor as high as 1032 can be realized and the value of Q factor is dependent on thickness of the buffer layer. In addition, in contrast to 1D HCGs, symmetric nature of 2D HCG makes it insensitive to the incident polarization. Such property may also be well suited for applications in amplifiers and optical filters.

Acknowledgments

This work was supported by Singapore Ministry of Education Academic Research Fund Tier 2, grant no. ARC 2/06 and A*Star SERC grant 072 101 0023.

References and links

1. A. E. Willner, “All mirrors are not created equal,” Nature Photon. 1, 87–88 (2007). [CrossRef]  

2. C. F. R. Mateus, M. C. Y. Huang, Y. Deng, A. R. Neureuther, and C. J. Chang-Hasnain, “Ultrabroadband Mirror Using Low-Index Cladded Subwavelength Grating,” IEEE Photon. Technol. Lett. 16, 518–520 (2004). [CrossRef]  

3. M. C. Y. Huang, Y. Zhou, and C. J. Chang-Hasnain, “A nanoelectromechanical tunable laser,” Nat. Photonics 2, 180 (2008). [CrossRef]  

4. Y. Ding and R. Magnusson, “Resonant leaky-mode spectral-band engineering and device applications,” Opt. Express 12, 5661–5674 (2004). [CrossRef]   [PubMed]  

5. A. Hardy, D. F. Welch, and W. Streifer, “Analysis of second-order gratings,“ IEEE J. Quantum Electron. 25, 2096–2105 (1989). [CrossRef]  

6. Y. Ding and R. Magnusson, “Band gaps and leaky-wave effects in resonant photonic-crystal waveguides,” Opt. Express 15, 680–694 (2007). [CrossRef]   [PubMed]  

7. R. Magnusson and M. Shokooh-Saremi, “Physical basis for wideband resonant reflectors,” Opt. Express 16, 3456–3462 (2008). [CrossRef]   [PubMed]  

8. A. Taflove and S. C. Hagness, Computational Electrodynamics: the finite difference time-domain method (Artech House, Boston, 2005).

9. D. Rosenblatt, A. Sharon, and A. A. Friesem, “Resonant grating waveguide structures,” IEEE J. Quantum Electron. 33, 2038–2059 (1997). [CrossRef]  

10. S. S. Wang and R. Magnusson, “Theory and applications of guided-mode resonance filters,” Appl. Opt. 32, 2606–2613 (1993). [CrossRef]   [PubMed]  

11. T. Kobayashi, Y. Kanamori, and K. Hane, “Surface laser emission from solid polymer dye in a guided mode resonant grating filter structure,” Appl. Phys. Lett. 87, 151106 (2005). [CrossRef]  

12. D. Sauder, A. Minassian, and M. J. Damzen, “Laser operation at 1.3 μm of 2at% doped crystalline Nd:YAG in a bounce geometry,“ Opt. Express 15, 3230–3235 (2007). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. Schematic of a 2D HCG. Thickness of grating layer, t g, was set to 460 nm. The refractive indices of grating layer, buffer layers and substrate were set to n H = 3.48, n L = 1.45 and n sub = 3.48 respectively.
Fig. 2.
Fig. 2. Plot of (a) reflection and (b) transmission spectra of 1D and 2D HCGs. The HCGs have same period (∧ = 700 nm) and thickness of gratings (t g = 460 nm) and were supported by an infinite thick buffer layer (i.e., t buf = ∞).
Fig. 3.
Fig. 3. Reflection spectra of 2D HCGs with varying thickness of buffer layer (i.e., t buf varies between 775 and 925 nm).
Fig. 4.
Fig. 4. Plot of the reflection spectra of the HCGs versus normalized wavelength (i.e., normalized to the wavelength of the reflection dip). For the reflection dip at 1.343 μm, the dimensions of 2D HCG were set to ∧ = 700 nm, fill factor = 65%, t g = 460 nm and t buf = 825 nm. For a reflection dip at 1.55 and 2 μm, the corresponding dimensions of 2D HCG were scaled by a factor of 1.157 and 1.489 respectively
Fig. 5.
Fig. 5. Surface emission spectra of 2D HCG lasers with t buf of 725 nm and infinite thick. The inset shows the reflectivity and 10% linewidth (defined as the width of the dip spectrum at a level 10% above the reflection dip, see Fig. 3) of the reflection dip, as well as the Q factor of the 2D HCGs versus t buf.

Equations (3)

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

tan(κitg)=κi(γi+δi) / (κi2γiδi) (TEmode)
tan(κitg)=ng2κi(nL2γi+δi) / (nL2κi2ng4γiδi) (TMmode)
ε(f)=εo+εlorωo2/(ωo2+22πf(2πf)2)
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.