SO(3) Symmetric potential
粒子在具有 SO(3) 球对称性的势中运动的定态 Schrodinger 方程为
\begin{align} \left( -\frac{\hbar^2}{2\mu} \nabla^2 + V(r) \right)\psi(\vec{r}) = E \psi(\vec{r}) \end{align}其中 \(\nabla^2\) 为
\begin{align} \nabla^2 = \frac{\partial^2}{\partial r^2} + \frac{2}{r} \frac{\partial}{\partial r} - \frac{1}{r^2} \frac{\hat{L}^2}{\hbar^2} \end{align}SO(3) 群有 \(2l + 1\) 维的不可约表示. 由不可约表示的维度, 可得出它的能级至少有 \(2l + 1\) 重的简并. \(Y_{l,m}(\theta, \phi)\) 给出了 \(SO(3)\) 群的所有不可约表示. 波函数可 以写成
\begin{align} \psi(\vec{r}) = f_l(r) Y_{l,m} (\theta, \phi) \end{align}其中 \(Y_{l,m}(\theta, \phi)\) 满足
\begin{align} \hat{\Pi} Y_{l,m}(\theta, \phi) = ( - 1)^l Y_{l,m}(\theta, \phi) \end{align}代回原来的方程有
\begin{align} \left[ -\frac{\hbar^2}{2\mu}\left( \frac{\mathrm{d}^2}{\mathrm{d}r^2} + \frac{2}{r} \frac{\mathrm{d}}{\mathrm{d}r} \right) + \frac{l(l + 1)}{r^2} \frac{\hbar^2}{2\mu} + V(r) \right] f_l (r) = E f_l (r) \end{align}如果做代换
\begin{align} \psi(\vec{r}) = \frac{\phi_l(r)}{r} Y_{l,m}(\theta, \phi) \end{align}那么径向的方程就等价成了一个一维问题
\begin{align} \left( -\frac{\hbar^2}{2\mu}\frac{\mathrm{d}^2}{\mathrm{d}r^2} + \frac{l(l + 1)}{r^2} \frac{\hbar^2}{2\mu} + V(r) \right) \phi_l(r) = E \phi_l(r) \end{align}不用说, 径向方程是正交归一的, 也可以验证
\begin{align} \langle \psi_{l,m} | \psi_{l', m'} \rangle = \delta_{l,l'} \delta_{m, m'} \int_0^{\infty} r^2 \mathrm{d}r\cdot \frac{\phi_l^{ * }(r)}{r} \frac{\phi_{l'}(r)}{r} = \delta_{l,l'} \delta_{m, m'} \int_0^{\infty} \mathrm{d}r\cdot \phi_l^{ * }(r) \phi_{l'}(r) \end{align}Boundary Condition
假设 \(\phi\) 在 \(r \to 0\) 时的行为为 \(r^{\alpha}\) . 如果 \(V(r)\) 趋于 \(0\) 的速度比 \(r^{-2}\) 要慢, 那么将 \(\phi_l (r) \propto r^{\alpha}\) 代回径向方程就可以得到
\begin{align} -\frac{\mathrm{d}^2}{\mathrm{d}r^2} + l(l + 1) r^{\alpha - 2} = k^2 r^{\alpha} \\ \left[ - \alpha(\alpha - 1) + l(l + 1)\right] r^{\alpha - 2} = k^2r^a \end{align}等式右边相对于左边在 \(r\to 0\) 时可以忽略, 可得 \(\alpha = l + 1\) 或者 \(\alpha = - l\) . 但 $α = -l $ 时, \(\psi\propto r^{-l-1}\) 会使得波函数发散, 无法归一化, 因此
\begin{align} \phi_l(0) = 0 \quad\mathrm{for}\quad \mathrm{all} \quad l\\ \phi_l(r) \propto r^{l+1} \quad \mathrm{for} \quad r \to 0 \end{align}Radially Symmetric Harmonic Oscillator
此时的势能为
\begin{align} V(r) = \frac{\mu}{2} \omega^2 r^2 \end{align}由于 \(r \to \infty\) 时, \(V(r) \to \infty\) , 所以解只存在 bound states
\begin{align} E_{n, l} = \left( 2n + l +\frac{3}{2} \right)\hbar\omega \end{align}波函数是 generalized Laguerre \(L_n^{\alpha}(x)\) 多项式的形式
\begin{align} \phi_{n, l} = 2(\sqrt{\pi}\beta)^{-1/2} \left[ \frac{2^{n + l} n!}{(2n + 2l + 1)!!} \right] \left( \frac{r}{\beta} \right)^{l+1} L_n^{l + 1/2} \left(\frac{r^2}{\beta^2}\right) e^{-\frac{r^2}{2\beta^2}} \end{align}其中
\begin{align} \beta = \sqrt{\frac{\hbar}{\mu\omega}} \end{align}Short ranged
短程力, 也就是说当 \(r > r_0\) 时, \(V(r) = 0\)
Bound states
当 \(E < 0\) 时, 在 \(r > r_0\) 处的解为 modified Bessel functions, 舍去发散的解, 只 剩下
\begin{align} \phi_-(r) = \sqrt{\kappa r} K_{l+1/2}(\kappa r) \end{align}如果在 \(r > r_0\) 处为 Coulomb 势, 也就是
\begin{align} V(r) = - \frac{C}{r}, \quad r > r_0 \end{align}那么它在 \(r > r_0\) 处的解是 Whittaker functions
\begin{align} \phi_-(r) = W_{\gamma, l+1/2}(2 \kappa r) \end{align}其中 \(\gamma = \frac{\mu C}{\hbar^2\kappa}\)
Unbound states
当 \(E > 0\) 时, 在 \(r > r_0\) 处的解有两个, 是 Spherical Bessel 和 Neumann functions
\begin{align} \phi_s(r) = k r j_l(kr) \\ \phi_c(r) = kr n_l(kr) \end{align}它们的渐近行为为
\begin{align} \phi_s(r) = \sin \left(kr - \frac{l \pi}{2} \right) \left[ 1 + \mathcal{O}\left( \frac{1}{r} \right) \right] \\ \phi_c(r) = \cos \left(kr - \frac{l \pi}{2} \right) \left[ 1 + \mathcal{O}\left( \frac{1}{r} \right) \right] \end{align}可以写成相移的形式
\begin{align} \phi(r) \propto \sin \left( kr - \frac{l\pi}{2} + \delta_l \right) \end{align}相称表征了波函数相对于 regular 解的 shift. 如果在 \(r > r_0\) 处为 Coulomb 势, 也就是
\begin{align} V(r) = - \frac{C}{r}, \quad r > r_0 \end{align}那么对应的结果为 regular Coulomb function 和 irregular Coulomb function
\begin{align} \phi_s(r) = F_l(\eta, kr) \\ \phi_c(r) = G_l(\eta, kr), \quad r > r_0 \end{align}其中
\begin{align} \eta = - \frac{\mu C}{\hbar^2 k} \end{align}渐近行为为
\begin{align} F_l(\eta, kr) \to \sin\left( kr - \eta \ln 2kr - \frac{l\pi}{2} + \sigma_l \right),\quad \mathrm{for} \quad r\to\infty \\ G_l(\eta, kr) \to \cos\left( kr - \eta \ln 2kr - \frac{l\pi}{2} + \sigma_l \right),\quad \mathrm{for} \quad r\to\infty \end{align}类似地, 也可以写成相移的形式.
Reference
- Harald Siegfried Friedrich, Theoretical Atomic Physics-Springer (2005) Chap 1.2, 1.3
后记
昨天返校, 今天终于把从上学期就想整理的东西草草整理了.