Next Article in Journal
Potential Efficient Separation of Oil from Bilgewater and Kitchen Wastewater by Fractional Freezing Process
Previous Article in Journal
Effect of Microstructure and Mechanical Properties of Al5083 Alloy Processed by ECAP at Room Temperature and High Temperature
Previous Article in Special Issue
Engineered Layer-Stacked Interfaces Inside Aurivillius-Type Layered Oxides Enables Superior Ferroelectric Property
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Gap-Size-Dependent Effective Phase Transition in Metasurfaces of Closed-Ring Resonators

1
Department of Mechanical Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
2
Department of Optical Engineering, Kongju National University, Cheonan 31080, Korea
3
Department of Future Convergence Engineering, Kongju National University, Cheonan 31080, Korea
4
Institute of Application and Fusion for Light, Kongju National University, Cheonan 31080, Korea
*
Author to whom correspondence should be addressed.
Crystals 2021, 11(6), 684; https://doi.org/10.3390/cryst11060684
Submission received: 13 April 2021 / Revised: 4 June 2021 / Accepted: 9 June 2021 / Published: 14 June 2021
(This article belongs to the Special Issue Advanced Research on Heterogeneous Materials)

Abstract

:
We theoretically investigate a metal-to-insulator transition in artificial two-dimensional (2D) crystals (i.e., metasurfaces) of tightly coupled closed-ring resonators. Strong interaction between unit resonators in the metasurfaces yields the effective permittivity highly dependent on the lattice spacing of unit resonators. Through our rigorous theory, we provide a closed form of effective permittivity of the metasurface and reveal that the permittivity possesses a Lorentzian-type resonant behavior, implying that the transition of the effective permittivity can arise when the lattice spacing passes a critical value.

1. Introduction

Metal-to-insulator transition is an essential physical phenomenon in crystal solids as it results in drastic changes in electrical and optical responses of materials [1,2,3]. In artificial crystals, such as three-dimensional metamaterials and two-dimensional metasurfaces, the notion of transition can be extended in the context of an effective medium theory. It was shown that artificial assemblies of subwavelength unit resonators can exhibit either metal or insulator phases [4,5] during light–matter interaction. Hyperbolic metasurfaces [6,7,8] that utilize both metal and insulator phases for separate principal axes are now at the core of various applications [9,10,11,12], realizing unprecedented optical properties. Moreover, in high-density metasurfaces, the transition is shown to be controllable by the coupling strength between unit resonators [13,14].
Likewise, owing to its importance, the metal-to-insulator transition in metasurfaces has been investigated in various structures. However, despite such extensive studies, detailed physical accounts (i.e., physical interpretations of the effective phase and the transition condition) have not been intensively made yet. More importantly, understanding of the underlying physics of the phase transition is still desired to capture an intuitive picture of the transition and for further development of relevant applications. In this paper, we consider metasurfaces with tightly coupled resonators and show that the metal-to-insulator transition arises in metasurfaces depending on the strength of coupling. In particular, we make a rigorous analysis on a model metasurface composed of closed-ring unit resonators and provide accurate and physically intuitive criteria for the transition. When resonators are capacitively coupled to each other with a deep subwavelength lattice spacing, we find that the metal-to-insulator transition occurs when the lattice spacing exceeds a critical value. We also clarify the underlying physical mechanism of the phase transition in terms of light transmission in metasurfaces through two competing transmission channels.

2. Theoretical Model

The phase of a metasurface can be expressed in terms of its effective permittivity and permeability. These effective parameters are governed by the interaction of light with constituent resonator elements. For tightly coupled resonators, the conventional effective medium approaches used with dilute metasurfaces, such as the Maxwell Garnett theory, do not apply [15]. Here, we make a rigorous coupled-mode theory calculation to derive effective parameters for dense metasurfaces. Consider an exemplary metasurface composed of subwavelength rectangular ring resonators, as shown in Figure 1. Resonators are spaced periodically with gap sizes gx and gy in the x- and y-directions.
To simplify the problem, we assume that the rectangular rings are made of a perfect electric conductor (PEC) with side lengths a, b (a < b) smaller than the wavelength λ of the incident light. For tightly coupled resonators with narrow gaps gx, gyλ, nearby resonators can interact with each other through charges induced by incident light [16]. These capacitive interactions open a highly efficient transmission channel through the gap (channel c(2)) mediated by the funneling process [17,18]. Together with the channel c(1) through the open rectangular hole region, the gap channel c(2) determines the phase of the metasurface. To determine the transition analytically, we adopt the coupled-mode theory [19,20]. For an x-polarized electric field incident normally onto the metasurface (region II, −h/2 < z < h/2) as illustrated in Figure 1, magnetic fields in regions I (z < −h/2) and III (z > h/2) can be expanded in terms of quantized diffracted waves of orders m and n:
H I = Z 0 1 [ y ^ e i k 0 ( z + h 2 ) + m , n = R m n e i χ m n ( z + h 2 ) + i ϕ m n ] H I I I = Z 0 1 m , n = T m n e i χ m n ( z h 2 ) + i ϕ m n ,
where Z0 ≡ (μ0/ε0)1/2 is the free-space impedance, φmnαm(xa/2) + βn(yb/2), αm ≡ 2πm/dx, βn ≡ 2πn/dy, χ m n ( k 0 2 α m 2 β n 2 ) 1 / 2 , and k0 ≡ 2π/λ. Inside the metasurface region II, we treat the rectangular hole region and gap region separately. In the gap region, we note that the x-polarized electric field barely couples to the gap region parallel to the x-axis (−w < x < a + w, b + w < y < b + w + gy) due to the strong cut-off condition. Subsequently, we neglect electromagnetic (EM) waves inside this gap region and consider only the remaining gap c(2) (a + wxa + w + gx). In the rectangular region c(1) (0 ≤ xa, 0 ≤ yb), we make a single-mode approximation and consider only the fundamental waveguide mode [20]. This allows for writing the EM waves inside the metasurface as
E I I ( 1 ) = ( sin ( π y b ) ( A 1 e i k 1 z + B 1 e i k 1 z ) , 0 , 0 ) ,   H I I ( 2 ) = ( 0 , Z 0 1 ( A 2 e i k 0 z + B 0 e i k 0 z ) , 0 ) ,
with k 1 ( k 0 2 π 2 / b 2 ) 1 / 2 . The remaining field components can be obtained from Maxwell’s equations. By applying continuity boundary conditions of tangential EM fields at two interfaces z = ±h/2, all the undetermined coefficients Aj, Bj, T m n ( T x , m n , T y , m n , T z , m n ) , and R m n ( R x , m n , R y , m n , R z , m n ) can be fixed. Specifically, one can readily show that the continuity of Ey at two interfaces (z = ±h/2) and the divergencelessness of the magnetic field give constraints such that [20]
R x , m n = α m β n α m 2 + χ m n 2 R y , m n ,   R z , m n = β n χ m n α m 2 + χ m n 2 R y , m n , T x , m n = α m β n α m 2 + χ m n 2 T y , m n ,   T z , m n = β n χ m n α m 2 + χ m n 2 T y , m n .
Then, the continuity conditions of Ex at two interfaces can be simplified as
1 m , n = ( k 0 χ m n α m 2 + χ m n 2 R y , m n e i ϕ m n ) = sin ( π y b ) ( A 1 e i k 1 h / 2 + B 1 e i k 1 h / 2 ) ( on   c ( 1 ) ) = A 2 e i k 0 h / 2 + B 2 e i k 0 h / 2 ( on   c ( 2 ) ) = 0 ( otherwise ) ,
and
m , n = ( k 0 χ m n α m 2 + χ m n 2 T y , m n e i ϕ m n ) = sin ( π y b ) ( A 1 e i k 1 h / 2 + B 1 e i k 1 h / 2 ) ( on   c ( 1 ) ) = A 2 e i k 0 h / 2 + B 2 e i k 0 h / 2 ( on   c ( 2 ) ) = 0   ( otherwise ) .
Now, the inverse Fourier transform of Equations (4) and (5) gives rise to
δ m 0 δ n 0 k 0 χ m n α m 2 + χ m n 2 R y , m n = Q m n [ A 1 e i k 1 h / 2 + B 1 e i k 1 h / 2 ] + Q ˜ m n [ A 2 e i k 0 h / 2 + B 2 e i k 0 h / 2 ] , k 0 χ m n α m 2 + χ m n 2 T y , m n = Q m n [ A 1 e i k 1 h / 2 + B 1 e i k 1 h / 2 ] + Q ˜ m n [ A 2 e i k 0 h / 2 + B 2 e i k 0 h / 2 ] .
Here, δmn is the Kronecker delta. Q m n and Q ˜ m n are vector dependencies between external diffracted waves and internal waveguide modes in two channels, and are defined as
Q m n 1 d x d y 0 a 0 b d x d y sin ( π y b ) e i ϕ m n , Q ˜ m n 1 d x d y a + w a + w + g x 0 d y d x d y e i ϕ m n .
Equation (6) gives the relationship between the waveguide modes in the metasurface and the diffracted waves in the free space. In particular, we note that the zero-th order transmission and reflection coefficients can be written as
T y , 00 = 2 a b d x d y π [ A 1 e i k 1 h / 2 + B 1 e i k 1 h / 2 ] + g x d x [ A 2 e i k 0 h / 2 + B 2 e i k 0 h / 2 ] , R y , 00 = 1 2 a b d x d y π [ A 1 e i k 1 h / 2 + B 1 e i k 1 h / 2 ] g x d x [ A 2 e i k 0 h / 2 + B 2 e i k 0 h / 2 ] .
The results so far are obtained from the continuity of the tangential components of the electric fields at two interfaces. In order to fix the undetermined coefficients A1,2 and B1,2, the same procedure can be applied to the continuity of the magnetic fields. After a straightforward calculation, one can find that A1,2 and B1,2 can be obtained from a matrix equation C = M−1G, where a 4 × 4 matrix M is defined by
M ( e i k 1 h 2 ( W 1 , 1 + k 1 k 0 ) e i k 1 h 2 ( W 1 , 1 k 1 k 0 ) e i k 0 h 2 W 1 , 2 e i k 0 h 2 W 1 , 2 e i k 1 h 2 ( W 1 , 1 k 1 k 0 ) e i k 1 h 2 ( W 1 , 1 + k 1 k 0 ) e i k 0 h 2 W 1 , 2 e i k 0 h 2 W 1 , 2 e i k 1 h 2 W 2 , 1 e i k 1 h 2 W 2 , 1 e i k 0 h 2 ( W 2 , 2 + 1 ) e i k 0 h 2 ( W 2 , 2 1 ) e i k 1 h 2 W 2 , 1 e i k 1 h 2 W 2 , 1 e i k 0 h 2 ( W 2 , 2 1 ) e i k 0 h 2 ( W 2 , 2 + 1 ) ) ,
with two column vectors CT ≡ (A1, B1, A2, B2) and GT ≡ (π/8, 0, 2, 0), where T denotes the transpose. Wj,l is the light-channel coupling factor given by
W j , l 2 δ j 2 S ( c ( j ) ) d x d y m , n = α m 2 + χ m n 2 k 0 χ m n ( c ( j ) d x d y u y ( j ) e i ϕ m n ) ( c ( l ) d x d y u y ( l ) e i ϕ m n ) .
Here, S(c(j)) is the cross-sectional area of the channel c(j) in the xy plane, and u y ( 1 ) sin ( π y / b ) and u y ( 2 ) 1 are the waveguide mode profile. When the ring resonators are in a deep subwavelength scale, the coupling factors can be approximately written as
W 1 , 1 8 a b π 2 d x d y + i 4 a d y 2 k 0 π b 3 d x , W 2 , 2 g x ( 1 d x + i k 0 π Ω ) , W 1 , 2 4 g x π d x i g x 2 4 k 0 a π 2 Ω ˜ , W 2 , 1 2 a b π d x d y i g x 2 k 0 b π 2 d y Ω ˜ ,
where
Ω ln ( 2 g x π d x ) 3 2 , Ω ˜ [ ln ( 2 π ( w + g x ) d x ) + 2 w g x ln ( 2 π ( w + g x ) d x ) + w 2 g x 2 ln ( w + g x w ) 3 w g x 3 2 ] .
The deep subwavelength thickness of the metasurface allows for approximating Equation (9). By expanding up to linear terms of k0h and k1h, we finally arrive at simplified zero-th order transmission and reflection coefficients that are necessary to define the effective parameters:
T y , 00 F + i σ , R y , 00 1 F + i σ , F 1 π η [ 2 a b d x d y π ( π W 1 , 2 + 4 W 2 , 2 2 i k 0 h ) + g x d x ( π W 1 , 1 i π k 1 2 h 2 k 0 4 W 2 , 1 ) ] , σ k 0 h ( 4 a b d x d y π 2 + g x 2 d x ) , η W 1 , 1 W 2 , 2 i W 1 , 1 k 0 h 2 i W 2 , 2 k 1 2 h 2 k 0 W 1 , 2 W 2 , 1 .

3. Retrieval of Effective Indices of the Metasurface from Transmission and Reflection Coefficients

The main idea of defining the effective indices is to replace the inhomogeneous structures by unstructured homogeneous media by keeping the scattering properties unchanged. Since the near and far fields in unstructured media are exactly the same, the homogenization is usually done by matching the far-field scattering properties, which is valid only for the structures having subwavelength periodicity. The retrieval methods of effective indices are introduced in various literatures. Here, we follow the scheme present in [21]. In terms of transmission and reflection coefficients, T and R, the effective refractive index and impedance can be written as [21].
n e f f = 1 k 0 h arccos ( 1 + T 2 R 2 2 T ) , Z e f f = ( 1 R 2 ) 2 T 2 ( 1 + R 2 ) 2 T 2 .
Since we assume that the metasurface is made of PEC and the lattice constant is smaller than the incident wavelength, the system is lossless, and only zero-th order diffracting waves can reach the far field. We note that this constraint allows a further simplification of Equation (14). One can readily show that the transmission and reflection coefficients in the complex plane are orthogonal to each other when the media are lossless:
Re ( T ) Re ( R ) + Im ( T ) Im ( R ) = 0 .
Together with energy conservation, | T | 2 + | R | 2 = 1 , one can find that the relation ( 1 + T 2 R 2 ) / 2 T = Re ( T ) / | T | 2 holds. Consequently, Equation (14) can be reduced as
n e f f = 1 k 0 h arccos ( Re ( T ) | T | 2 ) , Z e f f = Re ( R ) | R | 2 Re ( R ) + | R | 2 .
Replacing T and R in Equation (16) with Ty,00 and Ry,00 yields the effective index and impedance of the metasurface. For instance, let us examine the refractive index. By putting Equation (13) into | T y , 00 | 2 + | R y , 00 | 2 = 1 , one can find that | T y , 00 | 2 can be approximated as Re ( F ) + 2 Im ( F ) σ . This gives a significantly simplified form of the refractive index:
n e f f = 1 k 0 h arccos ( Re ( T y , 00 ) | T y , 00 | 2 ) 1 k 0 h arccos ( Re ( F ) Re ( F ) + 2 Im ( F ) σ ) = 1 k 0 h arccos ( 1 2 Im ( F ) σ Re ( F ) + 2 Im ( F ) σ ) = 1 k 0 h arccos ( 1 2 Re ( η ) σ Im ( η ) + 2 Re ( η ) σ ) .
In the last step in Equation (17), we used the fact that Re ( π W 1 , 2 + 4 W 2 , 2 ) = 0 and Re ( π W 1 , 1 4 W 2 , 1 ) = 0 from Equations (11) and (13). Note that σ is real-valued and very small when h , g x a , b . Equations (16) and (17) clearly show that there is a singularity at T y , 00 = 0 = Re ( F ) + 2 Im ( F ) σ . One can find that this singularity corresponds to the transition point at which the effective refractive indices change from real (imaginary) to imaginary (real), and that the absolute value of Im ( F ) σ / ( Re ( F ) + 2 Im ( F ) σ ) is always smaller than 1 except when it is very close to the singularity. This allows for approximating Equation (17) such that
n e f f 1 k 0 h 4 Re ( η ) σ Im ( η ) + 2 Re ( η ) σ ( 1 + O ( σ 2 ) + ) .
The validity of this approximation is shown in Figure 2. However, note that Equation (18) holds the same singularity with Equation (17), implying that the transition condition in Equation (17) is still valid in Equation (18). By putting all geometric parameters in Equations (11) and (13) into Equation (18) and keeping up to the linear terms of g x / Γ and h / Γ , where Γ is one of all other geometric parameters much larger than gx and h, we arrive at
n e f f 8 a d y 2 h k 0 2 d x b 3 π Ω π h / 2 g x π h / 2 g x Ω π 2 d y 3 / k 0 2 b 4 d x .
Depending on structural parameters, this index is either pure imaginary or real, thereby representing the metal and insulator phases, respectively. When structural parameters are varied continuously, we observe that the change of refractive index from pure imaginary to real, or metal-to-insulator transition, occurs when the denominator passes zero, or when the gap size gx exceeds the critical value gt, which is determined by
h π g t 2 ln ( 2 π g t d x ) = λ 2 d y 3 4 d x b 4 3 .
This is a condition for the phase transition with a critical lattice spacing. We point out that the same condition can be derived from the vanishing transmission, Ty,00 = 0. In the transition condition in Equation (20), gt is determined only implicitly. For convenience, we re-express the transition condition in terms of the critical frequency ωt for given gap sizes gxy such that
ω t π d y c 0 b 2 d y 2 d x ( h π 2 g x ln ( 2 g x π d x ) + 3 2 ) 1 / 2
where ω t = 2 π c 0 / λ t , c0 is the speed of light in a vacuum, and λt is the gap-size-dependent transition wavelength. Then, Equation (19) can be simplified in terms of transition frequency:
n e f f 8 a d y 2 c 0 2 d x b 3 h π 1 ω t 2 ω 0 2 .
Note that Equation (22) is valid when g x , h d x , d y < λ . In the same manner, explicit calculation also shows that the effective impedance in Equation (16) can be approximated as
Z e f f 8 a h b 5 c 0 2 π 3 d x d y 4 ω t 2 ω 2 * ,
and effective permittivity and permeability of the metasurface by
ε e f f x x c 0 2 π d y 3 h b 4 1 ω t 2 ω 2 , μ e f f x x 8 a b π 2 d x d y .
We note that the effective permittivity possesses a Lorentzian-type resonant behavior. We also note that similar expressions for ε e f f y y and μ e f f y y can be obtained with x, y interchanged.
Figure 3 shows the wavelength- and gap-dependent effective permittivity of the metasurface. One can clearly see that, for a fixed wavelength, the phase transition can be controlled by varying gap size at a deep subwavelength scale. We point out that our theory successfully resolves the extraordinarily high refractive index in [4]. We also note that the effective permeability is independent of gap size agreeing with the result in [22].

4. Discussion

4.1. The Origin of the Metal-to-Insulator Transition: Macroscopic Channel Competition

So far, we discussed the gap-size-dependent metal-to-insulator transition based on the far-field transmission and reflection coefficients. One of the main advantages of using the coupled-mode theory is that the theory provides the full descriptions of the electromagnetic field distribution near the metasurface, allowing understanding of how the transition is related to the macroscopic near-field features. Specifically, one can rewrite the far-field transmission coefficient in Equation (8) as
T y , 00 = ( c ( 1 ) d x d y + c ( 2 ) d x d y ) E x I I ( x , y ) d x d y .
We previously discussed that the transition condition can be also obtained by letting Ty,00 = 0. This means that, from Equation (25), the phase transition occurs when the transmissions of the Ex fields in two channels cancel each other out. To be specific, we first evaluate the integral over the c(1) channel, averaging Ex on a unit lattice. This gives rise to
E x ( 1 ) ¯ c ( 1 ) E x I I d x d y d x d y | z = 0 2 a b ω 2 π d x d y × ( ω 2 ω t 2 ) i α ( ω 2 ω t 2 ) 2 + α 2 ,
where α π d y 3 ω 2 / 2 b 4 k 0 > 0 . Comparing with Equation (24), we find that the sign of Re ( E x ( 1 ) ¯ ) is opposite to that of the effective permittivity and vanishes at the transition frequency, whereas Im ( E x ( 1 ) ¯ ) does not vanish. For the second channel c(2), one can find that an explicit calculation yields the same trend except for the overall sign change. Therefore, the metal-to-insulator transition condition can be restated in terms of lattice-averaged electric field components by
Re ( E x ( 1 ) ¯ ) = Re ( E x ( 2 ) ¯ ) = 0 , Im ( E x ( 1 ) ¯ ) + Im ( E x ( 2 ) ¯ ) = 0 .
On the other hand, the magnetic field components in each channel can be written as H y I I ( 1 ) 4 Z 0 1 sin ( π y / b ) / π and H y I I ( 2 ) Z 0 1 with negligible imaginary. Consequently, the total energy flux through the metasurface, equivalent to the z-component of the time-averaged Poynting vector S Re ( E × H * ) / 2 , can be simplified as S z Re ( E x ) H y / 2 . This gives rise to an intuitive criterion for the metal and insulator phases of a metasurface through
S z ( 1 ) > 0 , S z ( 2 ) < 0 , ε e f f < 0   for   metal , S z ( 1 ) < 0 , S z ( 2 ) > 0 , ε e f f > 0   for   insulator .
This is one of our main observations. The effective phase becomes metal when incident light transmits through the evanescent channel c(1); otherwise, it becomes insulator. This means that the effective phase follows the characteristic of the channel through which incident light passes. Therefore, the phase transition is a result of the competition between those light channels to admit the incident light.

4.2. Test of the Validity of the Effective Parameters

In order to check the validity of our theory, we examined light transmission through multilayered metasurfaces, as shown in Figure 4. The metasurfaces, samples 1 and 2, are made of rectangular rings as depicted in Figure 1 with a = b, dx = dy = d, and gxgy. Sample 2 is the 90-degree rotated version of sample 1. Both samples 1 and 2 are designed to exhibit strongly anisotropic optical responses in x- and y-directions and to have opposite signs in the effective permittivities in a particular spectral range (the pink-colored regions in Figure 4). The spacing between each layer is kept as d/2π to allow two neighboring metasurfaces to interact with each other via only zero-th order waves. Note that the momentum of the first-order diffracted wave is given by χ 01 = χ 10 = ( k 0 2 ( 2 π / d ) 2 ) 1 / 2 . Therefore, for a deep subwavelength periodicity d, the longest penetration depth for the diffracted wave is approximately given by ( k 0 2 ( 2 π / d ) 2 ) 1 / 2 d / 2 π . This spacing is required to use the same effective permittivities in both single- and multilayered systems.
Shown in Figure 4b,c are the transmission spectra of x- and y-polarized incident lights, respectively, through the 1-2-1 layered system. For the numerical calculations by implementing inhomogeneous 2D metasurfaces, we used the finite-difference time-domain (FDTD) method. In the colored spectral range, our system effectively becomes a metal–insulator–metal (M–I–M) or I–M–I multilayered system depending on the polarization. In Figure 4b, we can see a peak in the colored region that corresponds to resonant tunneling through two metallic barriers, and another peak at a longer wavelength that originates from the Fabry–Perot resonance. With an incident light rotated as shown in Figure 4c, two peaks appear in the colored regions. These peaks also correspond to resonant tunneling, which appear in negative-permittivity media as demonstrated in [23]. For comparison, we calculated the transmission spectra through multilayered homogeneous slabs possessing corresponding effective permittivities. We note that all transmission spectra are in good agreement.

5. Conclusions

In conclusion, we theoretically demonstrated metal-to-insulator transition in metasurfaces that can be controlled by the lattice spacing of the unit resonators. Through analytic calculations on the model metasurface, we rigorously defined the critical gap size gt at which the transition occurs. Additionally, by seeking the macroscopic details of the electromagnetic waves in the proximity of the metasurface, we defined quantitative criteria for the phase transition. Our theory demonstrates that drastic changes in optical responses can emerge from the interactions between unit resonators, an observation that is essential for the engineering of metasurfaces with tightly coupled unit resonators for deformable, reconfigurable, and functionalized metamaterials aiming for ultra-high-sensitive optical responses and dynamic manipulation of their effective properties.

Author Contributions

Conceptualization, J.-H.K.; methodology, S.L. and J.-H.K.; investigation, S.L. and J.-H.K.; data curation, S.L. and J.-H.K.; writing—original draft preparation, S.L. and J.-H.K.; writing—review and editing, S.L. and J.-H.K.; supervision, J.-H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a research grant from Kongju National University in 2020 and by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (grant nos. NRF-2020R1F1A1076886 and NRF-2021R1A2C2012617). S.L. was supported by an NRF grant funded by the Korean government (MSIT) (grant no. NRF-2020R1C1C1012138).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the corresponding author on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zylbersztejn, A.; Mott, N.F. Metal-insulator transition in vanadium dioxide. Phys. Rev. B 1975, 11, 4383. [Google Scholar] [CrossRef] [Green Version]
  2. Kibis, O.V. Metal-insulator transition in graphene induced by circularly polarized photons. Phys. Rev. B Condens. Matter Mater. Phys. 2010, 81, 165433. [Google Scholar] [CrossRef] [Green Version]
  3. Seo, M.; Kyoung, J.; Park, H.; Koo, S.; Kim, H.S.; Bernien, H.; Kim, B.J.; Choe, J.H.; Ahn, Y.H.; Kim, H.T.; et al. Active terahertz nanoantennas based on VO2 phase transition. Nano Lett. 2010, 10, 2064–2068. [Google Scholar] [CrossRef] [Green Version]
  4. Choi, M.; Lee, S.H.; Kim, Y.; Kang, S.B.; Shin, J.; Kwak, M.H.; Kang, K.Y.; Lee, Y.H.; Park, N.; Min, B. A terahertz metamaterial with unnaturally high refractive index. Nature 2011, 470, 369–373. [Google Scholar] [CrossRef] [PubMed]
  5. Takano, K.; Shibuya, K.; Akiyama, K.; Nagashima, T.; Miyamaru, F.; Hangyo, M. A metal-to-insulator transition in cut-wire-grid metamaterials in the terahertz region. J. Appl. Phys. 2010, 107, 024907. [Google Scholar] [CrossRef]
  6. Smith, D.R.; Schurig, D. Electromagnetic Wave Propagation in Media with Indefinite Permittivity and Permeability Tensors. Phys. Rev. Lett. 2003, 90, 077405. [Google Scholar] [CrossRef] [PubMed]
  7. Lee, S.; Kim, W.T.; Kang, J.H.; Kang, B.J.; Rotermund, F.; Park, Q.H. Single-Layer Metasurfaces as Spectrally Tunable Terahertz Half- A nd Quarter-Waveplates. ACS Appl. Mater. Interfaces 2019, 11, 7655–7660. [Google Scholar] [CrossRef] [PubMed]
  8. Luan, P.G. Effective electrodynamics theory for the hyperbolic metamaterial consisting of metal–dielectric layers. Crystals 2020, 10, 863. [Google Scholar] [CrossRef]
  9. Biehs, S.A.; Tschikin, M.; Ben-Abdallah, P. Hyperbolic metamaterials as an analog of a blackbody in the near field. Phys. Rev. Lett. 2012, 109, 104301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Salandrino, A.; Engheta, N. Far-field subdiffraction optical microscopy using metamaterial crystals: Theory and simulations. Phys. Rev. B Condens. Matter Mater. Phys. 2006, 74, 075103. [Google Scholar] [CrossRef] [Green Version]
  11. Slobozhanyuk, A.P.; Ginzburg, P.; Powell, D.A.; Iorsh, I.; Shalin, A.S.; Segovia, P.; Krasavin, A.V.; Wurtz, G.A.; Podolskiy, V.A.; Belov, P.A.; et al. Purcell effect in hyperbolic metamaterial resonators. Phys. Rev. B Condens. Matter Mater. Phys. 2015, 92, 195127. [Google Scholar] [CrossRef] [Green Version]
  12. Maas, R.; Parsons, J.; Engheta, N.; Polman, A. Experimental realization of an epsilon-near-zero metamaterial at visible wavelengths. Nat. Photonics 2013, 7, 907–912. [Google Scholar] [CrossRef]
  13. Decker, M.; Feth, N.; Soukoulis, C.M.; Linden, S.; Wegener, M. Retarded long-range interaction in split-ring-resonator square arrays. Phys. Rev. B Condens. Matter Mater. Phys. 2011, 84, 085416. [Google Scholar] [CrossRef] [Green Version]
  14. Kang, J.; Lee, S.; Kang, B.J.; Kim, W.T.; Rotermund, F.; Park, Q.-H. Anomalous Wavelength Scaling of Tightly-Coupled Terahertz Metasurfaces. ACS Appl. Mater. Interfaces 2018, 10, 19331–19335. [Google Scholar] [CrossRef] [PubMed]
  15. Yoo, S.J.; Park, Q.-H. Effective permittivity for resonant plasmonic nanoparticle systems via dressed polarizability. Opt. Express 2012, 20, 16480–16489. [Google Scholar] [CrossRef]
  16. Kim, K.-H.; Jung, G.-H.; Lee, S.-J.; Park, H.-G.; Park, Q.-H. Ultrathin Capacitive Metasurfaces for Strong Electric Response. Adv. Opt. Mater. 2016, 4, 1501–1506. [Google Scholar] [CrossRef]
  17. Kang, J.H.; Kim, D.S.; Park, Q.H. Local capacitor model for plasmonic electric field enhancement. Phys. Rev. Lett. 2009, 102, 093906. [Google Scholar] [CrossRef]
  18. Seo, M.A.; Park, H.R.; Koo, S.M.; Park, D.J.; Kang, J.H.; Suwal, O.K.; Choi, S.S.; Planken, P.C.M.; Park, G.S.; Park, N.K.; et al. Terahertz field enhancement by a metallic nano slit operating beyond the skin-depth limit. Nat. Photonics 2009, 3, 152–156. [Google Scholar] [CrossRef]
  19. Kang, J.-H.; Wang, S.; Shi, Z.; Zhao, W.; Yablonovitch, E.; Wang, F. Goos-Hänchen Shift and Even–Odd Peak Oscillations in Edge-Reflections of Surface Polaritons in Atomically Thin Crystals. Nano Lett. 2017, 17, 1768–1774. [Google Scholar] [CrossRef]
  20. Kang, J.H.; Choe, J.-H.; Kim, D.-S.; Park, Q.-H. Substrate effect on aperture resonances in a thin metal film. Opt. Express 2009, 17, 15652–15658. [Google Scholar] [CrossRef]
  21. Smith, D.R.; Vier, D.C.; Koschny, T.; Soukoulis, C.M. Electromagnetic parameter retrieval from inhomogeneous metamaterials. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2005, 71, 036617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Pendry, J.B.; Martín-Moreno, L.; Garcia-Vidal, F.J. Mimicking surface plasmons with structured surfaces. Science 2004, 305, 847–848. [Google Scholar] [CrossRef] [PubMed]
  23. Zhou, L.; Wen, W.; Chan, C.T.; Sheng, P. Electromagnetic-wave tunneling through negative-permittivity media with high magnetic fields. Phys. Rev. Lett. 2005, 94, 243905. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of a metasurface consisting of rectangular closed-ring resonators of perfect electric conductors. The incident is an x-polarized plane wave impinging upon the metasurface normally. c(1) and c(2) denote the two dominant light channels in the metasurface.
Figure 1. Schematic of a metasurface consisting of rectangular closed-ring resonators of perfect electric conductors. The incident is an x-polarized plane wave impinging upon the metasurface normally. c(1) and c(2) denote the two dominant light channels in the metasurface.
Crystals 11 00684 g001
Figure 2. Comparison of exact and approximated effective refractive indices calculated from Equations (14) and (19), respectively. Here, we set a = b, dx = dy = 1.021a, gx = gy = 0.001a, and h = 0.01a.
Figure 2. Comparison of exact and approximated effective refractive indices calculated from Equations (14) and (19), respectively. Here, we set a = b, dx = dy = 1.021a, gx = gy = 0.001a, and h = 0.01a.
Crystals 11 00684 g002
Figure 3. (a) Analytically calculated wavelength- and gap-dependent effective permittivity of a metasurface of a = b, w = 0.01a, h = 0.01a. (b) The upper panel shows a cross-cut spectrum of effective permittivity (purple) and permeability (blue) at gx = 0.005a. The cut line is shown in (a) as a dotted line. The bottom panel shows corresponding effective refractive indices.
Figure 3. (a) Analytically calculated wavelength- and gap-dependent effective permittivity of a metasurface of a = b, w = 0.01a, h = 0.01a. (b) The upper panel shows a cross-cut spectrum of effective permittivity (purple) and permeability (blue) at gx = 0.005a. The cut line is shown in (a) as a dotted line. The bottom panel shows corresponding effective refractive indices.
Crystals 11 00684 g003
Figure 4. (a) Schematic of two metasurfaces, 1 and 2, and corresponding effective permittivities where gx = 0.0025a, gy = 0.0188a, h = 0.03a, and d = 1.0375a. (b) FDTD-calculated transmission spectrum (solid line) through the 1-2-1 layer with x-polarized incident light and (c) with y-polarized incident light. Dotted lines are the spectrum obtained from the transfer matrix scheme applied to multilayered homogeneous slabs possessing effective permittivities corresponding to (a).
Figure 4. (a) Schematic of two metasurfaces, 1 and 2, and corresponding effective permittivities where gx = 0.0025a, gy = 0.0188a, h = 0.03a, and d = 1.0375a. (b) FDTD-calculated transmission spectrum (solid line) through the 1-2-1 layer with x-polarized incident light and (c) with y-polarized incident light. Dotted lines are the spectrum obtained from the transfer matrix scheme applied to multilayered homogeneous slabs possessing effective permittivities corresponding to (a).
Crystals 11 00684 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, S.; Kang, J.-H. Gap-Size-Dependent Effective Phase Transition in Metasurfaces of Closed-Ring Resonators. Crystals 2021, 11, 684. https://doi.org/10.3390/cryst11060684

AMA Style

Lee S, Kang J-H. Gap-Size-Dependent Effective Phase Transition in Metasurfaces of Closed-Ring Resonators. Crystals. 2021; 11(6):684. https://doi.org/10.3390/cryst11060684

Chicago/Turabian Style

Lee, Seojoo, and Ji-Hun Kang. 2021. "Gap-Size-Dependent Effective Phase Transition in Metasurfaces of Closed-Ring Resonators" Crystals 11, no. 6: 684. https://doi.org/10.3390/cryst11060684

APA Style

Lee, S., & Kang, J. -H. (2021). Gap-Size-Dependent Effective Phase Transition in Metasurfaces of Closed-Ring Resonators. Crystals, 11(6), 684. https://doi.org/10.3390/cryst11060684

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop