简介这份资源面向光学、磁光材料与物理仿真方向的学习者和研究人员围绕掺铈钇铁石榴石Ce:YIG这一典型磁光晶体提供一维磁光透射、反射与法拉第旋转效应的数值模拟代码。压缩包内共1个文件为MATLAB脚本.m格式整体约1KB可直接运行以计算不同磁场与波长条件下晶体的透射系数、反射系数及法拉第旋转角帮助理解偏振光在磁光介质中的传播规律。目前已有438人学习下载说明该方向具备一定关注度。通过这份脚本读者可快速复现Ce:YIG晶体的磁光响应曲线掌握透射反射系数的建模思路并在此基础上调整参数、拓展至磁光隔离器或调制器的性能评估适合作为磁光效应入门仿真与课程实验的参考素材。1. Ce_YIG 磁光晶体从透射谱到法拉第旋转一条能跑通的表征链路如果你手头有一片 Ce 掺杂 YIGCe:YIG薄膜或者正在调研一维磁光晶体的透射与法拉第旋转大概率会遇到一个很实际的问题同一片样品透射谱测出来挺漂亮法拉第旋转角却怎么都对不上文献值。这不是玄学而是磁光表征里最常见的翻车点——透射和偏振旋转本来就是两套耦合但不等价的测量逻辑。Ce_YIG 这类铈掺杂钇铁石榴石因为 Ce³⁺ 的强自旋-轨道耦合在 1550 nm 通信波段有远高于纯 YIG 的法拉第旋转能力同时吸收也明显上升所以透射和反射数据必须一起看。这篇笔记面向做磁光薄膜、集成光隔离器、一维磁光子晶体的实验和仿真人员把透射谱、法拉第旋转角、椭偏率这几件事的测量与建模链路讲清楚让你能照着搭一套可复现的流程而不是只停留在“测了个谱”。2. 先搞懂 Ce_YIG 的磁光响应到底由什么决定2.1 透射、反射、法拉第旋转三者的物理耦合关系在磁光介质中介电张量可以写成对角项 ε 和非对角项 ε_xy 的形式非对角项正比于磁化强度 M。法拉第旋转角 θ_F 在薄样品近似下正比于 ε_xy 的实部乘以传播距离而透射率 T 同时受 ε 的虚部吸收和界面反射影响。这意味着你测到的透射下降可能来自吸收增大也可能来自反射增强两者对法拉第旋转的贡献完全不同。Ce 掺杂会同时抬高 ε_xy 和吸收系数所以不能只看透射高低判断磁光性能好坏。常见做法是同时采集透射谱 T(λ) 和反射谱 R(λ)用 T R A 1 先做一次能量守恒校验。如果 T R 明显小于 1 且随波长变化剧烈说明吸收或散射不可忽略后续拟合必须把消光系数 k 作为自由参数放进去而不是固定成文献值。2.2 Ce 掺杂浓度与一维磁光晶体周期数的选型逻辑Ce:YIG 的 Ce 掺杂量通常在 x 0.51.5 之间掺杂越高法拉第旋转越大但晶格失配和吸收也越严重。一维磁光晶体比如 YIG/SiO₂ 多层膜的周期数 N 决定了透射带内的相位积累和磁光增强因子。N 太少法拉第旋转积累不够N 太多吸收和界面散射会把透射压到探测器噪声以下。我一般会先做一组 N 4、8、12 的对照样品固定总厚度只改周期数测透射带边和法拉第旋转谱。如果 N 12 时透射在 1550 nm 已经低于 1%那再增加周期数就没有工程意义应该回头调单层厚度或换更低吸收的缓冲层。2.3 测量前必须固定的三个系统参数光源偏振消光比至少 1000:1否则法拉第旋转的零点会漂。样品法线对准入射角偏差 0.5° 就能让反射谱出现干涉条纹偏移。磁场方向与光传播方向平行度法拉第旋转对磁场角度是余弦依赖偏 5° 损失约 0.4% 信号但如果是反射式测量误差会翻倍。提示先把这三个参数用标准样品比如纯 YIG 或已知旋转角的 TGG标定一遍再上 Ce_YIG否则后面拟合出来的 ε_xy 没有可比性。3. 透射谱与法拉第旋转的实测流程从光路到数据3.1 搭建透射-反射同步测量光路的最小命令级步骤下面是一套我常用的透射/反射同步采集流程用 Python 做数据记录和初步处理。硬件上需要可调谐激光器、偏振控制器、分束器、两个光电探测器和一个电磁铁。代码只负责采集和存盘不涉及具体仪器驱动你可以把read_power()替换成自己的 SDK 调用。import numpy as np import time # 波长扫描范围1450-1650 nm步进 0.5 nm wavelengths np.arange(1450, 1650.5, 0.5) T_data [] R_data [] for wl in wavelengths: set_laser_wavelength(wl) # 替换为你的激光器控制函数 time.sleep(0.05) # 等激光稳定 p_in read_power(input) # 入射光功率 p_t read_power(trans) # 透射光功率 p_r read_power(refl) # 反射光功率 T_data.append(p_t / p_in) R_data.append(p_r / p_in) np.savetxt(ce_yig_T_R.csv, np.column_stack([wavelengths, T_data, R_data]), delimiter,, headerwavelength_nm,T,R, comments)逻辑说明每个波长点先读入射功率做归一化避免激光器功率漂移影响透射率。参数说明步进 0.5 nm 是为了捕捉 Ce:YIG 在 1550 nm 附近的吸收边细节如果只关心带内旋转可以放宽到 2 nm。time.sleep(0.05)是给激光器波长切换留稳定时间具体值看你的激光器手册。3.2 法拉第旋转角的偏振调制测量与符号约定法拉第旋转不能直接用两个偏振片夹样品测因为样品本身的线性双折射会混进去。常见做法是用光弹调制器PEM加锁相放大器或者用旋转检偏器法做二次谐波检测。下面是一个旋转检偏器法的简化数据处理代码假设你已经采到了一系列检偏角下的光强。import numpy as np # 检偏器角度度和对应光强 angles np.arange(0, 360, 5) intensities np.array([...]) # 你的实测光强 # 用 cos(2θ) 拟合提取旋转角 from scipy.optimize import curve_fit def fit_func(theta, I0, theta0, offset): return I0 * np.cos(2 * np.deg2rad(theta - theta0)) offset popt, _ curve_fit(fit_func, angles, intensities) theta0 popt[1] # 这个就是偏振主轴方向 # 法拉第旋转角 theta0(加磁场) - theta0(零磁场)逻辑说明旋转检偏器法通过拟合光强随角度的 cos(2θ) 变化来提取偏振主轴比直接找极小值点更抗噪声。参数说明theta0是拟合出的主轴角加磁场和零磁场各测一次差值就是法拉第旋转角。注意符号约定如果磁场方向与光传播方向相同Ce:YIG 的旋转通常是负的具体看你的坐标系定义建议在论文里明确写清楚。3.3 从 T/R 数据反推吸收系数与磁光优值拿到 T 和 R 之后可以用下面的近似公式反推吸收系数 α 和磁光优值 FOM θ_F / α。对于弱吸收薄膜import numpy as np # 假设已知薄膜厚度 dcmT 和 R 是小数 d 500e-7 # 500 nm T np.array(T_data) R np.array(R_data) # 忽略干涉的近似A 1 - T - R A 1 - T - R alpha -np.log(T) / d # 粗略吸收系数cm^-1 # 磁光优值需要法拉第旋转角 theta_F单位度 theta_F 200 # 示例值度 FOM theta_F / alpha逻辑说明A 1 - T - R是最粗的能量守恒校验如果 A 出现负值说明反射测量有系统误差。参数说明d必须用实际薄膜厚度不能用法拉第旋转测量中的光程代替。alpha这里忽略了反射修正适合吸收较大的 Ce:YIG如果吸收很小需要用更完整的传输矩阵法。4. 用传输矩阵法把一维磁光晶体的透射和旋转算准4.1 磁光传输矩阵的构建与四个必调参数一维磁光晶体的透射和法拉第旋转可以用 4×4 传输矩阵精确计算。对每一层需要输入介电常数对角项 ε、非对角项 ε_xy、厚度 d、入射角 θ。下面是一个单层磁光介质的 4×4 矩阵构建代码。import numpy as np def mag_optical_layer(eps, eps_xy, d, theta, wavelength): 构建单层磁光介质的 4x4 传输矩阵 eps: 对角介电常数 eps_xy: 非对角项决定法拉第旋转 d: 厚度单位 m theta: 入射角弧度 wavelength: 波长单位 m k0 2 * np.pi / wavelength n np.sqrt(eps) # 简化正入射情况斜入射需要更完整的 Berreman 矩阵 delta k0 * n * d # 这里只展示正入射的简化形式 M np.eye(4, dtypecomplex) M[0,0] np.cos(delta) M[0,1] 1j * np.sin(delta) / n M[1,0] 1j * n * np.sin(delta) M[1,1] np.cos(delta) # 磁光耦合项 M[0,2] -eps_xy * np.sin(delta) / (n * eps) M[2,0] eps_xy * np.sin(delta) / (n * eps) return M逻辑说明正入射下磁光耦合主要体现在 x 和 y 偏振的交叉项。参数说明eps_xy是决定法拉第旋转大小的关键Ce:YIG 在 1550 nm 典型值在 0.010.05 之间具体取决于 Ce 掺杂量和磁化状态。d用米wavelength也用米避免单位混用。4.2 多层堆叠与透射率、法拉第旋转角的提取把每层的矩阵乘起来得到总矩阵然后从总矩阵的透射系数提取透射率和偏振旋转。def stack_transmission(layers, wavelength, theta0): layers: 列表每个元素是 (eps, eps_xy, d) M_total np.eye(4, dtypecomplex) for eps, eps_xy, d in layers: M_total M_total mag_optical_layer(eps, eps_xy, d, theta, wavelength) # 提取透射系数简化只看前两个模式 t_xx 1 / M_total[0,0] t_xy -M_total[0,2] / M_total[0,0] T np.abs(t_xx)**2 np.abs(t_xy)**2 # 法拉第旋转角 theta_F np.arctan2(np.abs(t_xy), np.abs(t_xx)) return T, np.rad2deg(theta_F)逻辑说明t_xx是共偏振透射t_xy是交叉偏振透射法拉第旋转角由两者的比值决定。参数说明theta默认正入射如果要做斜入射需要换成完整的 Berreman 4×4 矩阵这里不展开。layers列表里每一层的eps_xy符号要一致否则多层之间的磁光耦合会相互抵消。4.3 拟合实测数据时 ε_xy 和吸收系数的边界条件拟合时最容易翻车的是把 ε_xy 和吸收系数同时放开导致多解。我一般会固定吸收系数从 T/R 数据反推的值只拟合 ε_xy 的实部和虚部。如果透射带内出现法诺共振说明一维磁光晶体的周期结构引起了模式耦合这时候需要把层厚也作为弱拟合参数但步长不要超过 2 nm。注意Ce:YIG 的 ε_xy 对温度敏感室温到 80°C 之间可能变化 10% 以上。如果你在变温环境下测记得把温度作为协变量记录。5. 避坑与排查Ce_YIG 磁光表征里最容易翻车的五件事5.1 透射谱出现负吸收反射测量没校准现象算出来的 A 1 - T - R 出现负值或者吸收系数为负。原因反射测量没有扣除背景反射或者分束器引入了额外损耗。解决用已知反射率的标准镜比如金镜标定反射通道把分束器损耗单独测出来做修正。5.2 法拉第旋转角随波长振荡干涉效应没扣除现象θ_F 随波长出现周期性振荡周期和薄膜厚度对应。原因薄膜前后表面的多次反射形成了法布里-珀罗干涉偏振旋转被调制。解决要么用楔形样品或折射率匹配液消除背面反射要么在拟合时把干涉项一起建模不要直接对振荡曲线取平均。5.3 磁场加反向时旋转角不反号剩磁或磁光饱和现象正反向磁场下法拉第旋转角差不多没有反号。原因电磁铁剩磁太大或者样品已经磁光饱和外加磁场不足以翻转磁化。解决先用高斯计测样品位置的剩磁如果剩磁超过 10 Oe需要加退磁流程如果是饱和增大磁场或换更薄的样品。5.4 一维磁光晶体透射带边偏移层厚误差累积现象仿真透射带边和实测差 10 nm 以上。原因多层膜的每层厚度误差在堆叠后累积尤其是 SiO₂ 和 YIG 的溅射速率漂移。解决每层镀完用椭偏仪测一次厚度不要只信工艺参数如果已经镀完用透射谱反推总光学厚度做修正。5.5 拟合出的 ε_xy 比文献大一个量级单位或符号搞错现象拟合得到的 ε_xy 是 0.5 而不是 0.05。原因厚度单位用了 nm 但公式里按 m 算或者把度当成弧度代入了。解决把所有输入统一到 SI 单位角度全部转弧度拟合前先用一组已知文献值做正向验证。6. 进阶技巧用透射极值点快速估算法拉第旋转的符号和量级如果你不想每次都跑完整拟合可以用透射谱的极值点位置来快速估算法拉第旋转。一维磁光晶体在带边附近透射极小值对应的波长会随磁场移动移动量 Δλ 和 ε_xy 近似成正比。下面是一个快速估算的代码片段。import numpy as np from scipy.signal import find_peaks # 假设有两组透射谱零磁场和加磁场 wl np.array(wavelengths) T_zero np.array(T_data_zero) T_field np.array(T_data_field) # 找极小值点 peaks_zero, _ find_peaks(-T_zero) peaks_field, _ find_peaks(-T_field) # 取带边附近第一个极小值 idx_zero peaks_zero[0] idx_field peaks_field[0] delta_lambda wl[idx_field] - wl[idx_zero] # 经验系数 k 需要用自己的样品标定一次 k 0.5 # 示例值单位 nm^-1 theta_F_est k * delta_lambda print(f估算法拉第旋转角: {theta_F_est:.2f} 度)逻辑说明透射极小值移动量 Δλ 和磁光耦合强度在带边附近近似线性系数 k 需要用一组已知旋转角的样品标定。参数说明find_peaks(-T)是找透射极小值因为find_peaks默认找极大值。k的典型值在 0.30.8 nm^-1 之间取决于一维磁光晶体的周期结构。这个方法只能给量级和符号不能替代完整拟合但适合在镀膜后快速筛选样品。我自己的习惯是每批样品先测透射谱用极值点法筛一遍符号不对或量级差太多的直接重镀不浪费机时做法拉第旋转的精细测量。这套流程帮我省过至少两周的无效拟合时间。希望帮到你。本文还有配套的精品资源点击获取
