简介：这份资源面向光学、磁光材料与物理仿真方向的学习者和研究人员，围绕掺铈钇铁石榴石（Ce:YIG）这一典型磁光晶体，提供一维磁光透射、反射与法拉第旋转效应的数值模拟代码。压缩包内共1个文件，为MATLAB脚本（.m格式），整体约1KB，可直接运行以计算不同磁场与波长条件下晶体的透射系数、反射系数及法拉第旋转角，帮助理解偏振光在磁光介质中的传播规律。目前已有438人学习下载，说明该方向具备一定关注度。通过这份脚本，读者可快速复现Ce:YIG晶体的磁光响应曲线，掌握透射反射系数的建模思路，并在此基础上调整参数、拓展至磁光隔离器或调制器的性能评估，适合作为磁光效应入门仿真与课程实验的参考素材。
1. Ce_YIG 磁光晶体：从透射谱到法拉第旋转，一条能跑通的表征链路
如果你手头有一片 Ce 掺杂 YIG（Ce: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.5～1.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=',', header='wavelength_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 # 假设已知薄膜厚度 d（cm），T 和 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, dtype=complex) 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.01～0.05 之间，具体取决于 Ce 掺杂量和磁化状态。d用米，wavelength也用米，避免单位混用。
4.2 多层堆叠与透射率、法拉第旋转角的提取
把每层的矩阵乘起来得到总矩阵，然后从总矩阵的透射系数提取透射率和偏振旋转。
def stack_transmission(layers, wavelength, theta=0): """ layers: 列表，每个元素是 (eps, eps_xy, d) """ M_total = np.eye(4, dtype=complex) 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.3～0.8 nm^-1 之间，取决于一维磁光晶体的周期结构。这个方法只能给量级和符号，不能替代完整拟合，但适合在镀膜后快速筛选样品。
我自己的习惯是：每批样品先测透射谱，用极值点法筛一遍，符号不对或量级差太多的直接重镀，不浪费机时做法拉第旋转的精细测量。这套流程帮我省过至少两周的无效拟合时间。希望帮到你。
本文还有配套的精品资源，点击获取
