1. Introduction
Metasurfaces provide a way to manipulate light by precisely arranging sub-wavelength structures and fine-tuning the parameters of these structures on a flat substrate [
1,
2,
3]. By sophisticatedly designing the 2D-arrays, metasurfaces show the capability of controlling the phase, amplitude, and polarization of the light [
4,
5], resulting in a wide range of applications, such as imaging [
6], holography [
7], quantum optics [
8], spectrometry [
9], and structured light [
10]. One important application of interest is the highly reflective mirrors, which often use all-dielectric nanomaterials taking their advantages of the low-loss property and enhanced resonance [
11]. Such devices are also referred to as
perfect reflectors because they can achieve near-unity reflectivity at certain wavelengths by properly selecting the geometries of resonant nano-structures [
12,
13,
14,
15] and precisely arranging their positions [
16]. For example, Slovick
et al. has presented a structure of square-latticed Si cubes on
substrate with over 99.999% reflectivity and less than 0.001% absorptivity at
in simulation [
12]. Based on the principles of this work, cylindrical silicon resonators were experimentally demonstrated to achieve a reflectance of 99.9% at 1450
[
13]. Further, large-scale dielectric reflectors with frustums in hexagonal lattices applying the same material were fabricated and demonstrated to show 99.7% reflectivity at 1530
[
14].
To better guide the design of the highly reflective metasuraces, one can mainly follow the principles of guided-mode resonance (GMR) [
17,
18,
19]. The GMR gives a physical picture of how the incident light diffracts and propagates in a grating waveguide slab. By optimizing the parameters of the grating layer, it is possible to select the diffracted light at certain mode to achieve the constructive interference [
17], which comes from the electromagnetic resonances within the nanostructures of the grating [
20,
21,
22]. In order to have a deeper insight of the physics, the effective medium theory [
12,
13,
14,
23,
24,
25,
26] and the multipole decomposition method [
27,
28,
29] are often used to help the analysis. Following the effective medium theory, one can first calculate the effective permittivity and permeability via the scattering-parameters (S-parameter) [
12,
25], and correspondingly acquire the effective impedance and refractive index, which directly show the properties of the grating layer. On the other hand, electromagnetic resonances in the grating layer result in the scattered fields of electric and magnetic multipole. Through the multipole decomposition, the scattering cross-sections (SCSs) of different resonant modes can be quantitatively calculated, which reveals the contributions of various electromagnetic dipoles, quadrupoles, or higher-order multipoles [
28].
In recent years, metasurface reflectors have also been demonstrated to show a high laser induced damage threshold (LIDT) under high-power lasers, using single-crystal diamonds (SCDs) as the material [
15]. Because such materials have excellent optical and thermal properties: relatively high refractive index (
), wide bandgap (
), and most importantly, the highest thermal conductivity ( 2200
/
/
at room temperature) [
30]. Atikian
et al. fabricated a diamond-based metasurface by angled-etching nanostructures into "golf tee"-shaped columns, achieving over 98% reflectivity at 1070
and damage-free operation under 10
continuous-wave laser with the spot diameter at 750
. [
15]. Another advantage of SCDs is that they have extremely low absorption at a broad range of light spectrum, possessing the potential for shorter wavelengths, e.g. visible or even ultraviolet light [
30]. Therefore, the diamond-based metasurfaces have the potential to be used for high-power alkali lasers [
31], whose typical wavelength is at 795
. Moreover, the etching (vertical and angled) of diamond pillars has been achieved in the study by Atikian et al [
15], demonstrating the potential of more diamond-based metasurface designs in fabrication.
In this paper, we present the design of a diamond-based metasurface, which can achieve near-unity reflectivity at 795 . First, we utilize GMR as design principles to select the configuration and the substrate material of desired, and analyse the scattering properties of a single SCD resonator by multipole decomposition. Then, we select two configuration of cylinders in square lattice and frustums in hexagonal lattice, and by optimizing the geometries of cylinders in square lattice in simulation we can maximize reflectivity while minimizing absorption. Next, we adjust the geometries of the square-latticed cylinder to obtain dual-band reflectivity for certain applications. For the designs above, analyses of electromagnetic fields, effective medium properties and geometrical parameter sweepings are conducted to enhance the understanding of the optimized structures.
2. The Principles for Designing the Perfect Reflector
2.1. Guided-Mode Resonance
As mentioned before, GMR is a fundamental principle used to design metasurfaces with high reflectivity [
18].
Figure 1 (a) presents the schematic of a dual-layer GMR, a grating layer (
layer) and a waveguide layer (
layer), with a normal incident light at the intensity of
I and its reflected and transmitted light (in solid blue). The incident light is first diffracted by the grating layer, and is going to be coupled into the waveguide layer (in dashed blue) when meeting the resonance condition. Here, GMR occurs because the incident light and the waveguide mode constructively interference, resulting in the reflectivity
R approaching near-unity and the transmission
T near-zero. Otherwise, the light would transport directly through the two functional layers as well as the substrate (
layer) (in solid gray). There is another way of realizing GMR by using a a single-layer configuration, as illustrated in
Figure 1 (b). The grating layer can also act as the waveguide as long as the higher refractive index of the grating
, the lower index
and optical thickness
d of the grating structure are carefully optimized [
19,
32,
33]. Here,
represents the effective refractive index of the same layer of the grating and the waveguide. The single-layer GMR configuration can simplify the metasurface structure and facilitate dissipating heat [
15,
19].
The metasurfaces, commonly seen as a semi-finite-thickness slab, can be treated as an effective homogeneous medium, with the effective complex permittivity
and the effective permeability
expressed as
and
respectively [
34]. Following the effective medium theory, the perfect reflection can be achieved while the Eq. (
1a) and Eq. (1b) are both satisfied.
The first condition in Eq. (
1a) requires the real parts of the permittivity and permeability to have opposite signs, which is easily satisfied at the resonate wavelength of electric or magnetic resonances. Although the second condition in Eq. (1b) appears more strict, it can be satisfied at certain conditons. Particularly, it is more easily to be satisfied in all-dielectric materials with low loss, since the imaginary part of
and
are approaching zero [
12,
13]. Therefore, we can verify whether our design is indeed a perfect reflector judging by the two conditions above, of which the permittivity
and permeability
can be extracted from the S-parameter via simulation [
23,
24,
25,
26].
As mentioned before, the SCDs are the perfect candidate as the grating layer for metasurface reflector that can endure high-power lasers, but the materials for the substrate should also be thoroughly considered in order to meet our purpose. We propose the single-layer GMR configuration, using SCDs as the resonant layer with the effective index denoted as
. Here, the high refractive index equals to that of the SCDs, i.e.
, and the low refractive index is that of the air (
). And in this paper we have chosen the fused silica (silicon dioxide,
) as the substrate mainly for the two reasons below. First, its refractive index is lower than that of SCDs, meeting the conditions of the GMR theory. And second, silica has relatively high LIDT due to its extremely low thermal expansion coefficient, making it possible for high-power-laser applications [
35]. Last but not least, the diamond coating on
substrates can be achieved using Microwave Plasma Chemical Vapor Deposition technique (MPCVD) [
36].
2.2. The Multipole Decomposition
To fully understand the behavior of light scattering and reflection in these resonant structures, it is essential to decompose the scattering field into multipole components for the complexity of electromagnetic field [
29]. Multipole decomposition is a crucial method for analyzing the complex electromagnetic fields of scattering, which can separate the contributions of electric dipole (ED), magnetic dipole (MD), electric quadrupole (EQ), magnetic quadrupole (MQ), and higher-order multipoles. It effectively simplifies the description of the field, facilitating the analysis of the field distribution, which helps guide the design of highly reflective metasurfaces [
29,
37,
38].
The electric far-field scattered by an electric and a magnetic dipole (scalar projection) can be written as Eq. (
2a) and Eq. (2b) [
37].
Here,
(
) represents the amplitudes of the fields scattered by electric (magnetic) dipoles. The phase factor
accounts for a phase shift observed during dipole excitation at resonance. The function
is a Heaviside step function, i.e.
for
and
for
, which describes the antisymmetric behavior of the magnetic dipole resonance [
37]. The functions
and
present the properties shown as below.
The incident field
interact with the scattered field
or
, ending up with interference. When
or
, the destructive interference occurs in the forward (transmitted) direction, and a standing wave is generated in the backward (reflected) direction. Particularly, when only ED and MQ contribute to the scattering field, the node of the standing wave is exactly at the surface of the reflector, exhibiting the same behavior of a perfect electric conductor (PEC), and therefore it can also be defined as the
generalized electric mirror. In contrast, when the contribution purely comes from MD and EQ, an anti-node (instead of a node) should form at the surface of the mirror, which is referred to as the
generalized magnetic mirrors [
22,
37].
Here, we start from the simplest cylindrical geometry structure of diamond metasurfaces, seen in
Figure 2 (a). According to [
39], the resonance wavelengths of ED and MD in a cylindrical resonator depend on the aspect ratio (
), which is defined as
, where
H and
D represents the height and diameter of the cylinder respectively. Here, the simulations are conducted by 3D finite-difference time-domain solvers (FDTD, Lumerical Solutions, Inc.) and a total field scattered field (TFSF) source is used to analyze the scattering properties. The SCSs along with the wavelength
can be numerically calculated, and the peaks are shown at certain wavelengths, representing electric or magnetic resonance. By analysing the wavelength difference of the peaks, an optimized geometry is determined with
and
, where
P (pitch) represents the center-to-center distance between adjacent diamond columns.
Figure 2 (b) illustrates the SCS for a diamond cylinder resonator without substrate, with
ranging from 0.8 to 1.6. The calculation of the SCS relies on Hinamoto’s work [
27], and it requires the results of electric fields and refractive indices, which varies with the spatial position (
x, y, z) and the light frequency
f. Shown in
Figure 2 (b), the resonant modes happen at the certain wavelength peaks (see
Table 1 in the supplementary materials), forming either generalized electric mirror (generated by ED+MQ) or generalized magnetic mirror (generated by MD+EQ).
To further investigate the collective response of the diamond resonators, the reflective spectra at square-latticed configuration have been calculated, as shown in
Figure 2 (c). Corresponding to the SCSs, these illustrate that the yellow points of high-reflectivity are mainly the results of MD+EQ, while the peaks in green are the results of both ED+MQ and MD+EQ. However, for
and
, there are sharp drops of the spectra at
around 2.1. This is because the scattering field of different poles interference destructively [
40], resulting in the energy localized in the resonators without emitting electromagnetic radiation, which have a high Q-factor [
41].
3. Results and Analysis
3.1. Cylindrical Resonators in Square Lattices
Based on the guidelines above, we can begin to design our all-dielectric metasurfaces with high reflectivity by using diamond resonators and silica substrates.
Figure 3 (a) shows the schematic of diamond cylindrical resonators in square lattices. Simulations were performed at vertical incidence for p-polarized planwave at the wavelength between 700
and 890
.
Figure 3 (b) presents the efficiencies of reflectivity, transmission, and absorption along with the wavelength, at the optimized parameters of
,
and
. It is clearly seen that at 795
the reflectivity exceeds 99.99% and absorbance is below 0.001%, behaving as a perfect reflector. The electromagnetic fields at 795
were calculated, as displayed in
Figure 3 (c). It shows the electric field
(left-top) and magnetic field
(left-bottom) in the xy-plane, revealing the simultaneous presence of both ED and MD resonances respectively. The right diagram in
Figure 3 (c) presents the
field in the xz-plane, and a half-wave loss in the backward propagation direction is observed here [
22]. Therefore, ED is the main contribution of the scattering field, where the metasurface exhibits the behaviour of an electric mirror. Additionally, the field enhancement is predominantly localized within the diamond resonators. This reveals that most of the electromagnetic energy is confined in the diamond layer, rather than being dissipated into the silica substrate, which can minimize thermal damage that may affect the metasurface performance.
Then we have calculated the effective medium properties of the functional layer by S-parameters in FDTD solver. The results are shown in
Figure 4, with the effective permittivity (
), magnetic permeability (
), refractive index (
) and impedance (
).
Figure 4 (a,b) shows that for wavelength near 795
,
and
, so the condition for perfect reflection in Eq. (
1a) is satisfied. Meanwhile, the imaginary parts
and
change the signs at the resonance wavelength 795
, so that
and
are purely real; this meets the condition in Eq. (1b). As shown in
Figure 4 (c), the imaginary part of the refractive index (
) is proved to be maximized to prevent evanescent tunneling across the slab while the real part of the refractive index (
) equals to zero. In
Figure 4 (d) we also note that the real part of impedance approaches zero at 795
, which is corresponding to a reflectivity of unity.
To analyze the impact of structural errors on the reflection spectrum, we conducted parameter sweeping for the three structural parameters (
H,
D and
P) at the range of 700
to 890
, and the reflective spectra results are shown in
Figure 5. As shown in
Figure 5 (a), the reflectivity is higher than 98 % when
H is between 360
and 470
, providing high tolerance for the thickness of diamond film. It is advantageous for fabrication, because the highly precise control of cylindrical height is costing, as the realization of the film requires the heteroepitaxial growth of SCDs on a silica substrate. In contrast, the high-reflectivity peak redshifts to a significant degree with the growth of
P and a given pitch
P supports only a single resonance mode as seen in
Figure 5 (b).
Figure 5 (c) shows that variations in cylinder diameter
D result in multiple reflective peaks, each representing a different resonance mode. During the etching process in fabrication,
P can be controlled with relatively high precision, while the cylinder diameter
D is easily influenced by the vertical precision of etching processes. Therefore, more precise control of
D is needed.
3.2. Frustums in Hexagonal Lattices
In order to mimic the vertical errors of the cylindrical etching processes, the resonators often adopt a frustum-shape with a smaller top and larger base, rather than maintaining a perfectly cylindrical profile [
14]. On the other hand, the hexagonal lattice exhibiting six-fold symmetry (C6) have higher symmetry, resulting in the advantages of compactness, wide-angle response and broad photonic bandgap, compared with the square lattices, which possesses only four-fold symmetry (C4) [
16,
42,
43]. Consequently, we ultimately designed a metasurface reflector based on the hexagonal lattices of diamond frustums and the schematic is shown in
Figure 6 (a). Here,
P is the distance of adjacent frustums, while
H,
and
represent the height, the top diameter and bottom diameter of the frustum.
Simulations were performed at vertical incidence for p-polarized planwave at the wavelength between 700
and 900
.
Figure 6 (b) presents the efficiencies of reflectivity, transmission, and absorption along with the wavelength, at the optimized parameters of
,
,
and
. It is clearly seen that at 795
the reflectivity exceeds 99.99% and absorbance is below 0.001%, behaving as a perfect reflector.
To further analyze the impact of structural errors on the reflection spectrum, we conducted parameter sweeping for the structural parameters
H and
(see Eq. 4b) at the range of 700
to 890
and the reflective spectra results are shown in
Figure 6 (c-d).
From
Figure 6 (c), the spectral position of resonance redshifts with
growing but also remains greater than 98 % at the range of 370 - 500
which has a better robustness to the cylinders in square lattices. However, as shown in
Figure 6 (d), when we make
(see Eq.
4a) constant, the spectrum changes little when
is less than 50 nm. For frustums, the errors of
and
come from the lack of precision in the electron beam lithography and etching process. But by slightly changing the input diameter of the cylinders, we can make
equal to the theoretical value of
D, in order to compensate the resonance shift induced by machining errors.
3.3. Dual-Band High Reflectivity
The metasurface can not only be designed for single-peak high reflectivity, but can also be optimized for dual-band high-reflectivity, holding the potential for applications that utilize lasers at two or more wavelengths, e.g. optical tweezers [
44]. Here we focuse on cylindrical diamond resonators in square lattices. After adjusting the geometries, we have achieved near-unity reflectance at the wavelengths of 744
and 828
when
,
and
, and the reflective spectrum is shown in
Figure 7(a) in solid blue line. To further analyse the dual-band reflector, we have computed the SCSs of the ED+MQ, MD+EQ and the combination of their total scattering, as shown in
Figure 7(a). It is observed that the first reflection peak at 744
primarily corresponds to an electric resonance, while the second reflection peak at 828
is mainly due to a magnetic resonance.
To explain the high reflectance at the two wavelengths from the perspectives of electric and magnetic mirrors, we conducted electromagnetic field analysis at the reflective peaks. At 744
, the electric field
closely resembles the reflective behaviour of PEC, because the electric field reflected at the surface experiences a half-wave loss, which is presented in
Figure 7(b). In contrast, at the 828
wavelength, the magnetic mirror does not exhibit a half-wave loss presented in
Figure 7(c).
4. Conclusion
In this paper, we have conducted the design and analysis of highly reflective all-dielectric metasurfaces based on SCD resonators via FDTD simulations. First, utilizing GMR theory and multipole decomposition, we have acquired the preliminary values of the geometrical parameters (
and
) in a single diamond cylinder resonator for high refection at 795
. Then, we have designed two single-peak highly reflective configurations. The first is based on cylindrical resonators in square lattice, and we have achieve near-unity reflectivity (greater than 99.99%) and minimal absorption (below 0.001%) after optimization (
,
and
), as the design behaves as a general electric mirror judging from the electromagnetic field. And considering the fabrication errors, we have swept the geometries of
H,
P and
D, and found that
H has strong robustness to the reflection spectrum, tolerating an error of approximately 100
. The second is the frustum-shaped resonators in hexagonal lattices (
,
,
and
) and it also achieved 99.99% reflectivity as well as good robustness of height
H. Moreover, we have found that the adjustment of
can compensate the errors of cylindrical
D. For the two designs above, the reflective spectra within our simulated wavelength range change slightly with the height
H of the columns. It is a good robustness because in fabrication the total thickness variation (TTV) is less than 30 nm (see
Figure 8). Finally, we have demonstrated a dual-band high-reflectivity metasurfaces at 744
and 828
based on cylinders in square lattices, which are mainly resulted by ED+MQ and MD+EQ resonances respectively.
This work paves the way for low-loss all-dielectric reflective metasurfaces in high-power laser applications. For example, our method can be extended to other wavelengths, particularly in the infrared range, where a broader bandwidth is possible. Moreover, by selecting more complex pillar structures, such as cross-shaped or H-shaped columns, more design dimensions are provided, allowing for dual-band or even multi-band frequency control [
16]. Using materials with higher refractive indices, large-angle oblique incidence for mirrors with high reflectivity can also be achieved [
14].
Author Contributions
Conceptualization, Z.X.(Zhongyang Xing) and Z.X.(Zhongjie Xu); Formal Analysis, J.Z.; Methodology, J.L.; Software, J.L.; Validation, Z.X.(Zhongyang Xing); Data curation, J.L.; writing—original draft preparation, Z.X.(Zhongyang Xing) and J.L.; writing—review and editing, Z.X.(Zhongyang Xing) and J.L.; supervision, X.C., J.Z. and Z.X.(Zhongjie Xu); funding acquisition, J.Z. and Z.X.(Zhongyang Xing).All authors have read and agreed to the published version of the manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (12204541), the Science and Technology Innovation Program of Hunan Province (2021RC3083), and the High-level Talents Programs of National University of Defense Technology.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Data will be made available on request.
Acknowledgments
We thank Prof. Yueqiang Hu and Dr. Shaozhen Lou from Hunan University, Research Fellow. Mingcheng Panmai from Nanyang Technological University and Xinye Liao from National University of Defense Technology for the useful discussions.
Conflicts of Interest
The authors declare no conflicts of interest.
Abbreviations
The following abbreviations are used in this manuscript:
GMR |
Guided Mode Resonance |
FDTD |
Finite-Difference Time-Domain |
SCD |
Single-Crystal Diamond |
ED |
Electric Dipole |
MD |
Magnetic Dipole |
EQ |
Electric Quadrupole |
MQ |
Magnetic Quadrupole |
TFSF |
Total Field Scattered Field |
AR |
Aspect Ratio |
SCS |
Scattering Cross-Section |
PEC |
Perfect Electron Crystallograph |
References
- Zheludev, N.I.; Kivshar, Y.S. From Metamaterials to Metadevices. Nature Mater 2012, 11, 917–924. [CrossRef]
- Liu, Y.; Zhang, X. Metamaterials: A New Frontier of Science and Technology. Chem. Soc. Rev. 2011, 40, 2494. [CrossRef]
- Chen, W.T.; Zhu, A.Y.; Capasso, F. Flat Optics with Dispersion-Engineered Metasurfaces. Nat Rev Mater 2020, 5, 604–620. [CrossRef]
- Chen, H.T.; Taylor, A.; Nanfang, Y. A review of metasurfaces: physics and applications. Reports on Progress in Physics 2016, p. 076401.
- Quan, Y.; Qin, G.; Linsen, C.; Yi, Z.; Yuhang, Y.; Xun, C.; Shuming, W.; Shining, Z.; Zhenlin, W. Recent advanced applications of metasurfaces in multi-dimensions. Nanophotonics 2023, pp. 2295–2315.
- Fröch, J.E.; Huang, L.; Tanguy, Q.A.; Colburn, S.; Zhan, A.; Ravagli, A.; Seibel, E.J.; Böhringer, K.F.; Majumdar, A. Real time full-color imaging in a Meta-optical fiber endoscope. eLight 2023, p. 13.
- Ren, H.R.; Fang, X.Y.; Jang, J.; Burger, J.; Rho, J.; Maier, S.A. Complex-amplitude metasurface-based orbital angular momentum holography in momentum space. Nature Nanotechnology 2020, p. 948.
- Devlin, R.C.; Ambrosio, A.; Rubin, N.A.; Mueller, J.P.B.; Capasso, F. Arbitrary spin-to–orbital angular momentum conversion of light. Science 2017, 358, 896–901. [CrossRef]
- Zhang, Z.; Liu, Y.; Wang, Z.; Zhang, Y.; Guo, X.; Xiao, S.; Xu, K.; Song, Q. Folded Digital Meta-Lenses for on-Chip Spectrometer. Nano letters 2023, pp. 3459–3466.
- Ni, Y.; Chen, S.; Wang, Y.; Tan, Q.; Xiao, S.; Yang, Y. Metasurface for Structured Light Projection over 120° Field of View. Nano Letters 2020, pp. 6719–6724.
- Kuznetsov, A.I.; Miroshnichenko, A.E.; Brongersma, M.L.; Kivshar, Y.S.; Luk’yanchuk, B. Optically resonant dielectric nanostructures. Science 2016, 354, aag2472. [CrossRef]
- Slovick, B.; Yu, Z.G.; Berding, M.; Krishnamurthy, S. Perfect Dielectric-Metamaterial Reflector. Phys. Rev. B 2013, 88, 165116. [CrossRef]
- Moitra, P.; Slovick, B.A.; Gang Yu, Z.; Krishnamurthy, S.; Valentine, J. Experimental Demonstration of a Broadband All-Dielectric Metamaterial Perfect Reflector. Applied Physics Letters 2014, 104, 171102. [CrossRef]
- Moitra, P.; Slovick, B.A.; li, W.; Kravchencko, I.I.; Briggs, D.P.; Krishnamurthy, S.; Valentine, J. Large-Scale All-Dielectric Metamaterial Perfect Reflectors. ACS Photonics 2015, 2, 692–698. [CrossRef]
- Atikian, H.A.; Sinclair, N.; Latawiec, P.; Xiong, X.; Meesala, S.; Gauthier, S.; Wintz, D.; Randi, J.; Bernot, D.; DeFrances, S.; Thomas, J.; Roman, M.; Durrant, S.; Capasso, F.; Lončar, M. Diamond Mirrors for High-Power Continuous-Wave Lasers. Nat Commun 2022, 13, 2610. [CrossRef]
- Doiron, C.F.; Brener, I.; Cerjan, A. Realizing Symmetry-Guaranteed Pairs of Bound States in the Continuum in Metasurfaces. Nat Commun 2022, 13, 7534. [CrossRef]
- Magnusson, R.; Wang, S.S. New Principle for Optical Filters. Applied Physics Letters 1992, 61, 1022–1024. [CrossRef]
- Wei, CY (Wei, C.; Liu, SJ (Liu, S.; Deng, DG (Deng, D.; Shen, J (Shen, J.; Shao, JD (Shao, J.; Fan, ZX (Fan, Z. Electric field enhancement in guided-mode resonance filters. Optics letters 2006, pp. 1223–1225.
- Sang, T.; Wang, Z. Study on Guided-Mode Resonance Optical Devices; Tongji University Press, 2020.
- Lin, D.; Fan, P.; Hasman, E.; Brongersma, M.L. Dielectric Gradient Metasurface Optical Elements. Science 2014, 345, 298–302. [CrossRef]
- Liu, W.; Kivshar, Y.S. Generalized Kerker Effects in Nanophotonics and Meta-Optics [Invited]. Opt. Express 2018, 26, 13085. [CrossRef]
- Liu, W. Generalized Magnetic Mirrors. Phys. Rev. Lett. 2017, 119, 123902. [CrossRef]
- Smith, D.R.; Schultz, S.; Markoš, P.; Soukoulis, C.M. Determination of Effective Permittivity and Permeability of Metamaterials from Reflection and Transmission Coefficients. Phys. Rev. B 2002, 65, 195104. [CrossRef]
- Smith, D.R.; Schurig, D. Electromagnetic Wave Propagation in Media with Indefinite Permittivity and Permeability Tensors. PHYSICAL REVIEW LETTERS 2003, 90. [CrossRef]
- Smith, D.R.; Vier, D.C.; Koschny, Th.; Soukoulis, C.M. Electromagnetic Parameter Retrieval from Inhomogeneous Metamaterials. Phys. Rev. E 2005, 71, 036617. [CrossRef]
- Arslanagic, S.; Hansen, T.V.; Mortensen, N.A.; Gregersen, A.H.; Sigmund, O.; Ziolkowski, R.W.; Breinbjerg, O. A Review of the Scattering-Parameter Extraction Method with Clarification of Ambiguity Issues in Relation to Metamaterial Homogenization. IEEE Antennas Propag. Mag. 2013, 55, 91–106. [CrossRef]
- Hinamoto, T.; Fujii, M. MENP: An Open-Source MATLAB Implementation of Multipole Expansion for Nanophotonics. OSA Continuum 2021, 4, 1640. [CrossRef]
- Alaee, R.; Rockstuhl, C.; Fernandez-Corbaton, I. An Electromagnetic Multipole Expansion beyond the Long-Wavelength Approximation. Optics Communications 2018, 407, 17–21. [CrossRef]
- Terekhov, P.D.; Babicheva, V.E.; Baryshnikova, K.V.; Shalin, A.S.; Karabchevsky, A.; Evlyukhin, A.B. Multipole Analysis of Dielectric Metasurfaces Composed of Nonspherical Nanoparticles and Lattice Invisibility Effect. Phys. Rev. B 2019, 99, 045424. [CrossRef]
- Holly, C.; Traub, M.; Hoffmann, D.; Widmann, C.; Brink, D.; Nebel, C.E.; Gotthardt, T.; Sözbir, M.C.; Wenzel, C. Monocrystalline CVD-diamond optics for high-power laser applications. SPIE LASE, 2016.
- Krupke, W.F. Diode Pumped Alkali Lasers (DPALs)—A Review (Rev1). Progress in Quantum Electronics 2012, 36, 4–28. [CrossRef]
- Magnusson, R.; Ding, Y. Spectral-band engineering with interacting resonant leaky modes in thin periodic films. Proceedings of the SPIE - The International Society for Optical Engineering 2004, pp. 119–129.
- Sang, T (Sang, T.; Wang, ZS (Wang, Z.; Wang, L (Wang, L.; Wu, YG (Wu, Y.; Chen, LY (Chen, L. Resonant excitation analysis of sub-wavelength dielectric grating. JOURNAL OF OPTICS A-PURE AND APPLIED OPTICS 2006, pp. 62–66.
- Lalanne, P.; Lemercier-Lelanne, D. Depth dependence of the effective properties of subwavelength gratings. Journal of the Optical Society of America 1997, pp. 450–458.
- Yu, X.; Yang, W.; Shen, C.; Tao, W.; Deng, X. Polarization beam combining by fused silica subwavelength grating. Optics Communications 2024, p. 130135.
- Li, Y.; Hao, X.; Dai, B.; Shu, G.; Zhao, J.; Zhang, S.; Liu, X.; Wang, W.; Liu, K.; Cao, W.; Yang, L.; Zhu, J.; Han, J. Research Progress on High Rate and High Quality Growth of MPCVD Single Crystal Diamond. Journal of Synthetic Crystals 2020, pp. 979–989.
- Kruk, S.; Kivshar, Y. Tailoring Transmission and Reflection with Metasurfaces. In Dielectric Metamaterials; Elsevier, 2020; pp. 145–174.
- Savinov, V.; Fedotov, V.A.; Zheludev, N.I. Toroidal Dipolar Excitation and Macroscopic Electromagnetic Properties of Metamaterials. Physical Review B 2014, 89, 205112. [CrossRef]
- Campione, S.; Liu, S.; Basilio, L.I.; Warne, L.K.; Langston, W.L.; Luk, T.S.; Wendt, J.R.; Reno, J.L.; Keeler, G.A.; Brener, I.; Sinclair, M.B. Broken Symmetry Dielectric Resonators for High Quality Factor Fano Metasurfaces. ACS Photonics 2016, 3, 2362–2367. [CrossRef]
- Miroshnichenko, A.E.; Evlyukhin, A.B.; Yu, Y.F.; Bakker, R.M.; Chipouline, A.; Kuznetsov, A.I.; Luk’yanchuk, B.; Chichkov, B.N.; Kivshar, Y.S. Nonradiating Anapole Modes in Dielectric Nanoparticles. Nature Communications 2015, 6, 8069. [CrossRef]
- Hsu, C.W.; Zhen, B.; Stone, A.D.; Joannopoulos, J.D.; Soljačić, M. Bound States in the Continuum. Nature Reviews Materials 2016, 1, 16048. [CrossRef]
- Vaidya, S.; Benalcazar, W.A.; Cerjan, A.; Rechtsman, M.C. Point-Defect-Localized Bound States in the Continuum in Photonic Crystals and Structured Fibers. Phys. Rev. Lett. 2021, 127, 023605. [CrossRef]
- Wang, X.; Wang, J.; Zhao, X.; Shi, L.; Zi, J. Realizing Tunable Evolution of Bound States in the Continuum and Circularly Polarized Points by Symmetry Breaking. ACS Photonics 2023, 10, 2316–2322. [CrossRef]
- Bustamante, C.J.; Chemla, Y.R.; Liu, S.; Wang, M.D. Optical Tweezers in Single-Molecule Biophysics. Nat Rev Methods Primers 2021, 1, 25. [CrossRef]
Figure 1.
Dual-layer and single-layer GMR designs for high-reflectivity metasurfaces. (a) Schematic illustration of the dual-layer GMR design, where is the effective refractive index of grating layer, is the refractive index of waveguide layer, and is for the substrate. The blue and gray lines show the pathway of resonance and non-resonance relatively. Here, is the grating period and d represents the thickness of waveguide layer. I, R, T are the intensity of incident, reflected and transmitted light. (b) The single-layer GMR structure where the grating layer function as waveguide layer as well. Here, is the higher refractive index and is the lower.
Figure 1.
Dual-layer and single-layer GMR designs for high-reflectivity metasurfaces. (a) Schematic illustration of the dual-layer GMR design, where is the effective refractive index of grating layer, is the refractive index of waveguide layer, and is for the substrate. The blue and gray lines show the pathway of resonance and non-resonance relatively. Here, is the grating period and d represents the thickness of waveguide layer. I, R, T are the intensity of incident, reflected and transmitted light. (b) The single-layer GMR structure where the grating layer function as waveguide layer as well. Here, is the higher refractive index and is the lower.
Figure 2.
SCSs and reflectance spectrum of diamond resonators in square lattices by 3D FDTD solvers. (a) Schematic of the FDTD simulation setup for a single diamond resonator embedded in air, using a total field scattered field (TFSF) source to calculate the SCSs. (b) SCS of generalized electric mirror (ED+MQ) and magnetic mirror (MD+EQ) for . (c) Reflective spectrum of diamond resonators arrays in square lattices for .
Figure 2.
SCSs and reflectance spectrum of diamond resonators in square lattices by 3D FDTD solvers. (a) Schematic of the FDTD simulation setup for a single diamond resonator embedded in air, using a total field scattered field (TFSF) source to calculate the SCSs. (b) SCS of generalized electric mirror (ED+MQ) and magnetic mirror (MD+EQ) for . (c) Reflective spectrum of diamond resonators arrays in square lattices for .
Figure 3.
(a) Simulation schematic of diamond cylindrical resonators arranged in a square lattice, where Dis the diameter of cylinder while H is the height, and P represents the distance of adjecent cylinders. (b) Reflectivity, transmission, and absorption plotted on a logarithmic scale with , and . R is reflectance, T stands for transmittance and A is absorption. (c) The electric field and magnetic field in the xy-plane, and the field in the xz-plane. Colors indicate the electric field amplitude.
Figure 3.
(a) Simulation schematic of diamond cylindrical resonators arranged in a square lattice, where Dis the diameter of cylinder while H is the height, and P represents the distance of adjecent cylinders. (b) Reflectivity, transmission, and absorption plotted on a logarithmic scale with , and . R is reflectance, T stands for transmittance and A is absorption. (c) The electric field and magnetic field in the xy-plane, and the field in the xz-plane. Colors indicate the electric field amplitude.
Figure 4.
Effective medium properties of diamond metasurface mirror. (a) Effective permittivity () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of . (b) Effective permeability () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of . (c) Effective refractive index () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of . (d) Effective impedance () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of .
Figure 4.
Effective medium properties of diamond metasurface mirror. (a) Effective permittivity () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of . (b) Effective permeability () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of . (c) Effective refractive index () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of . (d) Effective impedance () as a function of wavelength. Re() represents the real part of and Im() represents the imaginary part of .
Figure 5.
Sweeping map of cylinder diamond resonantors in square lattices at the wavelength between 700 nm and 890 nm. (a) Cylinder height from 200 nm-600 nm. Colors indicates the reflectivity. (b) Pitch in the square lattice from 310 nm-570 nm. (c) Cylinder diameter from 240 nm-490 nm.
Figure 5.
Sweeping map of cylinder diamond resonantors in square lattices at the wavelength between 700 nm and 890 nm. (a) Cylinder height from 200 nm-600 nm. Colors indicates the reflectivity. (b) Pitch in the square lattice from 310 nm-570 nm. (c) Cylinder diameter from 240 nm-490 nm.
Figure 6.
Frustum-shaped diamond resonantors in hexagonal lattices. (a) Simulation schematic of diamond frustum-shaped resonators arranged in a hexagonal lattice.Here, H is the height, and P represents the distance of adjacent frustums, while represents the top diameter of the frustum and is the bottom. (b) Reflectivity, transmission, and absorption plotted on a logarithmic scale with , , and . (c) Sweeping map of frustum height, showing that the design remains robust within the height range of 380 nm to 480 nm. Colors indicates the reflectivity. (d) Sweeping map of (half the gap between the top and bottom diameters of the frustum), showing how increasing the diameter affects the reflection spectrum.
Figure 6.
Frustum-shaped diamond resonantors in hexagonal lattices. (a) Simulation schematic of diamond frustum-shaped resonators arranged in a hexagonal lattice.Here, H is the height, and P represents the distance of adjacent frustums, while represents the top diameter of the frustum and is the bottom. (b) Reflectivity, transmission, and absorption plotted on a logarithmic scale with , , and . (c) Sweeping map of frustum height, showing that the design remains robust within the height range of 380 nm to 480 nm. Colors indicates the reflectivity. (d) Sweeping map of (half the gap between the top and bottom diameters of the frustum), showing how increasing the diameter affects the reflection spectrum.
Figure 7.
Cylindrical diamond resonantors in square lattices and electric fields at two high-reflectivity wavelength, compared with PEC. (a) Reflectivity spectrum and SCS of ED+MQ and MD+EQ resonances, showing dual-band reflectivity peaks at 744 nm and 828 nm. (b) The electric field field in the xz-plane at 744 compared with PEC, showing a half-wave phase loss, indicating the behavior of an generalized electric mirror. (c) The electric field field in the xz-plane at 828 compared with PEC, showing a half-wave phase loss, where no half-wave phase loss is observed, demonstrating the characteristics of a generalized magnetic mirror.
Figure 7.
Cylindrical diamond resonantors in square lattices and electric fields at two high-reflectivity wavelength, compared with PEC. (a) Reflectivity spectrum and SCS of ED+MQ and MD+EQ resonances, showing dual-band reflectivity peaks at 744 nm and 828 nm. (b) The electric field field in the xz-plane at 744 compared with PEC, showing a half-wave phase loss, indicating the behavior of an generalized electric mirror. (c) The electric field field in the xz-plane at 828 compared with PEC, showing a half-wave phase loss, where no half-wave phase loss is observed, demonstrating the characteristics of a generalized magnetic mirror.
Figure 8.
The AFM test results for the four test points over a scan area of 10 mm × 10 mm. Measurements were taken at two locations: the center and the edge of the film on two separate samples, providing a total of four data points.
Figure 8.
The AFM test results for the four test points over a scan area of 10 mm × 10 mm. Measurements were taken at two locations: the center and the edge of the film on two separate samples, providing a total of four data points.
Table 1.
SCSs of the two peaks for different AR values and multipoles
Table 1.
SCSs of the two peaks for different AR values and multipoles
|
|
Peaks 1 |
Peaks 2 |
Multipole |
AR |
|
SCS () |
|
SCS () |
ED+MQ |
0.8 |
2.0761 |
12.679 |
2.4374 |
2.5640 |
|
1.0 |
2.0991 |
9.5938 |
1.5586 |
2.6994 |
|
1.2 |
2.0871 |
132.10 |
2.6196 |
3.3584 |
|
1.4 |
2.1361 |
1761.5 |
2.6456 |
4.5045 |
|
1.6 |
2.2652 |
605.26 |
2.6607 |
6.3820 |
MD+EQ |
0.8 |
2.0951 |
9.9177 |
2.3654 |
13.946 |
|
1.0 |
2.1552 |
9.7153 |
2.4765 |
11.432 |
|
1.2 |
2.0881 |
18.442 |
2.5025 |
11.549 |
|
1.4 |
2.1361 |
618.95 |
2.5365 |
13.930 |
|
1.6 |
2.1652 |
282.97 |
2.5816 |
14.364 |
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).