在实际光学加工中将联合优化得到的连续相位分布转化为可制造的超表面Metasurface结构是一个涉及离散化、单元库映射、制造约束集成和工艺补偿的关键步骤。以下是完整的转化流程、核心方法及程序实现。一、从连续相位到可制造结构的核心转化流程步骤输入输出关键操作与目标1. 相位量化连续相位分布 $\phi(x,y) \in [0, 2\pi)$$N$ 级离散相位分布 $\phi_q(x,y) \in {\frac{2\pi k}{N}}$将连续值映射到有限的制造级别减少工艺复杂度。2. 单元库映射离散相位值 $\phi_q$ 与目标波长 $\lambda$、入射角 $\theta$纳米结构几何参数如柱体直径 $D$、高度 $H$、材料 $n$通过预仿真的“相位-结构”查找表LUT或代理模型为每个相位值找到对应的物理结构。3. 制造约束嵌入理想结构参数符合工艺规范的结构参数如最小特征尺寸、最大深宽比、最小间隙应用设计规则检查DRC调整结构以避免制造失败。4. 邻近效应修正初始结构布局经过光学邻近校正OPC或反演优化的结构补偿光刻/刻蚀过程中的衍射、散射等效应确保加工后结构接近设计目标。5. 文件导出校正后的结构数据制造文件如GDSII、OASIS转换为标准格式供电子束光刻EBL或深紫外光刻DUV使用。二、关键技术实现与程序示例2.1 相位量化与离散化在联合优化中相位参数通常是连续的。制造时需将其量化为有限的离散级别如 2、4、8 级相位这通常由制造工艺如灰度光刻的精度或多层刻蚀的层数决定。import torch import numpy as np def quantize_phase(continuous_phase, levels, methodnearest): 将连续相位量化为离散级别。 参数: continuous_phase: 连续相位张量范围 [0, 2π) 或 [-π, π) levels: 离散级别数如 2, 4, 8 method: 量化方法nearest最近邻或 stochastic随机舍入用于训练中 返回: quantized_phase: 量化后的相位 phase_indices: 对应的级别索引用于单元库查找 # 确保相位在 [0, 2π) 范围内 phase continuous_phase % (2 * np.pi) # 计算每个级别对应的相位值 phase_values torch.linspace(0, 2*np.pi, levels1)[:-1] # 0, 2π/levels, 4π/levels, ... if method nearest: # 计算每个像素到所有级别的距离取最近 # 扩展维度以便广播计算 phase_expanded phase.unsqueeze(-1) # (..., levels) distances torch.abs(phase_expanded - phase_values) indices torch.argmin(distances, dim-1) quantized phase_values[indices] elif method stochastic: # 随机舍入在训练中引入量化噪声以增强鲁棒性 # 将连续相位映射到级别索引的连续值 scaled phase / (2*np.pi) * levels lower torch.floor(scaled) upper torch.ceil(scaled) # 随机选择上界或下界 prob scaled - lower rand torch.rand_like(prob) indices torch.where(rand prob, upper, lower).long() % levels quantized phase_values[indices] else: raise ValueError(method 必须是 nearest 或 stochastic) return quantized, indices # 示例将联合优化得到的相位分布量化为4级 continuous_phase torch.rand(256, 256) * 2 * np.pi # 模拟优化后的连续相位 quantized_phase, level_indices quantize_phase(continuous_phase, levels4, methodnearest) print(f连续相位范围: [{continuous_phase.min():.3f}, {continuous_phase.max():.3f}]) print(f量化后相位值集合: {torch.unique(quantized_phase)}) print(f级别索引形状: {level_indices.shape})2.2 单元库映射从相位到纳米结构超表面通常由亚波长纳米结构如硅或氮化硅纳米柱组成。每个离散相位值对应一个特定的结构几何参数。通常通过电磁仿真如FDTD、RCWA预先建立“相位-结构”查找表LUT。class MetasurfaceUnitCellLibrary: 超表面单元库存储相位与结构参数的映射关系。 参考来源 中基于相位梯度逆向设计的超表面架构。 def __init__(self, materialSiN, wavelength532e-9, substrateSiO2): self.material material # 纳米柱材料 self.wavelength wavelength self.substrate substrate self.library {} # 键相位值弧度值结构参数字典 def build_library_from_simulation(self, phase_levels, height_range, diameter_range): 通过仿真构建单元库实际中应使用FDTD/RCWA仿真数据。 这里用解析模型近似。 参数: phase_levels: 离散相位级别数 height_range: (min, max) 纳米柱高度范围米 diameter_range: (min, max) 纳米柱直径范围米 # 为每个相位级别生成结构参数 target_phases torch.linspace(0, 2*np.pi, phase_levels1)[:-1] for phase in target_phases: # 简化的解析模型相位延迟 φ (2π/λ) * (n_eff - n_sub) * H # 其中 n_eff 是纳米柱的有效折射率与直径D相关 # 实际中应使用严格的电磁仿真 H height_range[0] (phase / (2*np.pi)) * (height_range[1] - height_range[0]) # 假设有效折射率随直径线性变化简化模型 n_eff_min, n_eff_max 1.5, 2.5 # 硅氮化物在可见光波段的近似 D diameter_range[0] (phase / (2*np.pi)) * (diameter_range[1] - diameter_range[0]) self.library[phase.item()] { height: H.item(), diameter: D.item(), material: self.material, phase: phase.item() } def map_phase_to_structure(self, quantized_phase): 将量化后的相位分布映射为结构参数分布。 参数: quantized_phase: 量化相位张量值在单元库的键集中 返回: height_map: 高度分布图米 diameter_map: 直径分布图米 height_map torch.zeros_like(quantized_phase) diameter_map torch.zeros_like(quantized_phase) # 获取所有唯一的相位值 unique_phases torch.unique(quantized_phase) for phase_val in unique_phases: phase_key phase_val.item() if phase_key in self.library: params self.library[phase_key] mask (quantized_phase phase_val) height_map[mask] params[height] diameter_map[mask] params[diameter] else: # 如果相位值不在库中使用最近邻插值 available_phases torch.tensor(list(self.library.keys())) nearest_phase available_phases[torch.argmin(torch.abs(available_phases - phase_key))] params self.library[nearest_phase.item()] mask (quantized_phase phase_val) height_map[mask] params[height] diameter_map[mask] params[diameter] print(f警告: 相位 {phase_key:.3f} 不在库中使用 {nearest_phase.item():.3f} 替代) return height_map, diameter_map # 示例构建单元库并映射 unit_cell_lib MetasurfaceUnitCellLibrary(materialSiN, wavelength532e-9) unit_cell_lib.build_library_from_simulation( phase_levels8, height_range(200e-9, 600e-9), # 200nm 到 600nm diameter_range(100e-9, 300e-9) # 100nm 到 300nm ) # 将量化相位映射为结构参数 height_map, diameter_map unit_cell_lib.map_phase_to_structure(quantized_phase) print(f高度图范围: [{height_map.min():.2e}, {height_map.max():.2e}] 米) print(f直径图范围: [{diameter_map.min():.2e}, {diameter_map.max():.2e}] 米)2.3 制造约束嵌入与设计规则检查DRC实际制造工艺存在限制如最小特征尺寸、最大深宽比高度/直径、最小结构间隙等。需要在映射后应用这些约束。def apply_manufacturing_constraints(height_map, diameter_map, constraints): 应用制造约束调整结构参数以满足工艺要求。 参数: height_map: 原始高度分布 diameter_map: 原始直径分布 constraints: 约束条件字典包含: min_diameter: 最小直径米 max_aspect_ratio: 最大深宽比高度/直径 min_gap: 相邻结构最小间隙米 diameter_quantization_step: 直径量化步长米可选 返回: height_map_adj: 调整后的高度分布 diameter_map_adj: 调整后的直径分布 height_map_adj height_map.clone() diameter_map_adj diameter_map.clone() # 1. 最小直径约束 min_d constraints.get(min_diameter, 80e-9) # 例如80nm diameter_map_adj torch.where(diameter_map_adj min_d, min_d, diameter_map_adj) print(f应用最小直径约束: {min_d:.2e} 米) # 2. 最大深宽比约束 max_ar constraints.get(max_aspect_ratio, 4.0) # 例如高度/直径 ≤ 4 # 计算当前深宽比 aspect_ratio height_map_adj / (diameter_map_adj 1e-12) # 如果超过限制调整高度 mask_ar aspect_ratio max_ar if mask_any(mask_ar): new_heights diameter_map_adj[mask_ar] * max_ar height_map_adj[mask_ar] new_heights print(f调整了 {mask_ar.sum().item()} 个单元的高度以满足深宽比约束) # 3. 直径量化如果工艺要求离散直径 quant_step constraints.get(diameter_quantization_step, None) if quant_step is not None: # 将直径量化到最近的步长倍数 diameter_map_adj torch.round(diameter_map_adj / quant_step) * quant_step print(f直径量化步长: {quant_step:.2e} 米) # 4. 最小间隙约束需要布局信息这里简化为直径调整 # 实际中需要更复杂的几何检查可能涉及布局重排 min_gap constraints.get(min_gap, 20e-9) # 20nm # 简化处理确保直径间隙不超过单元周期 # 假设单元周期为P则直径应满足 D ≤ P - min_gap # 这里假设周期为400nm period 400e-9 max_diameter period - min_gap diameter_map_adj torch.where(diameter_map_adj max_diameter, max_diameter, diameter_map_adj) return height_map_adj, diameter_map_adj # 应用约束 constraints { min_diameter: 100e-9, # 100nm 最小直径 max_aspect_ratio: 5.0, # 深宽比 ≤ 5 min_gap: 30e-9, # 30nm 最小间隙 diameter_quantization_step: 10e-9 # 直径10nm步长量化 } height_map_adj, diameter_map_adj apply_manufacturing_constraints( height_map, diameter_map, constraints ) print(f调整后高度范围: [{height_map_adj.min():.2e}, {height_map_adj.max():.2e}] 米) print(f调整后直径范围: [{diameter_map_adj.min():.2e}, {diameter_map_adj.max():.2e}] 米)2.4 邻近效应修正OPC在纳米尺度光刻中由于衍射和散射加工出的结构会与设计图形有偏差。光学邻近校正OPC通过预变形设计图形来补偿这种效应。def simple_opc_correction(diameter_map, wavelength13.5e-9, na0.33, sigma0.5): 简化的光学邻近校正OPC模型。 实际OPC需要严格的光刻仿真和逆优化这里使用基于规则的边缘偏移作为示例。 参数: diameter_map: 设计直径分布 wavelength: 光刻波长例如EUV为13.5nm na: 数值孔径 sigma: 部分相干因子 返回: diameter_map_corrected: 校正后的直径分布 # 简化模型根据特征尺寸计算需要的偏移量 # 经验规则小尺寸特征需要更大的补偿 diameter_map_corrected diameter_map.clone() # 计算光刻系统的分辨率极限瑞利准则 resolution 0.61 * wavelength / na # 对于接近或低于分辨率的特征进行补偿 # 这里使用一个简化的线性补偿模型 min_d diameter_map.min() max_d diameter_map.max() # 补偿量尺寸越小补偿越大 # 实际中应基于严格的光刻仿真曲线 for i in range(diameter_map.shape[0]): for j in range(diameter_map.shape[1]): d diameter_map[i, j].item() # 如果特征尺寸小于2倍分辨率进行补偿 if d 2 * resolution: # 补偿量随尺寸减小而增加 compensation resolution * (1 - d / (2 * resolution)) * 0.3 diameter_map_corrected[i, j] d compensation print(fOPC校正完成平均补偿量: {(diameter_map_corrected - diameter_map).abs().mean():.2e} 米) return diameter_map_corrected # 应用简化的OPC校正 diameter_map_opc simple_opc_correction(diameter_map_adj, wavelength13.5e-9, na0.33)2.5 GDSII文件生成将最终的结构布局导出为标准制造格式如GDSII。import gdspy import numpy as np def export_to_gdsii(diameter_map, height_map, period, filenamemetasurface.gds): 将超表面结构导出为GDSII文件。 每个单元由一个圆柱体纳米柱表示。 参数: diameter_map: 直径分布米 height_map: 高度分布米- 在GDSII中通常用层表示而不是3D高度 period: 单元周期米 filename: 输出文件名 # 创建库和单元 lib gdspy.GdsLibrary() cell lib.new_cell(METASURFACE) # 设置图层根据工艺层定义 # 例如层1用于纳米柱图案层2用于标记等 layer_pillar 1 layer_text 2 # 转换单位通常GDSII使用微米或纳米 scale 1e9 # 转换为纳米 # 遍历所有单元 rows, cols diameter_map.shape for i in range(rows): for j in range(cols): # 计算单元中心位置 x_center j * period * scale y_center (rows - 1 - i) * period * scale # GDSII坐标系Y轴向上 # 直径转换为纳米 diameter diameter_map[i, j].item() * scale # 创建圆形表示纳米柱 # 实际纳米柱可能是圆柱但在GDSII中表示为圆盘俯视图 radius diameter / 2 # 如果直径大于0创建图形 if radius 0: circle gdspy.Round( (x_center, y_center), radius, number_of_points64, # 多边形近似圆的边数 layerlayer_pillar ) cell.add(circle) # 添加标记和注释 # 添加边界框 boundary gdspy.Rectangle( (0, 0), (cols * period * scale, rows * period * scale), layerlayer_text ) cell.add(boundary) # 添加文本标签 label gdspy.Text( fMetasurface {rows}x{cols}, period{period*1e9:.1f}nm, 10, # 字体大小 (10, -10), # 位置 layerlayer_text ) cell.add(label) # 保存文件 lib.write_gds(filename) print(fGDSII文件已保存: {filename}) # 可选显示预览 try: gdspy.LayoutViewer(lib) except: print(无法显示图形预览但文件已生成) # 导出为GDSII period 400e-9 # 单元周期400nm export_to_gdsii(diameter_map_opc, height_map_adj, period, metasurface_design.gds)2.6 完整转化流程封装将上述步骤整合为一个完整的转化流水线。class MetasurfaceManufacturingConverter: 将联合优化得到的相位分布转化为可制造结构的完整流水线。 参考来源中超表面设计流程与中深度学习辅助的逆向设计。 def __init__(self, phase_levels8, wavelength532e-9, materialSiN): self.phase_levels phase_levels self.wavelength wavelength self.material material self.unit_cell_lib MetasurfaceUnitCellLibrary(material, wavelength) def convert(self, continuous_phase, constraintsNone, apply_opcTrue): 主转换函数。 参数: continuous_phase: 联合优化得到的连续相位分布张量 constraints: 制造约束字典 apply_opc: 是否应用光学邻近校正 返回: design_data: 包含所有设计数据的字典 # 步骤1: 相位量化 print(步骤1: 相位量化...) quantized_phase, level_indices quantize_phase( continuous_phase, levelsself.phase_levels, methodnearest ) # 步骤2: 构建或加载单元库 print(步骤2: 单元库映射...) if not self.unit_cell_lib.library: # 如果库为空构建一个实际中应从仿真数据加载 self.unit_cell_lib.build_library_from_simulation( phase_levelsself.phase_levels, height_range(200e-9, 600e-9), diameter_range(100e-9, 300e-9) ) # 步骤3: 相位到结构映射 height_map, diameter_map self.unit_cell_lib.map_phase_to_structure(quantized_phase) # 步骤4: 应用制造约束 print(步骤4: 应用制造约束...) if constraints is None: constraints { min_diameter: 80e-9, max_aspect_ratio: 5.0, min_gap: 20e-9, diameter_quantization_step: 5e-9 } height_map_adj, diameter_map_adj apply_manufacturing_constraints( height_map, diameter_map, constraints ) # 步骤5: 光学邻近校正可选 if apply_opc: print(步骤5: 应用光学邻近校正...) diameter_map_final simple_opc_correction(diameter_map_adj) else: diameter_map_final diameter_map_adj # 收集所有数据 design_data { continuous_phase: continuous_phase, quantized_phase: quantized_phase, level_indices: level_indices, height_map: height_map_adj, diameter_map: diameter_map_final, constraints_applied: constraints, unit_cell_library: self.unit_cell_lib.library } print(转换完成!) return design_data def export_design(self, design_data, output_dir./manufacturing_data/): 导出所有设计文件。 import os os.makedirs(output_dir, exist_okTrue) # 1. 保存相位和结构数据为npy文件 np.save(os.path.join(output_dir, quantized_phase.npy), design_data[quantized_phase].cpu().numpy()) np.save(os.path.join(output_dir, height_map.npy), design_data[height_map].cpu().numpy()) np.save(os.path.join(output_dir, diameter_map.npy), design_data[diameter_map].cpu().numpy()) # 2. 导出GDSII文件 period 400e-9 # 假设单元周期 export_to_gdsii( design_data[diameter_map], design_data[height_map], period, os.path.join(output_dir, metasurface_layout.gds) ) # 3. 生成制造工艺说明书 self.generate_process_spec(design_data, output_dir) print(f所有设计文件已导出到: {output_dir}) def generate_process_spec(self, design_data, output_dir): 生成制造工艺说明书。 spec_content f 超表面制造工艺说明书 1. 基本参数 - 设计波长: {self.wavelength*1e9:.1f} nm - 材料: {self.material} - 相位级别数: {self.phase_levels} - 单元总数: {np.prod(design_data[diameter_map].shape)} - 阵列尺寸: {design_data[diameter_map].shape[0]} × {design_data[diameter_map].shape[1]} 2. 结构参数范围 - 纳米柱高度: {design_data[height_map].min()*1e9:.1f} - {design_data[height_map].max()*1e9:.1f} nm - 纳米柱直径: {design_data[diameter_map].min()*1e9:.1f} - {design_data[diameter_map].max()*1e9:.1f} nm - 单元周期: 400 nm (假设) 3. 制造约束 - 最小特征尺寸: {design_data[constraints_applied].get(min_diameter, N/A)*1e9:.1f} nm - 最大深宽比: {design_data[constraints_applied].get(max_aspect_ratio, N/A)} - 最小间隙: {design_data[constraints_applied].get(min_gap, N/A)*1e9:.1f} nm 4. 工艺步骤建议 a) 基底准备: 清洗硅片生长{self.material}薄膜至最大高度 b) 电子束光刻: 使用HSQ或PMMA抗蚀剂剂量根据特征尺寸调整 c) 显影: 使用TMAH或MIBK:IPA溶液 d) 刻蚀: ICP-RIE刻蚀{self.material}选择比 10:1 e) 去胶: O2等离子体去胶 5. 质量检查项目 - SEM检查特征尺寸和形状 - AFM检查高度均匀性 - 光学显微镜检查缺陷 - 白光干涉仪检查相位分布 6. 文件清单 - metasurface_layout.gds: GDSII布局文件 - quantized_phase.npy: 量化相位分布 - height_map.npy: 高度分布 - diameter_map.npy: 直径分布 生成时间: {datetime.datetime.now().strftime(%Y-%m-%d %H:%M:%S)} with open(os.path.join(output_dir, process_specification.txt), w) as f: f.write(spec_content) # 使用示例 if __name__ __main__: import datetime # 假设从联合优化得到连续相位分布 optimized_phase torch.randn(512, 512) * 0.5 np.pi # 模拟优化结果 # 创建转换器 converter MetasurfaceManufacturingConverter( phase_levels8, wavelength532e-9, materialSiN ) # 定义制造约束 manufacturing_constraints { min_diameter: 100e-9, # 100nm最小直径 max_aspect_ratio: 4.0, # 深宽比≤4 min_gap: 30e-9, # 30nm最小间隙 diameter_quantization_step: 10e-9 # 10nm直径步长 } # 执行转换 design_data converter.convert( optimized_phase, constraintsmanufacturing_constraints, apply_opcTrue ) # 导出设计文件 converter.export_design(design_data, ./metasurface_design/)三、先进制造技术考虑3.1 多层与三维超表面结构对于更复杂的相位分布可能需要多层或三维超表面结构。class MultiLayerMetasurfaceConverter: 处理多层超表面结构的转换。 参考来源中基于深度学习的生成式设计方法。 def __init__(self, num_layers3, materials(SiN, SiO2, SiN)): self.num_layers num_layers self.materials materials self.layer_converters [ MetasurfaceManufacturingConverter(phase_levels4, materialmat) for mat in materials ] def convert_multilayer(self, continuous_phase_vol): 转换多层相位分布。 参数: continuous_phase_vol: 多层连续相位形状为 (num_layers, H, W) 返回: layer_designs: 每层的设计数据列表 layer_designs [] total_phase torch.zeros_like(continuous_phase_vol[0]) for i in range(self.num_layers): print(f处理第 {i1}/{self.num_layers} 层...) # 当前层的目标相位总相位的分配 # 简化均匀分配总相位 target_phase continuous_phase_vol[i] # 使用单层转换器 design self.layer_converters[i].convert( target_phase, constraints{min_diameter: 120e-9, max_aspect_ratio: 3.0}, apply_opcTrue ) layer_designs.append(design) total_phase design[quantized_phase] return layer_designs, total_phase3.2 制造误差建模与鲁棒性优化在实际制造前可以通过仿真预测制造误差的影响并在设计中加入鲁棒性约束。def model_manufacturing_errors(design, error_params): 模拟制造误差对性能的影响。 参数: design: 设计数据字典 error_params: 误差参数包括 diameter_error_std: 直径误差标准差 height_error_std: 高度误差标准差 misalignment_error: 层间对准误差 返回: degraded_performance: 性能评估指标 # 添加随机误差到直径和高度 diameter_map design[diameter_map] height_map design[height_map] # 直径误差正态分布 if diameter_error_std in error_params: diameter_error torch.randn_like(diameter_map) * error_params[diameter_error_std] diameter_with_error diameter_map diameter_error # 确保非负 diameter_with_error torch.clamp(diameter_with_error, min50e-9) else: diameter_with_error diameter_map # 高度误差 if height_error_std in error_params: height_error torch.randn_like(height_map) * error_params[height_error_std] height_with_error height_map height_error height_with_error torch.clamp(height_with_error, min100e-9) else: height_with_error height_map # 计算误差后的相位使用单元库的逆映射 # 这里简化假设相位与高度/直径有简单线性关系 phase_error torch.zeros_like(diameter_map) # 实际中应使用从结构到相位的精确映射通过插值单元库 for phase_val, params in design[unit_cell_library].items(): target_height params[height] target_diameter params[diameter] # 计算当前结构与目标结构的差异 height_diff height_with_error - target_height diameter_diff diameter_with_error - target_diameter # 简化的误差模型相位误差与尺寸误差成正比 phase_error (height_diff.abs() / target_height diameter_diff.abs() / target_diameter) * 0.1 # 评估性能下降例如Strehl ratio下降 # Strehl ratio ≈ exp(-(2π/λ * phase_error_std)^2) phase_error_std phase_error.std() wavelength 532e-9 strehl_ratio torch.exp(-(2*np.pi/wavelength * phase_error_std)**2) return { diameter_with_error: diameter_with_error, height_with_error: height_with_error, phase_error: phase_error, strehl_ratio: strehl_ratio, estimated_efficiency_loss: 1 - strehl_ratio.item() } # 误差分析示例 error_params { diameter_error_std: 5e-9, # 5nm直径误差 height_error_std: 10e-9, # 10nm高度误差 } error_analysis model_manufacturing_errors(design_data, error_params) print(f制造误差导致的斯特列尔比下降: {error_analysis[strehl_ratio].item():.4f}) print(f估计效率损失: {error_analysis[estimated_efficiency_loss]*100:.2f}%)四、总结与最佳实践将联合优化的相位分布转化为可制造的超表面结构需要系统性的流程相位离散化根据制造能力确定相位级别数平衡性能与工艺复杂度 。单元库映射建立精确的相位-结构映射关系可通过FDTD/RCWA仿真或深度学习代理模型实现 。制造约束集成在设计阶段就考虑最小特征尺寸、深宽比、材料限制等工艺约束 。邻近效应修正使用OPC或基于机器学习的方法补偿光刻过程中的失真。误差分析与鲁棒性设计通过蒙特卡洛仿真评估制造误差的影响并在优化中加入鲁棒性约束 。最佳实践建议设计-制造协同与工艺工程师密切合作了解产线的具体能力和限制。渐进式验证先制造小尺寸测试结构验证单元性能再逐步扩大阵列尺寸。多物理场仿真考虑热效应、机械应力等对光学性能的影响。自动化流程将上述流程集成为自动化脚本支持快速设计迭代。通过上述方法可以将联合优化得到的理论相位分布可靠地转化为实际可制造的超表面结构实现从数字设计到物理器件的闭环 。参考来源人工智能与光学系统的深度融合大模型在光学设计与成像中的应用~[超表面论文快讯-156]Advanced Science-36通道自旋与波长协同复用超表面全息-浦项科技大学Junsuk Rho团队玻璃AI飞秒激光雕刻的存算一体光学神经网络原理与应用深度学习在超构表面设计中的四大神经网络应用【信息科学与工程学】【物理/化学科学和工程技术】第八篇 光学02基于条件变分自编码器的信息超材料生成式设计原理、实现与应用