当前位置：首页 > news >正文

用Python和FDTD仿真，手把手教你理解超表面中的几何相位与传输相位

news 2026/5/26 11:53:44

用Python和FDTD仿真解析超表面中的几何相位与传输相位

当一束光穿过由硅纳米柱构成的超表面时，它的相位会发生什么变化？这个问题困扰着许多初入光子学领域的研究者。在实验室里，我们常常看到导师用复杂的数学公式解释相位调控原理，但真正的理解往往始于亲手建模仿真的那一刻。本文将带你用Python和FDTD工具，从电磁场模拟的角度直观感受几何相位与传输相位的本质区别。

1. 超表面相位调控基础：从理论到仿真准备

超表面作为二维形式的超材料，其核心在于通过亚波长结构单元调控电磁波的相位分布。在可见光和近红外波段，硅纳米结构因其高折射率和低损耗特性成为理想选择。我们以一个简单的硅纳米椭圆柱阵列为例，其典型尺寸为长轴300nm、短轴200nm、高度850nm，放置在二氧化硅基底上。

传输相位的本质是光波在不同等效折射率介质中传播产生的相位延迟。当线偏振光入射到各向异性纳米结构时，x和y偏振分量会经历不同的相位变化。这可以通过Python的meep库进行仿真验证：

import meep as mp import numpy as np # 设置硅纳米椭圆柱参数 a = 0.3 # 长轴(μm) b = 0.2 # 短轴(μm) h = 0.85 # 高度(μm)

几何相位则与结构旋转角度相关，仅对圆偏振光有效。Berry相位理论指出，当圆偏振光被旋转的结构散射时，会产生与旋转角度成比例的附加相位。这种相位与材料色散无关，为宽带调控提供了可能。

相位类型	适用偏振	调控机制	色散特性	实现方式
传输相位	线偏振	等效折射率	强色散	改变结构尺寸
几何相位	圆偏振	旋转角度	弱色散	旋转相同结构

2. FDTD仿真环境搭建与纳米结构建模

Lumerical FDTD是超表面仿真的黄金标准工具，但开源方案如MEEP同样强大。我们先配置仿真环境：

# 安装MEEP和依赖 pip install pymeep pymeep-extras

建立纳米椭圆柱模型需要明确定义材料属性和几何参数。硅在近红外波段的折射率约为3.5，二氧化硅为1.45。在FDTD中，我们通常设置：

# 材料定义 si = mp.Medium(index=3.5) sio2 = mp.Medium(index=1.45) # 仿真区域设置 cell_size = mp.Vector3(1,1,2) geometry = [mp.Block(size=mp.Vector3(mp.inf,mp.inf,0.5), center=mp.Vector3(0,0,-0.25), material=sio2), mp.Ellipsoid(size=mp.Vector3(a,b,h), material=si)]

关键仿真参数包括：

网格精度：通常设为λ/20
边界条件：PML吸收边界
光源：平面波，波长范围800-1600nm
监视器：透射场分析

提示：在椭圆柱旋转模拟中，建议采用参数化扫描方式，每次旋转后重新计算透射场，避免引入误差。

3. 偏振态与琼斯矩阵分析实战

理解相位调控必须掌握琼斯矩阵方法。对于线偏振入射光，超表面可表示为：

J = [t_x*exp(iφ_x) 0 0 t_y*exp(iφ_y)]

而当结构旋转θ角时，琼斯矩阵变为：

def rotated_jones(theta, tx, ty, phix, phiy): R = np.array([[np.cos(theta), -np.sin(theta)], [np.sin(theta), np.cos(theta)]]) J = np.diag([tx*np.exp(1j*phix), ty*np.exp(1j*phiy)]) return R @ J @ R.T

通过FDTD仿真提取透射场后，我们可以计算琼斯矩阵元素：

分别用x和y线偏振光入射
记录透射场的Ex和Ey分量
计算各矩阵元：
- t_x = |E_x(x-pol入射)|/|E_inc|
- φ_x = angle(E_x(x-pol入射))
- t_y = |E_y(y-pol入射)|/|E_inc|
- φ_y = angle(E_y(y-pol入射))

对于圆偏振光，几何相位表现为正交偏振分量的相位差。例如左旋圆偏振(LCP)入射时，右旋(RCP)分量将携带2θ的几何相位：

E_out = t*J_0*E_LCP + t*J_1*exp(i2θ)*E_RCP

4. 仿真结果可视化与相位特性对比

完成仿真后，我们需要提取并分析数据。以下Python代码演示了如何计算和绘制相位分布：

import matplotlib.pyplot as plt # 提取传输相位 phix = np.angle(Ex_xpol) # x偏振入射时的x分量相位 phiy = np.angle(Ey_ypol) # y偏振入射时的y分量相位 # 计算几何相位 theta = np.linspace(0, np.pi, 10) # 旋转角度 phi_geo = [] # 存储几何相位 for t in theta: # 运行旋转后的仿真 E_out = simulate_rotated_structure(t) phi_geo.append(np.angle(E_out['rcp']))

典型结果分析应包括：