matlab模拟塞曼图谱,塞曼效应以及能级的计算

简单总结了一下在原子结构的基础上的Zeeman效应。很久之前就知道，这次算是复习，顺便计算一下。仅限于LS耦合，并且假设核自旋是0. 后面可能会考虑上核自旋非零的原子。

弱场下原子的Zeeman效应

把电子运动视为经典带电小球的圆周运动，按照电磁学，得到的磁矩为$\boldsymbol{\mu}=q/(2m_e)\boldsymbol{\ell}$，其中$q=-e$为电子电荷，$m_e$为电子质量。该磁矩还可以改写为$\boldsymbol{\mu}=q\hbar/(2m_e)\boldsymbol{\ell}/\hbar=-\mu_B\boldsymbol{\ell}/\hbar$ ，其中$\mu_B$为玻尔磁子，取正值。于是相互作用能量(选取外磁场沿$z$轴)为$H'=-\boldsymbol{B}\cdot\boldsymbol{\mu}=eB_z/(2m_e)\ell_z$。这是从经典图像中得到的结果。按照B.卡尼亚克《原子物理学(下册)》，直接分析静磁场中单电子的Schrodinger方程，可以得到两个附加项，一项为$H'=-\boldsymbol{B}\cdot\boldsymbol{\mu}=eB_z/(2m_e)\ell_z$，另一项正比于磁场平方。这两项分别为顺磁项和抗磁项。简单的计算表明，对于$1T$的磁场而言，第二项远小于第一项，因此略去。不论从哪里出发，得到的附加哈密顿是相同的。在多电子情形下即为 $$ H'=\frac{\mu_B}{\hbar}\boldsymbol{B}\cdot\sum_i \boldsymbol{\ell}_i $$ 进一步考虑各个电子的自旋，最终得到 $$ H'=\frac{\mu_B}{\hbar}\boldsymbol{B}\cdot\sum_i (\boldsymbol{\ell}i+g\text{s}\boldsymbol{s}i)=\frac{\mu_B}{\hbar}\vec{B}\cdot(\vec{L}+g_s\vec{S}) $$ 其中$g\text{s}$为电子自旋的Lande因子，值为2. 在Slater的Quantum Theory of Atomic Structure中，从磁场中单电子的Dirac方程出发，计算电流密度，分别得到电子轨道角动量贡献的磁矩(和这里一样)以及自旋贡献的磁矩，其中自旋磁矩含有值为2的Lande因子。也可以说电子轨道磁矩的Lande因子为1. 至此已经给出了磁场中哈密顿的附加项的来源，并说清了电子Lande因子的来源。

现在看LS耦合情形下，弱场Zeeman效应。弱场意味着磁场的附加哈密顿是所有微扰中最小的，因此最后考虑。如此考虑的必要条件是，磁场带来的精细结构态内部的能级分裂远小于不同精细结构态之间的能级间隔。LS耦合下一个精细结构子空间由$J^2$和$J_z$的共同本征态张成，它们有共同的$J$量子数，不同的$m_{J}$量子数，是简并的。限定在一个精细结构子空间内讨论，这样，基矢为${|E^0_J,J,m_J\rangle}$，其中$E^0_J$是该精细结构能级的简并能量，其实这个位置本来应该是包括$L,S$值在内的其他量子数。

由于Wigner-Eckart定理，两个同秩的不可约张量算符的各个分量在角动量基矢$|\alpha,J,m_J\rangle$上的矩阵元的比值相等，即(假设$k$秩) $$ \frac{\langle\alpha'J'm_J'|T^{(k)}_q|\alpha Jm_J\rangle}{\langle\alpha'J'm_J'|U^{(k)}_q|\alpha Jm_J\rangle}=C $$ 其中常数$C$和$q,m_J,m'J$无关。下面选取$U^{(k)}q$为一秩的总角动量算符$J_q$, 则有 $$ \langle\alpha'J'm_J'|\boldsymbol{A}|\alpha Jm_J\rangle=C\langle\alpha'J'm_J'|\boldsymbol{J}|\alpha Jm_J\rangle $$ 其中$\hat{\boldsymbol{A}}$是某矢量算符。因为