• Photonics Research
  • Vol. 6, Issue 5, A38 (2018)
Ziyao Feng, Jingwen Ma, Zejie Yu, and Xiankai Sun*
Author Affiliations
  • Department of Electronic Engineering, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong, China
  • show less
    DOI: 10.1364/PRJ.6.000A38 Cite this Article Set citation alerts
    Ziyao Feng, Jingwen Ma, Zejie Yu, Xiankai Sun. Circular Bragg lasers with radial PT symmetry: Design and analysis with a coupled-mode approach[J]. Photonics Research, 2018, 6(5): A38 Copy Citation Text show less

    Abstract

    Parity–time (PT) symmetry has been demonstrated in the frame of classic optics. Its applications in laser science have resulted in unconventional control and manipulation of resonant modes. PT-symmetric periodic circular Bragg lasers were previously proposed. Analyses with a transfer-matrix method have shown their superior properties of reduced threshold and enhanced modal discrimination between the radial modes. However, the properties of the azimuthal modes were not analyzed, which restricts further development of circular Bragg lasers. Here, we adopt the coupled-mode theory to design and analyze chirped circular Bragg lasers with radial PT symmetry. The new structures possess more versatile modal control with further enhanced modal discrimination between the azimuthal modes. We also analyze azimuthally modulated circular Bragg lasers with radial PT symmetry, which are shown to achieve even higher modal discrimination.

    1. INTRODUCTION

    Photonic ring and disk resonators have received significant attention because they are widely used for both fundamental research as well as practical applications in optical communications and signal processing [1]. A method to realize high optical quality factors is by confining light to a “defect” region defined by distributed Bragg reflectors. Such resonators are adopted for semiconductor lasers and have been studied with various techniques such as conformal mapping, coupled-wave approach, and transfer-matrix method [210]. Previously, much research attention was given to the radial modes, and the effort was focused on designing the structures for larger modal discrimination between modes of different radial orders. In order to achieve single-mode operation, it is equally important to achieve large discrimination between different azimuthal modes. The existing approaches to overcome the challenge include adoption of radially chirped Bragg gratings [7] and azimuthal modulation [8].

    “Parity–time (PT) symmetry” was first introduced by Bender et al. in quantum mechanics [11]. Due to the similarity between the Schrödinger equation and the Helmholtz equation, the concept of PT symmetry has been adopted in optics [1214]. The PT-symmetric and PT-broken phases are separated by the exceptional point, where the system’s eigenvalues coalesce. Interesting phenomena by mode management around the exceptional points include single-mode lasers [1520], reversing pump dependence [2123], orbital angular momentum lasers [24], and unidirectional reflection [25,26].

    In 2016, we proposed PT-symmetric circular Bragg lasers with radially periodic gratings. We employed a transfer-matrix method to calculate their reflection and transmission coefficients and analyzed their modal properties [16]. We found enhanced modal discrimination among the radial modes. However, the eigenwaves in circular Bragg lasers take the form of Hankel functions, whose zeros are not equidistantly distributed. In order to achieve perfect phase matching in the radial direction, the phase of the circular gratings of the laser structure has to follow that of the Hankel functions, rendering nonperiodicity (or chirping) in the radial direction. Here, we develop a coupled-mode approach for designing a PT-symmetric circular Bragg laser with the radially chirped gratings. We also investigate the effect of an additional azimuthal modulation of the refractive index and analyze the threshold conditions for different azimuthal modes. The rest of this paper is structured as follows. In Section 2, we focus on radially chirped circular Bragg lasers and derive a set of coupled equations. We compare the threshold gain of conventional and PT-symmetric radially chirped circular Bragg lasers. The PT-symmetric circular Bragg lasers possess the higher modal discrimination. Such lasers are also shown to be highly robust, as a small deviation away from the exceptional point does not affect the modal discrimination. In Section 3, we include an additional azimuthal modulation to the radially chirped PT-symmetric circular Bragg lasers. By simulating the structures with an finite-difference time-domain (FDTD) method, we find that the structures with an additional azimuthal modulation can effectively suppress other unwanted modes.

    2. PT-SYMMETRIC CIRCULAR BRAGG LASERS WITH RADIALLY CHIRPED GRATINGS

    We previously proposed a laser structure with its complex refractive index in the radial direction distributed periodically. We used a transfer-matrix method to analyze its threshold gain [16] and found that the threshold gain for different radial modes could have large difference. However, because circular Bragg lasers with radially periodic gratings cannot obtain large modal discrimination between the azimuthal modes [2], here we propose to apply the concept of PT symmetry to radially chirped gratings, as shown in Fig. 1(a). We suggest a structure that possesses the same dispersion relation as that of a planar PT-symmetric grating, where the refractive index distribution along the radial direction is shown in Fig. 1(c). We will compare modal properties of different types of circular Bragg lasers and show that the proposed structure performs better in controlling the azimuthal modes with a high modal discrimination.

    (a) Schematic of the PT-symmetric circular Bragg laser. (b) Schematic of the azimuthally modulated circular Bragg laser with radial PT symmetry. (c) Radial profile of the refractive index of the two types of lasers, showing the real (Δnr) and imaginary (Δni) refractive index modulations that satisfy the PT symmetry requirement.

    Figure 1.(a) Schematic of the PT-symmetric circular Bragg laser. (b) Schematic of the azimuthally modulated circular Bragg laser with radial PT symmetry. (c) Radial profile of the refractive index of the two types of lasers, showing the real (Δnr) and imaginary (Δni) refractive index modulations that satisfy the PT symmetry requirement.

    Let us consider a specific azimuthal mode with perturbation in the radial direction. Erdogan et al. derived the expression for different azimuthal modes of the TE polarization [2]. In contrast with Ref. [2], here we focus our analysis on the TM modes whose major field component is Ez. All the other electromagnetic field components of the TM modes can be obtained from Ez, which satisfies the Helmholtz equation Ez can generally be expressed as Ez=l=Ezlexp(ilφ), where Ezl=Al(r)Hl(1)(kr)+Bl(r)Hl(2)(kr) is the radial distribution of the electric field of the lth-order mode, l is the azimuthal mode number, r and φ are the radial and azimuthal coordinates, respectively, and k is the wavenumber, which becomes k0 in the vacuum and kav in the medium. The Hankel functions of the first (second) kind represent the radially outward (inward) waves, with Al and Bl being their amplitudes. With orthogonality of the function exp(ilφ), the equation for Hankel-phased circular Bragg resonators was derived by Scheuer et al. [8]: where Δϵ(r) can be expanded in terms of Hankel functions. The 0th-order Fourier expansion coefficient of Δϵ(r) determines how Al or Bl evolves due to self-contribution, and the ±1st-order coefficients, which contain exp[±2iphase(Hl(1)(kavr))] terms, determine the strength of coupling between Al and Bl. In Eq. (2), the derivatives of Al and Bl can be separated according to their spatial frequency in the radial direction. We suggest the ±1st-order perturbation of Δϵ(r) to be Δϵr(r)+iΔϵi(r) with where kde is the spatial frequency of the refractive index modulation. The phase shifts θ1 and θ2 are used for satisfying the PT symmetry requirements and for fine-tuning the laser wavelength such that kav=kde. We define δ as kavkde, which is the difference between the actual wavenumber and the designed spatial frequency. With the approximation Hl(1)(x+δx)Hl(1)(x)exp(iδx) and separation of the inward and outward waves, the coupled equations for the amplitudes are as follows:

    The above equations are typical coupled equations for chirped circular Bragg gratings. With κ=Δϵrk0/(4nav) and g=Δϵik0/(4nav) and expressing Al(r)=Rl(r)exp(iδr) and Bl(r)=Sl(r)exp(iδr), we obtain By expressing Rl(r)=D1exp(iγr)+D2exp(iγr) and similarly for Sl(r), we obtain the eigenvalues of the above equations γ2=κ2g2δ2+2iκgcos(θ1θ2). To ensure the same dispersion relation as that in a PT-symmetric planar Bragg grating, we should set θ1=θ2±π/2. Next, we focus our analysis on the 8th-order mode similar to that in Ref. [8].

    Let us suppose the refractive index is nav in the unperturbed region (r<r0) and is expressed as follows in the perturbed region (r0<r<rw): where sign(x) is 1 for positive x and 1 for negative x, and rw is the outer radius of the circular grating region. We also use N (the number of layers) to help describe the structure. Considering that Δn(r)=Δnr(r)+iΔni(r) and the perturbation is weak (Δnr, Δninav), we can approximate Δϵ(r) as 2nav[Δnr(r)+iΔni(r)] so that the Fourier expansion coefficients of Δϵ(r) respond linearly to Δn(r). The index profile in Eq. (6) results in the Fourier expansion Δϵ(r)=m=amexp[2miphase(H8(1)(kder))], where a0=0, a1=4nav[exp(iθ1)Δnr+iexp(iθ2)Δni]/π, a1=4nav[exp(iθ1)Δnr+iexp(iθ2)Δni]/π. In this case, κ=2k0Δnr/π and g=2k0Δni/π. We obtain the solutions for the 8th-order mode with Equation (8) reveals that, at the exceptional point (Δnr=Δni), if θ1=θ2π/2, then M12=0 and M11 is simplified as exp(iδL), yielding constant amplitude of the outward wave in the grating region. If θ1=θ2+π/2, then M21=0 and M22 is simplified as exp(iδL), yielding constant amplitude of the inward wave. In what follows, we focus on the case of θ=θ1=θ2π/2.

    Next, by applying the boundary conditions that the E field and its derivative must be continuous, and no inward wave should exist at r=rw, we obtain the reflection coefficient of the circular Bragg grating at r=r0, which is expressed as ρrexp(2iδr0) with ρr=M21/M22. We can further obtain the oscillation condition of circular Bragg lasers: ρrρ0exp(2iδr0+2αr0)=1. ρ0 denotes the reflection coefficient at r=0, which must be equal to 1, and α is the gain coefficient in the central region. In circular Bragg lasers, the central region (r<r0) provides the lasing modes with gain or loss to achieve the oscillation condition. Therefore, it is important to obtain the normalized threshold gain αr0=ln(|ρrρ0|)/2, where the resonant wavelength should satisfy that phase(ρr)+2δr0 be a multiple integer of 2π. For δ=0 and Δnr=Δni, the reflection coefficient at r=r0 can be simplified as i(κ+g)Lexp(iθ), which scales linearly with the length of the grating region. Additionally, the laser can lase at the desired frequency when θ is equal to π/2. It should be noted that the PT symmetry can be realized only for the specific order of Fourier expansion of the refractive index, while the refractive index does not have to follow Eq. (6).

    Now that we have designed radially chirped PT-symmetric circular Bragg lasers, we are ready to analyze the properties of their azimuthal modes. Due to different phase shifts in the approximation of even- and odd-order Hankel functions, the odd azimuthal modes are unlikely to be excited in lasers designed for even azimuthal modes [2]. By using the transfer-matrix method, which was verified by numerical solutions of the coupled equations, we calculated the normalized threshold gain (αr0) for even azimuthal modes for two types of circular Bragg lasers, namely, conventional radially chirped and PT-symmetric radially chirped, with the results plotted in Fig. 2(a). Here “conventional” refers to a structure with Δni=0 and Δnr=0.08, and “PT-symmetric” refers to a structure with Δni=Δnr=0.025, where a different Δnr was chosen to ensure a similar threshold gain for the targeted mode operating at the wavelength of 1.55 μm for both types of lasers. Other structural parameters include N=120, nav=3.4, and r0=1.164  μm.

    (a) Normalized threshold gain (αr0) for even-order azimuthal modes of PT-symmetric and conventional radially chirped circular Bragg lasers. Inset shows the E field profile of the targeted 8th-order azimuthal mode. (b) Effect of a small deviation away from the exceptional point on the normalized threshold gain values for the PT-symmetric radially chirped circular Bragg lasers.

    Figure 2.(a) Normalized threshold gain (αr0) for even-order azimuthal modes of PT-symmetric and conventional radially chirped circular Bragg lasers. Inset shows the E field profile of the targeted 8th-order azimuthal mode. (b) Effect of a small deviation away from the exceptional point on the normalized threshold gain values for the PT-symmetric radially chirped circular Bragg lasers.

    As shown in Fig. 2(a), in both types of lasers, the intended 8th-order azimuthal mode indeed possesses the lowest threshold gain and thus will be the lasing mode. However, the PT-symmetric laser structure achieves larger modal discrimination (i.e., the difference of threshold gain between the lasing mode and the other modes) than the conventional counterpart. The bottom line is that the PT-symmetric circular Bragg lasers can be designed to optimize the performance of a mode of any azimuthal order, with significantly enhanced modal discrimination compared with the conventional circular Bragg lasers.

    With the superior properties, the PT-symmetric circular Bragg lasers can be fabricated on a III–V wafer. Selective etching and metal deposition can modulate, respectively, the real (Δnr) and the imaginary (Δni) part of the refractive index. This method was employed to fabricate the PT-symmetric microring laser, without the problem of carrier diffusion to destroy the PT symmetry [15,24]. The effect of a small deviation away from the exceptional point was also studied, where the imaginary part (Δni) and the real part (Δnr) of the refractive index modulation are not equal. As shown in Fig. 2(b), increasing (decreasing) the imaginary modulation leads to an overall reduction (enhancement) of threshold gain for all modes, so the modal discrimination remains essentially unchanged.

    3. CIRCULAR BRAGG LASERS WITH RADIAL PT SYMMETRY AND AZIMUTHAL MODULATION

    In the previous section, a coupled-mode approach has been employed to analyze and design laser structures with modulation only in the radial direction. Because azimuthal modulation is usually employed to control azimuthal modes, in this section we propose to impose additional modulation in the azimuthal direction, as illustrated in Fig. 1(b). We will use the coupled-mode equations to analyze the modes and show the additional azimuthal modulation can also enhance the modal discrimination. To that end, we suggest the following refractive index profile: with where q denotes the number of periods in the azimuthal direction. The radial modulation in Eq. (9) is the same as that in Eq. (6). By following the approach in Ref. [8], under the weak perturbation condition, the responsible terms in the Fourier expansion of Δϵ(r,φ) are those of the 0th and 1st orders in the azimuthal direction and the 1st order in the radial direction. The azimuthal modulation defined by Θ(φ) actually introduces coupling between the lth- and (l±q)th-order azimuthal modes. From the previous work by Erdogan and Hall [2], we can analyze the structure with coupled-mode theory. For the PT-symmetric structure with refractive index defined in Eq. (9), the Fourier expansion of Δϵ(r) yields Δϵ(r,φ)=h=f=sf,he2hiqφexp[2fiphase(H8(1)(kder))] with s1,0=2nav(Δnr+Δni)/π, s1,1=s1,1=4nav(Δnr+Δni)/π2, s1,0=2nav(ΔnrΔni)/π, s1,1=s1,1=4nav(ΔnrΔni)/π2, and f(h) being the order of Fourier expansion in the radial (azimuthal) direction. These parameters determine the coupling coefficients between the different modes. It should be noted that the coupling strength relies not only on the coefficient but also on an overlap integral involving the modal profiles in the radial direction. Without loss of generality, we choose q=16, which leads to strong coupling between the 8th-order mode and the 8th-order mode, which have the same profile in the radial direction.

    We compared the reflectivity of the same azimuthal mode for two laser structures (without and with azimuthal modulation) by solving the coupled-mode equations numerically. The two structures share the same r0 (1.03 μm) and N (120). The modulation depth (Δnr=Δni) is set to be 0.025 and 0.03 to ensure the same reflectivity of the 8th-order mode for the laser with and without the azimuthal modulation, respectively. Assuming that the amplitudes of the outward wave of different modes are the same, we can obtain the two highest-reflectivity modes of the PT-symmetric laser without azimuthal modulation: 0.804 for the 1st-order mode and 0.904 for the 8th-order mode. The reflectivity of the 8th-order mode of the structure with azimuthal modulation remains 0.904, but that of the 1st-order mode reduces to 0.366. This suppression is attributed mainly to the mismatch of the modal profiles between the 1st-order mode and the 15th-order mode in the radial direction. As previously mentioned, such analysis is based on a 2D model. Some other effects of azimuthal modulation are not considered here. For example, the diffracted waves can be coupled with radiation modes, which results in reduced modal discrimination [27].

    We also simulated the two laser structures in Lumerical by an FDTD method [28]. Figure 3 plots the simulated laser output spectra with the modal field distribution for the resonant modes. As predicted by the coupled-mode theory, the 8th-order mode plays a major role, while the additional azimuthal modulation introduces significant mode suppression to the 1st-order mode. It is clear that the azimuthal modulation structures proposed in this section not only enjoy the benefit from the PT-symmetric radial gratings but also get further suppression to the unwanted modes from the deliberately designed azimuthal modulation.

    (a) Output spectra for two circular Bragg laser structures with radial PT symmetry. The upper corresponds to the one without azimuthal modulation, while the lower corresponds to the one with azimuthal modulation. (b) E field profiles of the modes at the three peaks in the laser output spectra. Modes (i) and (ii) from the upper spectrum are of the 8th and 1st azimuthal order, respectively. Mode (iii) from the lower spectrum is of the 8th azimuthal order.

    Figure 3.(a) Output spectra for two circular Bragg laser structures with radial PT symmetry. The upper corresponds to the one without azimuthal modulation, while the lower corresponds to the one with azimuthal modulation. (b) E field profiles of the modes at the three peaks in the laser output spectra. Modes (i) and (ii) from the upper spectrum are of the 8th and 1st azimuthal order, respectively. Mode (iii) from the lower spectrum is of the 8th azimuthal order.

    4. CONCLUSION

    We have adopted a coupled-mode approach to design and analyze circular Bragg lasers with radial PT symmetry. We analyzed and compared the threshold gain for different azimuthal modes of circular Bragg lasers with radially chirped gratings for both conventional and PT-symmetric structures. We have shown that the PT-symmetric radially chirped circular Bragg lasers have higher modal discrimination. Further suppression to unwanted modes is obtained by introducing an additional azimuthal modulation to the circular Bragg lasers with radial PT symmetry. With an intrinsic broadside circular aperture, such lasers will lend themselves to a variety of applications in integrated photonic and optoelectronic circuitry as well as fiber-optic communication.

    References

    [1] K. J. Vahala. Optical microcavities. Nature, 424, 839-846(2003).

    [2] T. Erdogan, D. G. Hall. Circularly symmetric distributed feedback semiconductor laser: an analysis. J. Appl. Phys., 68, 1435-1444(1990).

    [3] T. Erdogan, D. G. Hall. Circularly symmetric distributed feedback laser: coupled mode treatment of TE vector fields. IEEE J. Quantum Electron., 28, 612-623(1992).

    [4] C. Wu, T. Makino, S. I. Najafi, R. Maciejko, M. Svilans, J. Glinski, M. Fallahi. Threshold gain and threshold current analysis of circular grating DFB and DBR lasers. IEEE J. Quantum Electron., 29, 2596-2606(1993).

    [5] X. M. Gong, A. K. Chan, H. F. Taylor. Lateral mode discrimination in surface emitting DBR lasers with cylindrical symmetry. IEEE J. Quantum Electron., 30, 1212-1218(1994).

    [6] A. Shaw, B. Roycroft, J. Hegarty, D. Labilloy, H. Benisty, C. Weisbuch, T. F. Krauss, C. J. M. Smith, R. Stanley, R. Houdré, U. Oesterle. Lasing properties of disk microcavity based on a circular Bragg reflector. Appl. Phys. Lett., 75, 3051-3053(1999).

    [7] J. Scheuer. Radial Bragg lasers: optimal design for minimal threshold levels and enhanced mode discrimination. J. Opt. Soc. Am. B, 24, 2178-2184(2007).

    [8] J. Scheuer, A. Yariv. Coupled-waves approach to the design and analysis of Bragg and photonic crystal annular resonators. IEEE J. Quantum Electron., 39, 1555-1562(2003).

    [9] X. K. Sun, A. Yariv. Modal properties and modal control in vertically emitting annular Bragg lasers. Opt. Express, 15, 17323-17333(2007).

    [10] X. K. Sun, A. Yariv. Surface-emitting circular DFB, disk-, and ring-Bragg resonator lasers with chirped gratings. II: nonuniform pumping and far-field patterns. Opt. Express, 17, 1-6(2009).

    [11] C. M. Bender, S. Boettcher. Real spectra in non-Hermitian Hamiltonians having PT symmetry. Phys. Rev. Lett., 80, 5243-5246(1998).

    [12] A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, D. N. Christodoulides. Observation of PT-symmetry breaking in complex optical potentials. Phys. Rev. Lett., 103, 093902(2009).

    [13] K. G. Makris, R. El-Ganainy, D. N. Christodoulides, Z. H. Musslimani. Beam dynamics in PT symmetric optical lattices. Phys. Rev. Lett., 100, 103904(2008).

    [14] C. E. Ruter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, D. Kip. Observation of parity-time symmetry in optics. Nat. Phys., 6, 192-195(2010).

    [15] L. Feng, Z. J. Wong, R.-M. Ma, Y. Wang, X. Zhang. Single-mode laser by parity-time symmetry breaking. Science, 346, 972-975(2014).

    [16] J. H. Gu, X. Xi, J. W. Ma, Z. J. Yu, X. K. Sun. Parity-time-symmetric circular Bragg lasers: a proposal and analysis. Sci. Rep., 6, 37688(2016).

    [17] Y. Zhu, Y. Zhao, J. Fan, L. Zhu. Modal gain analysis of parity-time-symmetric distributed feedback lasers. IEEE J. Sel. Top. Quantum Electron., 22, 5-11(2016).

    [18] H. Hodaei, M.-A. Miri, M. Heinrich, D. N. Christodoulides, M. Khajavikhan. Parity-time-symmetric microring lasers. Science, 346, 975-978(2014).

    [19] M. Kulishov, B. Kress, R. Slavík. Resonant cavities based on Parity-Time-symmetric diffractive gratings. Opt. Express, 21, 9473-9483(2013).

    [20] H. Hodaei, A. U. Hassan, W. E. Hayenga, M. A. Miri, D. N. Christodoulides, M. Khajavikhan. Dark-state lasers: mode management using exceptional points. Opt. Lett., 41, 3049-3052(2016).

    [21] K.-H. Kim, M.-S. Hwang, H.-R. Kim, J.-H. Choi, Y.-S. No, H.-G. Park. Direct observation of exceptional points in coupled photonic-crystal lasers with asymmetric optical gains. Nat. Commun., 7, 13893(2016).

    [22] M. Brandstetter, M. Liertzer, C. Deutsch, P. Klang, J. Schöberl, H. E. Türeci, G. Strasser, K. Unterrainer, S. Rotter. Reversing the pump dependence of a laser at an exceptional point. Nat. Commun., 5, 4034(2014).

    [23] M. Liertzer, L. Ge, A. Cerjan, A. D. Stone, H. E. Türeci, S. Rotter. Pump-induced exceptional points in lasers. Phys. Rev. Lett., 108, 173901(2012).

    [24] P. Miao, Z. Zhang, J. Sun, W. Walasik, S. Longhi, N. M. Litchinitser, L. Feng. Orbital angular momentum microlaser. Science, 353, 464-467(2016).

    [25] M. Kulishov, J. M. Laniel, N. Bélanger, J. Azaña, D. V. Plant. Nonreciprocal waveguide Bragg gratings. Opt. Express, 13, 3068-3078(2005).

    [26] Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, D. N. Christodoulides. Unidirectional invisibility induced by PT-symmetric periodic structures. Phys. Rev. Lett., 106, 213901(2011).

    [27] M. Fujita, T. Baba. Proposal and finite-difference time-domain simulation of whispering gallery mode microgear cavity. IEEE J. Quantum Electron., 37, 1253-1258(2001).

    [28] .

    Ziyao Feng, Jingwen Ma, Zejie Yu, Xiankai Sun. Circular Bragg lasers with radial PT symmetry: Design and analysis with a coupled-mode approach[J]. Photonics Research, 2018, 6(5): A38
    Download Citation