• Photonics Research
  • Vol. 12, Issue 1, 123 (2024)
Zhi-Dan Lei1、†, Yi-Duo Xu1、†, Cheng Lei1、3, Yan Zhao1、2、4, and Du Wang1、*
Author Affiliations
  • 1The Institute of Technological Sciences, Wuhan University, Wuhan 430072, China
  • 2College of Materials Science and Engineering, Sichuan University, Chengdu 610065, China
  • 3e-mail: leicheng@whu.edu.cn
  • 4e-mail: yan2000@whu.edu.cn
  • show less
    DOI: 10.1364/PRJ.505991 Cite this Article Set citation alerts
    Zhi-Dan Lei, Yi-Duo Xu, Cheng Lei, Yan Zhao, Du Wang. Dynamic multifunctional metasurfaces: an inverse design deep learning approach[J]. Photonics Research, 2024, 12(1): 123 Copy Citation Text show less
    Illustration of multifunctional OMs and the material optical properties. (a) Scheme of structured Sb2Te3 OMs for optical data encryption. (b) Structure design and related parameters of meta-unit. The periods Px and Py of meta-cell are 300 nm. The MIM layers’ thicknesses h1, h2, and h3 are 30, 90, and 130 nm, respectively. (c) Distribution of the propagation phase with the size of meta-unit. (d) Distribution of PB phases with the rotation angle of meta-unit. (e), (f) Refractive index n and extinction coefficient k of Sb2Te3 in (e) amorphous and (f) crystalline states [58].
    Fig. 1. Illustration of multifunctional OMs and the material optical properties. (a) Scheme of structured Sb2Te3 OMs for optical data encryption. (b) Structure design and related parameters of meta-unit. The periods Px and Py of meta-cell are 300 nm. The MIM layers’ thicknesses h1, h2, and h3 are 30, 90, and 130 nm, respectively. (c) Distribution of the propagation phase with the size of meta-unit. (d) Distribution of PB phases with the rotation angle of meta-unit. (e), (f) Refractive index n and extinction coefficient k of Sb2Te3 in (e) amorphous and (f) crystalline states [58].
    Design process for the multifunctional OMs. (a) Proposed multifunctional OMs, (b) retrieving network model, (c) meta-unit library, and (d) predicting network model. These two deep learning network models are integrated into the GS iterative optimization algorithm, enabling bidirectional linkage between design objectives and OM geometric parameters.
    Fig. 2. Design process for the multifunctional OMs. (a) Proposed multifunctional OMs, (b) retrieving network model, (c) meta-unit library, and (d) predicting network model. These two deep learning network models are integrated into the GS iterative optimization algorithm, enabling bidirectional linkage between design objectives and OM geometric parameters.
    Schematic diagram of the deep learning model for a single meta-unit design. (a) The design parameters of the meta-unit and the crystalline phase organization of Sb2Te3 are amorphous. (b) Predicting model for reflectance spectra and (c) retrieving model of the required intensity and phase. The forward prediction from design parameters to reflection spectra is fixed, whereas the on-demand inverse design process is characterized by uncertainty, thus ensuring the multiplicity of design outcomes.
    Fig. 3. Schematic diagram of the deep learning model for a single meta-unit design. (a) The design parameters of the meta-unit and the crystalline phase organization of Sb2Te3 are amorphous. (b) Predicting model for reflectance spectra and (c) retrieving model of the required intensity and phase. The forward prediction from design parameters to reflection spectra is fixed, whereas the on-demand inverse design process is characterized by uncertainty, thus ensuring the multiplicity of design outcomes.
    Training results of the forward predicting network model. (a) Datasets for the forward deep learning network. The total data quantity of each cuboid is 23,328. (b) Distribution of datasets. (c) Correlation analysis of parameters. (d) Training and validation loss of X polarization, LCP, and RCP at wavelength of 520 nm.
    Fig. 4. Training results of the forward predicting network model. (a) Datasets for the forward deep learning network. The total data quantity of each cuboid is 23,328. (b) Distribution of datasets. (c) Correlation analysis of parameters. (d) Training and validation loss of X polarization, LCP, and RCP at wavelength of 520 nm.
    Training results of the retrieving model. (a) Distribution of datasets and (b) correspondence among the real, predicted, and generated values of X polarization, LCP, and RCP. The black plot represents the target values, the red plot represents the predicted values, and the blue plot represents the generated values. (c) Pearson correlations among the real, predicted, and generated values for each polarization.
    Fig. 5. Training results of the retrieving model. (a) Distribution of datasets and (b) correspondence among the real, predicted, and generated values of X polarization, LCP, and RCP. The black plot represents the target values, the red plot represents the predicted values, and the blue plot represents the generated values. (c) Pearson correlations among the real, predicted, and generated values for each polarization.
    Intensity and phase regeneration results of the retrieving model. Comparison images among the target phase, target image, reconstructed phase, and reconstructed image for each polarization state.
    Fig. 6. Intensity and phase regeneration results of the retrieving model. Comparison images among the target phase, target image, reconstructed phase, and reconstructed image for each polarization state.
    Target image and calculated results of multifunctional OMs under X polarization, LCP, and RCP. These calculated results are obtained through the combination of simulation results and the utilization of Fraunhofer diffraction integrals.
    Fig. 7. Target image and calculated results of multifunctional OMs under X polarization, LCP, and RCP. These calculated results are obtained through the combination of simulation results and the utilization of Fraunhofer diffraction integrals.
    BlocksLayersSize-InSize-Out
    ResBlockInput3×1
    ResBlockFirst linear3×1512×1
    ResBlockBasBlock1512×11024×1
    ResBlockBasBlock21024×12048×1
    ResBlockBasBlock32048×12048×1
    ResBlockBasBlock42048×11024×1
    ResBlockBasBlock51024×1512×1
    ResBlockBasBlock6512×1128×1
    ConvBlockResConv1128×13×128×512
    ConvBlockResConv23×128×5123×128×64
    GRUBlockGRU3×128×643×128
    Table 1. Detailed Configuration of the Predicting Model
    BlocksLayersSize-InSize-Out
    EncoderInput(3+6)×11000×1
    EncoderFirst linear1000×11000×1
    EncoderSecond linear1000×1800×1
    EncoderThird linear800×1400×1
    EncoderFourth linear400×1200×1
    EncoderFifth linear200×1128×1
    EncoderOutput128×1(9+9)×1
    DecoderInput(9+6)×1128×1
    DecoderFirst linear128×1200×1
    DecoderSecond linear200×1400×1
    DecoderThird linear400×1800×1
    DecoderFourth linear800×11000×1
    DecoderFifth linear1000×11000×1
    DecoderOutput1000×13 ×1
    Table 2. Detailed Configuration of the Retrieving Model
    Zhi-Dan Lei, Yi-Duo Xu, Cheng Lei, Yan Zhao, Du Wang. Dynamic multifunctional metasurfaces: an inverse design deep learning approach[J]. Photonics Research, 2024, 12(1): 123
    Download Citation