- Photonics Research
- Vol. 9, Issue 9, 1758 (2021)
Abstract
1. INTRODUCTION
Optical analog computing has attracted intensive attention recently because it holds the promise to tackle the issues of low speed, high power consumption, and complexity in traditional digital signal processing with electric solutions, especially for the scenarios that require real-time and high throughput [1–3]. As an efficient platform for realizing optical analog computing, artificially engineered photonic structures have been proposed and demonstrated to perform spatial differentiation, integration, convolution, even equation solving, and so on [4–9]. Of particular importance, differentiation is one of the most fundamental mathematical operations, and its realization has gained much research interest. Various principles have been exploited, including metasurface [10–16], Brewster effect [17,18], and surface plasmon polariton [19–22], among others [23–28]. By employing the optical analog differentiation on spatial functions, edge detection and further processing for images can be conducted in an efficient manner [29–33]. Therefore, optical analog computing can facilitate the development of image-processing technology together with electronic platforms.
Among differentiation operators, the Laplacian is a basic and important one that performs the isotropic second-order differentiation. It plays a necessary role in image processing for applications from medical imaging to object detection. It is also an essential operator in various fields of physical science and engineering, for example, wave phenomena of optics and vibration, diffusion process of carriers in semiconductors [34]. The realization of optical analog Laplace operator may offer a powerful tool for solving complex and large-scale mathematical problems. Nevertheless, few works have been reported because it is a challenge to realize the same second-order differentiation for every direction even in a 2D plane. Recently, a multilayer phase-shifted Bragg grating has been proposed [35]. The multiple reflections at different interfaces facilitate the optical transfer function (OTF) as a quadratic function of an in-plane wave vector, thus normally reflecting an incident Gaussian beam into a Laguerre–Gaussian beam. In order to minimize the bulky size of the devices, a photonic crystal slab system has been proposed [36]. The realization of Laplace operation relies on the excitation of a guided resonance by carefully tuning the photonic crystal size parameters and the distance between the photonic crystal slab and the uniform slab. A nearly perfect Laplacian and excellent edge detection of an image have been demonstrated, but the separated-slab system presents challenges for practical fabrication and applications. Meanwhile, it works only for a tiny incident angle and leads to a small numerical aperture (NA) of about 0.01, introducing rigorous requirements on the light path and hampering the improvement of spatial resolution to the wavelength scale.
In this paper, we introduce a single-layer dielectric metasurface that can perform Laplace operation nearly perfectly, the simple design of which can also release the restrictions of fabrications and applications. Most importantly, the maximum incident angle has been widely expanded, and the NA can approach 0.14, benefiting the practical setup and leading to a spatial resolution about 4 times the working wavelength. The proposed Laplace metasurface relies on the excitation of bound state in the continuum (BIC), which has been demonstrated with exotic optical properties [37–42]. The highly symmetric mode profile enables a nearly isotropic OTF of Laplace operation. We have demonstrated that such an OTF can not only produce the second-order differentiation results for a 1D spatial Gaussian function, but also output the correct results for a 2D Bessel function. In addition, we have also shown that the Laplace metasurface can be employed to realize edge detection of objects of concern in an image. The proposed Laplace metasurface can be tuned to work at different wavelengths in a transmission mode, thus benefiting the applications of optical computing, medical diagnostics, machine vision, and so on.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
2. PRINCIPLE AND DESIGN
The optical analog computing of Laplace operation for a 2D spatial function can be expressed in the Cartesian coordinate as
In order to satisfy the requirements of the OTF and the transmission properties, we have designed a novel Laplace metasurface, which is shown in Fig. 1(a). When passing through the Laplace metasurface, the input light field will be automatically processed, and the results of Laplace operation can be obtained just by recording the output light field in a simple and instantaneous way. The Laplace metasurface consists of a square lattice of a modified Si brick on a glass substrate, which is shown in Fig. 1(b). The high symmetry of square lattice and the defects at the four faces of Si brick ensure an isotropic OTF for in-plane directions. We have calculated the transmittance spectra of the Laplace metasurface at different incident angles and shown them in Fig. 1(c). It can be found that a resonance exhibits at the wavelength of approximately 740 nm when the light field is incident obliquely. The transmittance increases as the incident angle is increased. As a result, it provides the necessary condition to realize the Laplace operation.
Figure 1.(a) Illustration of a dielectric metasurface transforming an input 2D spatial function to another function as a Laplace operator; (b) unit cell of the dielectric metasurface. Left, 3D view of the unit cell. It consists of a Si brick (light blue color) with a thickness of
To unambiguously explain the dependence of the transmittance versus the incident angle for the Laplace metasurface, we have calculated the dispersion bands around 740 nm of the proposed dielectric metasurface, which are shown in Fig. 2(a). To capture the essential mechanism, we considered a lossless Si metasurface with a refractive index of 3.73 that was fitted from the experimental results [43]. Because the differentiation for the Laplace operator is isotropic, we consider both and directions of the irreducible Brillouin zone. It can be found that the middle band highlighted with red color is nearly the same along both directions at around the point, promising the isotropic property. The corresponding quality factors of the five bands are shown in Fig. 2(b). For the quasi-BIC band highlighted with red color, the quality factor approaches infinity at the point, indicating a BIC that cannot be accessed by the plane waves outside. Away from the point, the BIC becomes a quasi-BIC, and it can partially couple with the incident plane wave, as the corresponding quality factor is finite. Meanwhile, other bands feature lower quality factors, and they can also couple with the incident wave. The BIC at the point is protected by symmetry. This can be found by examining the mode patterns, which are shown in Fig. 2(c). For the normal component of the electric field of the mode shown at the top panel, it is even under the operation of 180° rotation around the normal direction () of the metasurface, while it is odd for the plane waves. In fact, the symmetry of the mode profile is necessary for achieving isotropic OTF. Furthermore, it can be found that the normal component dominates the electric field, as shown at the bottom panel of Fig. 2(c). Therefore, the quasi-BIC can be accessed by the -polarization plane waves with an incident angle. In addition, we have calculated the transmittance spectra as functions of incident angle along both and directions, which are shown in Figs. 2(d) and 2(e). It can be found that at the target wavelength of 740 nm, the quasi-BIC is excited gradually, with an increased transmittance as the incident angle increases. The resonating wavelength of the quasi-BIC shifts little, but the spectral bandwidth is extended, showing agreement with the calculated dispersion band structure and factors presented in Figs. 2(a) and 2(b).
Figure 2.(a) Dispersion bands around 740 nm of the proposed dielectric metasurface Laplace operator in both
3. RESULTS AND DISCUSSION
It is important to examine whether the transmittance is a quadratic function of the in-plane wave vector. Without loss of generality, we show the transmittance amplitude versus the incidence angle in Fig. 3(a) along the typical directions of and . One can notice that these two curves nearly overlap each other, intimating the possibility of achieving an isotropic OTF. To show the quadratic behavior of the OTF, we fitted the OTF in the direction, and a quadratic function with can be obtained, which is also plotted in Fig. 3(a) and considered as the ideal case. It is obvious that the transmittance curves for the Laplace metasurface are very close to the ideal case with a parabolic line shape within the incidence angle up to 8°, which is improved by nearly 1 order of magnitude compared to previous works and leads to an NA of about 0.14. Such a large NA can help to maintain more information on the incident light field and thus achieve a high spatial resolution approaching the wavelength scale. In addition, we have shown the corresponding phases in Fig. 3(b), and there is only a little difference between the results of Laplace metasurface and the ideal one. To further demonstrate the isotropic OTF, we have studied the transmittance of the Laplace metasurface for arbitrary in-plane wave vectors and the 2D OTF for an incident wave, which is shown in Fig. 3(c). The equal-transmittance contours are with a nearly circular shape, indicating good isotropic property. Moreover, the transmittance phase plotted in Fig. 3(d) shows that the Laplace metasurface can provide nearly the same phase for any in-plane wave vector, making the OTF greatly satisfy the definition of the Laplace operation. We have also studied the OTFs for the wave and polarization conversion, which are significantly smaller than that of the wave (see details in Appendix A). Because of the small values of the transmission for wave, it is hard to realize the Laplace operation for linearly polarized waves. Nevertheless, the OTF for unpolarized waves is similar to the case of the wave as an average result of both polarizations. In this case, the transmission drops a little compared to the case of -polarization. To increase the overall transmission efficiency, ensuring high symmetry of the structure might be a solution. Therefore, the OTF of the designed metasurface makes it possible to act as a practical Laplace operator for unpolarized waves. In addition, the isotropic OTF can be maintained even when the incidence angle is extended to 20°, at which the transmittance amplitude can reach about 0.76 and the phase is still close to (see details in Appendix B). As the excitation of quasi-BIC is general, similar OTFs can be obtained for various wavelengths. For example, we have designed another Laplace metasurface working at the wavelength of 1550 nm (see details in Appendix C). Therefore, the flexibility of the Laplace metasurface can be applied in various scenarios.
Figure 3.(a) Transmittance amplitudes of the Laplace metasurface as functions of the incidence angle along both
In order to verify whether the designed metasurface can output correct results of Laplace operation, we first investigate the second-order differentiation on a 1D spatial function, which is shown in Fig. 4(a). This spatial function consists of three Gaussian envelopes and can be regarded as the electric field profile of an input light beam. To calculate the output function, we first calculated the Fourier spectrum of the input function. We then applied the OTF to the Fourier spectrum. Finally, we conducted the inverse Fourier transform, and the output function can be expressed as . With this processing procedure, we have studied both cases of the ideal Laplacian OTF and the metasurface, which are shown in Figs. 4(b) and 4(c), respectively. It can be found that the output results from the metasurface faithfully show the spatial profile of the ideal one, indicating that the metasurface can correctly conduct the second-order differentiation. The tiny difference between the results of the ideal case and that of the metasurface can be attributed to the small deviation of the OTF of the metasurface from the ideal one. Furthermore, both OTFs in the and directions have been considered, and the results overlap each other very well, which can be clearly observed in Fig. 4(c). In addition, we have set the spatial resolution of for the input function, which is also nearly the maximum spatial resolution of the proposed metasurface [16]. Consequently, the metasurface can work at a resolution of wavelength-scale well, which is improved nearly 1 order of magnitude compared to previous work [36]. Therefore, the designed metasurface can indeed realize the second-order differentiation for a 1D spatial function.
Figure 4.(a) 1D spatial function with three Gaussian envelopes as the input light field; (b) output results after the operation of ideal second-order differentiation; (c) output results from the Laplace metasurface. Both OTFs along the
The Laplace operation not only includes the second-order differentiation but also requires that the differentiation is isotropic. In order to validate the Laplace operation of the designed metasurface, we have further investigated a 2D spatial function , where is a Bessel function of the first kind. Such a function is a possible solution of the Laplace equation in the polar coordinate. Furthermore, the results of the Laplace operation on this function can be obtained analytically (see details in Appendix D). Therefore, it is an appropriate example for validating the performance of the Laplace metasurface. Here we assumed that this spatial function was loaded on the electric field of an unpolarized light; the corresponding light-intensity profile is shown in Fig. 5(a). In the case of ideal Laplace operation, the output results can be derived analytically; its light-intensity profile is shown in Fig. 5(b). It should be noted that although the intensity distributions of both light fields are very similar, the phases between them are totally different (see details in Appendix E). For the Laplace metasurface, we have applied the similar processing procedure mentioned above to obtain the corresponding results, which are shown in Fig. 5(c). It is obvious that the result from metasurface is consistent with that of the ideal case, showing the effectiveness of Laplace operation with a metasurface. To further explain the accurate operation of the Laplace metasurface, we show the light-field intensity along the radial direction with both and in Figs. 5(d)–5(f), which correspond to the signals of input, output of the Laplace operating of the ideal case, and metasurface, respectively. The profiles are actually the cross-sectional views of the light-field patterns along the direction. One can find that the results of the ideal operation and Laplace metasurface have, indeed, a similar profile. In addition, we have also shown the results with varied from 0 to and in Figs. 5(g)–5(i). According to the analytical expression of the input function and its Laplacian results, the corresponding output should behave as , which is indeed correctly shown in the figures. Therefore, the proposed metasurface has demonstrated a valid Laplace operation.
Figure 5.(a) Input 2D spatial function for Laplace operation; (b) analytical solution of the input function; (c) output from the Laplace metasurface; (d)–(f) radial profiles of the corresponding functions along
One of the applications of Laplace operation is the edge detection of targets of concern in an image. Here, we also demonstrated that the proposed Laplace metasurface can be employed to recognize a traffic sign, which is critical for automatic driving. A general traffic sign was selected and is shown in Fig. 6(a). To extract the essential information, we first converted the false-color image to a gray-scale one, shown in Fig. 6(b). Such a gray-scale image can be considered as an intensity-modulated unpolarized light field that is impinged on the Laplace metasurface. From the results after ideal Laplace operation that are shown in Fig. 6(c), we can notice that the edges of both outer circles and the arrow can be clearly identified. Therefore, the necessary information can be safely retained, while the redundant one has been discarded, accelerating the information processing in the first step. We have also shown the output results from the Laplace metasurface in Fig. 6(d), where it can be found that similar edge detection can also be successfully realized. The good agreement between the results implies again that the designed metasurface can perform a Laplace operation with high quality. It should be noted that a monochromatic light beam is necessary for achieving excellent performance due to the narrow bandwidth of the BIC. In addition, it is also important to ensure the symmetry and the fabrication precision of the structure. Otherwise, the Laplace operation and edge detection will degrade. The edge detection for images can also be achieved at the wavelength of 1550 nm (see details in Appendix F).
Figure 6.(a) False-color image of a traffic sign; (b) corresponding gray-scale image as the input; (c) and (d) output image from the ideal Laplace operation and the Laplace metasurface corresponding to (b), respectively.
4. CONCLUSION
In summary, we have proposed and numerically demonstrated a Laplace metasurface that can perform the isotropic second-order spatial differentiation. The angle-resolved excitation of a BIC of the Laplace metasurface provides an OTF that is close to that for ideal Laplace operation at the wavelength of 740 nm. The maximum incident angle has been widely extended to about 8° compared to previous works, improving the spatial resolution. The mode pattern with high symmetry in the Laplace metasurface ensures the isotropic property of the Laplace operation. The Laplace operations of 1D Gaussian function and 2D spatial functions with the proposed Laplace metasurface have been successfully demonstrated, which agree well with that of the ideal Laplace operation. Furthermore, we have demonstrated that the Laplace metasurface can also be employed to realize the edge detection with high quality for images, which would find important applications not only in the fields of optical analog computing, but also in the fields of object recognition.
Acknowledgment
Acknowledgment. T. Feng acknowledges financial support of the Double Hundred Talents Program of Jinan University (DHT2020096).
APPENDIX A: OTFs FOR THE S WAVE AND POLARIZATION CONVERSION
We have calculated the OTFs for the incident wave with polarization and considered the case of polarization conversion, which are presented in Fig.
Figure 7.(a) OTF of the metasurface under the
APPENDIX B: ISOTROPY OF THE OTF
In context, we have shown that the Laplace operation could be achieved when the incidence angle is within 8°. However, even when the incidence angle approaches 20°, the OTF is still nearly isotropic and a transmittance of 0.76 is obtained at the expense of the deviation from the second-order differentiation. We have shown the OTF results at the wavelength of 740 nm in Fig.
Figure 8.(a) 2D transmittance amplitudes as functions of incidence angle under the
APPENDIX C: LAPLACE METASURFACES AT THE WAVELENGTH OF 1550??nm
The Laplace metasurface can be migrated to work at other wavelengths by adjusting the sizes of the structure. To demonstrate the possibility, another Laplace metasurface working at 1550 nm was designed in a similar way. The shape of the structure is the same as that in Fig.
Figure 9.(a) 2D transmittance of a Laplace metasurface at the wavelength of 1550 nm under the
APPENDIX D: ANALYTICAL RESULTS OF THE LAPLACE OPERATION ON THE 2D BESSEL FUNCTION
In context, we have examined the proposed Laplace metasurface with a 2D spatial function consisting of a Bessel function of the first kind and a cosine function as
APPENDIX E: RESULTS OF THE INPUT AND OUTPUT ELECTRIC FIELD DISTRIBUTIONS
We have examined the Laplace metasurface with an analytical function and shown the light-field intensity profiles as the output results. In fact, the light-field intensity cannot present the phase information, leading to a highly similar pattern as that of the input one. Nevertheless, the phases between the input and output functions are totally different. In order to show the difference, we have plotted the real part of the electric fields in Fig.
Figure 10.(a) Real part of the electric field of the input light field; (b) analytical solution of the Laplace operation on the input; (c) output from the Laplace metasurface; (d)–(f) radial profiles of the corresponding functions along the
APPENDIX F: EDGE DETECTION OF IMAGES AT THE WAVELENGTH OF 1550??nm
As the principle of Laplace metasurface is general, we demonstrated the edge detection of an image at the wavelength of 1550 nm with the results presented in Fig.
Figure 11.(a) Input image consisting of a QR code; (b) output image of the ideal Laplace operation; (c) output from the Laplace meatsurface. All images are the light-intensity profile; the pixel sizes are set as
References
[1] F. Zangeneh-Nejad, D. L. Sounas, A. Alù, R. Fleury. Analogue computing with metamaterials. Nat. Rev. Mater., 6, 207-225(2021).
[2] S. Abdollahramezani, O. Hemmatyar, A. Adibi. Meta-optics for spatial optical analog computing. Nanophotonics, 9, 4075-4095(2020).
[3] D. R. Solli, B. Jalali. Analog optical computing. Nat. Photonics, 9, 704-706(2015).
[4] A. Silva, F. Monticone, G. Castaldi, V. Galdi, A. Alù, N. Engheta. Performing mathematical operations with metamaterials. Science, 343, 160-163(2014).
[5] W. Zhang, C. Qu, X. Zhang. Solving constant-coefficient differential equations with dielectric metamaterials. J. Opt., 18, 075102(2016).
[6] N. M. Estakhri, B. Edwards, N. Engheta. Inverse-designed metastructures that solve equations. Science, 363, 1333-1338(2019).
[7] M. Camacho, B. Edwards, N. Engheta. A single inverse-designed photonic structure that performs parallel computing. Nat. Commun., 12, 1466(2021).
[8] X. Ding, Z. Wang, G. Hu, J. Liu, K. Zhang, H. Li, B. Ratni, S. N. Burokur, Q. Wu, J. Tan, C.-W. Qiu. Metasurface holographic image projection based on mathematical properties of Fourier transform. PhotoniX, 1, 16(2020).
[9] S. AbdollahRamezani, K. Arik, A. Khavasi, Z. Kavehvash. Analog computing using graphene-based metalines. Opt. Lett., 40, 5239-5242(2015).
[10] A. Pors, M. G. Nielsen, S. I. Bozhevolnyi. Analog computing using reflective plasmonic metasurfaces. Nano Lett., 15, 791-797(2015).
[11] A. Chizari, S. Abdollahramezani, M. V. Jamali, J. A. Salehi. Analog optical computing based on a dielectric meta-reflect array. Opt. Lett., 41, 3451-3454(2016).
[12] S. Abdollahramezani, A. Chizari, A. E. Dorche, M. V. Jamali, J. A. Salehi. Dielectric metasurfaces solve differential and integro-differential equations. Opt. Lett., 42, 1197-1200(2017).
[13] H. Kwon, D. Sounas, A. Cordaro, A. Polman, A. Alù. Nonlocal metasurfaces for optical signal processing. Phys. Rev. Lett., 121, 173004(2018).
[14] J. Zhou, H. Qian, C.-F. Chen, J. Zhao, G. Li, Q. Wu, H. Luo, S. Wen, Z. Liu. Optical edge detection based on high-efficiency dielectric metasurface. Proc. Natl. Acad. Sci. USA, 116, 11137-11140(2019).
[15] L. Wan, D. Pan, S. Yang, W. Zhang, A. A. Potapov, X. Wu, W. Liu, T. Feng, Z. Li. Optical analog computing of spatial differentiation and edge detection with dielectric metasurfaces. Opt. Lett., 45, 2070-2073(2020).
[16] H. Kwon, A. Cordaro, D. Sounas, A. Polman, A. Alù. Dual-polarization analog 2D image processing with nonlocal metasurfaces. ACS Photon., 7, 1799-1805(2020).
[17] A. Youssefi, F. Zangeneh-Nejad, S. Abdollahramezani, A. Khavasi. Analog computing by Brewster effect. Opt. Lett., 41, 3467-3470(2016).
[18] D. Xu, S. He, J. Zhou, S. Chen, S. Wen, H. Luo. Optical analog computing of two-dimensional spatial differentiation based on the Brewster effect. Opt. Lett., 45, 6867-6870(2020).
[19] T. Zhu, Y. Zhou, Y. Lou, H. Ye, M. Qiu, Z. Ruan, S. Fan. Plasmonic computing of spatial differentiation. Nat. Commun., 8, 15391(2017).
[20] W. Yang, X. Yu, J. Zhang, X. Deng. Plasmonic transmitted optical differentiator based on the subwavelength gold gratings. Opt. Lett., 45, 2295-2298(2020).
[21] W. Zhang, X. Zhang. Backscattering-immune computing of spatial differentiation by nonreciprocal plasmonics. Phys. Rev. Appl., 11, 054033(2019).
[22] Y. Lou, Y. Fang, Z. Ruan. Optical computation of divergence operation for vector fields. Phys. Rev. Appl., 14, 034013(2020).
[23] D. A. Bykov, L. L. Doskolovich, A. A. Morozov, V. V. Podlipnov, E. A. Bezus, P. Verma, V. A. Soifer. First-order optical spatial differentiator based on a guided-mode resonant grating. Opt. Express, 26, 10997-11006(2018).
[24] S. S. M. Khaleghi, P. Karimi, A. Khavasi. On-chip second-order spatial derivative of an optical beam by a periodic ridge. Opt. Express, 28, 26481-26491(2020).
[25] Y. Zhou, W. Wu, R. Chen, W. Chen, R. Chen, Y. Ma. Analog optical spatial differentiators based on dielectric metasurfaces. Adv. Opt. Mater., 8, 1901523(2019).
[26] T. Zhu, Y. Lou, Y. Zhou, J. Zhang, J. Huang, Y. Li, H. Luo, S. Wen, S. Zhu, Q. Gong, M. Qiu, Z. Ruan. Generalized spatial differentiation from spin Hall effect of light. Phys. Rev. Appl., 11, 034043(2019).
[27] D. Xu, S. He, J. Zhou, S. Chen, S. Wen, H. Luo. Goos-Hänchen effect enabled optical differential operation and image edge detection. Appl. Phys. Lett., 116, 211103(2020).
[28] T. Zhu, C. Guo, J. Huang, H. Wang, M. Orenstein, Z. Ruan, S. Fan. Topological optical differentiator. Nat. Commun., 12, 680(2021).
[29] Y. Zhou, H. Zheng, I. I. Kravchenko, J. Valentine. Flat optics for image differentiation. Nat. Photonics, 14, 316-323(2020).
[30] A. Cordaro, H. Kwon, D. Sounas, A. F. Koenderink, A. A. Alù, A. Polman. High-index dielectric metasurfaces performing mathematical operations. Nano Lett., 19, 8418-8423(2019).
[31] H. Wang, C. Guo, Z. Zhao, S. Fan. Compact incoherent image differentiation with nanophotonic structures. ACS Photon., 7, 338-343(2020).
[32] J. Zhou, S. Liu, H. Qian, Y. Li, H. Luo, S. Wen, Z. Zhou, G. Guo, B. Shi, Z. Liu. Metasurface enabled quantum edge detection. Sci. Adv., 6, eabc4385(2020).
[33] T. Zhu, J. Huang, Z. Ruan. Optical phase mining by adjustable spatial differentiator. Adv. Photon., 2, 016001(2020).
[34] M. L. Boas. Mathematical Methods in the Physical Sciences(2005).
[35] D. A. Bykov, L. L. Doskolovich, E. A. Bezus, V. A. Soifer. Optical computation of the Laplace operator using phase-shifted Bragg grating. Opt. Express, 22, 25084-25092(2014).
[36] C. Guo, M. Xiao, M. Minkov, Y. Shi, S. Fan. Photonic crystal slab Laplace operator for image differentiation. Optica, 5, 251-256(2018).
[37] C. W. Hsu, B. Zhen, A. D. Stone, J. D. Joannopoulos, M. Soljačić. Bound states in the continuum. Nat. Rev. Mater., 1, 16048(2016).
[38] C. Huang, C. Zhang, S. Xiao, Y. Wang, Y. Fan, Y. Liu, N. Zhang, G. Qu, H. Ji, J. Han, L. Ge, Y. Kivshar, Q. Song. Ultrafast control of vortex microlasers. Science, 367, 1018-1021(2020).
[39] X. Yin, J. Jin, M. Soljačić, C. Peng, B. Zhen. Observation of topologically enabled unidirectional guided resonances. Nature, 580, 467-471(2020).
[40] B. Wang, W. Liu, M. Zhao, J. Wang, Y. Zhang, A. Chen, F. Guan, X. Liu, L. Shi, J. Zi. Generating optical vortex beams by momentum-space polarization vortices centred at bound states in the continuum. Nat. Photonics, 14, 623-628(2020).
[41] Z. Liu, Y. Xu, Y. Lin, J. Xiang, T. Feng, Q. Cao, J. Li, S. Lan, J. Liu. High-
[42] K. Koshelev, S. Lepeshov, M. Liu, A. Bogdanov, Y. Kivshar. Asymmetric metasurfaces with high-
[43] E. Palik, G. Ghosh. Handbook of Optical Constants of Solids(1985).
Set citation alerts for the article
Please enter your email address