简介这份PDF文档围绕塞曼效应整理了经典理论、量子力学与相对论三种主流解释方法适合大学物理、光学或原子物理方向的师生及科研入门者使用。文档逐层梳理了Maxwell方程式与Lorentz力、Schrödinger方程式与Dirac方程式、狭义相对论与广义相对论等不同切入视角并统一推导出能级分裂公式ΔE±μB再通过对比三者的数学形式、物理假设与应用边界帮助读者跳出单一教材框架建立对磁场中原子光谱分裂机理的整体认识。全文共1个PDF文件压缩包大小1.31MB排版紧凑、内容集中便于课余或考前快速翻阅。资源已有111人学习/下载适合想通过横向比较掌握塞曼效应本质、并了解磁共振成像、光谱分析等实际应用的读者。1. 塞曼效应的三种解释法公式长得一样物理差一个 g 因子在互联网上检索“塞曼效应”相关 PDF最常见的教学版本会给出一个非常讨巧的结论经典理论、量子力学、相对论都能推出 ΔE ±μB。这个结论不是错而是省略了太多前提。经典解释里只讨论束缚电子在磁场中的回旋对应频率 eB/2m量子解释中多电子原子的能级分裂是 g_J m_J μ_B B朗德 g 因子可以等于 1、2、4/3、2/3到了狄拉克方程自旋磁矩的“2”自动出现精细结构和反常磁矩又对高 Z 或强场区域形成修正。三种方法对应三种计算模型各自有适用范围和失效边界。下面围绕一份网传对照讲义展开把三个模型的推导补齐并加入可直接运行的 Python 校验代码适合需要把物理讲义转成可复现计算脚本的工程师和教研人员。2. 经典理论解释从洛伦兹力到拉莫尔三重线先算清磁场能劈多少2.1 束缚电子在均匀磁场中的经典运动方程经典模型要把原子当作一个三维谐振子不能只写一句“由洛伦兹力推出 ΔE±μB”。取电子被势场束缚时恢复力为 -mω0²r外磁场 B 沿 z 方向电子电荷为 -e则运动方程是m r -mω0²r - e(v×B)。把 x 和 y 方向分量组合成复位移 u x iy在 B (0,0,B) 下可化为u i(eB/m)u ω0²u 0。设解为 u ~ exp(-iωt)特征方程是 ω² (eB/m)ω - ω0² 0。弱磁场下取正频率根得到ω± ω0 ± eB/(2m)。也就是说原来简并的圆偏振模式被拆成两个频率z 方向振动频率仍为 ω0。这三个本征模式对应谱线中心分量和左右两条谱线这是经典塞曼三重线的真正来源。这里 eB/(2m) 定义为拉莫尔频率 ωL它不是一个拟合参数而是由运动方程直接产生的特征量。这个推导在日常教学中经常被跳过导致读者只知道“经典解释得到正常塞曼效应”却不知道如何从运动方程得到三条线也不知道为什么两条 σ 分量是圆偏振。实际上用复位移把二维振动方程解耦是本问题最直观的数学工具。2.2 频移量级磁场多大才值得上实验设备拉莫尔频率与磁场强度成正比B 1 T 时ωL/(2π) eB/(4πm) ≈ 14.0 GHz。光学跃迁频率通常在 10¹⁴10¹⁶ Hz 量级所以塞曼分裂在频率轴上是窄谱线尺度上的微扰实验需要高分辨率干涉仪。工程上更常用波长偏移来估算Δλ λ²Δν/c。下面给出一组可直接对照的参数表。B(T)拉莫尔频率(GHz)500 nm 波长偏移(pm)典型用途0.11.4≈1.2教学演示需高分辨法珀1.014.0≈11.7常规磁场下的可测分裂5.070.0≈58.5精细结构不干扰时清晰分辨10.0140.0≈117.0偏向强场区需要检查帕邢-巴克效应表中最右列提醒你经典三重线只在磁场弱、能级结构简并且无自旋时严格成立。磁场一旦大到位移接近自旋轨道裂距就必须切换到量子模型。表中多普勒展宽没有画进去但实际对比时要把谱线线型叠加进去。2.3 用 Python 验证经典特征频率直接解运动方程再做傅里叶变换比套公式更能确认“±eB/2m”这个系数。为了数值稳定这里使用无量纲化令 ω01定义 alpha eB/(mω0)。对于电子洛伦兹力会产生速度耦合方程在无量纲时间下只有 alpha 一个参数。import numpy as np from scipy.integrate import solve_ivp # alpha eB/(m*omega0)这里取 0.1远小于 1对应弱磁场 alpha 0.1 def rhs(t, y): # 状态向量 [x, y, z, vx, vy, vz] x, y_, z, vx, vy, vz y # B 沿 z电子电荷为负洛伦兹力项为 -alpha*vy, alpha*vx ax -x - alpha*vy ay -y_ alpha*vx az -z return [vx, vy, vz, ax, ay, az] t_eval np.linspace(0, 2000, 20000) sol solve_ivp(rhs, [0, 2000], [1.0, 0.0, 0.0, 0.0, 0.0, 0.0], t_evalt_eval, rtol1e-9) freq np.fft.rfftfreq(len(t_eval), t_eval[1]-t_eval[0]) amp np.abs(np.fft.rfft(sol.y[0])) peaks np.argsort(amp)[-3:] print(主要特征频率:, np.sort(freq[peaks]))代码把 x 方向的振动轨迹做 FFT峰值会落在 1±alpha/2 附近。alpha 是磁场耦合强度与恢复力强度之比实际原子中这个量接近 10⁻⁵ 量级用 0.1 只是为了在数值图上明显显示两个峰。solve_ivp 的优势是可以同时保留线偏振与圆偏振的耦合信息而不必像解析方法那样先假设圆偏振基。2.4 经典模型的适用边界经典模型能够解释正常塞曼效应的三条线但它有一个内在缺陷电子自旋产生的固有磁矩完全没进模型。如果假设电子带着角动量做回旋运动经典回磁比必然是 e/(2m)对应 g1而实验和量子力学都要求电子的自旋回磁比 g≈2。多电子原子中轨道磁矩和自旋磁矩还要耦合经典谐振子没有办法描述朗德因子出现 2/3、4/3 这类分数的机制。因此经典解释的适用区域基本停留在自由电子回旋运动、等离子体磁光效应、以及 S0 的单重态跃迁。用来估算磁场下的谱线移动可以用来做高分辨原子光谱标定就不够。这是阅读“塞曼效应三种解释法”这类资料时首先要建立的边界意识。3. 量子力学解释泡利哈密顿量、朗德 g 因子和谱线合成3.1 磁场如何进入电子哈密顿量量子力学处理不能把洛伦兹力直接写成势能项而是要通过最小替换 p → p - qA 把磁矢势引入动能算符。取均匀磁场 B 沿 z 方向常用对称规范 A (-By/2, Bx/2, 0)展开动能项并加入自旋磁矩后非相对论泡利哈密顿量中的磁耦合部分为H μ_B B·(L 2S)/ħ (e²/8m)(B×r)²。第一项是顺磁塞曼项第二项是抗磁项。在常见磁场下抗磁项只是把所有能级整体移动一个零点几 GHz 量级的微小常数不会改变谱线结构所以谱线分裂主要看第一项。式子里 L 和 S 都是算符因此能级位移依赖于量子态。原文只说“由 Schrödinger 方程推出 ΔE±μB”严格讲只对单电子自旋态或极简单的 S0 态成立。真正要计算的是矩阵元ΔE μ_B B ⟨L_z 2S_z⟩/ħ。对于 LS 耦合的多电子原子这个平均值可以被朗德 g 因子吸收。3.2 朗德 g 因子的矢量模型推导在 LS 耦合下好量子数是 L、S、J、m_J。算符 L 2S 在 J 子空间内与 J 平行其投影写成L 2S g_J J。将等式两边与 J 做内积并用角动量代数可得到g_J 1 [J(J1) S(S1) - L(L1)] / [2J(J1)]。这时一阶塞曼能移为ΔE g_J m_J μ_B B。m_J 从 -J 到 J 共 2J1 个值因此每条能级裂成等间距的一组子能级。间距不是 μ_B B而是 g_J μ_B B。正常塞曼效应对应 g_J1反常塞曼效应绝大多数情况是 g_J≠1所以谱线成分数量比经典三重线多。这里有个容易踩的坑J0 时公式分母为零。实际上 J0 态没有角动量投影磁场线性项的平均值为零处理代码时需要做保护直接返回 1 只适用于一些特殊约定的占位。3.3 常见 LS 态 g 因子速查与代码实现手算 g 因子很机械适合做成函数。下表列出几个教科书态方便实现时验证。LSJ谱项g_J2J101/21/2²S₁/₂2211/21/2²P₁/₂2/3211/23/2²P₃/₂4/34111³P₁3/23112³P₂3/25213³D₃4/37下面是计算 g 因子和跃迁分量的 Python 代码。它把所有允许的 m 组合列出来并用选择定则过滤输出每个分量相对跃迁中心频率的位移。import numpy as np H 6.62607015e-34 MU_B 9.2740100783e-24 def lande_g(L, S, J): LS 耦合朗德 g 因子J0 时做占位处理 if J 0: return 1.0 return 1.0 (J*(J1) S*(S1) - L*(L1)) / (2.0*J*(J1)) def zeeman_spectrum(L1,S1,J1,L2,S2,J2,B): 返回上态 g、下态 g以及满足电偶极选择定则的跃迁频率位移 g1 lande_g(L1,S1,J1) g2 lande_g(L2,S2,J2) lines [] for m1 in np.arange(-J1, J11, 1): for m2 in np.arange(-J2, J21, 1): if abs(m1 - m2) 1: continue dE (g1*m1 - g2*m2) * MU_B * B dnu dE / H lines.append((m1, m2, dnu)) return g1, g2, lines g1, g2, lines zeeman_spectrum( L11, S10.5, J11.5, L20, S20.5, J20.5, B1.0 ) print(g1 , g1, g2 , g2) for line in lines[:8]: print(line)函数中 L1、S1、J1 是上能级量子数L2、S2、J2 是下能级量子数。输出 dnu 的单位是 Hz。选择定则过滤了 Δm 1 的组合保留 Δm0 和 ±1 三条链。你只要把上下态换成实际的谱项就能直接用这份输出画合成谱线。3.4 与经典解释的异同经典模型只有一个拉莫尔频率谱线结构固定为等间距三条量子模型则多了 g_J 和 m_J 两个自由度。当上态和下态 g 因子相同时子能级等距所有跃迁合并成三条线当上下态 g 因子不同不同 m 组合产生的频率位移不再等距谱线数量也会更多。实验上看到“多于三条”的塞曼分量说明经典模型已经失效需要用本节代码逐条列出。量子模型没有解释“为什么电子自旋磁矩是 2μ_B”这个“2”要到狄拉克方程里才自然出现。4. 相对论解释狄拉克方程给出自旋 g2但不要扯上广义相对论4.1 自由电子的 g 因子不是拟合值原子物理的常见表述是自旋磁矩等于 2μ_B S/ħ。这个 2 可以通过狄拉克方程在低能极限下得到。狄拉克方程本身是洛伦兹协变的把静磁场放进协变导数后重新构造正负能分量低能有效哈密顿量中会出现H_z eħ/(2m) σ·B / 2?严格写法是取自旋三分量的能量算符最后等价于自由电子的自旋磁矩算符为 g_s μ_B S_z/ħ其中 g_s2、S_z±ħ/2。所以自旋上态和自旋下态的能量差为ΔE g_s μ_B B 2μ_B B。单个子能级相对无磁场的位移为 ±μ_B B。原始资料把这一步也写成“ΔE±μB”不算全错但没讲清楚这个符号只描述自由电子自旋的子能级移动而两条子能级之间实际差 2μ_B B。计算微波共振频率时若用错系数结果会差一倍。电子自旋共振实验里B1 T 对应共振频率约 28 GHz而不是 14 GHz。4.2 相对论动能修正与精细结构相对论效应还体现在电子动能的色散关系上。将 E √(p²c² m²c⁴) 展开可以得到E ≈ mc² p²/(2m) - p⁴/(8m³c²) ...其中 p⁴ 项使能级略微下移形成氢原子的精细结构。自旋轨道耦合项同量级。对氢原子 2P₁/₂ 和 2P₃/₂精细结构裂距约为 10⁻⁵ eV换算成频率约 10 GHz正好和强场下的塞曼分裂同一尺度。这就是为什么在相对论修正不能被忽略的态上不能把朗德 g 公式当作恒定的参数。普通 LS 耦合计算把 g 因子当成纯角动量投影系数没有考虑相对论径向波函数的变形。高 Z 原子中轨道波函数向内收缩自旋轨道耦合和动能修正都会影响磁矩矩阵元单纯按 L、S、J 算出的 g 因子会偏离实验值。这时要使用 Breit-Pauli 哈密顿量做矩阵对角化或者在更高精度下用 Dirac-Coulomb 方法直接处理相对论波函数。4.3 强场区域与帕邢-巴克极限磁场强度增加后另一个判断标准是 μ_B B 与精细结构裂距 ΔE_fs 的相对大小。μ_B B 远小于 ΔE_fs 时J 是好量子数朗德 g 公式有效μ_B B 远大于 ΔE_fs 时L 和 S 分别绕磁场进动J 不再是好量子数能级重新按 m_L 和 m_S 排布这种现象叫帕邢-巴克效应。在帕邢-巴克区域朗德 g 公式失效但哈密顿量依然可以精确写出只是需要在 m_L 和 m_S 表象里对角化。对单电子原子能级位移近似为ΔE μ_B B (m_L 2m_S)。这是经典模型、朗德模型之外的第三个极限。实验上中等磁场往往处在两者之间最稳妥的方法是直接构造 Zeeman 加自旋轨道的矩阵通过数值对角化得到本征值而不是套用某一端的近似公式。4.4 数值核对相对论动能展开与 g 因子工程上常用 SymPy 展开相对论能量用来快速核对精细结构修正的作用面。下面代码输出动能级到 p⁴ 项。import sympy as sp m, c, p sp.symbols(m c p, positiveTrue) E sp.sqrt(p**2*c**2 m**2*c**4) print(sp.series(E, p, 0, 5))输出为 mc² p²/(2m) - p⁴/(8c²m³) O(p⁶)。第一项是静能第二项是非相对论动能第三项是相对论修正。这一项会让氢原子能级向低能方向移动移动幅度随主量子数 n 增大而减小。再核对自由电子的 g2。电子质量与电荷常量已知B1 T 时能量差为q_e 1.602176634e-19 hbar 1.054571817e-34 m_e 9.1093837015e-31 B 1.0 mu_B q_e * hbar / (2*m_e) delta_E 2 * mu_B * B print(delta_E / q_e, eV)结果约为 1.16×10⁻⁴ eV对应 28 GHz。把 g2 换成 g1得到的 14 GHz 不满足电子自旋共振实验这足以说明狄拉克方程中“2”的必要性。4.5 三种模型对照表模型核心对象能级位移关键参数适用边界经典谐振子 洛伦兹力束缚电子坐标 r(t)±ħωLω0、BS0 或自由电子谱线三重泡利/朗德量子力学角动量算符 L、S、Jg_J m_J μ_B BL、S、J、g_JLS 耦合B 小于精细结构裂距狄拉克方程/相对论修正旋量波函数(g_s m_s ...) μ_B Bg_s2、自旋轨道等高 Z、强场、覆盖反常磁矩修正表中后两行都能推出“±μB”形式的子能级位移但一个含 g_J一个含 g_s 与 m_s。使用资料时只要看公式前面有没有 g 因子和对应的量子数就能判断这段文字到底在讲哪种解释。5. 把三种解释法收敛到频率同一尺度一个方便自查的计算技巧拿到一个塞曼效应资料最实用的校验手段是先把所有模型换算成同一个频率尺度每特斯拉磁场对应的频率是多少。定义系数 k μ_B B/h单位是 GHz/T。这样经典模型的自旋子能级位移约 ±k B狄拉克自由电子自旋上下态能级差是 2k B朗德模型则是 g_J m_J k B。三者之间的差异全部集中在这个 k 前面的系数上。下面这段代码把朗德 g 因子、磁场和磁量子数统一包装成一个函数输出以 GHz 为单位的频率位移。它适合用于对比文献里的谱线图。def zeeman_ghz(g_factor, B, mJ): 把 g*mu_B*mJ*B 换算成 GHz muB_h 13.996245 # mu_B/h 单位 GHz/T return g_factor * muB_h * B * mJ # 对比自由电子 g2 与经典 g1 print(zeeman_ghz(2.0, 1.0, 0.5)) # 13.996 GHz 子能级下移/上移 print(zeeman_ghz(1.0, 1.0, 0.5)) # 6.998 GHz经典回旋频率减半 print(zeeman_ghz(4/3, 1.0, 0.5)) # 反常塞曼态 ²P3/2 的值这里 muB_h 是约化普朗克常量取法定单位后的结果。使用该函数时注意 mJ 必须是能级的磁量子数投影不能拿上下态混合使用。例如 ²P₃/₂ 的 mJ1/2 和 ²S₁/₂ 的 mJ1/2 属于不同能级它们之间没有直接位移替代关系必须分别计算上下态差异后相减。把这种方法固定为工作流能快速识别网传资料中的公式错位。一个典型的错误是经典拉莫尔频率和量子朗德频率混用最后实验谱线半宽对不上。你在整理讲义或做课程设计时只要把每个模型输出统一标为“GHz 1 T”对照表就可以直接复用到任意磁场强度。这个技巧同样适合 MRI 和谱学分析里快速估算谱线分离度不需要每次重新翻教材。本文还有配套的精品资源点击获取