1. Introduction
In recent years, infrared laser guidance has employed multi-wavelength lasers working together [
1]. Under the irradiation of multi-wavelength lasers, targets need to reduce their surface reflection, so as to achieve “laser stealth”. At present, there are two main ways to reduce reflection. One is to improve absorption with the help of absorbing materials; the other is to utilize light guide and/or light transmission effects. So far, a variety of laser-stealth materials have been developed. These materials include coated laser-absorbing materials, nano-absorbing materials, semiconductor compounds, and spectral conversion materials [
2]. It is difficult to achieve laser stealth based on light guide and/or light transmission effects. Therefore, great effort has been devoted to research related to absorbing materials.
Recently, metamaterials have opened up new directions for polarization detection, sensing [
3,
4] and infrared laser-absorbing materials. Various metamaterials, such as gradient metasurfaces [
5], typed porous metasurfaces [
6], and circular hole array metasurfaces [
7], have been proposed. Typically, metamaterials consist of large numbers of artificial “atomic” microstructures. Since Pendry et al. [
8,
9] proposed electromagnetic stealth metamaterials and perfect lenses, stealth materials based on metamaterials have emerged one after another. For example, Cheng et al. [
10] presented a non-resonant “I”-shaped metamaterial and a resonant ELC metamaterial, which demonstrate 99% absorption at microwave frequencies. Nie et al. [
11] proposed a two-layer asymmetric split-ring metamaterial that exhibits efficient multi-band absorption. Kadir et al. [
12] achieved ~90% absorption in the entire long-wave infrared (LWIR) by burying two parallel metasurfaces in an amorphous silicon layer. However, there are few studies on laser-absorbing metamaterials. Many metamaterial absorbers suffer from the disadvantages of a single absorption band [
13], polarization sensitivity [
14], and only being suitable for normal incidence [
15].
In particular, 1.55 μm lasers have been applied to achieve precise guidance [
16], which presents new challenges in the study of laser-absorbing metamaterials. This is because 1.55 μm lasers usually have wavelength tuning capabilities. In the face of wavelength tuning, if the metamaterial covering the target surface cannot maintain high absorption within the tuning range of the incident laser, the target will barely survive [
17]. To address this problem, some researchers have investigated the use of Helmholtz resonators [
18]. However, Helmholtz resonators operating at light frequencies are difficult to manufacture.
In the study of solar energy harvesting, researchers have found that two-dimensional gratings can produce high absorption of visible light [
19]. Inspired by this, we design and investigate a dual-band laser absorber based on Si-Ni metasurface. It can efficiently absorb the 1.55 μm and 1.064 μm lasers, providing a way of achieving laser stealth. For a recent 1.55 μm laser with a wavelength tuning range of 60 nm [
20], the metasurface can maintain an absorptivity greater than 93%. In addition, the metasurface has a low emissivity in the mid-infrared range, showing good thermal stealth performance under night conditions [
21]. To understand the underlying mechanism, electromagnetic field distributions, energy loss distributions, and effective impedance are investigated. Moreover, considering the fabrication error and practical applications, the effects of geometric parameters, polarization of incidence, and incident angle on the absorption performance are discussed.
2. Structure Design
Figure 1 shows the proposed metasurface based on Si-Ni microstructures. It is a two-dimensional periodic structure consisting of three functional layers. In each unit cell of the metasurface, there is a continuous silicon film between a silicon square and a nickel layer, as shown in
Figure 1a. The period of the metasurface is
, and the side length of the square is
, as shown in
Figure 1b. The thicknesses of the square, the continuous silicon film, and the nickel film are
,
, and
, respectively, as shown in
Figure 1c. We chose silicon because of its high refractive index, which can improve tolerance to oblique incidence. Nickel is use due to its relatively high loss at the wavelengths of interest, which can enhance absorption. The nickel film must be thick enough to prevent the penetration of incident light. In our structure, it has a thickness of 100 nm (i.e.,
nm). Please note that there is a substrate (not shown) under the nickel film. The substrate (such as a polished silicon wafer) only acts as a support. It has no effect on the performance of the metasurface.
When the metasurface is illuminated by a laser, surface plasmon polaritons (SPPs) could be excited. According to the wavevector matching theory [
22,
23], whether or not to excite the SPP is controlled by:
where
and
are respectively the in-plane wavevector of the SPP and the incident light,
is the grating vector,
m is a positive integer. The thickness of the silicon square is small, so only (+1) and (−1) SPP modes will be excited. Assuming that the polarization of the incident light is along the
-direction, Equation (1) can be expressed as:
where
is the effective refractive index of the SPP mode,
λ is the free-space wavelength, and
α is the incident angle. It is similar if the polarization of the incident light is along the
-direction.
Since silicon has a high refractive index, localized plasmon resonance may also be excited [
24,
25]. According to the impedance matching theory [
26], when the radiative and resistive losses of the eigenmode equal to each other, critical coupling can occur. In this situation, if the incident light excites the eigenmode due to the radiation coupling, the incident energy will be completely converted into the ohmic loss, leading to perfect absorption. The absorption rate
can be expressed as:
where
Z(
λ) is the effective impedance of the metasurface.
To optimize the structural parameters and achieve near-perfect dual-band laser absorption, finite-difference time-domain (FDTD) simulations were carried out. The optical parameters of silicon and nickel were taken from the literature [
27].
Figure 2 plots the absorbance, reflectivity and transmittance of the metasurface under normal incidence. The geometric parameters are as follows:
nm,
nm,
nm,
nm,
nm. As can be seen in
Figure 2, the transmittance is almost zero in the investigated spectral range. The reflectance has two dips centered at wavelengths of 1.064 μm and 1.55 μm, respectively. The absorptivity of the metasurface satisfies
, where
and
are the reflectance and transmittance, respectively. Therefore, the absorptivity is ~97.8% and ~94.4% at 1.064 μm and 1.55 μm wavelengths, respectively. The results indicate that the metasurface can achieve high absorption of dual-band guided laser. Especially for the 1.55 μm band, the metasurface can still maintain >90% absorbance even when the wavelength tuning range of the laser reaches ±50 nm. It implies that our absorber has the ability to counteract the wavelength tuning.
3. Absorption Mechanism
To understand the absorption mechanism, the distributions of electric field
and magnetic field
at the resonance wavelengths were studied. The results are shown in
Figure 3. In
Figure 3a,c, we observe that the electric field is mainly confined in the gap between adjacent squares. The magnetic field mainly distributes in the continuous silicon layer, and is relatively strong below the gap. This is because the free electrons in the nickel layer form an oscillating current. In other words, propagating plasmonic resonance is excited. However, there is a portion of the magnetic field below the square, implying some degree of localized plasmonic resonance. Therefore, the absorption at the 1.064 μm wavelength can be attributed to a mixed mode, which is dominated by propagating plasmonic resonance and supplemented by localized plasmonic resonance. In
Figure 3b,d, it can be seen that the electric field is mainly associated with the square. Meanwhile, the magnetic field is strongly concentrated below the square. These phenomena are consistent with localized plasmonic resonance. However, a small portion of the magnetic field leaks out. Therefore, the absorption at the 1.55 μm wavelength is also a mixed mode, which is dominated by localized plasmonic resonance and supplemented by propagating plasmonic resonance.
Then, the permittivity, permeability, effective refractive index and effective impedance of the metasurface were retrieved [
28], as shown in
Figure 4. In
Figure 4a,c, it can be seen that Im(
) and Im(
) are positive and large at wavelengths of 1.064 μm and 1.55 μm, indicating magnetic resonance. Meanwhile, Im(
) is positive and large at the two wavelengths, indicating high loss. According to
, we can obtain the effective impedance of the metasurface, as shown in
Figure 4d. For the wavelength of 1.064 μm, Re(
) is relatively large and close to 1, and Im(
) is approximately 0. For the wavelength of 1.55 μm, Re(
) is close to 1, and Im(
) is a relative maximum. These indicate that the impedance of the metasurface is matched to the free space at wavelengths of 1.064 μm and 1.55 μm, resulting in reduced reflections. Therefore, at 1.064 μm and 1.55 μm wavelengths, the loss increases and the reflection decreases, leading to an increase in absorption. However, since
does not reach 1, there will be little reflection, and the absorption does not reach 100%.
Because the structure in each unit cell is symmetric, absorption properties will not depend on the polarization of the incidence [
29]. Assuming a TM-polarized (the polarization is along the
-direction) incidence,
Figure 5 shows the energy dissipation in the nickel film. A cross-section close to the upper surface of the nickel film was studied, as shown in
Figure 5a. In
Figure 5b, for the wavelength of 1.064 μm, the energy loss occurs mainly below the gap between adjacent squares. A small portion of energy dissipates below the square. This is due to the oscillation of free electrons in the nickel film, which is consistent with the electromagnetic field distributions in
Figure 3. Similarly, in
Figure 5c, for the wavelength of 1.55 μm, the energy loss occurs mainly below the square. This is also consistent with the electromagnetic field distributions in
Figure 3.
4. Discussion
The absorber is easily fabricated by micro-nanotechnology. The following method can be used to prepare the absorber. First, by means of vacuum evaporation, a layer of nickel and a silicon film can be sequentially coated on a substrate. Then, a layer of E-beam resist is spin-coated on top of the silicon film. Using E-beam lithography, the desired squares can be patterned. After developing and fixing, a sample with periodic square holes will be obtained. Next, an additional layer of silicon is evaporated on the sample. Removing the E-beam resist, we can obtain the proposed absorber.
To analyze the influence of manufacturing errors, the absorption spectra were simulated by adjusting the geometric parameters
,
,
, and
(see
Figure 1). Since the proposed metasurface is polarization independent, the absorption performance is discussed by taking a TM-polarized incidence as an example. The thickness of the nickel film (i.e.,
) is not discussed.
Figure 6 shows the calculated absorption when a single geometric parameter is adjusted with other parameters unchanged. For ease of discussion, the peaks around the wavelength of 1.064 μm are marked as
, whereas the peaks around the wavelength of 1.55 μm are marked as
.
Figure 6a shows the absorption spectra as a function of
. A change in
implies a change in the gap, which affects the near-field coupling between adjacent squares. When
756 nm, the absorption at
decreases due to the increased in the duty cycle of the top grating. Owing to the small gap, the near-field coupling becomes strong, and the absorption peaks degenerate. In this situation, the impedance is not perfectly matched to the free space. When
decreases from 956 nm to 856 nm,
does not change significantly, whereas the absorption at
increases significantly. It is indicated that an appropriate period will enhance the resonance.
Figure 6b shows the absorption spectra as a function of
. It can be observed that changing
has little effect on the absorption at
. When
increases from 50 nm to 70 nm,
exhibits a slight red shift, indicating that increasing
can tune the position of
to a certain extent.
Figure 6c shows the absorption spectra as a function of
. When
26 nm, the two peaks are blue-shifted to 0.96 μm and 1.46 μm, respectively. This is because the coupling between the upper and lower layers is enhanced. Meanwhile, the impedance is not well matched to the free space. When
36 nm, the absorption is high. The reason for this is that the coupling between the upper and lower layers is weakened (the impedance matching is better). However, when
= 46 nm, the two peaks are red-shifted to 1.31 μm and 1.66 μm, respectively. In this situation, high absorption at the 1.064 μm wavelength cannot be obtained.
Figure 6d shows the absorption spectra as a function of
. As
increases, the absorption at
increases, whereas the absorption at
decreases. When
600 nm, the two peaks tend to degenerate. This is because the two resonances tend to influence each other as
increases.
In practical applications, the response to the angle of incidence needs to be considered.
Figure 7 shows the spectral response of the metasurface at different angles of incidence. As the incident angle increases, the absorption rate decreases to a certain extent. For the peak at the wavelength of 1.064 μm, its absorption rate remains ≥85% at incident angles up to 60 degrees. For the peak at the wavelength of 1.55 μm, its spectral location exhibits a blue shift at oblique incidence. When the incident angle is as high as 60 degrees, the absorption at the 1.55 μm wavelength is still ≥75%. The results show that the metasurface is tolerant to the incident angle. In particular, at large incident angles, if the 1.55 μm laser performs wavelength tuning (±50 nm), the metasurface can still have an absorptivity of ≥70%.
Moreover, the response to the polarization of incidence is also an important property of the metasurface.
Figure 8 illustrates the absorption spectra at different polarization angles. In
Figure 8, “0°” represents that the polarization of the incidence is along the
-direction; and “90°” represents that the polarization of the incidence is along the
-direction. As can be observed, when the polarization varies, absorption spectra of the metasurface remain almost unchanged. These results suggest that the metasurface is polarization insensitive.
It is well known that passive infrared detection (detecting the thermal radiation of a target) is a common and important detection method under night conditions.
Figure 9 plots the reflection of the metasurface in the mid-infrared. It is obvious that, in addition to high absorption of the 1.064 μm and 1.55 μm lasers, the reflectivity is ≥80% at wavelengths greater than 2.8 μm. As mentioned above, the thick nickel layer cannot transmit (infrared) light. Therefore, the observed high reflectivity is associated with low absorptivity for thermal infrared. According to Kirchhoff’s law, the low absorptivity means a low thermal emissivity. For a given temperature, the thermal infrared radiation intensity of a target coated with our metasurface is relatively small. Therefore, under night conditions, the proposed metasurface is useful for circumventing passive infrared detection.