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

Design and optimization of a passive PT-symmetric grating with asymmetric reflection and diffraction

Open Access Open Access

Abstract

In recent years, notions drawn from non-Hermitian physics and parity-time (PT) symmetry have raised considerable attention in photonics, enabling various novel structures with entirely new and unexpected features. Here we propose, design, and optimize a compact passive PT-symmetric grating to achieve asymmetric reflection and diffraction based on a silicon-on-insulator (SOI) platform. The structure is composed of two sets of interleaved tailored gratings, which are all well-defined on the top of a silicon waveguide. Without additional loss or gain materials, the effective index and the scattering loss of the waveguide mode are modulated by the structure design. To our knowledge, it is the first time that the scattering loss arising from grating elements is regarded as an efficient way to realize PT-symmetric structure. The complicated multi-parameter optimization process of the proposed PT-symmetric grating is completed by using the particle swarm optimization (PSO) algorithm. In the simulation, asymmetric reflection with high contrast ratio is realized. We also find that the waveguide-to-free-space diffraction from one side of the structure is significantly suppressed, leading to asymmetric diffraction. Moreover, we investigate the fabrication tolerance of the proposed PT-symmetric grating. Our work provides a new perspective for exploring and creating complicated on-chip PT-symmetric devices.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Optical losses are considered one of the major challenges in integrated optics, as they affect the performance of many optical devices. Therefore, considerable research efforts have been made to eliminate or mitigate underired optical losses [1]. In recent years, non-Hermitian notions in optics have opened a new avenue for exploiting optical losses, which have led to a number of fascinating optical devices and phenomena. In particular, the concept of PT symmetry stems from the discussion of operators in quantum mechanics, which offers a new approach to exploiting optical losses through tailored structural design [25]. In 2007, D. N. Christodoulides et al. introduced PT symmetry into optical systems and formulated a coupled-mode theory (CMT) for PT-symmetric optical waveguide couplers. Thanks to the equivalence between the Schrödinger equation in quantum mechanics and the wave equation in optics, complex PT-symmetric structures can be realized by spatial modulation of the refractive index with balanced gain/loss [6].

The progress in the field of non-Hermitian physics creates opportunities to study novel photonic integrated devices by using gain and loss as new degrees of freedom [720]. On the experimental front, PT symmetry has been investigated in various photonic integrated platforms including both silicon, organic and InP materials. Numerous optical structures associated with PT symmetry have been demonstrated, such as compact dimer structures. The first passive and the first active PT-symmetric optical systems based on coupled waveguides were experimentally demonstrated, which were used to realize optical analogies of spontaneous PT symmetry breaking in complex index potentials [7,8]. PT symmetry has also been realized with active-passive-coupled microresonators, where tunable mode splitting and optical isolation were observed [11,17]. Further theoretical and experimental studies of PT symmetry reveal many other interesting optical phenomena [2131]. An innovative phenomenon named unidirectional invisibility was experimentally demonstrated at the exceptional point in PT-symmetric grating, where reflection from one direction is almost suppressed totally [3235]. Such a PT-symmetric grating was further employed to achieve unidirectional single-mode orbital angular momentum (OAM) vortex lasing based on an InP-based microring resonator [36]. Despite impressive operation performance, the realization of on-chip PT symmetry remains challenging. It is worth stressing that it is still particularly difficult to integrate and control nanoscale optical gain and loss with high resolution in the fabrication process. This is because the optical gain and loss are limited to certain materials and are not generally compatible with all platforms. For example, the optical gain can be implemented through III–V semiconductors via optical or electrical pumping [37], and optical loss usually arises from material absorption, such as chromium (Cr). Therefore, the heterogeneous integration of different semiconductor materials will significantly increase the complexity of the fabrication process. Therefore, alternative designs of passive PT-symmetric structures with reduced complexity are highly desirable.

In this paper, we propose, design, and optimize a compact passive PT-symmetric grating structure based on a silicon photonic platform. The proposed PT-symmetric structure consists of two sets of interleaved tailored gratings, which are all well-defined on the top of a silicon waveguide. The scattering loss arising from the grating elements is engineered for the design of PT-symmetric structure, which provides a new degree of freedom for the design of non-Hermitian structures. By adjusting the period, the etching depth, and the widths of the gratings, the complex effective index of the waveguide mode could be modulated. In the optimization process, the particle swarm optimization (PSO) method is used for such a multi-parameter PT-symmetric grating structure. In the simulation, both unidirectional reflection and asymmetric waveguide-to-free-space diffraction with high contrast ratio are verified, which indicates the successful theoretical design of PT-symmetric structure. In our design, without the integration of additional gain or loss materials, the fabrication complexity could be significantly simplified. Our work may provide a new perspective for exploring the advantages and applications of chip-scale non-Hermitian photonic devices.

2. Principle, design and optimization

Before the notion of PT symmetry was introduced into optics, the complex grating structure has been studied. In 2005, Mykola Kulishov et al. analyzed such kind of grating systematically, and revealed its striking characteristic named unidirectional reflectionless. A complex grating with the modulated dielectric permittivity might be used to break the space-time reversibility between the mode interaction [35]. This leads to an asymmetrical behavior in this mode coupling process. Specifically, the complex grating will induce and amplify coupling from the forward-propagating mode into the backward-propagating mode when the light is launched from one end of the grating but the same is not true when the light is launched from the opposite end. Figure 1(a) shows the detail of an ideal PT-symmetric grating, which has complex modulation of the refractive index in the propagation. In Fig. 1(a), the process of unidirectional reflectionless is illustrated, where the forward reflection is significantly larger than the backward reflection. In the bottom of Fig. 1(a), the red line shows the real part nR of complex refractive index modulation, and the blue line represents the imaginary part nI (gain/loss). The ideal dielectric permittivity modulation can be expressed as follows:

$$\Delta \varepsilon = \cos (\frac{{2\pi }}{\Lambda }x) - i\delta \sin (\frac{{2\pi }}{\Lambda }x),$$
where PT-symmetric grating is actually a periodic longitudinal perturbation of the complex exponential form with simultaneous both index and gain/loss modulation. However, as shown in Fig. 1(a), nanoscale accurate complex index modulation with both gain and loss is quite challenging, which has not yet been experimentally demonstrated.

 figure: Fig. 1.

Fig. 1. (a) Diagram illustrating the process of unidirectional reflectionless and schematic showing ideal PT-symmetric grating with complex refractive index modulation. (b) A quasi-PT-symmetric grating structure with asymmetric reflection.

Download Full Size | PDF

To realize an equivalent complex index modulation, Feng et al. designed a passive grating structure by using additional Si and loss materials germanium (Ge)/chrome (Cr) bilayer combo structures on top of the waveguide to mimic the microscopic parity-time modulation on a macroscopic scale. As shown in Fig. 1(b), two sets of interleaved passive gratings are used to form a quasi-PT-symmetric grating. The dielectric permittivity modulation can be written as:

$$\Delta \varepsilon = C\left\{ {\begin{array}{{cc}} {\cos (qx)}&{4n\pi /q + 5\pi /2q \le x \le 4n\pi /q + 7\pi /2q}\\ { - i\delta \sin (qx)}&{4n\pi /q + \pi /q \le x \le 4n\pi /q + 2\pi /q}\\ 0&{\textrm{other else}} \end{array}} \right..$$
where q = 4π/Λ. When the constant C is 1, the first and second terms in Eq. (2) have a fixed phase difference, which represents the real index modulation and the loss modulation, respectively. The Fourier expansion coefficients can be further obtained by:
$$\begin{array}{l} \textrm{ }{A_0} = \frac{1}{\Lambda }\int_0^{\frac{{4\pi }}{q}} {\Delta \varepsilon dx} ,\textrm{ }\\ {A_m} = \frac{1}{\Lambda }\int_0^{\frac{{4\pi }}{q}} {\Delta \varepsilon exp ( - i\frac{{2\pi m}}{\Lambda }x)dx} ,\textrm{ }\\ {A_{ - m}} = \frac{1}{\Lambda }\int_0^{\frac{{4\pi }}{q}} {\Delta \varepsilon exp (i\frac{{2\pi m}}{\Lambda }x)dx} \textrm{. } \end{array}$$

Substituting Eq. (2) into Eq. (3), one can obtain the values as follows:

$$\begin{array}{l} {A_0} = i\frac{\delta }{{2\pi }},\\ {A_1} = \frac{1}{{3\pi }}(\delta - i(\delta + \sqrt 2 )),\\ {A_{ - 1}} = \frac{1}{{3\pi }}( - \delta + i(\sqrt 2 - \delta )),\\ {A_2} = \frac{{1 - \delta }}{8},\\ {A_{ - 2}} = \frac{{1 + \delta }}{8}, \end{array}$$
where the Fourier expansion coefficients A2 and A-2 determine the mode coupling coefficients between forward and backward fundamental modes. Obviously, when δ is set to be ±1, one of the cross-coupling coefficients becomes zero, and asymmetric reflection introduced by asymmetric coupling will be realized [33]. Meanwhile, the Fourier expansion coefficients A1 and A-1 are also different, leading to asymmetric vertical diffraction.

Here, we consider that the constant C in Eq. (2) is set to be exp(-iπ/4), thus, the dielectric permittivity modulation is expressed by:

$$\Delta \varepsilon = \left\{ {\begin{array}{{cc}} {(1 - i)\cos (qx)}&{4n\pi /q + 5\pi /2q \le x \le 4n\pi /q + 7\pi /2q}\\ {( - 1 - i)\delta \sin (qx)}&{4n\pi /q + \pi /q \le x \le 4n\pi /q + 2\pi /q}\\ 0&{\textrm{other else}} \end{array}}, \right.$$
where two sets of gratings are complex modulated with both refractive index and loss, which is different from the index distribution in Eq. (2). The first term in Eq. (5) represents both refractive index decreasing and loss modulation, and the second term indicates the index increasing and loss modulation.

Based on the above analysis, we propose a new guided-mode modulation method by meticulous structural design and algorithm optimization. Without additional loss materials, the scattering loss could be considered as an efficient way to form a PT-symmetric structure. As shown in Fig. 2, a passive PT-symmetric grating structure based on silicon platform is designed, which is composed of two sets of interleaved tailored gratings. One grating is formed with additional silicon on the top of the waveguide,which can lead to the mode effective index increasing and the scattering loss. Another grating with periodic trenching etching can decrease the mode index and also enhance the scattering loss. Two complex gratings are consistent with the modulations in Eq. (5). To simplify the design, rectangular-function modulated optical potentials using additional silicon and trench etching are employed, which could achieve the same unidirectional effect with sinusoidal-function modulation [33].

 figure: Fig. 2.

Fig. 2. (a) Structure of the proposed passive PT-symmetric grating. (b) Schematic cross section of the grating, where h is the height of the silicon waveguide, h1 and h2 are the heights of two sets of grating respectively. T, T1, T2 and T3 represent the period of PT-symmetric grating, the width of the additional silicon, the width of the trenches and the distance of two sets of grating, respectively.

Download Full Size | PDF

As shown in Fig. 2, the studied passive PT-symmetric grating is a 500-nm-wide and 220-nm-thick silicon waveguide with periodic etching. The waveguide supports a fundamental mode (quasi-TE mode) with an effective index of neff =2.37 around the wavelength of 1550 nm. As depicted in Fig. 2(a), two sets of grating are interleaved on the top of the silicon waveguide. The height h of the silicon waveguide is 220 nm. The period T is chosen to be 654 nm for the central wavelength of 1550 nm, corresponding to the second Bragg order. The number of grating units is set to 30. Considering that the mode effective index and the scattering loss could be greatly affected by the grating size, which should be meticulously designed. However, the complex effective index of the waveguide mode is actually difficult to calculate accurately due to microscopic structure perturbation. In order to obtain optimal structure parameters, the proposed PT-symmetric grating is optimized by PSO method. The distance between the additional silicon and the trench etching is T/8 (i.e. T3 = T/8) in Fig. 2(b). The widths (T1, T2) and heights (h1, h2) of the grating elements needs further optimization.

Recently, many iterative optimization algorithms have been used to optimize structural parameters of functional photonic integrated devices, such as PSO, genetic algorithm (GA), and direct binary search (DBS) algorithm [3841]. It is known that PSO has the characteristics of fast convergence and simple operation. Here, we use PSO algorithm to optimize the proposed multi-parameter PT-symmetric structure, which has simpler rules than genetic algorithm (GA) and is more suitable for multi-parameter continuous optimization than direct binary search (DBS) algorithm. The entire optimization procedure is shown in Fig. 3. A flow chart of PSO process is provided in Fig. 3(a). Initializing generation means randomly generating 20 grating structures with T1, T2, h1, h2 to serve as 20 populations of the first generation. The figure-of-merit (FOM) of the device designed by PSO is defined as:

$$FOM = 1 - 0.5(1 - {t_F}) - 0.5{t_B},$$
where tF and tB are the reflectivities when light is injected into the designed structure from forward and backward directions. In the simulation, the maximum values of tF and tB over an operating bandwidth of 60 nm are chosen for the FOM calculation. For an ideal PT-symmetric grating, the FOM is 1. In the initial step of the optimization process, particle’s velocity and position in the 4-dimensional parameter space are randomly selected. In PSO algorithm, x means particle’s position and v means particle’s velocity. We choose h2 as x1, h1 as x2, T2 as x3, and T1 as x4. Then v1 and v2 are randomly selected from -10 to 10 when x1 is randomly selected from 20 nm to 200 nm, and x2 is randomly selected from 240 nm to 420 nm, while v3 and v4 are randomly selected from -4 to 4 when x3 and x4 are randomly selected from 80 nm to 220 nm.

 figure: Fig. 3.

Fig. 3. (a) Flow chart for the optimization process. (b) Figure of merit at each iteration for the device with 20 populations. The dots are the individual values, and the blue line of the inset shows the maximum FOM values for each iteration. The optimized device has a maximum FOM value of 0.8845.

Download Full Size | PDF

We simulate the proposed device by using three-dimensional finite difference time domain (3D-FDTD) method. The simulation mesh size is 44 nm×44 nm×10 nm, and the power cutoff condition is 0.00001% of initial power. In the simulation, the structural parameters are adjusted and optimized by PSO algorithm. When the FOM increases, new populations would be generated and used for FOM calculation. Otherwise, the parameters go back to last state. When the FOM exhibits no great improvement after 10 iterations, the whole optimization will end. Figure 3(b) shows the convergence to a FOM of 0.8845, and the inset shows the details. It can be seen that the FOM growth trend slows down after 5 iterations. The final optimization parameters are obtained when 30 iterations are finished. The width T1 and depth h1 of additional silicon are 130 nm and 330 nm, the width T2 and depth h2 of trench etching are 190 nm and 190 nm.

3. Simulation results

3.1 Asymmetric reflection

In this part, we first simulate two gratings separately. Due to the left-right symmetry of the conventional grating structures, the transmission and reflection spectra are the same on both sides. The transmission and reflection spectra of the grating with trench etching are plotted in Fig. 4(a). One can see that the maximum reflectivity is about 0.4 around the wavelength of 1543 nm. In Fig. 4(b), the simulation results for the grating with addition silicon are presented. The center wavelength moves to 1567 nm, and the maximum reflectivity is about 0.35. It is worth noting that trench etching and additional silicon would decrease and increase the waveguide effective index, respectively, and thus shift the resonance peak.

 figure: Fig. 4.

Fig. 4. (a)(b) Numerically calculated transmission and reflection spectra for conventional gratings with trench etching and additional silicon, respectively. (c) Simulated reflection spectra of the designed PT-symmetric grating in forward and backward directions. (d) Calculated contrast ratio of reflectivities, showing high contrast ratios over the studied wavelength range from 1520 nm to 1580 nm.

Download Full Size | PDF

Different from the uniform gratings with the same reflection spectra on both sides, when two complex gratings are combined to form a passive PT-symmetric grating, the structural symmetry is broken and asymmetric reflection could be realized. Figure 4(c) shows the simulated results of the optimized PT-symmetric grating. One can see that the reflectivities are significantly distinguished in the forward and backward directions with high extinction ratio in the studied wavelength range from 1520 nm to 1580 nm. The resonance peak of the forward reflection spectrum is around 1550 nm owing to the balanced index modulation introduced by both trench etching and additional silicon. In such a passive system, optical loss is no longer compensated without optical gain. The corresponding S-matrix is not unitary and a similar exceptional point still exists, which causes unidirectional invisibility. It is evident that the designed on-chip PT-symmetric grating successfully mimics the unidirectional effect inherently associated with the exceptional point. The complex refractive index modulation essentially creates the asymmetric coupling coefficients, thus leading to nontrivial asymmetric reflection. The contrast ratio between the forward and backward reflection is defined as C = (RF-RB)/(RF+RB), where RF and RB represent the forward and backward reflectivities, respectively. The calculated contrast ratio is plotted in Fig. 4(d), showing that the reflection in the backward direction is significantly suppressed. At the wavelength of 1551 nm, the forward reflectivity is 0.779, while the backward reflectivity is 0.009, and the contrast ratio is 0.977. Compared with previously reported device [33], the reflectivity and contrast ratio of the optimized PT-symmetric grating have been significantly improved.

The asymmetric reflection from the designed PT-symmetric grating can also be visualized from the light fields propagating inside the waveguide. We simulate the transmitted light fields inside the waveguide grating at 1550 nm when light beams are launched from both sides. As shown in Fig. 5, the forward propagating light beam and its reflection form a strong interference, whereas the reflection is barely seen with backward incidence. These results show that the designed on-chip waveguide grating successfully mimics the unidirectional effect inherently associated with the exceptional point.

 figure: Fig. 5.

Fig. 5. (a)(b) Simulated transmitted light fields at 1550 nm when light beams are launched from both sides. The incidences are set at boundaries of the designed PT-symmetric grating along + x/-x in the forward/backward directions.

Download Full Size | PDF

3.2 Asymmetric diffraction

As shown in Eq. (4), when δ is close to 1, the second-order Fourier expansion coefficient A-2 is much larger than A2, which indicates asymmetric in-plane reflection. Meanwhile, the first-order Fourier expansion coefficient A-1 is smaller than A1, thus leading to asymmetric out-of-plane vertical diffraction. The simulated diffraction efficiencies with forward and backward incidences are plotted in Fig. 6(a), where the vertical diffraction efficiency with forward incidence is much smaller than that of the backward direction around the wavelength of 1550 nm. The contrast ratio between the forward and backward diffractions is calculated and plotted in Fig. 6(b), and the maximum value is 0.75 around 1550 nm. The simulated electric field distributions of out-of-plane emissions at 1550 nm are shown in Figs. 6(c) and 6(d). Since the transmitted guided mode is the fundamental transverse electric (TE) mode, the diffracted light fields are polarized in y-direction. As a comparison, it can be seen that the electric field amplitude of the forward diffraction pattern in Fig. 6(c) is smaller than that of the backward diffraction in Fig. 6(d), which agrees well with the results given in Fig. 6(a). Therefore, both asymmetric in-plane reflection and chip-to-free-space diffraction can be realized based on the proposed PT-symmetric grating.

 figure: Fig. 6.

Fig. 6. (a) Simulated diffraction spectra of the designed PT-symmetric grating with forward and backward incidences. (b) Calculated contrast ratio of the diffraction between forward and backward propagation. (c)(d) Electric field distributions of out-of-plane emissions with forward and backward incidences at 1550 nm.

Download Full Size | PDF

4. Fabrication tolerance

Since it is difficult to control the etching of the microscopic structures during fabrication, we further analyze the effect of the structure parameter deviations on device performance, including the widths (T1, T2) and heights (h1, h2) of the grating elements. In Figs. 7(a) and 7(b), the forward reflection spectrum shifts slowly when the height and width of the additional silicon vary in the range of ±30 nm. The corresponding contrast ratios are presented in Figs. 7(c) and 7(d). One can see that the average values of the contrast ratios are larger than 0.85, which indicates that the contrast ratio is not very sensitive to the parameter deviations of the additional silicon. In Fig. 8(a), when the trench depth varies in the range of ±20 nm, the forward reflection spectrum shifts drastically and the maximum reflectivity also varies greatly. In Fig. 8(b), the forward reflection spectrum shifts slowly when the trench width varies in the range of ±30 nm. The corresponding contrast ratios of reflectivities are presented in Figs. 8(c) and 8(d). When the deviation Δh2 of etching depth is 20 nm, the average value of the contrast ratio shown in Fig. 7(c) is reduced to 0.6. When the deviation ΔT2 of etching width is ±30 nm, the average value of the contrast ratio is still more than 0.9 in Fig. 7(d). Therefore, by comparison, it can be found that the device performance is very sensitive to the trench depth, and the corresponding fabrication tolerance is about ±20 nm.

 figure: Fig. 7.

Fig. 7. (a)(b) Simulated forward reflection spectra of the device by changing the height and width of the additional silicon, respectively. (c)(d) Calculated contrast ratios of reflectivities of the device by changing the height and width of the additional silicon, respectively.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. (a)(b) Simulated forward reflection spectra of the device by changing the height and width of the trenches, respectively. (c)(d) Calculated contrast ratios of reflectivities of the device by changing the height and width of the trenches, respectively.

Download Full Size | PDF

5. Discussion

Compared with previous reports, we have made some improvements in structure design and performance. First, we introduce the scattering loss to PT-symmetric optical system by tailoring microscopic structure design. Thus, the absorption materials could be avoided, and the fabrication process could be simplified. Second, we use the PSO algorithm to accelerate the complicated multi-parameter optimization process and improve device performance, which may open up a door to develop more superior on-chip non-Hermitian photonic integrated devices for wide optical applications. Last, we further reveal that asymmetric diffraction can be also achieved. Thanks to its intrinsic asymmetric properties, the proposed passive PT-symmetric grating could be further developed and employed for functional integrated devices such as grating couplers and unidirectional lasers. Specifically, the proposed PT-symmetric grating could help reduce reflection from grating couplers because the backward reflection is greatly suppressed. Moreover, the proposed device could be elegantly exploited to establish unidirectional single-mode lasing in microring resonators, which indicates a likely perspective for fostering a new generation of active photonic devices [36]. Furthermore, the scattering loss induced by gratings could be also used to explore more complicated non-Hermitian photonic devices, such as anti-PT symmetry and high-order exceptional point [4244].

6. Conclusion

Optical loss is usually undesirable for optical systems and micro-scale optical devices. Fortunately, the concepts from non-Hermitian and PT symmetry provide a new perspective for the design of optical devices by employing various optical losses. We propose, design and optimize a compact passive silicon PT-symmetric grating to achieve asymmetric reflection and diffraction by using PSO algorithm. By combining two sets of interleaved tailored gratings, equivalent PT-symmetric grating structure is designed and optimized. In the simulation, asymmetric reflection and diffraction with high contrast ratio have been realized. It is believed that the proposed passive PT-symmetric structure would be further developed and employed for more complicated functional integrated devices.

Funding

National Key Research and Development Program of China (2018YFE0201800); China Postdoctoral Science Foundation (2021M690390); National Natural Science Foundation of China (62071042, 62105028).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. R. Wolf, I. Breunig, H. Zappe, and K. Buse, “Scattering-loss reduction of ridge waveguides by sidewall polishing,” Opt. Express 26(16), 19815–19820 (2018). [CrossRef]  

2. T. Kottos, “Optical physics: Broken symmetry makes light work,” Nat. Phys. 6(3), 166–167 (2010). [CrossRef]  

3. C. M. Bender and S. Boettcher, “Real spectra in non-Hermitian Hamiltonians having PT symmetry,” Phys. Rev. Lett. 80(24), 5243–5246 (1998). [CrossRef]  

4. C. M. Bender, G. V. Dunne, and P. N. Meisinger, “Complex periodic potentials with real band spectra,” Phys. Lett. A 252(5), 272–276 (1999). [CrossRef]  

5. C. M. Bender, “Making sense of non-Hermitian Hamiltonians,” Rep. Prog. Phys. 70(6), 947–1018 (2007). [CrossRef]  

6. R. El-Ganainy, K. G. Makris, D. N. Christodoulides, and Z. H. Musslimani, “Theory of coupled optical PT-symmetric structures,” Opt. Lett. 32(17), 2632–2634 (2007). [CrossRef]  

7. A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. VolatierRavat, and V. Aimez, “Observation of PT-symmetry breaking in complex optical potentials,” Phys. Rev. Lett. 103(9), 093902 (2009). [CrossRef]  

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

9. C. Hang, D. A. Zezyulin, G. Huang, V. V. Konotop, and B. A. Malomed, “Tunable nonlinear double-core PT-symmetric waveguides,” Opt. Lett. 39(18), 5387–5390 (2014). [CrossRef]  

10. S. Phang, A. Vukovic, H. Susanto, T. M. Benson, and P. Sewell, “Impact of dispersive and saturable gain/loss on bistability of nonlinear parity-time Bragg gratings,” Opt. Lett. 39(9), 2603–2606 (2014). [CrossRef]  

11. L. Chang, X. Jiang, S. Hua, C. Yang, J. Wen, L. Jiang, G. Li, G. Wang, and M. Xiao, “Parity Time Symmetry and Variable Optical Isolation in Active Passive-Coupled Microresonators,” Nat. Photonics 8(7), 524–529 (2014). [CrossRef]  

12. S. Phang, A. Vukovic, T. M. Benson, H. Susanto, and P. Sewell, “A versatile all-optical parity-time signal processing device using a Bragg grating induced using positive and negative Kerr-nonlinearity,” Opt. Quantum Electron. 47(1), 37–47 (2015). [CrossRef]  

13. T. Wasak, P. Szankowski, V. V. Konotop, and M. Trippenbach, “Four-wave mixing in a parity-time (PT)-symmetric coupler,” Opt. Lett. 40(22), 5291–5294 (2015). [CrossRef]  

14. R. El-Ganainy, J. I. Dadap, and R. M. Osgood, “Optical parametric amplification via non-Hermitian phase matching,” Opt. Lett. 40(21), 5086–5089 (2015). [CrossRef]  

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

16. H. Hodaei, M.-A. Miri, M. Heinrich, D. N. Christodoulides, and M. Khajavikhan, “Parity-time-symmetric microring lasers,” Science 346(6212), 975–978 (2014). [CrossRef]  

17. B. Peng, S. K. özdemir, F. Lei, F. Monifi, M. Gianfreda, G. L. Long, S. Fan, F. Nori, C. M. Bender, and L. Yang, “Parity-time-symmetric whispering-gallery microcavities,” Nat. Phys. 10(5), 394–398 (2014). [CrossRef]  

18. F. Nazari, N. Bender, H. Ramezani, M. K. Moravvej-Farshi, D. N. Christodoulides, and T. Kottos, “Optical isolation via PT-symmetric nonlinear Fano resonances,” Opt. Express 22(8), 9574–9584 (2014). [CrossRef]  

19. Y. Yan and N. C. Giebink, “Passive PT symmetry in organic composite films via complex refractive index modulation,” Adv. Opt. Mater. 2(5), 423–427 (2014). [CrossRef]  

20. Y. Jia, Y. Yan, S. V. Kesava, E. D. Gomez, and N. C. Giebink, “Passive parity-time symmetry in organic thin film waveguides,” ACS Photonics 2(2), 319–325 (2015). [CrossRef]  

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

22. S Longhi, “PT-symmetric laser absorber,” Phys. Rev. A: At., Mol., Opt. Phys. 82(3), 031801 (2010). [CrossRef]  

23. H. Benisty, A. Degiron, A. Lupu, A. De Lustrac, S. Ch’enais, S. Forget, M. Besbes, G. Barbillon, A. Bruyant, S. Blaize, and G. Lérondel, “Implementation of PT Symmetric devices using plasmonics: Principle and Applications,” Opt. Express 19(19), 18004–18019 (2011). [CrossRef]  

24. A. Lupu, H. Benisty, and A. Derigon, “Switching using PT symmetry in plasmonic systems: positive role of the losses,” Opt. Express 21(18), 21651–21668 (2013). [CrossRef]  

25. K. G. Makris, R. El-Ganainy, D. N. Christodoulides, and Z. H. Musslimani, “PT-symmetric optical lattices,” Phys. Rev. A 81(6), 063807 (2010). [CrossRef]  

26. A. Regensburger, C. Bersch, M. A. Miri, G. Onishchukov, D. N. Christodoulides, and U. Peschel, “Parity-time Synthetic Photonic Lattices,” Nature 488(7410), 167–171 (2012). [CrossRef]  

27. H. Jing, S. K. Özdemir, X. Y. Lü, J. Zhang, L. Yang, and F. Nori, “PT-symmetric phonon laser,” Phys. Rev. Lett. 113(5), 053604 (2014). [CrossRef]  

28. Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. N. Christodoulides, “Optical solitons in PT periodic potentials,” Phys. Rev. Lett. 100(3), 030402 (2008). [CrossRef]  

29. C. M. Bender, M. Gianfreda, Ş. K. Özdemir, B. Peng, and L. Yang, “Twofold transition in PT-symmetric coupled oscillators,” Phys. Rev. A 88(6), 062111 (2013). [CrossRef]  

30. A. Szameit, Y. V. Kartashov, M. Heinrich, F. Dreisow, T. Pertsch, S. Nolte, A. Tunnermann, F. Lederer, V. A. Vysloukh, and L. Torner, “Observation of two-dimensional defect surface solitons,” Opt. Lett. 34(6), 797–799 (2009). [CrossRef]  

31. A. A. Sukhorukov, Z. Xu, and Y. S. Kivshar, “Nonlinear suppression of time reversals in PT-symmetric optical couplers,” Phys. Rev. A 82(4), 043818 (2010). [CrossRef]  

32. Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, and D. N. Christodoulides, “Unidirectional Invisibility Induced by PT-Symmetric Periodic Structures,” Phys. Rev. Lett. 106(21), 213901 (2011). [CrossRef]  

33. L. Feng, Y.-L. Xu, W. S. Fegadolli, M.-H. Lu, J. E. B. Oliveira, V. R. Almeida, Y.-F. Chen, and A. Scherer, “Experimental Demonstration of a Unidirectional Reflectionless Parity-Time Metamaterial at Optical Frequencies,” Nat. Mater. 12(2), 108–113 (2013). [CrossRef]  

34. S. Zheng and J. Wang, “Design of Asymmetrical Silicon Waveguide Grating by Introducing the Scattering Loss,” Asia Communications and Photonics Conference, Optical Society of America, paper ASu2A-2 (2015).

35. M. Kulishov, J. M. Laniel, N. Belanger, J. Azaña, and D. Plant, “Nonreciprocal waveguide Bragg gratings,” Opt. Express 13(8), 3068–3078 (2005). [CrossRef]  

36. P. Miao, Z. Zhang, J. Sun, W. Walasik, S. Longhi, N. M. Litchinitser, and L. Feng, “Orbital angular momentum microlaser,” Science 353(6298), 464–467 (2016). [CrossRef]  

37. Z. Gao, S. T. Fryslie, B. J. Thompson, P. S. Carney, and K. D. Choquette, “Parity-time symmetry in coherently coupled vertical cavity laser array,” Optica 4(3), 323–329 (2017). [CrossRef]  

38. Y. Zhang, S. Yang, A. E. J. Lim, G. Q. Lo, C. Galland, T. Baehr-Jones, and M. Hochberg, “A compact and low loss Y-junction for submicron silicon waveguide,” Opt. Express 21(1), 1310–1316 (2013). [CrossRef]  

39. W. Chang, L. Lu, X. Ren, D. Li, Z. Pan, M. Cheng, D. Liu, and M. Zhang, “Ultracompact dual-mode waveguide crossing based on subwavelength multimode-interference couplers,” Photonics Res. 6(7), 660–665 (2018). [CrossRef]  

40. J. C. Mak, C. Sideris, J. Jeong, A. Hajimiri, and J. K. Poon, “Binary particle swarm optimized 2×2 power splitters in a standard foundry silicon photonic platform,” Opt. Lett. 41(16), 3868–3871 (2016). [CrossRef]  

41. Z. Yu, H. Cui, and X. Sun, “Genetically optimized on-chip wideband ultracompact reflectors and Fabry–Perot cavities,” Photonics Res. 5(6), B15–B19 (2017). [CrossRef]  

42. H. Hodaei, A. U. Hassan, S. Wittek, H. Garcia-Gracia, R. El-Ganainy, D. N. Christodoulides, and M. Khajavikhan, “Enhanced sensitivity at higher-order exceptional points,” Nature 548(7666), 187–191 (2017). [CrossRef]  

43. H. Fan, J. Chen, Z. Zhao, J. Wen, and Y. P. Huang, “Antiparity-time symmetry in passive nanophotonics,” ACS Photonics 7(11), 3035–3041 (2020). [CrossRef]  

44. Y. Wei, H. Zhou, Y. Chen, Y. Ding, J. Dong, and X. Zhang, “Anti-parity-time symmetry enabled on-chip chiral polarizer,” Photonics Res. 10(1), 76–83 (2022). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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 (8)

Fig. 1.
Fig. 1. (a) Diagram illustrating the process of unidirectional reflectionless and schematic showing ideal PT-symmetric grating with complex refractive index modulation. (b) A quasi-PT-symmetric grating structure with asymmetric reflection.
Fig. 2.
Fig. 2. (a) Structure of the proposed passive PT-symmetric grating. (b) Schematic cross section of the grating, where h is the height of the silicon waveguide, h1 and h2 are the heights of two sets of grating respectively. T, T1, T2 and T3 represent the period of PT-symmetric grating, the width of the additional silicon, the width of the trenches and the distance of two sets of grating, respectively.
Fig. 3.
Fig. 3. (a) Flow chart for the optimization process. (b) Figure of merit at each iteration for the device with 20 populations. The dots are the individual values, and the blue line of the inset shows the maximum FOM values for each iteration. The optimized device has a maximum FOM value of 0.8845.
Fig. 4.
Fig. 4. (a)(b) Numerically calculated transmission and reflection spectra for conventional gratings with trench etching and additional silicon, respectively. (c) Simulated reflection spectra of the designed PT-symmetric grating in forward and backward directions. (d) Calculated contrast ratio of reflectivities, showing high contrast ratios over the studied wavelength range from 1520 nm to 1580 nm.
Fig. 5.
Fig. 5. (a)(b) Simulated transmitted light fields at 1550 nm when light beams are launched from both sides. The incidences are set at boundaries of the designed PT-symmetric grating along + x/-x in the forward/backward directions.
Fig. 6.
Fig. 6. (a) Simulated diffraction spectra of the designed PT-symmetric grating with forward and backward incidences. (b) Calculated contrast ratio of the diffraction between forward and backward propagation. (c)(d) Electric field distributions of out-of-plane emissions with forward and backward incidences at 1550 nm.
Fig. 7.
Fig. 7. (a)(b) Simulated forward reflection spectra of the device by changing the height and width of the additional silicon, respectively. (c)(d) Calculated contrast ratios of reflectivities of the device by changing the height and width of the additional silicon, respectively.
Fig. 8.
Fig. 8. (a)(b) Simulated forward reflection spectra of the device by changing the height and width of the trenches, respectively. (c)(d) Calculated contrast ratios of reflectivities of the device by changing the height and width of the trenches, respectively.

Equations (6)

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

Δ ε = cos ( 2 π Λ x ) i δ sin ( 2 π Λ x ) ,
Δ ε = C { cos ( q x ) 4 n π / q + 5 π / 2 q x 4 n π / q + 7 π / 2 q i δ sin ( q x ) 4 n π / q + π / q x 4 n π / q + 2 π / q 0 other else .
  A 0 = 1 Λ 0 4 π q Δ ε d x ,   A m = 1 Λ 0 4 π q Δ ε e x p ( i 2 π m Λ x ) d x ,   A m = 1 Λ 0 4 π q Δ ε e x p ( i 2 π m Λ x ) d x
A 0 = i δ 2 π , A 1 = 1 3 π ( δ i ( δ + 2 ) ) , A 1 = 1 3 π ( δ + i ( 2 δ ) ) , A 2 = 1 δ 8 , A 2 = 1 + δ 8 ,
Δ ε = { ( 1 i ) cos ( q x ) 4 n π / q + 5 π / 2 q x 4 n π / q + 7 π / 2 q ( 1 i ) δ sin ( q x ) 4 n π / q + π / q x 4 n π / q + 2 π / q 0 other else ,
F O M = 1 0.5 ( 1 t F ) 0.5 t B ,
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.