量子简谐振子(级数法)

                     

贡献者: JierPeter; addis

  • 本文处于草稿阶段。
预备知识 简谐振子(升降算符)

   量子力学中,简谐振子的定态薛定谔方程为

\begin{equation} -\frac{\hbar^2}{2m} \frac{\mathrm{d}^{2}{\psi}}{\mathrm{d}{x}^{2}} + \frac12 m \omega^2 x^2\psi = E\psi~, \end{equation}
能级为
\begin{equation} E_n = \left(\frac12 + n \right) \hbar \omega ~. \end{equation}
归一化的束缚态波函数为
\begin{equation} \psi_n (x) = \frac{1}{\sqrt{2^n n!}} \left(\frac{\alpha^2}{\pi } \right) ^{1/4} H_n(u) \mathrm{e} ^{-u^2/2}~. \end{equation}
其中
\begin{equation} \alpha \equiv \sqrt{\frac{m\omega}{\hbar }}, \qquad u \equiv \alpha x~, \end{equation}
$H_n(u)$ 是 Hermite 多项式

图
图 1:简谐振子的前几个束缚态波函数(式 3 ,$\alpha = 1$)

   前 4 个波函数分别为(注意函数的奇偶性与角标的奇偶性相同)

\begin{equation} \begin{aligned} \psi_0(x) &= \left(\frac{\alpha^2}{\pi} \right) ^{1/4} \mathrm{e} ^{-u^2/2}~, & \psi_2(x) &= \left(\frac{\alpha^2}{\pi } \right) ^{1/4} \frac{1}{\sqrt 2 } (2u^2 - 1) \mathrm{e} ^{-u^2/2}~,\\ \psi_1(x) &= \left(\frac{\alpha ^2}{\pi } \right) ^{1/4} \sqrt2u \mathrm{e} ^{-u^2/2}~, \quad & \psi_3(x) &= \left(\frac{\alpha ^2}{\pi } \right) ^{1/4} \frac{1}{\sqrt 3} u(2u^2 - 3) \mathrm{e} ^{-u^2/2}~. \end{aligned} \end{equation}

   Matlab 代码如(使用原子单位制

代码 1:psi_SHO.m
function psi = psi_SHO(m, w, n, x)
alpha = sqrt(m*w); u = alpha * x;
psi = 1/sqrt(2^n*factorial(n)) *...
    (alpha^2/pi)^0.25 * hermiteH(n, u) .* exp(-u.^2/2);
end

1. 推导(生成函数法)

   为了方便,记 $\varepsilon \equiv 2mE/\hbar^2$,再使用 $\hbar=1$ 的单位制,于是式 1 化为

\begin{equation} \frac{\mathrm{d}^2}{ \,\mathrm{d}{y} ^2}\psi(x) + (\varepsilon-m^2\omega^2x^2)\psi(x)=0~, \end{equation}
如果记 $\psi(x)=h(x) \exp\left(-m\omega x^2/2\right) $1,那么代入式 6 后得
\begin{equation} h''(x)-2m\omega xh'(x)+(\varepsilon-m\omega)h=0~. \end{equation}
只要能解出式 7 即可。

生成函数

   考虑生成函数 $g(x, t)\equiv \exp\left(-t^2+2tx\right) $,用它来定义一系列 Hermite 多项式 $H_n$:

\begin{equation} g(x, t) = \exp\left(-t^2+2tx\right) = \sum_{n=0}^\infty H_n(x)\frac{t^n}{n!}~, \end{equation}
\begin{equation} \frac{\partial^n}{\partial t^n}g(x, t)\mid_{t=0} = H_n(x)~, \end{equation}
容易得
\begin{equation} H_0(x) = g(x, 0) = 1~. \end{equation}

   对式 8 求关于 $x$ 的偏导数得:

\begin{equation} \begin{aligned} \sum_{n=0}^\infty H_n(x)\frac{2t^{n+1}}{n!}&=2t \exp\left(-t^2+2tx\right) \\ &=\frac{\partial}{\partial x} \exp\left(-t^2+2tx\right) \\ &=\frac{\partial}{\partial x}\sum_{n=0}^\infty H_n(x)\frac{t^n}{n!}\\ &=\sum_{n=0}^\infty H'_n(x)\frac{t^n}{n!}~. \end{aligned} \end{equation}
比较式 11 两端后得
\begin{equation} H'_n(x) = 2nH_{n-1}(x)~. \end{equation}

   由式 12 式 10 ,再考虑到 $H_n(x)$ 的常数项是 $H_n(0)=\frac{\partial^n}{\partial t^n}g(0, t)\mid_{t=0}=\frac{\partial^n}{\partial t^n} \exp\left(-t^2\right) \mid_{t=0}$,可以推知所有 Hermite 多项式:

\begin{equation} H_0(x)=1\implies H'_1(x)=2\implies H_1(x)=2x~, \end{equation}
\begin{equation} H_1(x)=2x\implies H'_2(x)=8x\implies H_2(x)=4x^2-2~, \end{equation}
\begin{equation} H_2(x)=4x^2-2\implies H'_3(x)=24x^2-12\implies H_3(x)=8x^3-12x~. \end{equation}
以此类推。

证明 $H_n$ 满足 $h$ 的方程

   对式 8 求关于 $t$ 的偏导数得:

\begin{equation} \begin{aligned} (2x-2t)\sum_{n=0}^\infty H_n(x)\frac{t^n}{n!}&=(2x-2t) \exp\left(-t^2+2tx\right) \\ &=\frac{\partial}{\partial t} \exp\left(-t^2+2tx\right) \\ &=\frac{\partial}{\partial t}\sum_{n=0}^\infty H_n(x)\frac{t^n}{n!}\\ &=\sum_{n=0}^\infty H_n(x)\frac{t^{n-1}}{(n-1)!}~. \end{aligned} \end{equation}
比较式 16 两端后得
\begin{equation} \frac{2xH_n(x)}{n!}-\frac{2H_{n-1}(x)}{(n-1)!}=\frac{H_{n+1}}{n!}(x~, \end{equation}
整理得
\begin{equation} H_{n+1}(x) = 2xH_n(x)-2nH_{n-1}(x)~, \end{equation}
等价于
\begin{equation} H_{n}(x) = 2xH_{n-1}(x)-2(n-1)H_{n-2}(x)~. \end{equation}

   将式 19 式 12 组合后得

\begin{equation} \begin{aligned} H_n''(x)&=2nH'_{n-1}(x)\\ &=2n\cdot 2(n-1)H_{n-2}(x)\\ &=2n\cdot(2xH_{n-1}(x)-H_n(x))\qquad\text{这一步代入式 19 }\\ &=2xH'_n(x)-2nH_n(x)~, \end{aligned} \end{equation}
整理一下式 20
\begin{equation} H_n''-2xH_n'-2nH_n=0~. \end{equation}

   式 21 式 7 形式完全相同,系数的差异可以通过合适的变量代换或单位选择来消除。将 $\hbar$ 添回式 6 ,进行合适的变量代换(式 4 ),即可从式 7 得到式 3


1. ^ 指数上的符号是为了让波函数在无穷远处趋于 0。


致读者: 小时百科一直以来坚持所有内容免费,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 10 元,我们一个星期内就能脱离亏损, 并保证在接下来的一整年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利