Home High accuracy inverse design of reconfigurable metasurfaces with transmission-reflection-integrated achromatic functionalities
Article Open Access

High accuracy inverse design of reconfigurable metasurfaces with transmission-reflection-integrated achromatic functionalities

  • Xiao-Qiang Jiang ORCID logo , Wen-Hui Fan ORCID logo EMAIL logo , Xu Chen ORCID logo , Lv-Rong Zhao , Chong Qin , Hui Yan , Qi Wu and Pei Ju
Published/Copyright: March 25, 2025
Become an author with De Gruyter Brill

Abstract

Artificial intelligence algorithms based on deep neural network (DNN) have become an effective tool for conceiving metasurfaces recently. However, the complex and sharp resonances of metasurfaces will tremendously increase the training difficulty of DNNs with non-negligible prediction errors, which hinders their development in designing multifunctional metasurfaces. To overcome the obstacles, the interaction mechanisms between meta-atoms and terahertz (THz) waves via multipole decomposition are investigated to establish a high-quality dataset, which can decrease the complexity of DNN and improve the prediction accuracy. Meanwhile, transfer learning is also employed to reduce the large quantity of training data required by the DNN. Accordingly, two broadband and transmission-reflection-integrated reconfigurable metasurfaces for focused vortex beam generation are inversely designed by counter propagating the DNN with fraction error less than 10−4. The results indicate that transmission-reflection-integrated achromatic performances are well achieved in the frequency range of 0.7–1.3 THz, which have the average focusing efficiency and mode purity higher than 48 % and 92 %, respectively. Moreover, transmission-reflection-integrated achromatic THz imaging and edge detection can also be realized by the metasurfaces. This work provides a high accuracy inverse design method for conceiving multifunctional meta-devices, which may promise further progress for the on-chip THz imaging systems.

1 Introduction

Terahertz (THz) waves, with frequency range of 0.1–10 THz in the electromagnetic (EM) spectrum, have demonstrated numerous promising applications that can greatly reshape future researches and industrial fields [1], [2], [3]. As for the further exploitation of THz waves, it can be realized by manipulating their wavefronts to obtain desired special beams [4], [5]. Therefore, THz waves carrying additional EM characteristics can greatly extend their application boundaries. For instance, THz vortex beams carrying orbital angular momentum (OAM) make them extraordinarily attractive in terabit wireless communications due to the mutually orthogonal eigenstates of OAM [6], [7]. Besides, the helical phase-fronts and hollow properties of THz vortex beams can realize super-resolution imaging and edge detection [8], [9]. Specifically, the spiral phase contrast (SPC) imaging via OAM is of great importance in edge detection since it can reconstruct the edge information for both amplitude and phase objects, which may have great potential in machine vision and remote sensing [10], [11]. Until now, plentiful methods have been implemented for generating optical vortex beams, including spiral phase plates and spatial light modulators [7], [12]. However, these conventional devices confront great challenges in THz range since natural materials barely have EM responses to THz waves, which severely impedes the further development of THz technology.

Due to the powerful capabilities for manipulating EM waves, metasurfaces may have considerable potential to realize revolutionary advancements that overcome the constraints of natural materials on the development of THz technology [13], [14], [15]. By employing metasurfaces, the generation of THz vortex beams and their applications have been developed rapidly [5], [16]. However, these devices still confront limited function, which can only be operated in either transmission mode (T-mode) or reflection mode (R-mode), greatly wasting space resources [17]. To realize broadband achromatic focusing in both T-mode and R-mode, multi-layer metallic metasurfaces with frequency- and polarization-multiplexing capabilities were proposed [18], [19]. Nevertheless, the focusing efficiencies of these transmission-reflection-integrated metasurfaces are only about 20 % due to the complex structures and strict operating requirements. Although monolayer dielectric metasurfaces were conducted for transmission-reflection-integrated optical imaging, they still confront inefficiency and limited bandwidth in the implementations [20], [21]. In contrast with passive counterparts, reconfigurable metasurfaces with dynamic properties can enable more flourishing performances, becoming a major frontier in further development of metasurfaces [22], [23]. Since vanadium dioxide (VO2) can span five orders of magnitude changes in its conductivity during the insulator-to-metal transition (IMT) [24], the transmission-reflection-integrated operation can be accomplished by employing VO2 instead of complicated multi-layer structures.

In order to accomplish transmission-reflection-integrated achromatic metasurfaces, each meta-atom is expected to satisfy the phase requirements in both T-mode and R-mode with relatively high efficiency. The complicated demands and in-depth connections are difficult to be recognized by limited human brainpower, which becomes a common issue in engineering sophisticated meta-devices [25]. Recently, deep learning has been employed to solve the bottlenecks [26], [27]. However, current deep neural networks (DNNs) are primarily focused on amplitude-modulation metasurfaces due to their relatively simple EM responses [28], [29], [30]. Since both amplitude and phase responses should be considered in the design of phase-modulation metasurfaces, their training difficulty is much harder than that of amplitude-modulation metasurfaces [31]. Besides, the complicated resonances will also cause non-negligible prediction errors. For instance, previous research for phase prediction suffered from test error up to 16° due to the sharp peaks [32]. This defect is still one of the most severe obstacles to design phase-modulation metasurfaces by DNNs [33], [34]. Although several methods were proposed by topological optimization or other advanced algorithms, their performances are highly depended on algorithm itself and require a large amount of training data and testing time. More importantly, the EM responses of meta-atoms barely receive enough attention, which are the foundation for designing metasurfaces. In general, the inverse design of broadband metasurfaces with advanced functionalities is still a huge challenge due to the complex spectral responses and phase requirements.

In this work, a high accuracy inverse design method is proposed to accomplish transmission-reflection-integrated achromatic reconfigurable metasurface (TRARM) for the first time to our knowledge. To satisfy complicated phase requirements, DNN and transfer learning are employed in the design scheme. Moreover, the fundamental mechanisms of interactions between meta-atoms and THz waves are investigated via multipole decomposition to exclude undesirable resonances. Thus, a high-quality dataset can be established to decrease the complexity of DNN and improve its prediction accuracy. To validate the inverse design method, two focused vortex beam (FVB) generators with different parameters are constructed and marked as TRARM-Ⅰ and TRARM-Ⅱ, respectively. Numerical results indicate that chromatic aberrations in both R-mode (0.7–1.0 THz) and T-mode (1.0–1.3 THz) can be well corrected. The maximum deviation ratios of TRARM-Ⅰ are 5.14 % (R-mode) and 4.30 % (T-mode), while that of TRARM-Ⅱ are 7.73 % (R-mode) and 4.32 % (T-mode). The average focusing efficiencies of these two metasurfaces are all higher than 48 %, and the OAM mode purity of TRARM-Ⅱ is higher than 92 %. Moreover, transmission-reflection-integrated achromatic THz imaging and edge detection can be also achieved by the metasurfaces. The presented high accuracy inverse design method has great significance in conceiving multifunctional metasurfaces.

2 Design principle for TRARM

The operation schematics of TRARM are illustrated in Figure 1(a), which can generate FVB and correct chromatic aberration in both T-mode (insulator phase, I-VO2) and R-mode (metallic phase, M-VO2), and it is a reversible transition. The polarization state of incident THz waves are left-circularly polarization (LCP), and that of generated FVBs are right-circularly polarization (RCP). The meta-atoms with lattice constant p, VO2 thickness h v , pillar height h, length a, and width b are shown in Figure 1(b), which consists of a high-resistance silicon (Si) block on silicon dioxide (SiO2) hexagonal substrate and VO2 is introduced as a bottom plate. The design principle of TRARM is shown in Figure 1(c), composed of three phase profiles of hyperbolic phase, spiral phase, and propagation phase, which are responsible for converging THz waves, generating vortex beam, and correcting chromatic aberration, respectively.

Figure 1: 
Schematics of (a) TRARM, (b) meta-atoms and (c) phase profiles. The yellow and green circles represent the polarization states of incident and outgoing THz waves, respectively.
Figure 1:

Schematics of (a) TRARM, (b) meta-atoms and (c) phase profiles. The yellow and green circles represent the polarization states of incident and outgoing THz waves, respectively.

Benefitting from extraordinary capabilities of metasurfaces in manipulating EM waves, the hyperbolic phase and spiral phase can be super-imposed on a monolayer metasurface [35]:

(1) φb(x,y,f)=2πfc(x2+y2+Fl2Fl)+larctan(yx)φb(x,y,f)=2πfc(x2+y2+Fl2Fl)+larctan(yx)

where f represents the operating frequency, c denotes the light speed in vacuum, l is the topological charge of OAM, and F l is the focal length. To correct chromatic aberration in the broadband frequency range, an additional phase compensation should be introduced [36]. Nevertheless, current broadband design method is employed to compensate the phase difference between f and f min (minimum frequency) [37], which is not applicable to design the transmission-reflection-integrated achromatic metasurfaces. Therefore, the design method should be further improved for designing TRARM, and the modified phase profiles are expressed as:

(2) φr(x,y,fr)=φbr(x,y,fmaxr)+Δφr(x,y,fr)φr(x,y,fr)=φbr(x,y,fmaxr)+Δφr(x,y,fr)

(3) φt(x,y,ft)=φbt(x,y,fmint)+Δφt(x,y,ft)φt(x,y,ft)=φbt(x,y,fmint)+Δφt(x,y,ft)

where the footnote max and min are maximum and minimum frequency, and the footnote t and r represent T-mode and R-mode of TRARM. And the frequency range of f r and f t are f r ∈ {f minr, f maxr} and f t ∈ {f mint, f maxt}, respectively. Therefore, the phase compensation Δφ can be described as:

(4) Δφr(x,y,fr)=[2πc(x2+y2+F2lrFlr)](fmaxrfr)Δφr(x,y,fr)=[2πc(x2+y2+F2lrFlr)](fmaxrfr)

(5) Δφt(x,y,ft)=[2πc(x2+y2+Flt2Flt)](ftfmint)Δφt(x,y,ft)=[2πc(x2+y2+Flt2Flt)](ftfmint)

In Eqs. (2) and (3), φ b is considered as basic phase for THz FVB generation and it can be satisfied by Pancharatnam–Berry (PB) phase for both modes. Meanwhile, Δφ is responsible for correcting chromatic aberration by exploiting propagation phase. The illustration of phase compensation for T-mode and R-mode is depicted in Figure 2(a), where f maxr and f mint are both set as 1 THz with identical focal length. And the propagation phase is employed to compensate the phase differences between f r (f t ) and f maxr (f mint) for R-mode (T-mode), respectively.

Figure 2: 
Phase requirements. (a) Illustration of phase compensation. Phase compensations for TRARM-Ⅰ in (b) R-mode and (c) T-mode. Phase compensations for TRARM-Ⅱ in (d) R-mode and (e) T-mode.
Figure 2:

Phase requirements. (a) Illustration of phase compensation. Phase compensations for TRARM-Ⅰ in (b) R-mode and (c) T-mode. Phase compensations for TRARM-Ⅱ in (d) R-mode and (e) T-mode.

In this work, the focal length and diameter of TRARM-Ⅰ are set as F lt-Ⅰ = F lr-Ⅰ = 2.80 mm and D 1 = 3.38 mm, and that of TRARM-Ⅱ are F lt-Ⅱ = F lr- = 1.80 mm and D 2 = 2.80 mm. Meanwhile, the operating frequency ranges are 0.7–1.0 THz for R-mode and 1.0–1.3 THz for T-mode, respectively. The selected parameters and frequency range are explained in Supplementary Section 1. Accordingly, the phase compensations for TRARM-Ⅰ and TRARM-Ⅱ are depicted in Figure 2(b)–(e). Since f maxr and f mint are both set as 1 THz with identical focal length, the meta-atoms can share the same rotation angles at each position. Meanwhile, the phase compensation for TRARM in T-mode and R-mode should also be satisfied simultaneously by the meta-atoms, which is one of the major difficulties for designing TRARM. Fortunately, this obstruction can be readily solved by employing the proposed inverse design method.

3 High accuracy inverse design method

To resolve the defects of current DNNs in designing multifunctional phase-modulation metasurfaces, a high accuracy inverse design method via physical analyses, deep learning, and transfer learning is presented. To begin with, the multipole decomposition is performed to investigate the fundamental mechanisms of interactions between meta-atoms and THz waves, which can provide a high-quality dataset to significantly improve the prediction accuracy. Subsequently, DNN and transfer learning are both employed to predict the EM responses of meta-atoms and inversely design the TRARM.

3.1 Investigation on resonance mechanism

According to Eqs. (2)(5), the phase responses of meta-atoms should be linearly proportional to operating frequency, indicating that resonances are not allowed to appear in the EM spectra. Moreover, the basic phase profile for R-mode is set at maximum frequency of 1.0 THz, requiring a large range of phase compensation to cover the full 2π phase range at minimum frequency of 0.7 THz. Therefore, the p and h of meta-atoms should be as large as possible to provide sufficient phase compensation. However, meta-atoms with large p and h can also induce undesired resonances to deteriorate their performances [38]. Therefore, the underlying physical mechanisms of EM responses from meta-atoms should be investigated to delicately design the p and h, which are expected to satisfy the phase requirements and avoid undesired resonances in the operating frequency range.

Since identical meta-atoms will have larger effective refractive index if they operate in a higher frequency range [39], the undesired resonances for R-mode (0.7–1.0 THz) will not appear if the amplitude and phase spectra for T-mode (1.0–1.3 THz) are smooth curves. To avoid undesired resonances, the multipole decomposition (Supplementary Section 2) of meta-atoms operating in T-mode under different lattice constant is thoroughly investigated, which have geometric parameters of a = 72 μm, b = 25 μm, h v = 5 μm, and h = 245 μm. The polarization conversion ratio (PCR) and phase spectra varied with lattice constant of meta-atom are depicted in Figure 3(a)–(c) by red and blue curves, respectively. It can be clearly observed that sharp resonances are excited and they will gradually appear in the concerned frequency range (yellow parts) as the lattice constant increases. These sharp resonances are marked as points A1 (1.38 THz), B1 (1.49 THz), C1 (1.24 THz), D1 (1.38 THz), E1 (1.17 THz), F1 (1.28 THz), and G1 (1.37 THz), which are also painted with blue color in Figure 3(d)–(f). The sharp resonances will not only deteriorate the PCR efficiency but also introduce abrupt phase shift near the resonances according to Kramers–Kronig relation [40], which are verified by the points A2, B2, C2, D2, E2, F2, and G2. The inferior PCR efficiency will eventually decrease the focusing efficiency of constructed metasurface, and the abrupt phase shifts inevitably make horrible impacts on broadband operation. Moreover, transmission and reflection amplitude of the meta-atoms are also calculated for comprehensive understanding, as depicted in Supplementary Section 2.

Figure 3: 
Optical responses of meta-atoms. (a)–(c) PCR efficiency and phase responses of meta-atoms. (d)–(f) Scattering power of multipole decomposition. (g)–(m) Magnetic field distributions. Black boxes are the outline of Si blocks, purple circles are the rotation of magnetic fields, and the black dots and forks are the outward and inward plane.
Figure 3:

Optical responses of meta-atoms. (a)–(c) PCR efficiency and phase responses of meta-atoms. (d)–(f) Scattering power of multipole decomposition. (g)–(m) Magnetic field distributions. Black boxes are the outline of Si blocks, purple circles are the rotation of magnetic fields, and the black dots and forks are the outward and inward plane.

As depicted in Figure 3(d)–(f), the sharp peaks from A1, B1, C1, D1, E1, and F1 are primarily attributed to toroidal dipole (TD) resonances, and magnetic dipole (MD) resonance is responsible for G1. For instance, there are two magnetic fields with opposite circulating direction in yz plane at A1, as depicted in left part of Figure 3(g). Therefore, the circulating magnetic dipoles are excited in xy plane and eventually induce TD resonance, which accords with the theoretical results of multipole decomposition. Moreover, it can be noticed that the scattering power of magnetic quadrupole is also increased around A1, as shown in Figure 3(d). And its magnetic field is depicted in right part of Figure 3(g), which is a parasitical response from coherent oscillation between TD resonance and suppressed electric dipole resonance [41], [42]. From Figure 3(h)–(l), the circulating magnetic fields are appeared in xy plane for B1, C1, and E1, while the circulating magnetic fields of D1 and F1 are excited in xz plane, and those circulating magnetic fields are responsible for TD resonances similar to A1. Different from the others, G1 is induced by MD resonance as shown in Figure 3(f), which is verified by the magnetic fields depicted in Figure 3(m). Furthermore, the h of meta-atom also has evident impacts to its performances, as depicted in Supplementary Section 3. According to the investigation, the lattice constant and structural height of meta-atom are eventually set as p=75/3p=75/3 μm and h = 235 μm to suppress undesired resonances, which are thoroughly explained in Supplementary Section 4.

3.2 DNN construction and inverse design of TRARM

Firstly, T-mode operating meta-atoms are trained to learn the connections between meta-atoms and their EM responses. Subsequently, transfer learning, permitting trained modes to share their acquired knowledge and experience to improve the performances in other similar scenarios [43], [44], is introduced to decrease the samples from 16,915 data down to 5,039 data. It spent about one month on collecting 16,915 data for T-mode training, while it only took about one week on collecting 5,039 data for R-mode training. As depicted in Figure 4, forward prediction network (FPN) consists of six fully-connected hidden-layers with 200, 500, 1,000, 1,000, 300, and 61 neurons, respectively, and each output layer passes through a ReLU activation function before it is sent to next layer. Besides, the learning rate is set as 10−4, and the batch size is 256. Moreover, the weight decay is employed to alleviate the overfitting and ensure the accuracy of the network for test dataset [45].

Figure 4: 
The illustration of design method. Yellow circles are trained neurons, gray circles are transferred neurons, and green circles represent the neurons operating in R-mode.
Figure 4:

The illustration of design method. Yellow circles are trained neurons, gray circles are transferred neurons, and green circles represent the neurons operating in R-mode.

It should point out that the PCR efficiencies from established dataset (T-mode meta-atoms) range from 0.23 to 0.99 and the phase responses range from −400.65° to 179.99°, which are depicted in Supplementary Section 5. Therefore, the prediction of PCR efficiency by FPN is much easier than that of phase responses, and it can achieve high accuracy with only 1,000 iterations. The mean square error (MSE) of the PCR efficiency predicted by FPN decreases down to 7.41 × 10−9 with around 100 s for the whole training process, as shown in Figure 5(a). Moreover, the transmitted phase response predicted by FPN is set as 24,000 iterations since the large range of variation inevitably requires more iterations to acquire relatively high accuracy, which has the MSE of 2.46 × 10−1 and it takes about 65 min for the whole training process, as shown in Figure 5(b). After the direct learning, transfer learning is employed to copy several layers as initialized weights and biases for the training of reflected PCR efficiency and phase response, which have MSE of 8.68 × 10−8 (500 iterations) and 3.20 × 10−1 (10,000 iterations), respectively, as depicted in Figure 5(c)–(d). The MSE of reflected responses of various transferred layers are also investigated, as shown in Figure 5(e)–(f). Similarly, reflected PCR efficiencies are much easier for prediction than reflected phase responses due to the small variation range.

Figure 5: 
Performances of DNN and transfer learning. (a)–(b) Learning curves of T-mode by direct learning. (c)–(d) Learning curves of R-mode by transfer learning. (e)–(f) MSE of reflected responses of various transferred layers. Red circles are optimal MSE for transfer learning, and blue circles are MSE for direct learning.
Figure 5:

Performances of DNN and transfer learning. (a)–(b) Learning curves of T-mode by direct learning. (c)–(d) Learning curves of R-mode by transfer learning. (e)–(f) MSE of reflected responses of various transferred layers. Red circles are optimal MSE for transfer learning, and blue circles are MSE for direct learning.

Moreover, the fractional errors of transmitted PCR efficiency and phase response by direct learning are 5.47 × 10−6 and 3.08 × 10−4, respectively. And that of reflected PCR efficiency and phase response employing transfer learning are 2.24 × 10−5 and 2.34 × 10−4, respectively. Such high accuracy demonstrates that the fundamental physical investigations are extremely essential to provide high-quality dataset for DNN training and it can offer a solid foundation to inversely design multifunctional metasurfaces. The comparison between predicted results by FPN and simulated results from finite-integrated method (FIM) are shown in Supplementary Section 6, indicating that the EM responses of various meta-atoms can be precisely predicted by FPN, and it only take a few milliseconds for each meta-atom.

Eventually, the trained FPN can be regarded as an EM simulator to replace conventional simulation methods. In the inverse design, a closed loop design network is employed to satisfy the various phase requirements [31], and its network schematic as well as the training curve of inverse network are depicted in Supplementary Section 7. The inverse design of TRARM can be considered as several steps described as follows. (Ⅰ) The objective metasurface is decomposed into multiple pixels following design principle discussed in Section 2. (Ⅱ) For each pixel corresponding to a target phase, the trained FPN can be activated to predict the EM responses of the meta-atoms, and their parameters are not restricted to the established datasets. (Ⅲ) The output parameters (a and b) obtained by the counter-propagation of DNN and multi-regression will be fed back to the FPN. It is supposed to provide a new design to further minimize the divergences between predicted results and target phase profiles. (Ⅳ) The previous step will be constantly implemented until the minimum error is achieved, and then output a and b. Moreover, the meta-atom with higher PCR efficiency will be selected if there are several structures with identical phase profile. (V) The designed metasurface can be eventually acquired by properly arranging the output meta-atoms.

It is also important to note that the proposed high accuracy inverse design method can also be employed to efficiently design the metasurfaces with similar functionalities operating in other frequency range (e.g., visible, infrared, and so on) by appropriately engineering the feature size of corresponding meta-atoms. Moreover, the established high-quality dataset can be regarded as training data and reused in similar cases.

4 Results and discussion of TRARM

4.1 The performances of TRARM-Ⅰ

The phase profiles and configuration of TRARM-Ⅰ, which can be well compensated by the output meta-atoms from the DNN, are shown in Supplementary Section 8. The intensity profiles of reflected and transmitted THz waves in propagation plane of TRARM-Ⅰ are depicted in Figure 6(a) and (b). It is clear that the focal lengths of TRARM-Ⅰ are barely deviated in both modes, indicating that TRARM-Ⅰ is capable of eliminating chromatic aberration in broadband range of 0.7–1.0 THz (R-mode) and 1.0–1.3 THz (T-mode). Moreover, the focal lengths are plotted in Figure 6(c), where the blue and red lines correspond to TRARM-Ⅰ under working temperature of 23 °C and 80 °C, respectively. The average focal length (F la ) of TRARM-Ⅰ in R-mode is 2.39 mm and that of T-mode is 2.78 mm, indicating that their corresponding numerical aperture (NA) are 0.58 and 0.52, respectively. The maximum deviation ratios of focal length are 5.14 % (R-mode) and 4.30 % (T-mode). The maximum deviation ratio is described as 100 % × (F lmaxF lmin)/F la , where F lmax and F lmin indicate the maximum and minimum values of focal length, respectively.

Figure 6: 
Achromatic performances of TRARM-Ⅰ. (a)–(b) Intensities of TRARM-Ⅰ for R-mode and T-mode in x–z plane. (c) Corresponding focal length. (d)–(e) Focal plane of TRARM-Ⅰ for R-mode and T-mode. (f)–(g) Corresponding values of FWHM and focusing efficiency.
Figure 6:

Achromatic performances of TRARM-Ⅰ. (a)–(b) Intensities of TRARM-Ⅰ for R-mode and T-mode in xz plane. (c) Corresponding focal length. (d)–(e) Focal plane of TRARM-Ⅰ for R-mode and T-mode. (f)–(g) Corresponding values of FWHM and focusing efficiency.

To comprehensively characterize the performances of converged THz waves, the point spread functions (PSFs) of TRARM-Ⅰ are also calculated and depicted in Figure 6(d) and (e). The bright focal spots without any diffusions can be evidently observed in both modes, which further confirms the transmission-reflection-integrated achromatic functionalities of TRARM-Ⅰ. Their intensities are exhibiting a typical Gaussian-type distributions along x-direction, as depicted in Supplementary Section 9. The focusing efficiency, which is one of the most crucial criteria to evaluate the performances of achromatic metasurfaces, is also considered and it is defined as 100 % × I f /I in, where I f and I in represent optical power of focal spot (with a radius of 3 × FWHM, Full Width at Half Maximum) and incident THz waves, respectively [46]. The values of FWHM under R-mode and T-mode are depicted in Figure 6(f). As shown in Figure 6(g), the focusing efficiency of R-mode ranges from 39.27 % at 1.00 THz to 61.47 % at 0.85 THz (blue dotted-line), and that of TRARM-Ⅰ in T-mode ranges from 48.87 % at 1.30 THz to 75.52 % at 1.03 THz (red dotted-line). And the average focusing efficiencies of R-mode and T-mode are 48.79 % and 66.45 %, respectively.

The convolution calculations are also performed to demonstrate the capabilities of transmission-reflection-integrated THz imaging based on TRARM-Ⅰ, described as o(x, y) = g(x, y) ⊗ p(x, y), where ⊗ indicates the convolution operation, o(x, y) is the output function at imaging plane, g(x, y) is the complex amplitude of the object, and p(x, y) represents the PSF [10]. In this case, the objects for R-mode and T-mode are chosen as the logo of our institution (grayscale, 1,596 × 1,680 pixels) and the character of “THz” (522 × 522 pixels), respectively. The calculation results are shown in Figure 7, the transmission-reflection-integrated achromatic THz imaging can be well accomplished by TRARM-Ⅰ and the dispersions are barely observed from imaged picture. Moreover, their peak signal-to-noise ratio are discussed in Supplementary Section 10.

Figure 7: 
The transmission-reflection-integrated achromatic THz imaging by TRARM-Ⅰ.
Figure 7:

The transmission-reflection-integrated achromatic THz imaging by TRARM-Ⅰ.

4.2 The performances of TRARM-Ⅱ

The phase profiles and configuration of TRARM-Ⅱ are depicted in Supplementary Section 11, and it can be perfectly matched by output meta-atoms from the DNN. It can be clearly observed that the focal length of all the sampled frequency points ranging from 0.7–1.0 THz (R-mode) to 1.0–1.3 THz (T-mode) are nearly constant in Figure 8(a) and (b), implying that TRARM-Ⅱ can correct chromatic aberration in the concerned frequency range and generate FVB in both T-mode and R-mode. To quantitatively elucidate the achromatic performances, the focal lengths of TRARM-Ⅱ are plotted in Figure 8(c). The F la of TRARM-Ⅱ in R-mode is 1.54 mm and that of T-mode is 1.79 mm, indicating that the NA of two operation modes is 0.67 and 0.62, respectively. And the maximum deviation ratios of focal length are 7.73 % (R-mode) and 4.32 % (T-mode), respectively.

Figure 8: 
Achromatic performances of TRARM-Ⅱ. (a)–(b) Intensities of TRARM-Ⅱ for R-mode and T-mode in x–z plane. (c) Corresponding focal length. (d)–(e) Focal plane of TRARM-Ⅱ for R-mode and T-mode. (f)–(g) Corresponding values of FWHM and focusing efficiency.
Figure 8:

Achromatic performances of TRARM-Ⅱ. (a)–(b) Intensities of TRARM-Ⅱ for R-mode and T-mode in xz plane. (c) Corresponding focal length. (d)–(e) Focal plane of TRARM-Ⅱ for R-mode and T-mode. (f)–(g) Corresponding values of FWHM and focusing efficiency.

After the exploration of focal length, intensity distributions of focal plane can be subsequently acquired and depicted in Figure 8(d) and (e). Different from the Gaussian-type focal points from TRARM-Ⅰ, the shapes of focal points from TRARM-Ⅱ are similar to doughnuts. The doughnut-shaped distributions manifest the existence of FVB accompanied with a null-amplitude region at focal plane center, which can be also verified by their normalized intensity distributions along x-direction (Supplementary Section 12). In addition, the focusing efficiency of generated FVB are also calculated for evaluation, defined as the ratio of optical power inside doughnut ring to that of incident waves [47]. Figure 8(f) shows the radii of doughnut rings of R-mode and T-mode. As shown in Figure 8(g), the focusing efficiencies of TRARM-Ⅱ in R-mode range from 40.18 % at 1.00 THz to 62.26 % at 0.91 THz (blue dotted-line), and that of T-mode range from 36.89 % at 1.30 THz to 71.34 % at 1.09 THz (red dotted-line). Their average focusing efficiencies are 55.14 % (R-mode) and 56.60 % (T-mode), respectively.

Moreover, the phase profiles of focal points are shown in Figure 9(a) and (b). It can be clearly observed that the phase distributions of focal points from TRARM-Ⅱ are all carrying an OAM with l = −2 in both R-mode (0.7–1.0 THz) and T-mode (1.0–1.3 THz) according to the spiral branches. Moreover, the calculated values of mode purity are all higher than 92 % at l = −2 with average value of 97.02 %. Although there are several parasitic phase noises appeared at l = 0 or l = +2, they can be neglected compared to the operation mode of l = −2.

Figure 9: 
Phase distributions and mode purity of TRARM-Ⅱ at focal point.
Figure 9:

Phase distributions and mode purity of TRARM-Ⅱ at focal point.

As one of the most promising optical methods for edge detection, SPC is a 4-f imaging system employing a vortex beam plate to filter the spatial-frequency signals [11]. Here, the transmission-reflection-integrated edge detections performed by TRARM-Ⅱ associated with l = −2 are investigated via convolution operation [48], which can facilitate the integration and miniaturization of SPC system. Same objects are employed for calculation. As shown in Figure 10, the edge information of the input images can be well accomplished in both R-mode and T-mode. Their root-mean-square errors are discussed in Supplementary Section 13. Results indicate that edge information from input images can be well recovered by TRARM-Ⅱ.

Figure 10: 
The transmission-reflection-integrated achromatic THz imaging by TRARM-Ⅱ.
Figure 10:

The transmission-reflection-integrated achromatic THz imaging by TRARM-Ⅱ.

It is worth to point out that the TRARM is composed of three phase profiles that greatly increase the design difficulty. And the improved phase design method discussed in Eqs. (2)(5) requires excellent optical responses of meta-atoms to provide sufficient phase compensation with relatively high reflection (0.7–1.0 THz) and transmission (1.0–1.3 THz). Therefore, DNN is employed to recognize the in-depth connections between geometric parameters and optical responses of metasurfaces, and eventually achieve transmission-reflection-integrated achromatic performances. Moreover, it can be noted that the focusing efficiency of the T-mode is higher than that of the R-mode for both TRARM-Ⅰ and TRARM-Ⅱ, which can be attributed to different optical losses of VO2 under different states (metallic phase, M-VO2; and insulator phase, I-VO2). In specific, the M-VO2 (80 °C) has relative dielectric constant of ε = –314.47 + i316.25 at 1.0 THz, while relative dielectric constant of I-VO2 (23 °C) is ε = 9 + i8.96 at 1.0 THz. It can be clearly observed that the imaginary part of M-VO2 is much higher than that of I-VO2, indicating M-VO2 has stronger absorption of THz waves. Therefore, the focusing efficiency of the T-mode is higher than that of the R-mode in both TRARM-Ⅰ and TRARM-Ⅱ.

In general, different from the current state-of-the-art techniques that highly depend on algorithm itself and require a large amount of training data as well as testing time, the proposed high accuracy inverse design method in this paper concentrate more on the interaction mechanisms between meta-atoms and incident waves. The physical analyses can provide a high-quality dataset to dramatically decrease the complexity of DNN and significantly improve the prediction accuracy. Therefore, a simple deep learning network with relatively small amount data is enough for the design of advanced meta-devices instead of complicated algorithms. To elucidate this, the clear comparison between state-of-the-art techniques for the inverse design of phase-modulation metasurfaces and the unique benefits of our approach are discussed in Supplementary Section 14.

5 Conclusions

In summary, an efficient scheme with high accuracy is proposed to conceive TRARM, which can accomplish achromatic performances while generating FVB in both T-mode and R-mode. The interaction mechanisms of meta-atoms and THz waves are investigated to obtain high-quality dataset. Moreover, transfer learning is employed in design scheme to reduce the quantity requirements of training data by DNN. The training results indicate that the fractional errors of predicted results are less than 10−4 for predicting PCR efficiencies and phase responses from FPN. Based on the accurate design scheme and improved design principles, two transmission-reflection-integrated achromatic metasurfaces with different parameters are constructed accordingly. Furthermore, transmission-reflection-integrated achromatic THz imaging and edge detection are also accomplished by the metasurfaces via convolution operation. The results demonstrate that the physical analyses are extremely crucial for the inverse design method, which are able to tremendously simplify the DNNs and greatly improve the prediction accuracy. This work provides a high accuracy inverse design method for conceiving multifunctional meta-devices, which may promise further progress for on-chip THz imaging systems.

6 Methods

6.1 Calculation of meta-atoms

The FIM is implemented for simulation, and the periodic boundaries are employed in transversal directions and open boundaries are considered in longitudinal axis. The refractive index of Si and SiO2 are n Si = 3.48 and n SiO2 = 1.94 [49], [50]. Moreover, the temperature-dependent dielectric permittivity of VO2 is discussed in Supplementary Section 15. Moreover, the related references on VO2-based tunable THz metasurfaces have been discussed and compared in Table S4.

6.2 Data collection

After the theoretical consideration of p and h of meta-atoms, a dataset of T-mode operating meta-atoms is collected by employing FIM under LCP incident THz waves ranging from 1.0 THz to 1.3 THz with sampling frequency points of N = 61. The spectral responses of each meta-atom are determined by a and b, and their sweeping ranges are a ∈ [46.0 μm, 74.0 μm] and b ∈ [10.0 μm, 42.0 μm] with a total of 16,915 data (including 80 % training data and 20 % testing data). Since the substrate is a hexagonal structure, the meta-atoms should abide by certain order of 3ab75b<1,406.253ab75b<1,406.25 . The selected sweeping range and restricted condition are explained in Supplementary Section 16. In general, the FPN for T-mode meta-atoms has two input parameters (a, b) associated with the labels of transmitted PCR efficiency t (t 1, t 2, …, t 61) and phase profiles p t (p t1, p t2, …, p t61), and it has 61 outputs for each prediction. As for R-mode operating meta-atoms, the FPN also has two input parameters (a, b) with the labels of reflected PCR efficiency r (r 1, r 2, …, r 61) and phase profiles p r (p r1, p r2, …, p r61), and it also has 61 outputs for predicting each meta-atom of R-mode, which has a total of 5,039 data (including 80 % training data and 20 % testing data) in the identical sweeping range of a and b. In this case, PyTorch is employed for the experiment, operated on a system with a CPU of Intel Core i9-13900K and an NVIDIA GeForce GTX 4060Ti GPU. And it took about one and a half months for the training data collection, and it spent about an hour and a half for the four sets of data training.

6.3 Loss function and error

The loss function of MSE is defined to characterize the FPN, which can evaluate the squared differences between prediction results by FPN (S prediction) and simulated results by FIM (S FIM) [31]:

(6) MSE=1Ni=1,2,,N(SpredictionSFIM)2MSE=1Ni=1,2,,N(SpredictionSFIM)2

The fractional error (FE) is also employed to evaluate the differences between simulated results by FIM and prediction results by FPN [31]:

(7) FE=1Ni=1,2,,N(SpredictionSFIM)SFIMFE=1Ni=1,2,,N(SpredictionSFIM)SFIM

Moreover, the comprehensive details regarding the weight decay implementation, including its mathematical formulation within the loss function and the specific value are described in Supplementary Section 18. Besides, the thorough discussion and illustration of the transfer learning process are described in Supplementary Section 19.

6.4 Evaluation of TRARM

In order to validate the reliability of proposed inverse design method, two FVB generators with different parameters marked as TRARM-Ⅰ (l = 0) and TRARM-Ⅱ (l = −2) are rapidly constructed (the corresponding structural parameters of the meta-atoms are listed in Supplementary Section 17). The FIM and vectorial ray-based diffraction integral (VRBDI) are both employed to characterize their performances. To begin with, the FIM is utilized to calculate the near-field EM distributions of TRARM, where the open boundary is applied around the cubic simulation area. In addition, electric boundary is applied at the back of VO2 layer for R-mode. Subsequently, VRBDI is employed to further evaluate their far-field performances, which can simultaneously consider ray tracing and Fourier optics [51]. After the calculation of PSF simulated by FIM, the convolution operation can be performed to obtain THz imaging and edge detection. Moreover, in order to demonstrate the merits of the proposed method, the comparison of THz achromatic metasurfaces is depicted in Supplementary Section 20.

6.5 Mode purity of OAM

The mode purity is investigated to characterize the performances of FVB generated from TRARM-Ⅱ, which is performed by employing modal decomposition in Fourier transform [16]:

(8) φ(θ)=+[12π2π0φ(θ)eilθdθ]eilθφ(θ)=+12π2π0φ(θ)eilθdθeilθ

where φ(θ) is the phase samples around the doughnut-shaped annulus, and ei represents the harmonic item connected with OAM eigenstate.


Corresponding author: Wen-Hui Fan, State Key Laboratory of Transient Optics and Photonics, Xi’an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences, Xi’an 710119, P.R. China; University of Chinese Academy of Sciences, Beijing 100049, P.R. China; and Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan 030006, P.R. China, E-mail: 

Funding source: The CAS “Light of West China” Program

Award Identifier / Grant number: XAB2021YN07

Funding source: Shaanxi young science and technology star

Award Identifier / Grant number: 2022KJXX-97

Award Identifier / Grant number: 2023YFF0718602

Award Identifier / Grant number: 61675230, 61905276

  1. Research funding: National Natural Science Foundation of China (61675230, 61905276); National Key Research and Development Program of China (2023YFF0718602); The CAS “Light of West China” Program (XAB2021YN07); Shaanxi young science and technology star (2022KJXX-97).

  2. Author contribution: XQJ and LRZ conceived the idea. XQJ, XC, and PJ wrote the manuscript. WHF and XC revised the manuscript. CQ, HY, and QW participated in the discussions. All authors read and approved the final manuscript. WHF guided the entire work.

  3. Conflict of interest: Authors state no conflict of interest.

  4. Data availability: The datasets generated and analysed during the current study are available from the corresponding author upon reasonable request.

References

[1] T. L. Cocker, V. Jelic, R. Hillenbrand, and F. A. Hegmann, “Nanoscale terahertz scanning probe microscopy,” Nat. Photonics, vol. 15, no. 1, pp. 1–28, 2023. https://doi.org/10.1038/s41566-021-00835-6.Search in Google Scholar

[2] Y. E, L. Zhang, A. Tcypkin, S. Kozlov, C. Zhang, and X. C. Zhang, “Broadband THz sources from gases to liquids,” Ultrafast Sci., vol. 10, no. 2021, p. 9892763, 2021. https://doi.org/10.34133/2021/9892763.Search in Google Scholar

[3] R. Degl’Innocenti, H. Y. Lin, and M. Navarro-Cía, “Recent progress in terahertz metamaterial modulators,” Nanophotonics, vol. 11, no. 8, pp. 1485–1514, 2022. https://doi.org/10.1515/nanoph-2021-0803.Search in Google Scholar PubMed PubMed Central

[4] X. Feng, et al.., “Direct emission of focused terahertz vortex beams using indium-tin-oxide-based fresnel zone plates,” Adv. Opt. Mater., vol. 11, no. 1, p. 2201628, 2023. https://doi.org/10.1002/adom.202201628.Search in Google Scholar

[5] F. Y. Li, et al.., “All-dielectric terahertz metasurface for linearly-polarized multichannel transmission and superposition states of spherical and vortex waves,” Photonics Res., vol. 11, no. 3, pp. 485–493, 2023. https://doi.org/10.1364/prj.477381.Search in Google Scholar

[6] Y. H. Guo, et al.., “Spin-decoupled metasurface for simultaneous detection of spin and orbital angular momenta via momentum transformation,” Light Sci. Appl., vol. 10, no. 4, p. 63, 2021. https://doi.org/10.1038/s41377-021-00497-7.Search in Google Scholar PubMed PubMed Central

[7] X. Y. Li, et al.., “Monolithic spiral metalens for ultrahigh-capacity and single-shot sorting of full angular momentum state,” Adv. Funct. Mater., vol. 34, no. 7, p. 202311286, 2021. https://doi.org/10.1002/adfm.202311286.Search in Google Scholar

[8] J. Y. Xie, et al.., “Integrated terahertz vortex beam emitter for rotating target detection,” Adv. Photonics, vol. 5, no. 6, p. 066002, 2023. https://doi.org/10.1117/1.ap.5.6.066002.Search in Google Scholar

[9] Y. Zhu, et al.., “Geometric metasurface for polarization synthesis and multidimensional multiplexing of terahertz converged vortices,” Photonics Res., vol. 10, no. 6, pp. 1517–1532, 2022. https://doi.org/10.1364/prj.455459.Search in Google Scholar

[10] X. L. Zhu, et al.., “Inverse design of a spatial filter in edge enhanced imaging,” Opt. Lett., vol. 45, no. 9, pp. 2542–2545, 2020. https://doi.org/10.1364/ol.391429.Search in Google Scholar

[11] C. H. Chu, et al.., “Intelligent phase contrast meta-microscope system,” Nano Lett., vol. 23, no. 24, pp. 11630–11637, 2023. https://doi.org/10.1021/acs.nanolett.3c03484.Search in Google Scholar PubMed

[12] H. Zhao, X. K. Wang, S. T. Liu, and Y. Zhang, “Highly efficient vectorial field manipulation using a transmitted tri-layer metasurface in the terahertz band,” Opto-Electron. Adv., vol. 6, no. 2, p. 220012, 2023. https://doi.org/10.29026/oea.2023.220012.Search in Google Scholar

[13] J. Li, et al.., “Racemic dielectric metasurfaces for arbitrary terahertz polarization rotation and wavefront manipulation,” Opto-Electron. Adv., vol. 7, no. 10, p. 240075, 2024. https://doi.org/10.29026/oea.2024.240075.Search in Google Scholar

[14] L. Ke, et al.., “Tight focusing field of cylindrical vector beams based on cascaded low-refractive index metamaterials,” Nanophotonics, vol. 12, no. 18, pp. 3563–3578, 2023. https://doi.org/10.1515/nanoph-2023-0261.Search in Google Scholar PubMed PubMed Central

[15] J. X. Fan, et al.., “Hybrid bound states in the continuum in terahertz metasurfaces,” Opto-Electron. Sci., vol. 2, no. 4, p. 230006, 2023. https://doi.org/10.29026/oes.2023.230006.Search in Google Scholar

[16] X. Q. Jiang, W. H. Fan, C. Qin, and X. Chen, “Ultra-broadband polarization conversion metasurface with high transmission for efficient multi-functional wavefront manipulation in the terahertz range,” Nanomaterials, vol. 11, no. 11, 2021, https://doi.org/10.3390/nano11112895.Search in Google Scholar PubMed PubMed Central

[17] S. Sun, H. F. Ma, Y. T. Chen, and T. J. Cui, “Transmission-reflection-integrated metasurface with simultaneous amplitude and phase controls of circularly polarized waves in full space,” Laser Photonics Rev., vol. 18, no. 3, p. 2300945, 2024. https://doi.org/10.1002/lpor.202300945.Search in Google Scholar

[18] Z. J. Zhang, et al.., “Broadband achromatic transmission-reflection-integrated metasurface based on frequency multiplexing and dispersion engineering,” Adv. Opt. Mater., vol. 9, no. 7, p. 1801748, 2021. https://doi.org/10.1002/adom.202001736.Search in Google Scholar

[19] J. Y. Wang, et al.., “Transmission-reflection-integrated quadratic phase metasurface for multifunctional electromagnetic manipulation in full space,” Adv. Opt. Mater., vol. 10, no. 6, p. 2102111, 2022. https://doi.org/10.1002/adom.202102111.Search in Google Scholar

[20] G. Kim, et al.., “Metasurface-driven full-space structured light for three-dimensional imaging,” Nat. Commun., vol. 13, no. 1, p. 5920, 2022. https://doi.org/10.1038/s41467-022-32117-2.Search in Google Scholar PubMed PubMed Central

[21] J. Sung, G. Y. Lee, C. Choi, J. Hong, and B. Lee, “Polarization-dependent asymmetric transmission using a bifacial metasurface,” Nanoscale Horiz., vol. 5, no. 11, pp. 1487–1495, 2020. https://doi.org/10.1039/d0nh00319k.Search in Google Scholar PubMed

[22] H. Park, et al.., “Electrically tunable THz graphene metasurface wave retarders,” Nanophotonics, vol. 12, no. 3, pp. 2553–2562, 2023. https://doi.org/10.1515/nanoph-2022-0812.Search in Google Scholar PubMed PubMed Central

[23] T. Gu, H. J. Kim, C. Rivero-Baleine, and J. J. Hu, “Reconfigurable metasurfaces towards commercial success,” Nat. Photonics, vol. 17, no. 1, pp. 48–58, 2023. https://doi.org/10.1038/s41566-022-01099-4.Search in Google Scholar

[24] S. Wall, et al.., “Ultrafast disordering of vanadium dimers in photoexcited VO2,” Science, vol. 362, no. 6414, pp. 572–576, 2018. https://doi.org/10.1126/science.aau3873.Search in Google Scholar PubMed

[25] N. Wang, W. Yan, Y. R. Qu, S. Q. Ma, T. Li, and M. Qiu, “Intelligent designs in nanophotonics: from optimization towards inverse creation,” PhotoniX, vol. 2, no. 1, 2021, https://doi.org/10.1186/s43074-021-00044-y.Search in Google Scholar

[26] W. Ma, Z. C. Liu, Z. A. Kudyshev, A. Boltasseva, W. S. Cai, and Y. M. Liu, “Deep learning for the design of photonic structures,” Nat. Photonics, vol. 15, no. 2, pp. 77–90, 2021. https://doi.org/10.1038/s41566-020-0685-y.Search in Google Scholar

[27] R. X. Peng, S. M. Ren, J. Malof, and W. J. Padilla, “Transfer learning for metamaterial design and simulation,” Nanophotonics, vol. 13, no. 13, pp. 2323–2334, 2024. https://doi.org/10.1515/nanoph-2023-0691.Search in Google Scholar PubMed PubMed Central

[28] W. Ma, F. Cheng, Y. H. Xu, Q. L. Wen, and Y. M. Liu, “Probabilistic representation and inverse design of metamaterials based on a deep generative model with semi-supervised learning strategy,” Adv. Mater., vol. 31, no. 35, p. 1901111, 2019. https://doi.org/10.1002/adma.201901111.Search in Google Scholar PubMed

[29] W. Ma, F. Cheng, and Y. M. Liu, “Deep-learning-enabled on-demand design of chiral metamaterials,” ACS Nano, vol. 12, no. 6, pp. 6326–6334, 2018. https://doi.org/10.1021/acsnano.8b03569.Search in Google Scholar PubMed

[30] F. Sterl, N. Strohfeldt, S. Both, E. Herkert, T. Weiss, and H. Giessen, “Design principles for sensitivity optimization in plasmonic hydrogen sensors,” ACS Sens., vol. 5, no. 4, pp. 917–927, 2020. https://doi.org/10.1021/acssensors.9b02436.Search in Google Scholar PubMed

[31] S. S. An, et al.., “A deep learning approach for objective-driven all-dielectric metasurface design,” ACS Photonics, vol. 6, no. 12, pp. 3196–3207, 2019. https://doi.org/10.1021/acsphotonics.9b00966.Search in Google Scholar

[32] D. J. Liu, Y. X. Tan, E. Khoram, and Z. F. Yu, “Training deep neural networks for the inverse design of nanophotonic structures,” ACS Photonics, vol. 5, no. 4, pp. 1365–1369, 2018. https://doi.org/10.1021/acsphotonics.7b01377.Search in Google Scholar

[33] F. L. Wang, S. Q. Zhao, Y. Z. Wen, J. B. Sun, and J. Zhou, “High efficiency visible achromatic metalens design via deep learning,” Adv. Opt. Mater., vol. 11, no. 18, 2023, https://doi.org/10.1002/adom.202300394.Search in Google Scholar

[34] I. Tanriover, W. Hadibrata, and K. Aydin, “Physics-based approach for a neural networks enabled design of all-dielectric metasurfaces,” ACS Photonics, vol. 7, no. 8, pp. 1957–1964, 2020. https://doi.org/10.1021/acsphotonics.0c00663.Search in Google Scholar

[35] K. Ou, et al.., “Mid-infrared polarization-controlled broadband achromatic metadevice,” Sci. Adv., vol. 6, no. 37, p. eabc0711, 2020. https://doi.org/10.1126/sciadv.abc0711.Search in Google Scholar PubMed PubMed Central

[36] S. M. Wang, et al.., “Broadband achromatic optical metasurface devices,” Nat. Commun., vol. 8, no. 1, p. 187, 2017. https://doi.org/10.1038/s41467-017-00166-7.Search in Google Scholar PubMed PubMed Central

[37] F. Aieta, M. A. Kats, P. Genevet, and F. Capasso, “Multiwavelength achromatic metasurfaces by dispersive phase compensation,” Science, vol. 347, no. 6228, pp. 1342–1345, 2015. https://doi.org/10.1126/science.aaa2494.Search in Google Scholar PubMed

[38] L. R. Zhao, X. Q. Jiang, C. X. Li, S. X. Gong, and W. X. Yu, “High-NA and broadband achromatic metalens for sub-diffraction focusing of long-wavelength infrared waves,” Results Phys., vol. 46, 2023, https://doi.org/10.1016/j.rinp.2023.106308.Search in Google Scholar

[39] S. Shrestha, A. C. Overvig, M. Lu, A. Stein, and N. F. Yu, “Broadband achromatic dielectric metalenses,” Light: Sci. Appl., vol. 7, no. 1, p. 85, 2018. https://doi.org/10.1038/s41377-018-0078-x.Search in Google Scholar PubMed PubMed Central

[40] Y. C. Zhao, et al.., “Dynamic photoinduced controlling of the large phase shift of terahertz waves via vanadium dioxide coupling nanostructures,” ACS Photonics, vol. 5, no. 8, pp. 3040–3050, 2018. https://doi.org/10.1021/acsphotonics.8b00276.Search in Google Scholar

[41] X. Q. Jiang, W. H. Fan, X. Chen, and H. Yan, “Ultrahigh-Q terahertz sensor based on simple all-dielectric metasurface with toroidal dipole resonance,” Appl. Phys. Express, vol. 14, no. 10, p. 102008, 2021. https://doi.org/10.35848/1882-0786/ac27b7.Search in Google Scholar

[42] P. C. Wu, et al.., “Optical anapole metamaterial,” ACS Nano, vol. 12, no. 2, pp. 1920–1927, 2018. https://doi.org/10.1021/acsnano.7b08828.Search in Google Scholar PubMed

[43] Z. X. Fan, et al.., “Transfer-learning-assisted inverse metasurface design for 30% data savings,” Phys. Rev. Appl., vol. 18, no. 2, p. 024022, 2022. https://doi.org/10.1103/physrevapplied.18.024022.Search in Google Scholar

[44] P. Z. Ramirez, A. Cardace, L. De Luigi, A. Tonioni, S. Salti, and L. D. Stefano, “Learning good features to transfer across tasks and domains,” IEEE Trans. Pattern Anal. Mach. Intell., vol. 45, no. 8, pp. 9981–9995, 2023. https://doi.org/10.1109/tpami.2023.3240316.Search in Google Scholar

[45] X. X. Jia, X. C. Feng, H. W. Yong, and D. Y. Meng, “Weight decay with tailored Adam on scale-invariant weights for better generalization,” IEEE Trans. Neural Netw. Learn. Syst., vol. 35, no. 5, pp. 6936–6947, 2024. https://doi.org/10.1109/tnnls.2022.3213536.Search in Google Scholar

[46] Y. J. Wang, et al.., “High-efficiency broadband achromatic metalens for near-IR biological imaging window,” Nat. Commun., vol. 12, no. 1, p. 5560, 2021. https://doi.org/10.1038/s41467-021-25797-9.Search in Google Scholar PubMed PubMed Central

[47] Y. F. Hu, et al.., “Dielectric metasurface zone plate for the generation of focusing vortex beams,” PhotoniX, vol. 2, no. 1, p. 10, 2021. https://doi.org/10.1186/s43074-021-00035-z.Search in Google Scholar

[48] C. Zhou, et al.., “Laplace differentiator based on metasurface with toroidal dipole resonance,” Adv. Funct. Mater., vol. 34, no. 19, p. 2313777, 2024. https://doi.org/10.1002/adfm.202313777.Search in Google Scholar

[49] S. Han, et al.., “All-dielectric active terahertz photonics driven by bound states in the continuum,” Adv. Mater., vol. 31, no. 37, p. 1901921, 2019. https://doi.org/10.1002/adma.201901921.Search in Google Scholar PubMed

[50] X. Q. Jiang, et al.., “Continuously varifocal metalens for broadband achromatic focusing of terahertz waves,” J. Sci. – Adv. Mater. Dev., vol. 8, no. 3, p. 100560, 2023. https://doi.org/10.1016/j.jsamd.2023.100560.Search in Google Scholar

[51] E. Otte, V. Bobkova, S. Trinschek, C. Rosales-Guzmán, and C. Denz, “Single-shot all-digital approach for measuring the orbital angular momentum spectrum of light,” APL Photonics, vol. 7, no. 8, p. 086150, 2022. https://doi.org/10.1063/5.0086536.Search in Google Scholar


Supplementary Material

This article contains supplementary material (https://doi.org/10.1515/nanoph-2024-0680).


Received: 2024-11-26
Accepted: 2025-03-03
Published Online: 2025-03-25

© 2025 the author(s), published by De Gruyter, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 9.5.2025 from https://www.degruyterbrill.com/document/doi/10.1515/nanoph-2024-0680/html
Scroll to top button