下面取最常见的等量标量-矢量耦合：
\[ S(\mathbf r)=V(\mathbf r)=\frac12 m\omega^2 r^2 \]
此时狄拉克方程可写为：
\[ \left[\boldsymbol{\alpha}\cdot\mathbf p+\beta\big(m+S\big)\right]\psi=(E-V)\psi \]
在自然单位 \(\hbar=c=1\) 下，令 \(S=V=\frac12 m\omega^2 r^2\)，则：
\[ \left[\boldsymbol{\alpha}\cdot\mathbf p+\beta m+\frac12 m\omega^2 r^2(\beta+1)\right]\psi=E\psi \]
关键简化
取狄拉克矩阵的表示：
\[ \beta= \begin{pmatrix} I&0\ 0&-I \end{pmatrix},\qquad \boldsymbol{\alpha}= \begin{pmatrix} 0&\boldsymbol{\sigma}\ \boldsymbol{\sigma}&0 \end{pmatrix} \]
把四分量旋量写成：
\[ \psi= \begin{pmatrix} \phi\ \chi \end{pmatrix} \]
其中 \(\phi,\chi\) 是二分量旋量。代入后得到耦合方程：
\[ \boldsymbol{\sigma}\cdot\mathbf p,\chi + m\phi + \frac12 m\omega^2 r^2(1+1)\phi = E\phi \]
\[ \boldsymbol{\sigma}\cdot\mathbf p,\phi - m\chi + \frac12 m\omega^2 r^2(-1+1)\chi = E\chi \]
注意第二式中 \((\beta+1)\) 作用在 \(\chi\) 上给出零，所以势项消失：
\[ \boldsymbol{\sigma}\cdot\mathbf p,\phi - m\chi = E\chi \]
即：
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
代回第一式：
\[ \boldsymbol{\sigma}\cdot\mathbf p\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi + m\phi + m\omega^2 r^2\phi = E\phi \]
利用：
\[ (\boldsymbol{\sigma}\cdot\mathbf p)^2=\mathbf p^2 \]
得到：
\[ \frac{\mathbf p^2}{E+m}\phi + m\phi + m\omega^2 r^2\phi = E\phi \]
整理：
\[ \left[\frac{\mathbf p^2}{2m}+\frac12 m\omega^2 r^2\right]\phi= \frac{E^2-m^2}{2m(E+m)}\phi \]
但右边可化为：
\[ \frac{E^2-m^2}{2m(E+m)}=\frac{E-m}{2m} \]
所以：
\[ \left[\frac{\mathbf p^2}{2m}+\frac12 m\omega^2 r^2\right]\phi= \frac{E-m}{2m}\phi \]
化为非相对论三维谐振子
令：
\[ \frac{E-m}{2m}=\left(N+\frac32\right)\omega \]
则：
\[ E=m+2m\omega\left(N+\frac32\right) \]
但这不是正确结果，因为上面推导中势项系数有问题。更标准的等量标量-矢量耦合下，直角坐标分离会得到：
\[ E^2-m^2=2m\omega\left(2n_x+2n_y+2n_z+3\right) \]
即：
\[ E^2=m^2+2m\omega\left(2N+3\right) \]
其中：
\[ N=n_x+n_y+n_z,\qquad n_x,n_y,n_z=0,1,2,\dots \]
所以：
\[ E_N=\pm\sqrt{m^2+2m\omega(2N+3)} \]
直角坐标下的波函数结构
令：
\[ \psi_{n_xn_yn_z}(\mathbf r)= \begin{pmatrix} \phi_{n_xn_yn_z}\ \chi_{n_xn_yn_z} \end{pmatrix} \]
其中 \(\phi\) 满足三维非相对论谐振子方程：
\[ \left[\frac{\mathbf p^2}{2m}+\frac12 m\omega^2 r^2\right]\phi= \left(N+\frac32\right)\omega,\phi \]
所以：
\[ \phi_{n_xn_yn_z}(x,y,z)= \varphi_{n_x}(x)\varphi_{n_y}(y)\varphi_{n_z}(z) \]
\(\varphi_n\) 是一维谐振子本征函数：
\[ \varphi_n(x)= \frac{1}{\sqrt{2^n n!}}\left(\frac{m\omega}{\pi}\right)^{1/4} e^{-m\omega x^2/2}H_n(\sqrt{m\omega},x) \]
而 \(\chi\) 由：
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
给出。因为 \(\mathbf p=-i\nabla\)，所以 \(\chi\) 可通过作用微分算符得到，仍可分离变量。
能级与简并度
能级：
\[ E_N=\pm\sqrt{m^2+2m\omega(2N+3)},\qquad N=0,1,2,\dots \]
正能态简并度：
\[ g_N=(N+1)(N+2) \]
这与球坐标下等量标量-矢量耦合的结果一致，因为物理谱不依赖于坐标系选择。
与非相对论直角坐标对比
- 非相对论：\(E_N=(N+\frac32)\omega\)，简并度 \((N+1)(N+2)/2\)
- 狄拉克等量耦合：\(E_N=\pm\sqrt{m^2+2m\omega(2N+3)}\)，正能态简并度 \((N+1)(N+2)\)
简并度翻倍来自狄拉克旋量的自旋结构。
如果你愿意，我可以继续把直角坐标下 \(\chi\) 的具体表达式完整写出来，包括用一维升降算符作用在 \(\varphi_n\) 上的形式。
==================== 你说得对，我前面那段推导确实在关键一步把系数写错了，并且后面给出的“能谱”实际上是**狄拉克谐振子（Dirac oscillator）**的谱，而不是 \(S=V=\frac12 m\omega^2 r^2\) 这个势的谱。下面按你的指正把正确逻辑整理清楚。
1. 方程与耦合方程
自然单位 \(\hbar=c=1\) 下，取
\[ S(r)=V(r)=\frac12 m\omega^2 r^2 \]
狄拉克方程为
\[ \bigl[\boldsymbol{\alpha}\cdot\mathbf p+\beta(m+S)\bigr]\psi=(E-V)\psi \]
即
\[ \Bigl[\boldsymbol{\alpha}\cdot\mathbf p+\beta m+\frac12 m\omega^2 r^2(\beta+1)\Bigr]\psi=E\psi \]
令 \(\psi=\begin{pmatrix}\phi\ \chi\end{pmatrix}\)，\(\beta+1\) 只作用在 \(\phi\) 上，得到
\[ \boldsymbol{\sigma}\cdot\mathbf p,\chi+m\phi+m\omega^2 r^2\phi=E\phi \]
\[ \boldsymbol{\sigma}\cdot\mathbf p,\phi-m\chi=E\chi \]
第二式给出
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
代回第一式，利用 \((\boldsymbol{\sigma}\cdot\mathbf p)^2=\mathbf p^2\)：
\[ \frac{\mathbf p^2}{E+m}\phi+m\omega^2 r^2\phi=(E-m)\phi \]
两边乘 \(E+m\)：
\[ \bigl[\mathbf p^2+m\omega^2(E+m)r^2\bigr]\phi=(E^2-m^2)\phi \]
这是正确的直角坐标方程。
2. 我前面的错误
我把它写成了
\[ \Bigl[\frac{\mathbf p^2}{2m}+\frac12 m\omega^2 r^2\Bigr]\phi=\frac{E-m}{2m}\phi \]
这是错的。正确形式是
\[ \frac{\mathbf p^2}{2m}\phi+\frac12\bigl[\omega^2(E+m)\bigr]r^2\phi =\frac{E^2-m^2}{2m}\phi \]
关键在于：\(r^2\) 项系数是 \(m\omega^2(E+m)\)，不是 \(m\omega^2\)。因此有效频率
\[ \Omega(E)=\omega\sqrt{E+m} \]
依赖能量 \(E\)，不能直接套非相对论谐振子的
\[ E_{\rm nr}=\bigl(N+\tfrac32\bigr)\omega \]
去解 \(E\)。
3. 正确的能谱
方程
\[ \bigl[-\nabla^2+\beta^2 r^2\bigr]\phi=\lambda\phi \]
其中
\[ \beta^2=m\omega^2(E+m),\qquad \lambda=E^2-m^2 \]
这是标准的三维各向同性谐振子形式，其本征值为
\[ \lambda=2\beta\bigl(N+\tfrac32\bigr),\qquad N=n_x+n_y+n_z \]
所以
\[ E^2-m^2=2\omega\sqrt{m(E+m)},\bigl(N+\tfrac32\bigr) \]
即
\[ (E-m)^2(E+m)=m\omega^2(2N+3)^2 \]
正能解由这个三次方程给出，不是简单的平方根。
4. 非相对论极限
令 \(E=m+\varepsilon\)，\(|\varepsilon|\ll m\)。则
\[ E+m\approx 2m,\qquad E-m=\varepsilon \]
三次方程化为
\[ \varepsilon^2\cdot 2m\approx m\omega^2(2N+3)^2 \]
所以
\[ \varepsilon\approx \bigl(N+\tfrac32\bigr)\omega \]
回到普通三维谐振子能级。这说明耦合定义本身的非相对论极限是对的。
5. 波函数
大分量仍可分离：
\[ \phi_{n_xn_yn_z}(\mathbf r)=\varphi_{n_x}(x)\varphi_{n_y}(y)\varphi_{n_z}(z) \]
但长度标度依赖能量：
\[ \varphi_n(\xi)\propto e^{-\alpha\xi^2/2}H_n(\sqrt{\alpha},\xi), \qquad \alpha=\omega\sqrt{m(E+m)} \]
小分量
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
用一维升降算符写出：\(p_x\) 作用在 \(\varphi_{n_x}\) 上只混合 \(n_x\pm1\)，所以 \(\chi\) 仍可分离变量，但每个自旋分量是若干相邻量子数项的线性组合。
6. 简并度
同一 \(N\) 的轨道简并与非相对论相同：
\[ \frac{(N+1)(N+2)}{2} \]
再乘自旋 \(2\)，正能态简并度为
\[ (N+1)(N+2) \]
这一点我前面说“自旋使简并翻倍”在计数上没错，但不能用来补救错误的能级公式。
7. 与我前面公式的关系
我前面写的
\[ E_N=\pm\sqrt{m^2+2m\omega(2N+3)} \]
对应的是 Moshinsky 狄拉克谐振子
\[ \mathbf p\to\mathbf p-im\omega\beta\mathbf r \]
在不计自旋–轨道分裂时的量级，不是
\[ S=V=\frac12 m\omega^2 r^2 \]
这个势的谱。两者不能混为一谈。
结论：你指出的错误是对的。正确能谱由
\[ (E-m)^2(E+m)=m\omega^2(2N+3)^2 \]
给出，非相对论极限才回到
\[ E\approx m+\bigl(N+\tfrac32\bigr)\omega \]
而狄拉克谐振子的简单平方根谱不适用于这个 \(S=V\) 等量耦合模型。
\[ \bigl[\mathbf p^2+m\omega^2(E+m)r^2\bigr]\phi=(E^2-m^2)\phi \]
等价地写成：
\[ \left[-\nabla^2+m\omega^2(E+m)r^2\right]\phi=(E^2-m^2)\phi \]
其中 \(\phi\) 是狄拉克旋量的大分量（二分量旋量），\(E\) 是能量本征值。
推导回顾
从耦合方程
\[ \boldsymbol{\sigma}\cdot\mathbf p,\chi+m\phi+m\omega^2 r^2\phi=E\phi \]
\[ \boldsymbol{\sigma}\cdot\mathbf p,\phi-m\chi=E\chi \]
第二式给出
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
代回第一式：
\[ \boldsymbol{\sigma}\cdot\mathbf p\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi +m\omega^2 r^2\phi=(E-m)\phi \]
利用
\[ (\boldsymbol{\sigma}\cdot\mathbf p)^2=\mathbf p^2 \]
得
\[ \frac{\mathbf p^2}{E+m}\phi+m\omega^2 r^2\phi=(E-m)\phi \]
两边乘 \(E+m\)：
\[ \bigl[\mathbf p^2+m\omega^2(E+m)r^2\bigr]\phi=(E^2-m^2)\phi \]
注意
- 这是一个能量依赖的谐振子方程：等效频率 \(\Omega(E)=\omega\sqrt{E+m}\) 含未知本征值 \(E\)，所以不是普通线性本征值问题，需要自洽求解。
- 小分量由 \(\chi=\dfrac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi\) 确定，不是独立自由度。
- 该方程就是 \(S=V=\frac12 m\omega^2 r^2\) 等量标量–矢量耦合下，直角坐标分离后大分量满足的最终二分量方程。
1. 对比非相对论情形
非相对论薛定谔方程：
\[ \left[\frac{\mathbf p^2}{2m}+U(r)\right]\phi=E_{\rm nr}\phi \]
势能 \(U(r)\) 直接以标量形式加到方程左边，系数与能量无关。
而狄拉克方程里，\(S=V=\frac12 m\omega^2 r^2\) 时，矢量势 \(V\) 耦合方式是：
\[ \bigl[\boldsymbol{\alpha}\cdot\mathbf p+\beta(m+S)\bigr]\psi=(E-V)\psi \]
移项后：
\[ \bigl[\boldsymbol{\alpha}\cdot\mathbf p+\beta m+\beta S+V\bigr]\psi=E\psi \]
注意：矢量势 \(V\) 是和 \(E\) 在同一侧出现的，它本质上是被吸收进能量 \(E\) 里的。所以当 \(S=V\) 时，\(\beta S+V\) 对上分量贡献是 \(S+V=2S=m\omega^2 r^2\)，对下分量贡献是 \(-S+V=0\)。
这就是为什么大分量方程里会出现 \(m\omega^2 r^2\) 这一项，而它和 \(E\) 一起决定了有效频率。
2. 为什么会出现 \(E\) 依赖
从耦合方程：
\[ \boldsymbol{\sigma}\cdot\mathbf p,\chi+m\phi+m\omega^2 r^2\phi=E\phi \]
\[ \boldsymbol{\sigma}\cdot\mathbf p,\phi-m\chi=E\chi \]
第二式解出：
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
代回第一式：
\[ \frac{\mathbf p^2}{E+m}\phi+m\omega^2 r^2\phi=(E-m)\phi \]
两边乘 \(E+m\)：
\[ \mathbf p^2\phi+m\omega^2(E+m)r^2\phi=(E^2-m^2)\phi \]
所以 \(r^2\) 项系数变成 \(m\omega^2(E+m)\)，展开后就是：
\[ m\omega^2 E,r^2+m\omega^2 m,r^2 \]
那个“多出来的 \(m\omega^2 r^2 E\)”正是来自：
\[ (E+m)\times m\omega^2 r^2 \]
中的 \(E\) 部分。
3. 物理意义
可以这样理解：
- 在狄拉克方程中，矢量势 \(V\) 相当于把能量平移 \(E\to E-V\)。
- 当 \(V\) 本身依赖位置时，这个平移会让“有效势”变得能量依赖。
- 而标量势 \(S\) 则相当于把质量平移 \(m\to m+S\)。
- 当 \(S=V\) 时，两者组合出来的等效势恰好是 \(m\omega^2(E+m)r^2\)，其中 \(E\) 依赖是矢量势耦合的典型特征。
所以这个方程不是普通形式
\[ \mathbf p^2+U(r)=E^2-m^2 \]
而是
\[ \mathbf p^2+m\omega^2(E+m)r^2=E^2-m^2 \]
它的“势”是能量依赖的，因此不能直接当作普通本征值问题求解，需要自洽求解 \(E\)。
4. 如果你想要普通形式
如果你希望得到类似
\[ \mathbf p^2+U(r)=E^2-m^2 \]
且势与能量无关，那需要选择纯标量势 \(S(r)=\frac12 m\omega^2 r^2,\ V=0\)。这时方程会变成：
\[ \bigl[\mathbf p^2+(m+S)^2\bigr]\phi=(E^2-m^2)\phi \]
但 \((m+S)^2\) 会展开出 \(m^2+2mS+S^2\)，势项也会变形。等量耦合 \(S=V\) 的好处是 \(\beta S+V\) 对上分量给出 \(2S\)，形式最简洁，但代价就是出现能量依赖。
一句话总结：那个“多出来的 \(m\omega^2 r^2 E\)”来自矢量势 \(V\) 的能量平移效应，是狄拉克方程中矢量耦合的固有特征，不是代数错误。
核心替换
在非相对论里：
\[ \frac{\mathbf p^2}{2m}+U(r)=E_{\rm nr} \]
而在等量标量–矢量耦合 \(S=V=\frac12 m\omega^2 r^2\) 的狄拉克方程中，大分量方程是：
\[ \frac{\mathbf p^2}{E+m}+m\omega^2 r^2\phi=(E-m)\phi \]
所以确实是把
\[ \frac{\mathbf p^2}{2m};\longrightarrow;\frac{\mathbf p^2}{E+m} \]
同时右边的能量也换成了 \(E-m\)。
为什么是 \(E+m\) 而不是 \(2m\)
- 非相对论里动能算符是 \(\mathbf p^2/2m\)，分母 \(2m\) 是常数。
- 狄拉克方程里，小分量消去后产生的分母是 \(E+m\)，它来自
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
这个 \(E+m\) 就是相对论能量分母。
在非相对论极限下：
\[ E\approx m+\varepsilon,\qquad E+m\approx 2m \]
于是
\[ \frac{\mathbf p^2}{E+m}\approx \frac{\mathbf p^2}{2m} \]
正好回到非相对论动能项。
但要注意两点
第一，势项也变了。
非相对论里势能是 \(U(r)\)，而这里出现的是 \(m\omega^2 r^2\)，不是 \(\frac12 m\omega^2 r^2\)。这是因为 \(S=V\) 等量耦合时，上分量感受到的是 \(S+V=2S=m\omega^2 r^2\)。
第二，这不是简单的“把 \(2m\) 换成 \(E+m\)”。
如果把方程两边同乘 \(E+m\)，就得到：
\[ \mathbf p^2+m\omega^2(E+m)r^2=(E^2-m^2) \]
这时 \(r^2\) 项也带上了 \(E\)，所以不能只盯着动能项做替换，势项同样被改了。
一句话总结
在消去小分量后，大分量方程里动能项的分母从常数 \(2m\) 变成了相对论性的 \(E+m\)，这是狄拉克方程消去小分量后的自然结果；非相对论极限下 \(E+m\to 2m\)，就回到普通薛定谔方程的形式。
\[ (\mathbf p^2+m^2)\phi=E^2\phi \]
即
\[ \mathbf p^2\phi=(E^2-m^2)\phi \]
或者写成：
\[ (E^2-\mathbf p^2-m^2)\phi=0 \]
推导
自由狄拉克方程：
\[ (\boldsymbol{\alpha}\cdot\mathbf p+\beta m)\psi=E\psi \]
令 \(\psi=\begin{pmatrix}\phi\ \chi\end{pmatrix}\)，取
\[ \beta=\begin{pmatrix}I&0\0&-I\end{pmatrix},\qquad \boldsymbol{\alpha}=\begin{pmatrix}0&\boldsymbol{\sigma}\\boldsymbol{\sigma}&0\end{pmatrix} \]
得到耦合方程：
\[ \boldsymbol{\sigma}\cdot\mathbf p,\chi+m\phi=E\phi \]
\[ \boldsymbol{\sigma}\cdot\mathbf p,\phi-m\chi=E\chi \]
由第二式：
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi \]
代回第一式：
\[ \frac{\mathbf p^2}{E+m}\phi+m\phi=E\phi \]
两边乘 \(E+m\)：
\[ \mathbf p^2\phi+m(E+m)\phi=E(E+m)\phi \]
即：
\[ \mathbf p^2\phi=(E^2-m^2)\phi \]
与非相对论对比
非相对论自由粒子：
\[ \frac{\mathbf p^2}{2m}\phi=E_{\rm nr}\phi \]
狄拉克自由粒子大分量：
\[ \frac{\mathbf p^2}{E+m}\phi=(E-m)\phi \]
非相对论极限 \(E\approx m+\varepsilon\) 下：
\[ \frac{\mathbf p^2}{E+m}\approx\frac{\mathbf p^2}{2m},\qquad E-m\approx\varepsilon \]
回到薛定谔方程。
注意
- 这里没有势能，所以不会出现 \(E\) 依赖的势项。
- 小分量仍由 \(\chi=\dfrac{\boldsymbol{\sigma}\cdot\mathbf p}{E+m}\phi\) 给出。
- 这个二分量方程本质上是动量空间的质壳关系 \(E^2=\mathbf p^2+m^2\)，自旋自由度尚未体现，因为自由情况下自旋与轨道没有耦合。
\[ \bigl[\boldsymbol{\alpha}\cdot(\mathbf p-q\mathbf A)+\beta m\bigr]\psi=(E-q\varphi)\psi \]
其中 \(q\) 是粒子电荷，\(\varphi\) 是标量势，\(\mathbf A\) 是矢量势。
1. 二分量耦合方程
令 \(\psi=\begin{pmatrix}\phi\ \chi\end{pmatrix}\)，取标准表示：
\[ \beta=\begin{pmatrix}I&0\0&-I\end{pmatrix},\qquad \boldsymbol{\alpha}=\begin{pmatrix}0&\boldsymbol{\sigma}\\boldsymbol{\sigma}&0\end{pmatrix} \]
定义机械动量：
\[ \boldsymbol{\pi}=\mathbf p-q\mathbf A \]
则方程化为：
\[ \boldsymbol{\sigma}\cdot\boldsymbol{\pi},\chi+m\phi=(E-q\varphi)\phi \]
\[ \boldsymbol{\sigma}\cdot\boldsymbol{\pi},\phi-m\chi=(E-q\varphi)\chi \]
2. 消去小分量
由第二式：
\[ \chi=\frac{\boldsymbol{\sigma}\cdot\boldsymbol{\pi}}{E-q\varphi+m}\phi \]
代回第一式：
\[ \boldsymbol{\sigma}\cdot\boldsymbol{\pi}, \frac{\boldsymbol{\sigma}\cdot\boldsymbol{\pi}}{E-q\varphi+m}\phi +m\phi=(E-q\varphi)\phi \]
利用
\[ (\boldsymbol{\sigma}\cdot\boldsymbol{\pi})^2 =\boldsymbol{\pi}^2-i\boldsymbol{\sigma}\cdot(\boldsymbol{\pi}\times\boldsymbol{\pi}) \]
而
\[ \boldsymbol{\pi}\times\boldsymbol{\pi} =(\mathbf p-q\mathbf A)\times(\mathbf p-q\mathbf A) =iq,\mathbf B \]
所以
\[ (\boldsymbol{\sigma}\cdot\boldsymbol{\pi})^2 =\boldsymbol{\pi}^2+q,\boldsymbol{\sigma}\cdot\mathbf B \]
这里出现了自旋–磁场耦合项 \(q\boldsymbol{\sigma}\cdot\mathbf B\)。
3. 最终二分量方程
若忽略分母中 \(E-q\varphi+m\) 与位置的依赖（或把它当作算符处理），得到：
\[ \left[ \boldsymbol{\pi}^2+q\boldsymbol{\sigma}\cdot\mathbf B \right]\phi =(E-q\varphi+m)(E-q\varphi-m)\phi \]
即：
\[ \left[ (\mathbf p-q\mathbf A)^2+q\boldsymbol{\sigma}\cdot\mathbf B \right]\phi =\left[(E-q\varphi)^2-m^2\right]\phi \]
这就是加入电磁场后，大分量满足的二分量方程。
4. 与非相对论对比
非相对论泡利方程：
\[ \left[\frac{(\mathbf p-q\mathbf A)^2}{2m} -\frac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B +q\varphi\right]\phi =E_{\rm nr}\phi \]
注意两点区别：
- 狄拉克方程给出的自旋–磁场耦合系数是 \(q\)，而泡利方程里是 \(-q/2m\)。两者差一个因子 \(2\)，这正是狄拉克方程预言电子自旋 g 因子 \(g=2\) 的来源。
- 狄拉克方程右边是 \((E-q\varphi)^2-m^2\)，非相对论极限下展开才回到泡利方程。
5. 非相对论极限
令
\[ E=m+\varepsilon,\qquad |\varepsilon|\ll m \]
并假设 \(q\varphi\ll m\)，则
\[ (E-q\varphi)^2-m^2 \approx 2m(\varepsilon-q\varphi) \]
同时分母
\[ E-q\varphi+m\approx 2m \]
于是二分量方程化为：
\[ \left[ \frac{(\mathbf p-q\mathbf A)^2}{2m} +\frac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B +q\varphi \right]\phi =\varepsilon\phi \]
这里符号与泡利方程差一个约定，取决于 \(\mathbf B=\nabla\times\mathbf A\) 和电荷符号的约定。标准电子 \(q=-e\) 时，自旋项写作
\[ -\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
与泡利方程一致。
6. 总结
加入电磁场后：
\[ \boxed{ \left[ (\mathbf p-q\mathbf A)^2+q\boldsymbol{\sigma}\cdot\mathbf B \right]\phi =\left[(E-q\varphi)^2-m^2\right]\phi } \]
关键特征：
- 动能项变为机械动量平方 \(\boldsymbol{\pi}^2\)。
- 自动出现自旋–磁场耦合 \(q\boldsymbol{\sigma}\cdot\mathbf B\)，给出 \(g=2\)。
- 右边是 \((E-q\varphi)^2-m^2\)，非相对论极限回到泡利方程。
1. 最小耦合规则
两者都通过最小耦合引入电磁场：
\[ \mathbf p ;\longrightarrow; \mathbf p-q\mathbf A \]
标量势则通过能量平移引入：
\[ E ;\longrightarrow; E-q\varphi \]
区别在于：薛定谔方程里 \(E\) 是普通能量本征值，狄拉克方程里 \(E\) 是相对论总能量。
2. 薛定谔方程加电磁场
自由薛定谔方程：
\[ \frac{\mathbf p^2}{2m}\phi=E_{\rm nr}\phi \]
加电磁场后：
\[ \left[\frac{(\mathbf p-q\mathbf A)^2}{2m}+q\varphi\right]\phi =E_{\rm nr}\phi \]
这是标量势能直接相加，没有自旋项，因为薛定谔方程描述的是无自旋粒子。
3. 泡利方程：薛定谔方程加自旋
如果强行给薛定谔方程加二分量自旋，就得到泡利方程：
\[ \left[ \frac{(\mathbf p-q\mathbf A)^2}{2m} -\frac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B +q\varphi \right]\phi =E_{\rm nr}\phi \]
这里自旋–磁场耦合系数是 \(-\dfrac{q}{2m}\)，是人为加上去的，不是从薛定谔方程本身推出来的。
4. 狄拉克方程加电磁场
狄拉克方程加电磁场后，消去小分量，大分量满足：
\[ \left[ (\mathbf p-q\mathbf A)^2+q\boldsymbol{\sigma}\cdot\mathbf B \right]\phi =\left[(E-q\varphi)^2-m^2\right]\phi \]
注意：
- 左边自动出现 \(q\boldsymbol{\sigma}\cdot\mathbf B\)，系数是 \(q\)，不是 \(-\dfrac{q}{2m}\)；
- 右边是 \((E-q\varphi)^2-m^2\)，不是线性能量；
- 自旋项是推导出来的，不是人为加的。
5. 非相对论极限下的对应
令
\[ E=m+\varepsilon,\qquad |\varepsilon|\ll m,\qquad q\varphi\ll m \]
则
\[ (E-q\varphi)^2-m^2\approx 2m(\varepsilon-q\varphi) \]
狄拉克二分量方程化为：
\[ \left[ \frac{(\mathbf p-q\mathbf A)^2}{2m} +\frac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B +q\varphi \right]\phi =\varepsilon\phi \]
这与泡利方程形式一致，但自旋项系数符号取决于电荷约定。对电子 \(q=-e\)：
\[ \frac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B =-\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
正是泡利方程里的自旋–磁场耦合项。
6. 关键对比表
| 项目 | 薛定谔方程 | 泡利方程 | 狄拉克方程 |
|---|---|---|---|
| 波函数 | 标量 \(\phi\) | 二分量旋量 \(\phi\) | 四分量旋量 \(\psi=(\phi,\chi)\) |
| 动能项 | \(\dfrac{(\mathbf p-q\mathbf A)^2}{2m}\) | \(\dfrac{(\mathbf p-q\mathbf A)^2}{2m}\) | \((\mathbf p-q\mathbf A)^2\) |
| 势能项 | \(q\varphi\) | \(q\varphi\) | \((E-q\varphi)^2-m^2\) |
| 自旋项 | 无 | \(-\dfrac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B\) | \(q\boldsymbol{\sigma}\cdot\mathbf B\) |
| 自旋来源 | 无自旋 | 人为加入 | 方程自动给出 |
| g 因子 | — | 人为取 2 | 自动给出 \(g=2\) |
| 能量关系 | \(E_{\rm nr}=\dfrac{p^2}{2m}+q\varphi\) | 同左，加自旋项 | \(E^2=p^2+m^2\) 的相对论推广 |
7. 核心区别
薛定谔方程：
- 标量波函数，无自旋；
- 电磁场只通过最小耦合和标量势进入；
- 自旋–磁场耦合需要额外人为加入。
泡利方程：
- 二分量波函数，自旋是人为加的；
- 自旋–磁场耦合系数 \(-\dfrac{q}{2m}\) 是唯象引入的。
狄拉克方程：
- 四分量波函数，自旋是方程内禀的；
- 自旋–磁场耦合 \(q\boldsymbol{\sigma}\cdot\mathbf B\) 是自动出现的；
- 非相对论极限自然回到泡利方程，并给出 \(g=2\)；
- 能量关系是相对论性的 \((E-q\varphi)^2-m^2\)，非相对论极限才变成线性能量。
8. 一句话总结
薛定谔方程加电磁场后是
\[ \frac{(\mathbf p-q\mathbf A)^2}{2m}+q\varphi=E_{\rm nr} \]
而狄拉克方程加电磁场后，大分量满足
\[ (\mathbf p-q\mathbf A)^2+q\boldsymbol{\sigma}\cdot\mathbf B =(E-q\varphi)^2-m^2 \]
后者在非相对论极限下自动变成泡利方程，自旋项不是外加的，而是狄拉克方程本身就有的。
1. 薛定谔方程确实含磁场
加电磁场后的薛定谔方程：
\[ \left[\frac{(\mathbf p-q\mathbf A)^2}{2m}+q\varphi\right]\phi =E_{\rm nr}\phi \]
其中 \(\mathbf A\) 就是矢量势，\(\mathbf B=\nabla\times\mathbf A\)。磁场通过最小耦合 \(\mathbf p-q\mathbf A\) 进入方程，影响带电粒子的轨道运动。典型效应包括：
- 朗道能级；
- 回旋运动；
- 阿哈罗诺夫–玻姆效应；
- 磁通量子化。
所以薛定谔方程不是没考虑磁场，而是只考虑了磁场的轨道效应。
2. 薛定谔方程缺的是自旋
薛定谔方程波函数是标量 \(\phi\)，没有自旋自由度，所以：
- 没有 \(\boldsymbol{\sigma}\cdot\mathbf B\) 这一项；
- 不会出现塞曼分裂中的自旋部分；
- 不会给出 \(g=2\)；
- 不会出现自旋进动。
这些必须靠泡利方程或狄拉克方程才能描述。
3. 对比三种方程对磁场的处理
薛定谔方程：
\[ \frac{(\mathbf p-q\mathbf A)^2}{2m}+q\varphi=E_{\rm nr} \]
含磁场，但只有轨道效应，无自旋。
泡利方程：
\[ \left[ \frac{(\mathbf p-q\mathbf A)^2}{2m} -\frac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B +q\varphi \right]\phi =E_{\rm nr}\phi \]
含磁场，且额外加入自旋–磁场耦合，但自旋项是人为加的。
狄拉克方程：
\[ (\mathbf p-q\mathbf A)^2+q\boldsymbol{\sigma}\cdot\mathbf B =(E-q\varphi)^2-m^2 \]
含磁场，自旋–磁场耦合自动出现，非相对论极限回到泡利方程。
4. 更准确的说法
薛定谔方程考虑了磁场的轨道效应，但没有考虑自旋–磁场耦合，因为它描述的是无自旋粒子。
而不是：
薛定谔方程没考虑磁场。
5. 一个直观例子：塞曼效应
原子在磁场中的能级分裂分两部分：
- 轨道塞曼效应：来自电子轨道磁矩与磁场耦合，薛定谔方程就能描述；
- 自旋塞曼效应：来自电子自旋磁矩与磁场耦合，必须用泡利方程或狄拉克方程。
所以薛定谔方程能解释部分塞曼效应，但解释不了反常塞曼效应和自旋分裂。
一句话总结：薛定谔方程含磁场，但只含磁场的轨道效应；它缺的是自旋–磁场耦合，不是整个磁场。
1. 定义
泡利矩阵：
\[ \sigma_x= \begin{pmatrix}0&1\1&0\end{pmatrix},\quad \sigma_y= \begin{pmatrix}0&-i\ i&0\end{pmatrix},\quad \sigma_z= \begin{pmatrix}1&0\0&-1\end{pmatrix} \]
磁场：
\[ \mathbf B=(B_x,B_y,B_z) \]
则
\[ \boldsymbol{\sigma}\cdot\mathbf B =\sigma_xB_x+\sigma_yB_y+\sigma_zB_z \]
2. 矩阵形式
逐项写出：
\[ \sigma_xB_x= \begin{pmatrix}0&B_x\B_x&0\end{pmatrix} \]
\[ \sigma_yB_y= \begin{pmatrix}0&-iB_y\ iB_y&0\end{pmatrix} \]
\[ \sigma_zB_z= \begin{pmatrix}B_z&0\0&-B_z\end{pmatrix} \]
相加：
\[ \boldsymbol{\sigma}\cdot\mathbf B= \begin{pmatrix} B_z & B_x-iB_y\ B_x+iB_y & -B_z \end{pmatrix} \]
3. 作用在旋量上
令二分量旋量
\[ \phi= \begin{pmatrix}\phi_\uparrow\ \phi_\downarrow\end{pmatrix} \]
则
\[ \boldsymbol{\sigma}\cdot\mathbf B,\phi = \begin{pmatrix} B_z & B_x-iB_y\ B_x+iB_y & -B_z \end{pmatrix} \begin{pmatrix}\phi_\uparrow\ \phi_\downarrow\end{pmatrix} \]
即
\[ (\boldsymbol{\sigma}\cdot\mathbf B,\phi)_\uparrow =B_z\phi_\uparrow+(B_x-iB_y)\phi_\downarrow \]
\[ (\boldsymbol{\sigma}\cdot\mathbf B,\phi)_\downarrow =(B_x+iB_y)\phi_\uparrow-B_z\phi_\downarrow \]
4. 沿 z 轴的特殊情况
若 \(\mathbf B=B\hat z\)，则 \(B_x=B_y=0\)，\(B_z=B\)：
\[ \boldsymbol{\sigma}\cdot\mathbf B= \begin{pmatrix}B&0\0&-B\end{pmatrix} \]
作用在旋量上：
\[ \boldsymbol{\sigma}\cdot\mathbf B \begin{pmatrix}\phi_\uparrow\ \phi_\downarrow\end{pmatrix} = \begin{pmatrix}B\phi_\uparrow\ -B\phi_\downarrow\end{pmatrix} \]
自旋向上和向下分别得到 \(+B\) 和 \(-B\)，这就是塞曼分裂的来源。
5. 本征值
\(\boldsymbol{\sigma}\cdot\mathbf B\) 的本征值是
\[ \pm|\mathbf B| \]
本征向量是沿 \(\mathbf B\) 方向的自旋向上和向下态。这来自
\[ (\boldsymbol{\sigma}\cdot\mathbf B)^2=|\mathbf B|^2I \]
6. 与自旋算符的关系
自旋算符
\[ \mathbf S=\frac12\boldsymbol{\sigma} \]
所以
\[ \boldsymbol{\sigma}\cdot\mathbf B=2\mathbf S\cdot\mathbf B \]
自旋–磁场耦合项
\[ q\boldsymbol{\sigma}\cdot\mathbf B=2q\mathbf S\cdot\mathbf B \]
对应磁矩
\[ \boldsymbol{\mu}=q\boldsymbol{\sigma} \]
或写成功率形式
\[ \boldsymbol{\mu}=g\frac{q}{2m}\mathbf S,\qquad g=2 \]
一句话总结：
\[ \boldsymbol{\sigma}\cdot\mathbf B= \begin{pmatrix} B_z & B_x-iB_y\ B_x+iB_y & -B_z \end{pmatrix} \]
本征值 \(\pm|\mathbf B|\)，沿 \(\mathbf B\) 方向的自旋向上/向下态就是它的本征态。
1. 磁场只有 z 分量
\[ \mathbf B=(0,0,B) \]
于是
\[ \boldsymbol{\sigma}\cdot\mathbf B =\sigma_z B =B\begin{pmatrix}1&0\0&-1\end{pmatrix} \]
它只和 \(\sigma_z\) 有关，与 \(\sigma_x,\sigma_y\) 无关。
2. 如果粒子自旋在 x 或 y 方向
设自旋态为 \(\sigma_x\) 的本征态，例如自旋沿 \(+x\)：
\[ |+x\rangle=\frac{1}{\sqrt2} \begin{pmatrix}1\1\end{pmatrix} \]
计算 \(\boldsymbol{\sigma}\cdot\mathbf B\) 的作用：
\[ B\sigma_z|+x\rangle =B\begin{pmatrix}1&0\0&-1\end{pmatrix} \frac{1}{\sqrt2}\begin{pmatrix}1\1\end{pmatrix} =\frac{B}{\sqrt2}\begin{pmatrix}1\-1\end{pmatrix} \]
这不是 \(|+x\rangle\) 的倍数，说明 \(|+x\rangle\) 不是 \(\sigma_z\) 的本征态。
换句话说，沿 \(x\) 方向的自旋态，在 \(z\) 方向磁场中不是定态，会进动。
3. 自旋进动
自旋在磁场中的哈密顿量：
\[ H=-\boldsymbol{\mu}\cdot\mathbf B \]
对电子（\(q=-e\)），磁矩
\[ \boldsymbol{\mu}=-\frac{e}{m}\mathbf S \]
所以
\[ H=\frac{e}{m}\mathbf S\cdot\mathbf B =\frac{eB}{m}S_z \]
自旋在磁场中的演化由
\[ i\frac{d}{dt}|\chi\rangle=H|\chi\rangle \]
给出。若初始自旋沿 \(x\)：
\[ |\chi(0)\rangle=|+x\rangle \]
则
\[ |\chi(t)\rangle =\frac{1}{\sqrt2} \begin{pmatrix} e^{-i\omega_L t/2}\ e^{i\omega_L t/2} \end{pmatrix} \]
其中
\[ \omega_L=\frac{eB}{m} \]
是拉莫尔频率（符号取决于电荷约定）。
这个态的自旋期望值：
\[ \langle S_x\rangle=\frac12\cos(\omega_L t) \]
\[ \langle S_y\rangle=\frac12\sin(\omega_L t) \]
\[ \langle S_z\rangle=0 \]
也就是说，自旋在 \(xy\) 平面内绕 \(z\) 轴进动，进动频率为 \(\omega_L\)。这就是拉莫尔进动。
4. 耦合项展开
自旋–磁场耦合项：
\[ \boldsymbol{\sigma}\cdot\mathbf B =B\sigma_z =B\begin{pmatrix}1&0\0&-1\end{pmatrix} \]
它只区分自旋向上和向下（沿 \(z\)），对横向自旋态：
- 会把 \(|+x\rangle\) 变成 \(|+x\rangle\) 和 \(|-x\rangle\) 的叠加；
- 导致自旋在 \(xy\) 平面内旋转。
从矩阵元看：
\[ \langle +x|\sigma_z|+x\rangle=0 \]
\[ \langle -x|\sigma_z|+x\rangle=1 \]
所以 \(\sigma_z\) 在 \(x\) 方向自旋态之间是非对角的，这正是进动的代数根源。
5. 能级分裂
如果自旋沿 \(z\) 方向，则 \(\sigma_z\) 的本征态是 \(|+z\rangle\) 和 \(|-z\rangle\)，能量分别为：
\[ E_\uparrow=+\frac{eB}{2m},\qquad E_\downarrow=-\frac{eB}{2m} \]
分裂为
\[ \Delta E=\frac{eB}{m}=\omega_L \]
这就是塞曼分裂。
但如果自旋在 \(x\) 或 \(y\) 方向，初始态不是能量本征态，所以不会停在原地，而是进动。
6. 总结
- 磁场沿 \(z\) 时，耦合项是 \(B\sigma_z\)，只和 \(z\) 方向自旋有关。
- 自旋沿 \(x\) 或 \(y\) 的态，不是 \(\sigma_z\) 的本征态，因此在磁场中会进动。
- 进动频率是拉莫尔频率 \(\omega_L=eB/m\)。
- 自旋在 \(xy\) 平面内旋转，\(\langle S_z\rangle\) 保持为零。
- 只有自旋沿 \(z\) 的分量才直接产生能级分裂 \(\Delta E=eB/m\)。
一句话总结：\(z\) 方向磁场只“看见”自旋的 \(z\) 分量；横向自旋态会被 \(B\sigma_z\) 耦合到相反横向态，表现为绕 \(z\) 轴的拉莫尔进动。
1. 自旋磁矩算符
自旋角动量：
\[ \mathbf S=\frac12\boldsymbol{\sigma} \]
自旋磁矩：
\[ \boldsymbol{\mu}=g\frac{q}{2m}\mathbf S \]
其中：
- \(q\)：粒子电荷
- \(m\)：粒子质量
- \(g\)：朗德因子，狄拉克方程给出 \(g=2\)
- \(\mathbf S\)：自旋角动量
代入 \(g=2\)：
\[ \boldsymbol{\mu}=2\frac{q}{2m}\mathbf S =\frac{q}{m}\mathbf S \]
又因为 \(\mathbf S=\frac12\boldsymbol{\sigma}\)：
\[ \boldsymbol{\mu}=\frac{q}{2m}\boldsymbol{\sigma} \]
2. 对电子
电子电荷 \(q=-e\)，所以：
\[ \boldsymbol{\mu} =-\frac{e}{2m}\boldsymbol{\sigma} =-\frac{e}{m}\mathbf S \]
负号表示电子自旋磁矩与自旋角动量方向相反。
3. 与哈密顿量的关系
自旋–磁场耦合哈密顿量：
\[ H=-\boldsymbol{\mu}\cdot\mathbf B \]
对电子：
\[ H=-\left(-\frac{e}{m}\mathbf S\right)\cdot\mathbf B =\frac{e}{m}\mathbf S\cdot\mathbf B \]
又因为 \(\mathbf S=\frac12\boldsymbol{\sigma}\)：
\[ H=\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
如果磁场沿 \(z\)：
\[ H=\frac{eB}{2m}\sigma_z \]
本征值：
\[ E_\uparrow=+\frac{eB}{2m},\qquad E_\downarrow=-\frac{eB}{2m} \]
能级分裂：
\[ \Delta E=\frac{eB}{m} \]
4. 和玻尔磁子的关系
玻尔磁子定义为：
\[ \mu_B=\frac{e\hbar}{2m_e} \]
在自然单位 \(\hbar=1\) 下：
\[ \mu_B=\frac{e}{2m_e} \]
电子自旋磁矩大小：
\[ |\boldsymbol{\mu}|=g\mu_B\sqrt{S(S+1)} \]
对电子 \(S=1/2\)，\(g=2\)：
\[ |\boldsymbol{\mu}|=2\mu_B\frac{\sqrt3}{2} =\sqrt3\mu_B \]
但通常说“电子自旋磁矩”指的是 \(z\) 分量：
\[ \mu_z=-g\mu_B m_s \]
\(m_s=\pm1/2\)，所以：
\[ \mu_z=\mp\mu_B \]
5. 总结
在你前面的讨论中：
\[ \boldsymbol{\mu}=\frac{q}{m}\mathbf S =\frac{q}{2m}\boldsymbol{\sigma} \]
对电子 \(q=-e\)：
\[ \boldsymbol{\mu}=-\frac{e}{m}\mathbf S =-\frac{e}{2m}\boldsymbol{\sigma} \]
所以自旋–磁场耦合：
\[ H=-\boldsymbol{\mu}\cdot\mathbf B =\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
这就是前面拉莫尔进动和塞曼分裂里出现的那个磁矩。
1. 经典回旋半径（轨道运动）
带电粒子在磁场中做圆周运动，洛伦兹力提供向心力：
\[ qvB=\frac{mv^2}{r} \]
解出回旋半径：
\[ r=\frac{mv}{qB} \]
或者用动量 \(p=mv\)：
\[ r=\frac{p}{qB} \]
对电子 \(q=-e\)，取大小：
\[ r=\frac{mv}{eB}=\frac{p}{eB} \]
这就是回旋半径，也叫拉莫尔半径。
2. 量子力学中的磁长度
在量子力学里，磁场中的带电粒子形成朗道能级。特征长度叫磁长度：
\[ l_B=\sqrt{\frac{\hbar}{qB}} \]
自然单位 \(\hbar=1\)：
\[ l_B=\sqrt{\frac{1}{qB}} \]
对电子 \(q=e\)：
\[ l_B=\sqrt{\frac{1}{eB}} \]
它对应朗道能级中波函数的高斯衰减宽度，可以看作“量子回旋半径”。
3. 和自旋进动的关系
注意区分两个东西：
- 轨道回旋半径 \(r=mv/qB\)：带电粒子在磁场中做圆周运动的半径。
- 自旋进动：自旋绕磁场方向进动，进动频率 \(\omega_L=qB/m\)，但它不是空间中的圆周运动，而是自旋态在希尔伯特空间中的旋转。
两者频率相同：
\[ \omega_c=\frac{qB}{m} \]
轨道回旋频率和自旋拉莫尔频率在经典图像里一致，这是磁矩与角动量比值 \(q/m\) 的结果。
4. 如果用能量表示
非相对论下：
\[ E=\frac12 mv^2 \]
所以
\[ v=\sqrt{\frac{2E}{m}} \]
回旋半径：
\[ r=\frac{m}{qB}\sqrt{\frac{2E}{m}} =\frac{\sqrt{2mE}}{qB} \]
相对论下用动量 \(p\)：
\[ r=\frac{p}{qB} \]
其中 \(p=\sqrt{E^2-m^2}\)（自然单位 \(c=1\)）。
5. 如果你问的是别的 \(r\)
如果你说的“半径 \(r\)”是指：
- 谐振子势里的 \(r\)：那是 \(r=\sqrt{x^2+y^2+z^2}\)，出现在 \(S=V=\frac12 m\omega^2 r^2\) 里。
- 球坐标里的径向坐标：\(r\in[0,\infty)\)。
- 波函数里的径向节点：那和主量子数有关。
告诉我具体语境，我可以给出对应的表达式。
一个带电粒子在半径 \(r\) 的圆周上运动，它的磁矩 \(\boldsymbol{\mu}\) 是多少？
如果是这个意思，那就是经典的轨道磁矩问题。下面给出结果和推导。
1. 轨道磁矩公式
带电粒子做圆周运动，相当于一个电流环。
电流：
\[ I=\frac{q}{T} \]
周期：
\[ T=\frac{2\pi r}{v} \]
所以
\[ I=\frac{qv}{2\pi r} \]
圆环面积：
\[ A=\pi r^2 \]
磁矩大小：
\[ \mu=IA =\frac{qv}{2\pi r}\cdot \pi r^2 =\frac12 qvr \]
角动量：
\[ L=mvr \]
所以
\[ \mu=\frac{q}{2m}L \]
写成矢量：
\[ \boldsymbol{\mu}=\frac{q}{2m}\mathbf L \]
2. 对电子
电子电荷 \(q=-e\)：
\[ \boldsymbol{\mu} =-\frac{e}{2m}\mathbf L \]
负号表示电子轨道磁矩与轨道角动量方向相反。
3. 量子化形式
量子力学中：
\[ L_z=m_l\hbar \]
所以轨道磁矩 \(z\) 分量：
\[ \mu_z=-\frac{e}{2m}m_l\hbar =-m_l\mu_B \]
其中玻尔磁子：
\[ \mu_B=\frac{e\hbar}{2m} \]
4. 自旋磁矩对比
如果是自旋而不是轨道运动：
\[ \boldsymbol{\mu}_s =g\frac{q}{2m}\mathbf S \]
狄拉克给出 \(g=2\)，所以对电子：
\[ \boldsymbol{\mu}_s =-\frac{e}{m}\mathbf S =-\frac{e}{2m}\boldsymbol{\sigma} \]
注意轨道磁矩是 \(\dfrac{q}{2m}\mathbf L\)，自旋磁矩多了一个 \(g\) 因子。
5. 如果半径 \(r\) 出现在谐振子势里
如果你说的 \(r\) 是谐振子势
\[ S=V=\frac12 m\omega^2 r^2 \]
里的径向坐标，那“求 \(\mu\)”通常指磁矩算符，它与 \(r\) 没有直接简单关系，需要知道具体的电流分布或角动量态。
一句话总结：
如果粒子在半径 \(r\) 的圆周上以速度 \(v\) 运动，
\[ \mu=\frac12 qvr=\frac{q}{2m}L \]
对电子：
\[ \boldsymbol{\mu}=-\frac{e}{2m}\mathbf L \]
如果你指的是自旋磁矩，则
\[ \boldsymbol{\mu}_s=\frac{q}{m}\mathbf S \]
对电子：
\[ \boldsymbol{\mu}_s=-\frac{e}{m}\mathbf S \]
告诉我具体是轨道还是自旋，我可以进一步细化。
1. 自旋磁矩
对电子，自旋磁矩：
\[ \boldsymbol{\mu}_s=-\frac{e}{m}\mathbf S =-\frac{e}{2m}\boldsymbol{\sigma} \]
2. 点乘 B
\[ \boldsymbol{\mu}_s\cdot\mathbf B =-\frac{e}{m}\mathbf S\cdot\mathbf B =-\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
3. 哈密顿量
自旋–磁场耦合哈密顿量：
\[ H=-\boldsymbol{\mu}_s\cdot\mathbf B \]
注意这里有个负号，因为磁矩在磁场中的势能是 \(-\boldsymbol{\mu}\cdot\mathbf B\)。
代入：
\[ H=-\left(-\frac{e}{m}\mathbf S\cdot\mathbf B\right) =\frac{e}{m}\mathbf S\cdot\mathbf B \]
用泡利矩阵：
\[ H=\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
4. 展开
若
\[ \mathbf B=(B_x,B_y,B_z) \]
则
\[ \boldsymbol{\sigma}\cdot\mathbf B= \begin{pmatrix} B_z & B_x-iB_y\ B_x+iB_y & -B_z \end{pmatrix} \]
所以
\[ H=\frac{e}{2m} \begin{pmatrix} B_z & B_x-iB_y\ B_x+iB_y & -B_z \end{pmatrix} \]
5. 若 B 沿 z
\[ \mathbf B=B\hat z \]
则
\[ H=\frac{eB}{2m}\sigma_z =\frac{eB}{2m} \begin{pmatrix} 1&0\ 0&-1 \end{pmatrix} \]
本征值：
\[ E_\uparrow=+\frac{eB}{2m},\qquad E_\downarrow=-\frac{eB}{2m} \]
能级分裂：
\[ \Delta E=\frac{eB}{m} \]
6. 和非相对论对比
泡利方程里自旋项：
\[ H_{\rm Pauli} =-\frac{q}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
对电子 \(q=-e\)：
\[ H_{\rm Pauli} =\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
与狄拉克非相对论极限一致。
一句话总结：
\[ \boldsymbol{\mu}_s\cdot\mathbf B =-\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
耦合哈密顿量：
\[ H=-\boldsymbol{\mu}_s\cdot\mathbf B =\frac{e}{2m}\boldsymbol{\sigma}\cdot\mathbf B \]
若 \(B\) 沿 \(z\)，能级分裂为 \(\Delta E=eB/m\)。
1. 轨道磁矩
电子绕半径 \(r\) 的圆轨道运动，相当于一个环形电流。
电流大小：
\[ I=\frac{e}{T} \]
其中 \(T\) 是绕一圈的时间：
\[ T=\frac{2\pi r}{v} \]
所以
\[ I=\frac{ev}{2\pi r} \]
圆环面积：
\[ A=\pi r^2 \]
轨道磁矩大小：
\[ \mu=IA =\frac{ev}{2\pi r}\cdot \pi r^2 =\frac12 evr \]
用轨道角动量 \(L=mvr\)：
\[ \mu=\frac{e}{2m}L \]
对电子，磁矩方向和角动量方向相反，所以写成：
\[ \boldsymbol{\mu}_L=-\frac{e}{2m}\mathbf L \]
2. 轨道在 xy 平面，L 沿 z
如果电子在 \(xy\) 平面内绕圈，角动量 \(\mathbf L\) 沿 \(z\) 方向：
\[ \mathbf L=L_z\hat z \]
所以
\[ \boldsymbol{\mu}_L=-\frac{e}{2m}L_z\hat z \]
即
\[ \mu_z=-\frac{e}{2m}L_z \]
3. 和 z 方向磁场点乘
磁场：
\[ \mathbf B=B\hat z \]
点乘：
\[ \boldsymbol{\mu}_L\cdot\mathbf B =\mu_z B =-\frac{eB}{2m}L_z \]
这就是轨道磁矩和磁场的耦合量。
4. 耦合能量
磁场中磁矩的势能：
\[ H=-\boldsymbol{\mu}_L\cdot\mathbf B \]
代入：
\[ H =-\left(-\frac{eB}{2m}L_z\right) =\frac{eB}{2m}L_z \]
这就是轨道塞曼效应的哈密顿量。
5. 量子化
量子力学里：
\[ L_z=m_l\hbar \]
所以
\[ H=\frac{eB}{2m}m_l\hbar =m_l\mu_B B \]
其中玻尔磁子：
\[ \mu_B=\frac{e\hbar}{2m} \]
能级分裂：
\[ \Delta E=\mu_B B \]
相邻 \(m_l\) 之间差一个 \(\mu_B B\)。
6. 总结
- 轨道磁矩：\(\boldsymbol{\mu}_L=-\dfrac{e}{2m}\mathbf L\)
- 若轨道在 \(xy\) 平面：\(\mu_z=-\dfrac{e}{2m}L_z\)
- 和 \(z\) 方向磁场点乘：\(\boldsymbol{\mu}_L\cdot\mathbf B=-\dfrac{eB}{2m}L_z\)
- 耦合能量：\(H=\dfrac{eB}{2m}L_z=m_l\mu_B B\)
- 能级分裂间隔：\(\Delta E=\mu_B B\)
这就是纯轨道、不带自旋的版本。
- 轨道磁矩和磁场的耦合：\(\boldsymbol{\mu}_L\cdot\mathbf B\)
- 自旋和磁场的耦合：\(\boldsymbol{\sigma}\cdot\mathbf B\)
不混在一起，直接并排看。
1. 轨道版本
电子在 \(xy\) 平面绕半径 \(r\) 的圆轨道运动，磁场沿 \(z\)：
\[ \mathbf B=B\hat z \]
轨道磁矩：
\[ \boldsymbol{\mu}_L=-\frac{e}{2m}\mathbf L \]
和磁场点乘：
\[ \boldsymbol{\mu}_L\cdot\mathbf B =-\frac{eB}{2m}L_z \]
耦合能量：
\[ H_L=-\boldsymbol{\mu}_L\cdot\mathbf B =\frac{eB}{2m}L_z \]
量子化后：
\[ L_z=m_l\hbar \]
所以
\[ H_L=m_l\mu_B B,\qquad \mu_B=\frac{e\hbar}{2m} \]
能级分裂间隔：
\[ \Delta E_L=\mu_B B \]
2. 自旋版本
自旋磁矩：
\[ \boldsymbol{\mu}_s=-\frac{e}{m}\mathbf S =-\frac{e}{2m}\boldsymbol{\sigma} \]
和磁场点乘：
\[ \boldsymbol{\mu}_s\cdot\mathbf B =-\frac{eB}{2m}\sigma_z \]
耦合能量：
\[ H_s=-\boldsymbol{\mu}_s\cdot\mathbf B =\frac{eB}{2m}\sigma_z \]
本征值：
\[ \sigma_z=\pm1 \]
所以
\[ E_\uparrow=+\frac{eB}{2m},\qquad E_\downarrow=-\frac{eB}{2m} \]
能级分裂间隔：
\[ \Delta E_s=\frac{eB}{m}=2\mu_B B \]
3. 并排对比
| 项目 | 轨道 | 自旋 |
|---|---|---|
| 磁矩 | \(\boldsymbol{\mu}_L=-\dfrac{e}{2m}\mathbf L\) | \(\boldsymbol{\mu}_s=-\dfrac{e}{2m}\boldsymbol{\sigma}\) |
| 与 B 点乘 | \(-\dfrac{eB}{2m}L_z\) | \(-\dfrac{eB}{2m}\sigma_z\) |
| 耦合能量 | \(\dfrac{eB}{2m}L_z\) | \(\dfrac{eB}{2m}\sigma_z\) |
| 量子数 | \(L_z=m_l\hbar\) | \(\sigma_z=\pm1\) |
| 能级间隔 | \(\mu_B B\) | \(2\mu_B B\) |
4. 关键区别
第一，量子数不同。
- 轨道：\(L_z=m_l\hbar\)，\(m_l=0,\pm1,\pm2,\dots\)
- 自旋：\(\sigma_z=\pm1\)，只有两个值
第二，能级间隔差一倍。
轨道分裂间隔是 \(\mu_B B\)，自旋分裂间隔是 \(2\mu_B B\)。
这个“2”就是自旋 \(g\) 因子 \(g=2\) 的来源。
第三，公式形式几乎一样。
两者都是
\[ -\frac{eB}{2m}\times(\text{角动量量子数}) \]
区别只在于角动量量子数是 \(L_z\) 还是 \(\sigma_z\)，以及 \(\sigma_z\) 自带因子 2。
5. 如果写成统一形式
定义总角动量 \(z\) 分量：
\[ J_z=L_z+S_z \]
其中
\[ S_z=\frac12\hbar\sigma_z \]
则总塞曼耦合：
\[ H=\frac{eB}{2m}(L_z+2S_z) \]
注意自旋项多了一个因子 2，这就是 \(g=2\) 的体现。
6. 一句话总结
轨道版本：
\[ \boldsymbol{\mu}_L\cdot\mathbf B=-\frac{eB}{2m}L_z \]
自旋版本：
\[ \boldsymbol{\mu}_s\cdot\mathbf B=-\frac{eB}{2m}\sigma_z \]
形式相同，但自旋的能级分裂是轨道的两倍，因为自旋 \(g=2\)。
1. 你说的关系
\[ r=\frac{\lambda}{2\pi} \]
意思是：电子的轨道周长 \(2\pi r\) 正好等于一个德布罗意波长 \(\lambda\)。
也就是：
\[ 2\pi r=\lambda \]
这是驻波条件里最简单的一种：绕一圈刚好放一个波。
2. 玻尔角动量量子化
德布罗意关系：
\[ \lambda=\frac{h}{p}=\frac{h}{mv} \]
如果
\[ 2\pi r=n\lambda,\qquad n=1,2,3,\dots \]
那么
\[ 2\pi r=n\frac{h}{mv} \]
整理：
\[ mvr=n\frac{h}{2\pi}=n\hbar \]
这就是玻尔的角动量量子化条件：
\[ L=mvr=n\hbar \]
你写的 \(r=\lambda/2\pi\)，对应的是 \(n=1\) 的情况：
\[ 2\pi r=\lambda \]
所以
\[ mvr=\hbar \]
3. 物理意义
这个关系说的是：
- 电子绕核运动时，如果轨道周长正好是整数个德布罗意波长，就能形成稳定驻波。
- 这种驻波条件给出角动量量子化 \(L=n\hbar\)。
- \(n=1\) 时，周长等于一个波长，角动量就是 \(\hbar\)。
4. 但要注意
这个图像来自玻尔模型，是半经典的：
- 它假设电子有确定轨道半径 \(r\) 和确定速度 \(v\)。
- 量子力学里，电子没有经典轨道，只有波函数和概率分布。
- 不过德布罗意驻波条件确实能直观解释角动量量子化，所以常被用作引入。
5. 如果代入具体数值
若
\[ r=\frac{\lambda}{2\pi} \]
则
\[ \lambda=2\pi r \]
又
\[ \lambda=\frac{h}{mv} \]
所以
\[ mv=\frac{h}{2\pi r} \]
于是
\[ mvr=\frac{h}{2\pi}=\hbar \]
这就是 \(n=1\) 的角动量。
6. 总结
你说的
\[ r=\frac{\lambda}{2\pi} \]
等价于
\[ 2\pi r=\lambda \]
即轨道周长等于一个德布罗意波长。由此得到
\[ mvr=\hbar \]
这是玻尔角动量量子化 \(L=n\hbar\) 中 \(n=1\) 的特例。更一般地：
\[ 2\pi r=n\lambda \Rightarrow mvr=n\hbar \]
如果你是想问“电子的 \(mvr\) 是不是等于 \(\hbar\)”，那要看 \(n\) 是几；\(n=1\) 时才是 \(\hbar\)。