Next Article in Journal
Theoretical Analysis of Thermal Distribution and Waveform Evolution in Pulsed Ytterbium-Doped Fiber Amplifier with Extra Feedback
Previous Article in Journal
Quantitative Endogenous Fluorescence Analysis Discriminates Glioblastoma Tumor through Visible and NIR Excitation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

An Efficient Method for Light Beaming from Subwavelength Slits Surrounded by Surface Gratings

School of Mathematics and Computational Science, Xiangtan University, Xiangtan 411105, China
Submission received: 8 March 2023 / Revised: 5 April 2023 / Accepted: 10 April 2023 / Published: 12 April 2023

Abstract

:
This paper investigates the directional beaming of metallic subwavelength slits surrounded by dielectric gratings. The design of the structure for light beaming was formulated as an optimization problem for the far-field angular transmission. A vertical mode expansion method was developed to solve the diffraction problem, which was then integrated into a genetic algorithm and an active set method to obtain the optimal result. By using the proposed method for a two-slit structure, we demonstrate that both the single- and dual-beaming effects can be efficiently achieved. Moreover, the beaming directions can be flexibly adjusted and precisely controlled.

1. Introduction

Since the discovery of extraordinary optical transmission through subwavelength hole arrays perforated in metal films [1], optical field diffraction by small apertures has been intensively studied from both the theoretical and application aspects. One of the interesting results is the light beaming from a subwavelength aperture with metal grating structures discovered by Lezec et al. [2]. Springing from the demonstration of this beaming effect, many efforts have been made to realize and optimize directional beaming. The most-popular models are the grating-based beaming structures, such as a single slit surrounded by metallic grooves or dielectric gratings [3,4,5,6,7,8,9], subwavelength slit arrays [10,11,12], and annular aperture array structures [13,14,15,16], to name a few. These methods generally rely on the excitation of surface plasmon polaritons (SPPs) on the corrugated metal surface and reradiated in a specific direction. In recent years, there have also been interesting methods for manipulating beams that involve the use of metasurfaces and photonic crystals, etc. [17,18,19,20,21,22].
The physical mechanism for light beaming from grating-based beaming structures has been thoroughly investigated [3,23,24,25,26,27]. These studies provide valuable information on choosing some key parameters for light beaming structures. However, the effects of many geometrical parameters are interlinked, and efficient numerical methods are still required to optimize the performance of light beaming structures in practical applications. This would facilitate the efficient and precise manipulation of light beams. General computational methods, such as the finite-difference time domain (FDTD) method and the finite-element method (FEM), have been used to study this problem [9,20]. However, they are time consuming, since a small grid or mesh size must be taken, which gives rise to a large indefinite linear system and is expensive to solve.
In this work, we developed a vertical mode expansion method (VMEM) for the light beaming structure, as shown in Figure 1. The structure can be separated into a number of regions, where each region is piecewise uniform along the vertical direction. In each uniform region, using a proper reference solution, the wave field can be expanded by eigenmodes of a one-dimensional differential operator. The expansion coefficients are solved from a linear system obtained by matching the wave field at the interfaces between neighboring regions. Since only the vertical direction is discretized, the method is extremely fast and efficient, which provides a powerful computational engine for the optimization of light beaming structures.
The aim of this paper was to demonstrate that beaming of light can be achieved by optimizing the far-field angular transmission over geometrical parameters that define the optical structures, and the beaming directions can be precisely controlled. In this work, the problem was formulated as a bound constrained optimization problem. The objective function, depending on practical purposes, was formulated as a functional of the far-field angular transmission. The problem was solved by the integration of the VMEM and optimization algorithms. More precisely, the objective function was evaluated using the VMEM, and the genetic algorithm and active set method were combined to obtain the optimal result.
The rest of this paper is organized as follows. In Section 2, we formulate the diffraction problem. In Section 3, we develop a vertical mode expansion method. In Section 4, the far-field angular transmission is defined and evaluated analytically by the VMEM. In Section 5, a two-slit beaming structure is considered, and both the single- and dual-beaming effects are realized to illustrate the effectiveness of the method we developed. The paper is concluded with some general remarks in Section 6.

2. Problem Formulation

We considered the diffraction of light beaming structures with metallic subwavelength slits surrounded by gratings in the output surface. The schematic approach of the structure with a single slit is shown in Figure 1.
In the Cartesian coordinate system, the structure is invariant in y and parallel to the x y plane. The lower boundary is placed on or slightly above z = 0 , and we assumed the half-plane R 2 = { ( x , z ) R 2 : z < 0 } is the free space. The metal slab was assumed to be thick enough that it is impenetrable. The structure is multi-layered, and it can be separated into a number of regions S l , l = 1 , 2 , , N ; in each region, the material is uniform in the x-direction. Denoting x 1 = and x N + 1 = + , then S l can be represented as
S l = { ( x , z ) : x l < x < x l + 1 , z R } , l = 1 , 2 , , N .
In each region S l , the materials vary only in the vertical direction; let ε and μ be the relative permittivity and the relative permeability, respectively, then
ε = ε ( l ) ( z ) , μ = μ ( l ) ( z ) , ( x , z ) S l .
Above the structure, we specified a TM polarized incident plane wave H y ( i ) = e i ( α x β z ) , where α = k 0 sin θ i n and β = k 0 cos θ i n , θ i n ( π / 2 , π / 2 ) is the angle of incidence with respect to the negative z-axis, and k 0 is the wavenumber of the free space. For TM polarization ( H y , E x and E z are the nonzero components of the electromagnetic fields), the total field H y satisfies the following Helmholtz equation:
x 1 ε H y x + z 1 ε H y z + k 0 2 μ H y = 0 .
To address the diffraction problem, we developed a new vertical mode expansion method, which provides a powerful tool for computing and optimizing plasmonic structures.

3. Vertical Mode Expansion Method

In this section, we developed a vertical mode expansion method for the problem formulated above. The light beaming structure is multi-layered and separated into a number of regions S l , l = 1 , 2 , , N . The basic idea of the vertical mode expansion method is to represent the total field in each region as the superposition of a reference field H y ( l ) and a sub-field H y ( s ) = H y H y ( l ) , i.e.,
H y = H y ( l ) + H y ( s ) .
The reference field H y ( l ) was chosen such that it is a solution excited by the incident wave on a hypothetical infinite 1D structure with the same profile as S l . Thus, H y ( l ) also satisfies the Helmholtz equation and can be written as
H y ( l ) = e i α x H ˜ y ( l ) ( z ) ,
where H ˜ y ( l ) ( z ) is a one-dimensional function of z and can be evaluated analytically; the details can be found in [28]. With the properly chosen reference field, the sub-field H y ( s ) is shown to be outgoing as z ± and satisfies the homogeneous Helmholtz equation in S l :
x 1 ε ( l ) H y ( s ) x + z 1 ε ( l ) H y ( s ) z + k 0 2 μ ( l ) H y ( s ) = 0 .
We then performed a separation of variables for H y ( s ) such that
H y ( s ) = ϕ ( l ) ( z ) V ( l ) ( x ) ,
then the function ϕ ( l ) and V ( l ) satisfy
ε ( l ) d d z 1 ε ( l ) d ϕ ( l ) d z + k 0 2 ε ( l ) μ ( l ) ϕ ( l ) = η ( l ) 2 ϕ ( l ) ,
d 2 V ( l ) d x 2 + η ( l ) 2 V ( l ) = 0 .
Since the vertical axis z is unbounded, the PML technique can be applied to truncate z to a finite interval. A PML corresponds to a complex coordinate stretching, which replaces z by a complex z ^ [29]. In practice, we only need to replace d z in (8) by s ( z ) d z for some complex function s ( z ) . With the truncation of z by PMLs, the eigenvalue problem (8) can be solved by a pseudospectral method as in [30,31], and a discrete sequence of eigenpairs { η j ( l ) , ϕ j ( z ) ( l ) } for j = 1 , 2 , 3 , can be obtained. Based on those vertical modes, the total field in region S l , l = 1 , 2 , , N can be represented as
H y = H y ( l ) + j = 1 ϕ j ( l ) V j ( l ) .
The other two components E x and E z can then be obtained from Maxwell’s equations, and we have
E x = E x ( l ) + 1 i k 0 ε ( l ) j = 1 d ϕ j ( l ) z V j ( l ) , E z = E z ( l ) 1 i k 0 ε ( l ) j = 1 ϕ j ( l ) d V j ( l ) d x ,
where E x ( l ) = e i α x E ˜ z ( z ) and E z ( l ) = e i α x E ˜ z ( z ) are components of the reference solution. The general solution of V j ( l ) in (9) is simply
V j ( l ) = a j l e i η j ( l ) ( x x l ) + b j l e i η j ( l ) ( x x l + 1 ) , l = 1 , , N ,
where a j l and b j l are the coefficients of the left-going and right-going eigenmodes in region S l . Noticing that, in region S ( 1 ) and S ( N ) , the fields are outgoing as x ± , thus a j 1 = b j N = 0 , j = 1 , 2 , 3 , . A linear system:
A u = c
can be established by matching the continuity condition of H y and E z on the vertical boundaries between the interfaces x l , l = 2 , 3 , , N . The unknown u is a vector for the coefficients a j l and b j l . In the fully discretized version, the variable z is discretized by N z points z n for 1 n N z . We obtain N z numerical modes, then the index j ranges from 1 to N z , and A is a 2 ( N 1 ) N z × 2 ( N 1 ) N z matrix. The right-hand side of (13) is related to the difference between the reference solutions in S l and S l 1 . Notice that the reference solutions for regions S l and S l 1 are different. Therefore, the right-hand side of (13) is nonzero in general. After the linear system is solved, the field components in all regions can be easily evaluated.

4. Far-Field Angular Transmission

For light beaming, it is important to calculate the far-field angular transmission, and we derived it in this section. As depicted in Figure 1, the total transmitted power T can be represented as the limit of the integral of Poynting vector S ( r ) over the half circle B r :
T = lim r B r S · r ^ d s = π / 2 π / 2 lim r r S · r ^ d θ ,
where r = | r | is the radius of the circle and r ^ = r / r is the unit vector from the origin to the point r . The far-field angular transmission can, thus, be defined as P ( θ ) = lim r r S · r ^ ; since the far fields are transverse to r ^ , the wave energy propagates only in the radial direction, and we have
P ( θ ) = 1 2 lim r r H y ( r ) 2 .
For any point r R 2 , the magnetic field H y ( r ) is related to the electric field E x on z = 0 by
H y ( r ) = k 0 2 H 0 ( 1 ) ( k 0 | r x x ^ | ) E x ( x , 0 ) d x ,
where H 0 ( 1 ) is the zeroth-order Hankel function of the first kind and x ^ is the unit vector of the x-axis [32]. The far-field can then be simplified by using the asymptotic expression for the Hankel function, and we have
H y ( r ) k 0 2 π r e i π 4 E x ( x , 0 ) e i k 0 sin θ x d x .
By substituting (17) into (15), the far-field angular transmission can, thus, be evaluated as
P ( θ ) = k 0 4 π E x ( x , 0 ) e i k 0 sin θ x d x 2 = k 0 4 π J ( θ ) 2 ,
where we denote the integral as J ( θ ) = E x ( x , 0 ) e i k 0 sin θ x d x . For slit–groove structures, i.e., slits surrounded by metallic grooves, the evaluation of J ( θ ) can be simplified if a PEC or a surface impedance boundary condition is applied, since by placing the bottom surface of the structure on z = 0 , the infinite integral interval can be reduced to the finite opening areas of the slits and grooves [3,33]. However, for a real metal structure or a structure with a grating protruding through the metal surface, the integral has to be evaluated on the infinite interval. In such cases, if the domain is truncated or discretized in the x-direction, the evaluation will be impossible or inefficient. One attractive advantage of the VMEM for light beaming is that J ( θ ) can be calculated analytically.
In the VMEM, the fields on z = 0 are given by the expansions (10) and (11), and the integral can be written as the summation over each region; we obtain
J ( θ ) = l = 1 N E x ˜ ( 0 ) I 1 l + 1 i k 0 ε 0 j = 1 d ϕ j ( l ) d z ( 0 ) I 2 , j l ,
where
I 1 l = x l x l + 1 e i α x e i k 0 sin θ x d x , I 2 , j l = x l x l + 1 V j ( l ) e i k 0 sin θ x d x .
By substituting the general solution of V j ( l ) (12) into (20), the integrals above can be evaluated analytically. More precisely, in region S l , l = 2 , 3 , , N 1 , we have
I 1 l = e i ( α k 0 sin θ ) x l + 1 e i ( α k 0 sin θ ) x l i ( α k 0 sin θ ) , I 2 , j l = a j l e i ( η j ( l ) k 0 sin θ ) δ l 1 i ( η j ( l ) k 0 sin θ ) e i k 0 sin θ x l + b j l e i ( η j ( l ) + k 0 sin θ ) δ l 1 i ( η j ( l ) + k 0 sin θ ) e i k 0 sin θ x l + 1 ,
where δ l = x l + 1 x l is the width of region S l . In the first and last regions, the widths of the regions are infinite, and the integrals should be considered separately. Since the structure is impenetrable in S 1 and S N , the reference fields vanish at z = 0 in these two regions, and we can simply let I 1 l = 0 , for l { 1 , N } . In addition, the fields in these two regions are outgoing as x ± ; thus, a j 1 = 0 and b j N = 0 , and this yields
I 2 , j 1 = b j 1 i e i k 0 sin θ x 2 η j ( 1 ) + k 0 sin θ , I 2 , j N = a j N i e i k 0 sin θ x N η j ( N ) k 0 sin θ .
After the linear system (13) is solved, J ( θ ) can then be evaluated. Consequently, the far-field angular transmission P ( θ ) is obtained from (18), and the total transmitted power is
T = k 0 4 π π / 2 π / 2 J ( θ ) 2 d θ .
Light beaming can be achieved by manipulating the far-field angular transmission, and we demonstrate this in the next section.

5. Optimization and Discussion

In this section, we considered the diffraction of a two-slit structure and demonstrate that light beaming can be achieved by manipulating the far-field angular transmission. Depending on the practical purposes, the objective functions can be formulated as functionals of the far-field angular transmission P ( θ ) . Light beaming can then be achieved by optimizing the objective function over incident angles and geometrical parameters that define the plasmonic structures, and the beaming angles can be precisely controlled.
We considered a structure with two slits surrounded by a few surface dielectric gratings, as depicted in Figure 2. The thickness of the metal slab is 0.3 μm, and the widths of the slits are d 1 and d 2 . We set d 1 + d 2 = 0.1 μm to maintain a constant total slit width. The slits are surrounded by dielectric gratings in the output surface, with a height of t = 0.08 μm. The structure has five, three, and five dielectric bars on the left, middle, and right part, respectively. The dielectric bars are periodically arranged in each part, with a fill factor of 0.5 and periods of p 1 , p 0 and p 2 . The offsets distances are denoted as l i ( i = 0 , 1 , 2 , 3 ) and are depicted in Figure 2. A TM polarized plane wave with wavelength λ = 0.532 μmm is incident on the top surface of the structure. At this wavelength, the refractive index of the metal slab and the dielectric bar are 0.13 + 3.19 i and 1.72, respectively. We adopted a Cartesian coordinate system with the bottom of the structure at z = 0 and the origin at the center of the middle part of the structure. For this structure, our purpose was to show that both the single- and dual-directional beaming effects can be efficiently achieved, and the beaming angles can be precisely controlled.
We first considered the single-beaming effect; it is expected to achieve highly efficient directional beaming with any desired angle by optimizing the structure parameters and incident angle. We define the normalized far-field angular transmission as P ( θ ) / T , which describes the far-field distribution, where T is the total transmission, as given in (23). For a given predefined beaming angle θ , we denote the parameters to be optimized as x = [ p 0 , p 1 , p 2 , l 0 , l 1 , l 2 , l 3 , d 1 , θ i n ] . The optimization problem can then be formulated as
min x Λ P ( θ , x ) / T ,
where Λ is the feasible set of x . In this work, Λ represents the bound constraints of the variables, and we set the search ranges as follows: p i ( i = 0 , 1 , 2 ) ( 0 , 2 λ ] , l i ( i = 0 , 1 , 2 , 3 ) ( 0 , 2 λ ] , d 1 ( 0 , 0.1 ) , and θ i n ( π / 2 , π / 2 ) . The optimization problem (24) was solved by the integration of the VMEM and efficient optimization algorithms. Specifically, the objective function was evaluated using the VMEM, and a genetic algorithm was employed to work out the optimal result to the first significant decimal place. This result was then used as an initial guess, and the final optimal result was obtained using the active set method [34].
For illustration purposes, we considered predefined beaming angles ranging from 0 to 30 in steps of 10 , using the proposed method, the final optimal results were obtained and are listed in Table 1. It was found that the periods of the gratings have some regularity. As the beaming angle θ increased, the period p 2 became larger, while p 0 and p 1 became smaller. Furthermore, for all the off-axis cases (i.e., θ 0 ), p 1 < p 2 . The periods of the gratings play a crucial role in determining the beaming angles.
In order to directly observe the beaming effect, the magnetic field intensity | H y | 2 and far-field angular transmission for the optimal cases are shown in Figure 3 and Figure 4. Figure 3a corresponds to the on-axis beaming case (i.e., θ = 0 ), while Figure 3b–d correspond to the off-axis beaming cases with θ = 10 , 20 , and 30 , respectively. The far-field angular spectrum in shown in Figure 4; it can be seen that the angles correspond to the peak intensities coinciding with the predefined angles; this means that the beaming angles are precisely controlled.
We next considered the dual-beaming effect. For light beaming of a two-slit plasmonic structure, a natural question to consider is whether dual-beams can be realized and well-separated. While a method for rendering bundle beams from multiple subwavelength slits was proposed in [12], it relies on the destructive interference effects of beams, and the resulting beaming angles can only be very close to the on-axis direction and cannot be controlled. In this paper, we propose a method for achieving separated beams by optimizing the far-field angular transmission. Suppose the predefined two beaming angles are θ 1 and θ 2 ; if the far-field angular transmission has two separated peaks at these two angles and the intensities are very similar to each other, then dual-separated beams can be realized. To achieve this, we formulated the optimization problem as follows:
min x Λ P ( θ 1 , x ) + P ( θ 2 , x ) P ( θ 1 + θ 2 2 , x ) T + γ P ( θ 2 , x ) P ( θ 1 , x ) ,
where the feasible set Λ and the parameter x to be optimized are the same as the single-beaming cases, and γ > 0 is a penalty parameter. The first term of the objective function aims to produce two separated peaks in the far-field angular transmission at θ 1 and θ 2 , while the second term penalizes differences in the peak intensities to ensure that they are similar. The optimization problem (25) can be solved using the same method as in the single-beaming cases.
To illustrate the results, we considered both symmetric and non-symmetric dual-beaming effects (i.e., θ 1 + θ 2 = 0 and θ 1 + θ 2 0 , respectively). The optimal results are listed in Table 2. To directly observe the beaming effect, we present the magnetic field intensity | H y | 2 and far-field angular transmission for the optimal cases in Figure 5 and Figure 6. Symmetric dual-beaming effects are shown in Figure 5a,b, with beaming angles of ( 5 , 5 ) and ( 10 , 10 ) , respectively. Non-symmetric dual-beaming effects are shown in Figure 5c,d, with beaming angles of ( 15 , 5 ) and ( 10 , 20 ) . It can be seen that the beaming angles can be flexibly adjusted to reside on the same or opposite side with respect to the on-axis direction. Figure 6 presents the far-field angular spectrum, in which two separated peaks are clearly displayed with almost identical intensity. Additionally, the positions of the peaks closely match the predefined angles, indicating precise control of the beaming angles.
Finally, we give some discussions on the physical mechanism of the directional beaming in the design of plasmonic structures. Firstly, the two slits were illuminated with TM polarized incident light from above, with the slit width much smaller than a wavelength. Consequently, only fundamental modes were considered to exist in the slits and the propagation constants were determined by the slit widths [35]. Secondly, the transmitted light through the slits diffracted into gratings, and under certain resonance conditions, surface plasmon polaritons (SPPs) were excited, with energy confined to the bottom of the dielectric bars. The SPPs were then radiated into free space, forming the directional beaming. In this study, it was evident that the incident angle and slit width were the parameters used to tune the optical phase retardation at the exits of the slits. The beaming angles are mainly determined by the periods of the gratings from the grating equation [35]. More explicitly, for a predefined beaming angle θ , the periods of the gratings are related to θ approximately as follows:
p 1 2 π k s p p + k 0 sin θ , p 2 2 π k s p p k 0 sin θ .
Here, k s p p denotes the wave number of the excited SPP. This explains the data in Table 1, which show that, as the beaming angle θ increases, p 1 decreases and p 2 increases, and p 1 < p 2 for θ > 0 . In fact, the SPP in the left part of the structure propagates left, with the direction of the radiated power determined by p 1 as in the first equation of (26). Similarly, the SPP in the right part of the structure propagates right, and the direction of the radiated power is determined by p 2 . When the left-going and right-going SPPs are radiated in the same direction, a single-beaming effect was observed. For the dual-beaming cases, however, the SPPs are radiated into two different directions.

6. Conclusions

In this paper, we investigated the phenomenon of light beaming in plasmonic structures comprised of a metallic slab with subwavelength slits surrounded by gratings. The design of the structure for light beaming was formulated as a bound constrained optimization problem. To solve the diffraction problem efficiently, we developed a vertical mode expansion method, which served as a direct solver for the optimization problem. The optimal result was obtained through the integration of the VMEM and efficient optimization algorithms, namely the genetic algorithm and active set method. Numerical experiments demonstrated that light beaming can be achieved by optimizing the far-field angular transmission over the incident angle and the geometrical parameters that define the structure. For a two-slit structure, we showed that both the single- and dual-beaming effects can be realized, and the beaming directions can be flexibly adjusted and precisely controlled. These results have some potential applications in directional antennas, optical sensing, data storage, and so on.

Funding

This work was supported by the Natural Science Foundation of Hunan Province (Grant No. 2019JJ50611) and by the Natural Science Foundation of China (Grant No. 11801484).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Ebbesen, T.W.; Lezec, H.J.; Ghaemi, H.; Thio, T.; Wolff, P.A. Extraordinary optical transmission through sub-wavelength hole arrays. Nature 1998, 391, 667–669. [Google Scholar] [CrossRef]
  2. Lezec, H.J.; Degiron, A.; Devaux, E.; Linke, R.; Martin-Moreno, L.; Garcia-Vidal, F.; Ebbesen, T. Beaming light from a subwavelength aperture. Science 2002, 297, 820–822. [Google Scholar] [CrossRef] [PubMed]
  3. Martin-Moreno, L.; Garcia-Vidal, F.; Lezec, H.; Degiron, A.; Ebbesen, T. Theory of highly directional emission from a single subwavelength aperture surrounded by surface corrugations. Phys. Rev. Lett. 2003, 90, 167401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Garcia-Vidal, F.; Martin-Moreno, L.; Lezec, H.; Ebbesen, T. Focusing light with a single subwavelength aperture flanked by surface corrugations. Appl. Phys. Lett. 2003, 83, 4500–4502. [Google Scholar] [CrossRef]
  5. Shi, H.; Wang, C.; Du, C.; Luo, X.; Dong, X.; Gao, H. Beam manipulating by metallic nano-slits with variant widths. Opt. Express 2005, 13, 6815–6820. [Google Scholar] [CrossRef] [PubMed]
  6. Kim, S.; Kim, H.; Lim, Y.; Lee, B. Off-axis directional beaming of optical field diffracted by a single subwavelength metal slit with asymmetric dielectric surface gratings. Appl. Phys. Lett. 2007, 90, 051113. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, Y.; Shi, H.; Wang, C.; Du, C.; Luo, X. Multiple directional beaming effect of metallic subwavelength slit surrounded by periodically corrugated grooves. Opt. Express 2008, 16, 4487–4493. [Google Scholar] [CrossRef]
  8. Kim, H.; Park, J.; Lee, B. Tunable directional beaming from subwavelength metal slits with metal–dielectric composite surface gratings. Opt. Lett. 2009, 34, 2569–2571. [Google Scholar] [CrossRef] [Green Version]
  9. Daniel, S.; Bawuah, P. Plasmonic implanted nanogrooves for optical beaming. Sci. Rep. 2019, 9, 391. [Google Scholar] [CrossRef] [Green Version]
  10. Xu, T.; Wang, C.; Du, C.; Luo, X. Plasmonic beam deflector. Opt. Express 2008, 16, 4753–4759. [Google Scholar] [CrossRef]
  11. Jiang, J.; Wu, H.; Jiang, L.; Li, X. Genetic optimization of double subwavelength metal slits surrounded by surface dielectric gratings for directional beaming manipulation. Opt. Commun. 2012, 285, 2201–2206. [Google Scholar] [CrossRef]
  12. Kim, S.; Lim, Y.; Park, J.; Lee, B. Bundle beaming from multiple subwavelength slits surrounded by dielectric surface gratings. J. Light. Technol. 2010, 28, 2023–2029. [Google Scholar]
  13. Zhang, Q.; Hu, P.; Liu, C. Realization of enhanced light directional beaming via a Bull’s eye structure composited with circular disk and conical tip. Opt. Commun. 2015, 339, 216–221. [Google Scholar] [CrossRef]
  14. Wang, J.; Zhou, W. Subwavelength beaming using depth-tuned annular nanostructures. J. Mod. Opt. 2009, 56, 919–926. [Google Scholar] [CrossRef]
  15. Zhang, M.; Du, J.; Shi, H.; Yin, S.; Xia, L.; Jia, B.; Gu, M.; Du, C. Three-dimensional nanoscale far-field focusing of radially polarized light by scattering the SPPs with an annular groove. Opt. Express 2010, 18, 14664–14670. [Google Scholar] [CrossRef]
  16. Caglayan, H.; Bulu, I.; Ozbay, E. Beaming of electromagnetic waves emitted through a subwavelength annular aperture. J. Opt. Soc. Am. B 2006, 23, 419–422. [Google Scholar] [CrossRef]
  17. Caglayan, H.; Bulu, I.; Ozbay, E. Off-axis directional beaming via photonic crystal surface modes. Appl. Phys. Lett. 2008, 92, 092114. [Google Scholar] [CrossRef] [Green Version]
  18. Kramper, P.; Agio, M.; Soukoulis, C.; Birner, A.; Müller, F.; Wehrspohn, R.; Gösele, U.; Sandoghdar, V. Highly directional emission from photonic crystal waveguides of subwavelength width. Phys. Rev. Lett. 2004, 92, 113903. [Google Scholar] [CrossRef] [Green Version]
  19. Hajian, H.; Ozbay, E.; Caglayan, H. Beaming and enhanced transmission through a subwavelength aperture via epsilon-near-zero media. Sci. Rep. 2017, 7, 4741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Zhu, H.; Yin, X.; Chen, L.; Li, X. Directional beaming of light from a subwavelength metal slit with phase-gradient metasurfaces. Sci. Rep. 2017, 7, 12098. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Deng, Z.L.; Cao, Y.; Li, X.; Wang, G.P. Multifunctional metasurface: From extraordinary optical transmission to extraordinary optical diffraction in a single structure. Photon. Res. 2018, 6, 443–450. [Google Scholar] [CrossRef]
  22. Fu, T.; Wang, D.; Yang, Z.; Deng, Z.l.; Liu, W. Steering Smith-Purcell radiation angle in a fixed frequency by the Fano-resonant metasurface. Opt. Express 2021, 29, 26983–26994. [Google Scholar] [CrossRef] [PubMed]
  23. Yu, L.B.; Lin, D.Z.; Chen, Y.C.; Chang, Y.C.; Huang, K.T.; Liaw, J.W.; Yeh, J.T.; Liu, J.M.; Yeh, C.S.; Lee, C.K. Physical origin of directional beaming emitted from a subwavelength slit. Phys. Rev. B 2005, 71, 041405. [Google Scholar] [CrossRef]
  24. Dai, W.; Soukoulis, C. Theoretical analysis of the surface wave along a metal-dielectric interface. Phys. Rev. B 2009, 80, 155407. [Google Scholar] [CrossRef] [Green Version]
  25. Burrow, G.M.; Gaylord, T.K. Multi-beam interference advances and applications: Nano-electronics, photonic crystals, metamaterials, subwavelength structures, optical trapping, and biomedical structures. Micromachines 2011, 2, 221–257. [Google Scholar] [CrossRef]
  26. Li, L.; Li, T.; Wang, S.; Zhu, S. Collimated plasmon beam: Nondiffracting versus linearly focused. Phys. Rev. Lett. 2013, 110, 046807. [Google Scholar] [CrossRef] [Green Version]
  27. Hong, J.S.; Chen, K.R. Light diffraction by a slit and grooves with a point source model based on wave dynamics. Phys. Rev. A 2017, 96, 043813. [Google Scholar] [CrossRef]
  28. Lu, X. An efficient method for the RCS reduction of cavities in the ground plane. Adv. Appl. Math. Mech. accepted.
  29. Chew, W.C.; Weedon, W.H. A 3D perfectly matched medium from modified Maxwell’s equations with stretched coordinates. Microw. Opt. Technol. Lett. 1994, 7, 599–604. [Google Scholar] [CrossRef]
  30. Song, D.; Yuan, L.; Lu, Y.Y. Fourier-matching pseudospectral modal method for diffraction gratings. J. Opt. Soc. Am. A 2011, 28, 613–620. [Google Scholar] [CrossRef] [Green Version]
  31. Lu, X.; Shi, H.; Lu, Y.Y. Vertical mode expansion method for transmission of light through a single circular hole in a slab. J. Opt. Soc. Am. A 2014, 31, 293–300. [Google Scholar] [CrossRef] [PubMed]
  32. Jin, J.M. The Finite Element Method in Electromagnetics; John Wiley & Sons: Hoboken, NJ, USA, 2015. [Google Scholar]
  33. Dai, W.; Soukoulis, C. Control of beaming angles via a subwavelength metallic slit surrounded by grooves. Phys. Rev. B 2010, 82, 045427. [Google Scholar] [CrossRef] [Green Version]
  34. Nocedal, J.; Wright, S.J. Numerical Optimization; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  35. Lee, B.; Kim, S.; Kim, H.; Lim, Y. The use of plasmonics in light beaming and focusing. Prog. Quantum Electron. 2010, 34, 47–87. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of a metal slit surrounded in the exit surface by dielectric gratings. A TM polarized plane wave impinges on the metal surface from above.
Figure 1. Schematic representation of a metal slit surrounded in the exit surface by dielectric gratings. A TM polarized plane wave impinges on the metal surface from above.
Photonics 10 00435 g001
Figure 2. Schematic of the two-slit beaming structure to be optimized. In a Cartesian coordinate system with coordinate ( x , z ) , the structure’s bottom is located at z = 0 , and the origin is at the center of the middle part of the structure.
Figure 2. Schematic of the two-slit beaming structure to be optimized. In a Cartesian coordinate system with coordinate ( x , z ) , the structure’s bottom is located at z = 0 , and the origin is at the center of the middle part of the structure.
Photonics 10 00435 g002
Figure 3. The magnetic field intensity | H y | 2 distributions for optimum single-beaming cases. The beaming angles are (a) θ = 0 , (b) θ = 10 , (c) θ = 20 , and (d) θ = 30 .
Figure 3. The magnetic field intensity | H y | 2 distributions for optimum single-beaming cases. The beaming angles are (a) θ = 0 , (b) θ = 10 , (c) θ = 20 , and (d) θ = 30 .
Photonics 10 00435 g003
Figure 4. Far-field angular spectrum for optimum single-beaming effects. The beaming angles are (a) θ = 0 , (b) θ = 10 , (c) θ = 20 , and (d) θ = 30 .
Figure 4. Far-field angular spectrum for optimum single-beaming effects. The beaming angles are (a) θ = 0 , (b) θ = 10 , (c) θ = 20 , and (d) θ = 30 .
Photonics 10 00435 g004
Figure 5. The magnetic field intensity | H y | 2 distributions for optimum dual-beaming cases. The beaming angles ( θ 1 , θ 2 ) are (a) ( 5 , 5 ) , (b) ( 10 , 10 ) , (c) ( 15 , 5 ) , and (d) ( 10 , 20 ) .
Figure 5. The magnetic field intensity | H y | 2 distributions for optimum dual-beaming cases. The beaming angles ( θ 1 , θ 2 ) are (a) ( 5 , 5 ) , (b) ( 10 , 10 ) , (c) ( 15 , 5 ) , and (d) ( 10 , 20 ) .
Photonics 10 00435 g005
Figure 6. Far-field angular spectrum for optimum dual-beaming effects. The beaming angles are (a) ( 5 , 5 ) , (b) ( 10 , 10 ) , (c) ( 15 , 5 ) , and (d) ( 10 , 20 ) .
Figure 6. Far-field angular spectrum for optimum dual-beaming effects. The beaming angles are (a) ( 5 , 5 ) , (b) ( 10 , 10 ) , (c) ( 15 , 5 ) , and (d) ( 10 , 20 ) .
Photonics 10 00435 g006
Table 1. The optimum structural parameters for single-beaming effect.
Table 1. The optimum structural parameters for single-beaming effect.
θ p 0 p 1 p 2 l 0 l 1 l 2 l 3 d 1 θ i n
   0 0.4270.4000.4000.2030.1570.1570.2030.0500 0
10 0.4180.3590.4390.1320.0270.1300.1710.0346 45 . 54
20 0.3430.3060.5100.2610.1520.1380.0040.0588 0 . 34
30 0.3120.2850.5920.2680.4660.2320.0000.0678 8 . 97
Table 2. The optimum structural parameters for dual-beaming effects.
Table 2. The optimum structural parameters for dual-beaming effects.
( θ 1 , θ 2 ) p 0 p 1 p 2 l 0 l 1 l 2 l 3 d 1 θ i n
( 5 , 5 ) 0.4040.3960.3960.4570.1480.1480.4570.0501 23 . 78
( 10 , 10 ) 0.4460.4510.4510.1480.0980.0980.1480.0498 10 . 37
( 15 , 5 ) 0.3260.4770.4110.4570.1250.1570.2120.0485 26 . 64
( 10 , 20 ) 0.1020.3320.4690.2090.0180.2240.9170.0580 22 . 43
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.

Share and Cite

MDPI and ACS Style

Lu, X. An Efficient Method for Light Beaming from Subwavelength Slits Surrounded by Surface Gratings. Photonics 2023, 10, 435. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics10040435

AMA Style

Lu X. An Efficient Method for Light Beaming from Subwavelength Slits Surrounded by Surface Gratings. Photonics. 2023; 10(4):435. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics10040435

Chicago/Turabian Style

Lu, Xun. 2023. "An Efficient Method for Light Beaming from Subwavelength Slits Surrounded by Surface Gratings" Photonics 10, no. 4: 435. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics10040435

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