简介这份资源是发表于《大学物理》的学术论文《磁各向异性介质中的平面电磁波》PDF适合电磁场理论、通信技术及物理学相关专业的研发人员和学生阅读参考。文中基于Maxwell方程系统推导了磁各向异性理想介质中平面电磁波的完整方程组梳理了电场、磁场、波矢量与能流方向之间的关系并引入对偶关系导出磁晶体中的“菲涅耳方程”进而分析平面波的结构、传播及其偏振特性。内容还涉及能量密度分布、相速度与光线速度的关系为理解磁各向异性介质中的电磁波行为提供了较完整的理论基础。资源共1个PDF文件压缩包大小139KB轻量便于下载目前已有83人学习。对于从事电磁波传播研究、各向异性介质建模及相关技术开发的人员这篇论文具有较强的专业参考价值。1. 磁各向异性介质中的平面电磁波为什么它不是换个折射率的事磁各向异性介质最常见的是外加磁场下的旋磁铁氧体以及各类磁化薄膜和多层膜。这里磁感应强度 B 与磁场强度 H 不再平行B μ₀μ̿ᵣHμ̿ᵣ 是一个含非对角元的 3×3 张量。平面波在这样的介质里传播时同一波矢量方向往往对应两支不同折射率的波偏振态也由介质张量决定而不是由入射条件随便指定。微波移相器、环行器、隔离器的厚度计算材料测量时的 S 参数反演都要先把这个“波矢到折射率”的关系写对。适合的读者是想在仿真中准确放入铁氧体张量、或正为各向异性反演参数发愁的工程师和研究人员。2. 磁各向异性介质的本构关系与特征方程先写对 Maxwell 方程组再谈波2.1 无源区平面波表示与时间约定考虑无源、线性、时不变介质频率域 Maxwell 方程组为∇×E −jωB∇×H jωD本构关系写成D ε₀ε̿ᵣEB μ₀μ̿ᵣH这里 ε̿ᵣ 和 μ̿ᵣ 都可以是张量。磁各向异性问题的“分量”在 μ̿ᵣ下文先假设介电常数为标量 εᵣ再在实现中保留张量形式。平面波解取 e^{j(ωt−k·r)} 形式。令 k nk₀s其中 s 为单位波矢方向k₀ ω√(μ₀ε₀)n 是待求的复折射率。把 ∇ 换成 −jk由两个旋度方程可以直接推出两个常用关系s×E Z₀μ̿ᵣH / ns×H −ε̿ᵣE / (nZ₀)其中 Z₀ 是自由空间波阻抗。消去 H得到关于电场强度 E 的特征方程n²s×(μ̿ᵣ⁻¹s×E) ε̿ᵣE 0写成矩阵形式M E (1/n²)EM −ε̿ᵣ⁻¹[s]×μ̿ᵣ⁻¹[s]×其中 [s]× 是叉乘矩阵0 −s_z s_y [s]× s_z 0 −s_x −s_y s_x 0时间约定必须统一。上面推导按 e^{jωt} 写如果仿真或者教材用 e^{−jωt}Polder 张量中的非对角元 jκ 要跟着变号否则解出来的旋向会反。2.2 Polder 张量从铁氧体出发的磁导率张量形式旋磁铁氧体在外加偏置磁场下磁导率张量可写成 Polder 形式。设偏置磁场沿 z 方向无损耗近似下μ jκ 0 μ̿ᵣ −jκ μ 0 0 0 1μ 和 κ 的频率依赖由旋磁共振决定μ 1 ω₀ωₘ/(ω₀²−ω²)κ ωωₘ/(ω₀²−ω²)其中 ω₀ μ₀γH₀ 是进动角频率ωₘ μ₀γMₛ 与饱和磁化强度相关。有损耗材料通常把 ω₀ 换成 ω₀ jαω 再代入得到的 μ、κ 都是复数。不同状态下的张量性质对求解方式影响很大见表 2-1。表Polder 张量状态与特征方程性质材料状态μ、κ 性质特征方程特征典型场景无损耗实 μ、实 κM 可取 Hermitiann² 为实数移相器、环行器设计频率附近的初步计算有损耗复 μ、复 κ复特征值n 为复数吸收材料、宽带测量接近共振μ 实部可能为负存在倏逝或禁带需要避开或利用共振吸收Polder 张量只在磁化方向为 z 时是上面这个形状。偏置场旋转时整个张量要经过坐标旋转常见做法是把外加磁场方向定义为局部坐标系的 z 轴算完电场后再转回全局坐标。2.3 特征方程的适用边界特征方程假设介质无空间色散即 μ̿ᵣ 只依赖频率、不依赖波矢。对于大多数宏观尺寸的铁氧体器件这个假设成立到了太赫兹或光学频段若涉及磁子-光子耦合的近场问题空间色散不能忽略本文的方程只能给出参考。另外方程只覆盖无源区域如果介质内有传导电流需要先把 J 项写回 Maxwell 方程组不能直接套用上节公式。3. 色散与偏振磁各向异性介质平面波的两类典型结果3.1 先用各向同性情况做校验把 μ̿ᵣ μI、ε̿ᵣ εᵣI 代入特征方程可以得到M (−1/εᵣ)[s]×μ⁻¹[s]×由于 μ⁻¹ 是标量[s]×[s]× ssᵀ − I特征值出现二重简并n² μεᵣ这就是最常见的各向同性折射率结果。这个退化的校验点很重要后续写代码时可以先用它确认符号和单位没有写反。3.2 纵向传播两支圆偏振与法拉第旋转波矢沿磁化方向即 s (0,0,1)。代入 Polder 张量后特征向量对应两个圆偏振E₊ ∝ (1, j, 0)ᵀE₋ ∝ (1, −j, 0)ᵀ两支波折射率分别n₊ √(εᵣ(μ κ))n₋ √(εᵣ(μ − κ))线偏振波进入纵向磁化介质后可以看成右旋和左旋圆偏振波的叠加。由于两支折射率不同经过距离 d 后叠加出的线偏振方向会旋转旋转角为θF k₀d·Re(n₊ − n₋)/2这就是 Faraday 旋转。它直观解释了为什么磁各向异性介质中不能只给一个折射率线偏振的“折射率”根本不存在只有圆偏振本征模的折射率存在。3.3 横向传播寻常波与异常波再看 s (1,0,0)即波矢垂直磁化方向。磁各向异性张量在这时候给出两支完全不同的模第一支电场沿磁化方向即 E ∝ (0,0,1)ᵀ与张量的非对角元无关n₁ √(εᵣμ_z) √εᵣ第二支电场在 xy 平面内n₂ √(εᵣ(μ²−κ²)/μ)这两支对应介质光学中很熟悉的寻常波和异常波。横向传播时会出现偏振面相关的折射率差异这和各向异性晶体中的双折射在现象上非常接近但来源是磁张量。表纵向与横向传播的解析对照εᵣ5μ1.2κ0.3传播方向本征模解析折射率数值纵向 s(0,0,1)右旋圆偏振√(5×(1.20.3))2.738纵向 s(0,0,1)左旋圆偏振√(5×(1.2−0.3))2.121横向 s(1,0,0)电场沿 z√(5×1)2.236横向 s(1,0,0)电场在 xy 面√(5×(1.2²−0.3²)/1.2)2.372这组解析值可以直接用于验证后文的数值代码。3.4 场矢量关系的易错点平面波条件下 s·D 0 和 s·B 0 始终成立但电场 E 不一定垂直于 s。在磁各向异性介质中E 与 s 的夹角由特征方程给出不能想当然地认为电磁波一定是横电波。计算坡印廷矢量时要直接用 E×H*不要先假设 E 和 H 都垂直于 s。4. 用特征值问题求解磁各向异性介质中的平面波一份可复现的 Python 实现4.1 最小实现输入张量、输出折射率和偏振把上一章的特征方程翻译成 NumPy 代码下面这个函数可以处理任意 3×3 磁导率张量和介电张量。import numpy as np def cross_mat(s): sx, sy, sz s return np.array([ [0.0, -sz, sy], [sz, 0.0, -sx], [-sy, sx, 0.0], ]) def plane_wave_modes(eps_r, mu_r, s): # 单位化传播方向 s s / np.linalg.norm(s) S cross_mat(s) # 特征方程: M E (1/n^2) E T S np.linalg.inv(mu_r) S M -np.linalg.inv(eps_r) T eigvals, eigvecs np.linalg.eig(M) n2 1.0 / eigvals n np.sqrt(n2) # 无损耗主分支; 有损耗需按传播方向统一取分支 # 由 E 恢复 H: B (n / c0) * s x E, H ... c0 2.99792458e8 mu0 4.0 * np.pi * 1e-7 modes [] for ni, ei in zip(n, eigvecs.T): Ei ei / np.linalg.norm(ei) B (ni / c0) * np.cross(s, Ei) Hi np.linalg.solve(mu_r * mu0, B) modes.append((ni, Ei, Hi)) return modes代码核心是先构造叉乘矩阵再求 M 的特征系统。特征值 λ 对应 1/n²所以取倒数后开方得到折射率。特征向量就是本征波的电场方向。H 由 B (n/c₀)s×E 恢复再代入本构关系解出 H。代码中把特征向量按列取出后转成行存储避免把电场和磁场放反。开方时使用 NumPy 的主分支无损耗介质没问题有损耗介质中 n 的实部和虚部符号需要结合波的传播方向手动修正否则会出现沿传播方向增幅的伪解。4.2 验证算例解析解对照用上一章的解析结果检验代码。写一个短脚本eps np.eye(3) * 5.0 mu_diag 1.2 kapp 0.3 mu_polder np.array([ [mu_diag, 1j * kapp, 0.0], [-1j * kapp, mu_diag, 0.0], [0.0, 0.0, 1.0], ]) # 各向同性校验 mu_iso np.eye(3) * 2.0 modes_iso plane_wave_modes(np.eye(3) * 4.0, mu_iso, np.array([1.0, 0.0, 0.0])) print([np.round(abs(m[0]), 3) for m in modes_iso]) # 应为 [2.828, 2.828,...] # 纵向传播 modes_z plane_wave_modes(eps, mu_polder, np.array([0.0, 0.0, 1.0])) for m in modes_z: print(np.round(abs(m[0]), 3)) # 2.121 / 2.738 # 横向传播 modes_x plane_wave_modes(eps, mu_polder, np.array([1.0, 0.0, 0.0])) for m in modes_x: print(np.round(abs(m[0]), 3)) # 2.236 / 2.372校验结果汇总见表 4-1。表数值计算与解析解的对照算例解析 n代码输出 n各向同性ε4μ22.8282.828纵向 Polderμ1.2κ0.32.121 / 2.7382.121 / 2.738横向 Polderμ1.2κ0.32.236 / 2.3722.236 / 2.372误差在浮点精度范围内。三个算例分别验证了退化本征值、圆偏振本征值和双折射结构覆盖了大多数工程场景。4.3 参数说明与常见误用使用这套代码时有几个隐含参数需要注意。首先是频率附在张量上不是独立变量。μ、κ 都是在特定频率下算出来的扫频时要在每个频率点重新构造 Polder 张量不能拿一个固定张量扫整个带宽。其次是材料损耗。把 ω₀ 换成 ω₀ jαω 后μ、κ 变成复数特征值也变成复数。开方时应该用复开方并对每个 n 检查传播方向的实部与衰减方向选取物理可行的一支。第三是张量的坐标关系。代码中输入的是全局坐标系下的 μ̿ᵣ。如果偏置磁场不在 z 轴需要先把 Polder 张量旋转到全局坐标常见做法是 Rμ̿ᵣRᵀ其中 R 是磁化方向与 z 轴之间的旋转矩阵。5. 界面与测量把磁各向异性平面波的结论用于实测5.1 界面上的相位匹配与四个本征模平面波遇到磁各向异性介质界面时不能只用一个斯涅尔定律。斜入射时切向波矢连续但反射波和透射波都要在当前介质内展开成本征模的叠加各向异性介质中的本征模一般不是纯 TE/TM反射系数需要在偏振之间耦合。常见做法是用 Berreman 4×4 传递矩阵把所有场的切向分量组成状态向量逐层相乘得到反射透射系数。若只是垂直入射单层板可以用层内两个本征模分别处理但入射偏振要与本征模匹配才能得到单一模的传播否则会出现模式混叠。5.2 用 Faraday 旋转加总相移反推 μ 和 κ纵向构型下测量线偏振波穿过厚度 d 的平板后得到两个量极化面旋转角 θF 和总相移 ψ。二者与本征折射率的关系是θF k₀d·Re(n₊ − n₋)/2ψ k₀d·Re(n₊ n₋)/2于是可以直接反解Re(n₊) (ψ θF)/(k₀d)Re(n₋) (ψ − θF)/(k₀d)已知介电常数 εᵣ 后得到 μ κ Re(n₊)²/εᵣμ − κ Re(n₋)²/εᵣ再解出 μ、κ。需要注意 ψ 包含 2π 相位模糊测量前要先用低损耗薄样品估计厚度电长度避免把 n₊ 和 n₋ 错位。5.3 实测前的三个验证习惯先用非磁性样板校准法拉第旋转为零确认系统没有残余旋光多频率点联合反演看 μ 和 κ 是否落在 Polder 模型的平滑曲线上如果单点频率反演结果跳动剧烈多半是相位模糊或张量坐标没旋转再用横向传播的寻常波折射率 √εᵣ 标定介电常数减少同时拟合多个参数的耦合。最后一个习惯在测量中特别实用μ 和 κ 对相位测量误差的敏感度不同κ 主要影响两支模之差先固定总相移减掉 εᵣ 贡献再反演 κ会稳定得多。本文还有配套的精品资源点击获取