单向辐射ugr模型 包含单向辐射电场模，上下表面辐射损耗，能带，q因字。
单向辐射UGR模型在光子学领域最近存在感越来越强。这玩意儿本质上是个光场调控系统，核心在于让电磁波像装了单向阀一样只能往特定方向传播。咱们今天拆解下它的电场分布特征和损耗机制，顺手用Python撸几个关键参数的计算。
先看电场模的分布规律。这类结构通常采用非对称光栅设计，上表面做亚波长级刻蚀，下表面保持平整。用有限差分法模拟时要注意边界条件设置：
import numpy as np from scipy.sparse import diags def solve_waveguide(dx=0.02, n_layers=200): epsilon = np.ones(n_layers) epsilon[:50] = 3.5**2 # 上表面高折射率层 epsilon[-50:] = 1.0**2 # 空气层 # 构建波动方程矩阵 main_diag = -2/dx**2 + (2*np.pi/1.55)**2 * epsilon off_diag = 1/dx**2 * np.ones(n_layers-1) A = diags([main_diag, off_diag, off_diag], [0, -1, 1]) # 求解本征模 eigvals, eigvecs = np.linalg.eigh(A.toarray()) return eigvals, eigvecs这段代码用五点差分格式离散了波动方程。注意第50层设置的突变界面，这会直接影响模式泄露。跑出来的本征值虚部直接对应辐射损耗——这个量在普通波导里可不存在。
单向辐射ugr模型 包含单向辐射电场模，上下表面辐射损耗，能带，q因字。
辐射Q因子计算是验证设计的关键。传统方法是算模式能量衰减：
def calc_Q_factor(eigenvalues): omega_real = np.real(eigenvalues[0]) # 基模频率 omega_imag = np.imag(eigenvalues[0]) # 衰减系数 return 0.5 * omega_real / abs(omega_imag)但单向辐射结构得考虑上下表面损耗不对称性。实际工程中会引入权重因子α修正公式：Q_corrected = Q/(1-α)，其中α≈0.7时模型最贴近实测数据。
能带分析更带劲——在k空间观察禁带特性时，会发现单向辐射对应的奇异点（Exceptional Point）。用平面波展开法算能带时要注意基矢选取：
% 非对称光栅的傅里叶系数计算 for n = -N:N eps_G(n+N+1) = 1/period * integral(@(x) eps_real(x).*exp(-1i*2*pi*n*x/period),0,period); end这里如果傅里叶级数截断阶数N不够，会漏掉高阶衍射对单向性的影响。建议至少取到N=15，否则禁带边缘会出现鬼影峰。
最后说个坑：仿真时网格尺寸如果大于结构特征尺寸的1/8，Q因子计算结果会严重失真。但网格太细又容易算到天荒地老，建议用自适应网格加密，重点区域局部细化。这平衡点的把握，大概就是理论和实操的距离吧。
