1. 非厄米超表面与异常点（EP）的物理基础
在传统光学系统中，厄米性（Hermiticity）往往被视为基本要求——这意味着系统能量守恒且时间反演对称。然而非厄米系统的出现彻底打破了这一认知框架。这类系统允许增益（gain）和损耗（loss）的存在，其哈密顿量不再满足H=H†的条件。这种非厄米特性催生出许多反直觉的物理现象，其中最引人注目的就是异常点（Exceptional Points, EP）。
异常点是参数空间中两个或多个本征态及其对应本征值发生简并的特殊点。与厄米系统中的简并点不同，EP处的本征态不仅本征值相同，本征向量也会坍缩到同一方向。这种独特的性质使得系统在EP附近表现出：
- 非互易的光学响应
- 偏振态的突变转换
- 对微扰的超高灵敏度
2021年Science报道的工作正是利用超表面（metasurface）这种二维人工结构，通过精心设计的非厄米单元实现了可控的EP。超表面相比传统体材料具有显著优势：
- 亚波长尺度单元可实现局域场增强
- 几何参数灵活调控偏振响应
- 平面结构与现有光电器件兼容
2. 超表面单元设计与FDTD建模要点
2.1 双共振元结构设计
实现EP的关键在于构造具有特定非对称性的传输矩阵。文献中采用的双共振元结构包含两个核心组件：
-
x偏振共振单元：
- 材料：金（Au）纳米棒
- 尺寸：100nm×30nm×50nm（长×宽×高）
- 功能：主要响应x方向线偏振光
-
y偏振耦合单元：
- 相同材料但旋转45度放置
- 尺寸调整为30nm×100nm×50nm
- 通过旋转引入与x偏振的耦合
lua复制-- Lumerical FDTD单元建模示例
set("x span", 300e-9); -- 单元周期300nm
set("y span", 300e-9);
-- x偏振纳米棒
addrect(
"name"="x_rod",
"material"="Au",
"x span"=100e-9,
"y span"=30e-9,
"z span"=50e-9,
"x"=0,
"y"=0
);
-- 倾斜45度的y偏振棒
addrect(
"name"="y_rod",
"material"="Au",
"x span"=30e-9,
"y span"=100e-9,
"z span"=50e-9,
"rotation z"=45,
"x"=50e-9,
"y"=50e-9
);
2.2 非厄米损耗的引入
为实现非厄米特性，需要在系统中引入可控损耗：
- 使用SiO2作为基底材料
- 在特定区域掺入吸收材料（如石墨烯量子点）
- 通过虚部介电常数量化损耗强度
lua复制-- 定义损耗材料
addmaterial(
"name"="SiO2_Absorber",
"Imaginary permittivity"=0.1 -- 损耗强度参数
);
-- 包裹纳米棒的损耗层
addrect(
"material"="SiO2_Absorber",
"x span"=150e-9,
"y span"=150e-9,
"z span"=60e-9
);
关键参数说明：
- 旋转角度：控制偏振耦合强度（通常40-50度最优）
- 损耗系数：虚部介电常数在0.05-0.3范围可调
- 单元间距：影响近场耦合效率
3. 传输矩阵提取与本征值分析
3.1 Jones矩阵的仿真获取
在FDTD中设置斜入射（如15度）的线偏振光源，通过频域监视器记录透射场。对x和y偏振分别仿真，提取透射场的复数振幅构成Jones矩阵：
| 入射偏振 | Ex透射 | Ey透射 |
|---|---|---|
| x偏振 | Txx | Tyx |
| y偏振 | Txy | Tyy |
完整的传输矩阵T可表示为：
[ T = \begin{bmatrix}
T_{xx} & T_{xy}\
T_{yx} & T_{yy}
\end{bmatrix} ]
3.2 本征值计算的MATLAB实现
将仿真数据导出至MATLAB进行特征分析：
matlab复制% 波长扫描范围（单位：nm）
lambda_range = 600:5:800;
num_points = length(lambda_range);
% 预分配存储空间
eigenvalues = zeros(2, num_points);
eigenvectors = zeros(2, 2, num_points);
for i = 1:num_points
% 从FDTD导出数据构建传输矩阵
T = [Txx(i), Txy(i);
Tyx(i), Tyy(i)];
% 归一化处理（消除相位不确定性）
T = T / sqrt(det(T));
% 特征分解
[V, D] = eig(T);
eigenvalues(:, i) = diag(D);
eigenvectors(:, :, i) = V;
end
% EP判定条件
diff = abs(eigenvalues(1,:) - eigenvalues(2,:));
[~, ep_index] = min(diff);
ep_wavelength = lambda_range(ep_index);
3.3 EP点的实验验证特征
当系统参数接近EP时，会出现以下典型现象：
- 本征值轨迹交叉：实部和虚部同时重合
- 本征态坍缩：两个本征向量趋同
- 偏振响应突变：如线偏振→圆偏振转换
图：本征值在复平面的运动轨迹，红色星号标记EP点
4. 实操中的关键问题与解决方案
4.1 FDTD仿真优化技巧
| 问题类型 | 解决方案 | 参数建议 |
|---|---|---|
| 收敛问题 | 增加网格精度 | 纳米棒边缘网格≤2nm |
| 数值噪声 | 使用共形网格 | Conformal Meshing等级≥1 |
| 耗时过长 | 采用对称边界 | 对于对称结构可节省50%时间 |
4.2 实验数据处理的注意事项
-
相位校准：
- 使用参考臂校正干涉相位
- 建议采用Mach-Zehnder干涉仪结构
-
噪声抑制：
matlab复制% 奇异值降噪示例 [U,S,V] = svd(T_matrix); S(3:end,3:end) = 0; % 保留前两个奇异值 T_denoised = U*S*V'; -
EP点精确定位：
- 采用双参数扫描（如波长+损耗系数）
- 使用最小二乘法拟合交叉区域
5. 应用前景与扩展方向
基于EP的非厄米超表面在以下领域展现独特优势：
-
超灵敏传感：
- 折射率灵敏度可达1000nm/RIU
- 比传统表面等离子体共振高1-2个数量级
-
偏振调控器件：
- 可编程偏振转换器
- 动态可调光学隔离器
-
量子模拟平台：
- 模拟PT对称系统
- 研究拓扑能带结构
实验组最新进展表明，通过引入相变材料（如GST），可实现EP的动态调控。这种可重构非厄米超表面的调制速度已达GHz级别，为下一代集成光子器件提供了新思路。
6. 复现研究的完整流程总结
-
结构设计阶段：
- 确定目标波长范围（如通信波段1550nm）
- 根据材料色散曲线优化几何参数
-
仿真验证阶段：
- 参数扫描（旋转角、损耗系数等）
- 收敛性测试（网格独立性验证）
-
数据分析阶段：
- 传输矩阵提取与归一化
- 本征值轨迹绘制与EP定位
-
实验验证阶段（可选）：
- 电子束光刻制备样品
- 椭偏仪测量偏振响应
对于初次尝试的研究者，建议先从可见光波段（如630-800nm）开始，该范围内金的光学性质稳定且易于加工。一个典型的项目周期约为2-3周，其中FDTD仿真耗时约占60%。