原子单位制

                     

贡献者: addis

预备知识 1 国际单位制,无单位的物理公式,玻尔原子模型

  1在量子力学的许多理论或数值计算中,选用原子单位(Hartree atomic unit)会更方便。需要特别注意的是,原子单位下的物理公式在使用时的习惯和其他单位制有所不同,例如会出现 $E = \omega$(能量等于角频率,见式 6 )这种看似不符合量纲分析的公式。要理解这种记号,这里推荐的方法是把使用原子单位的公式都理解成无单位的物理公式

   若声明使用原子单位,物理量的数值后面就不需要标记单位,例如 “氢原子基态能量为 $E \approx -0.5$”。有时候为了强调我们使用原子单位,我们会在数值后面加上 “a.u.”,如 “$E \approx 0.5 \,\mathrm{a.u.} $”。这里的 “a.u.” 可等效为 “单位 1”,类似弧度单位 “Rad”2

1. 推导

   原子单位下的物理量和公式都可以看作是无量纲的数字。沿用 “无单位的物理公式” 中的记号,长度、质量、时间、能量、一维波函数的转换常数分别记为 $\beta_x$,$\beta_m$,$\beta_t$,$\beta_E$,$\beta_\Psi$。

   在量子力学中,由于普朗克常数 $\hbar$ 大量出现,我们希望它在原子单位下等于 $1$,即 $\beta_L = \hbar$。注意 $\hbar$ 具有角动量量纲。令角动量公式 $ \boldsymbol{\mathbf{L}} = m \boldsymbol{\mathbf{r}} \boldsymbol\times \boldsymbol{\mathbf{v}} $ 在原子单位下依然成立,也就是

\begin{equation} \beta_L \boldsymbol{\mathbf{L}} _a = (\beta_m m_a) (\beta_x \boldsymbol{\mathbf{r}} _a) \boldsymbol\times \left(\frac{\beta_x}{\beta_t} \boldsymbol{\mathbf{v}} _a \right) ~. \end{equation}
可以化为 $ \boldsymbol{\mathbf{L}} _a = m_a \boldsymbol{\mathbf{r}} _a \boldsymbol\times \boldsymbol{\mathbf{v}} _a$,即
\begin{equation} \beta_t = \frac{\beta_m \beta_x^2}{\hbar}~. \end{equation}
这样,量子力学中所有公式中的 $\hbar$ 都可以看作是一个角动量常数并用数字 $1$ 代替了。

   再看能量,令动能公式 $E_k = mv^2/2$ 在原子单位中成立,那么同理要求

\begin{equation} \beta_E = \frac{\beta_m \beta_x^2}{\beta_t^2} = \frac{\hbar}{\beta_t}~. \end{equation}

例 1 光子的能量

   我们来看光子的能量公式

\begin{equation} E = \hbar \omega \Longrightarrow E_a = \frac{\hbar\beta_\omega}{\beta_E}\omega_a~. \end{equation}
为满足周期的定义($T_a = 2\pi/\omega_a$),必须令角频率的转换常数为
\begin{equation} \beta_\omega = \frac{1}{\beta_t}~, \end{equation}
再将上文的 $\beta_E$ 和 $\beta_t$ 代入可得原子单位下光子能量等于角频率
\begin{equation} E_a = \omega_a~. \end{equation}

2. 薛定谔方程

预备知识 2 薛定谔方程

例 2 无单位的薛定谔方程

   国际单位下的一维薛定谔方程为(式 3

\begin{equation} -\frac{\hbar^2}{2m} \frac{\partial^{2}{\Psi}}{\partial{x}^{2}} + V\Psi= \mathrm{i} \hbar \frac{\partial \Psi}{\partial t} ~, \end{equation}
令 $x = x_a\beta_x$,$m = m_a\beta_m$,$t = t_a\beta_t$,$V = V_a\beta_E$,$\Psi = \Psi_a \beta_\Psi$(注意 $x_a, m_a$ 等都是不带单位的)。 代入式 7 ,各项同除 $\beta_E\beta_\Psi$,得3
\begin{equation} - \left(\frac{\hbar^2}{\beta_m\beta_x^2\beta_E} \right) \frac{1}{2m_a} \frac{\partial^{2}{\Psi_a}}{\partial{x_a}^{2}} + V_a\Psi_a= \mathrm{i} \left(\frac{\hbar}{\beta_E\beta_t} \right) \frac{\partial \Psi_a}{\partial t_a} ~. \end{equation}
使用式 2 式 3 ,发现两个括号中的常数都是 1,于是无单位的薛定谔方程为
\begin{equation} -\frac{1}{2m_a} \frac{\partial^{2}{\Psi_a}}{\partial{x_a}^{2}} + V_a\Psi_a= \mathrm{i} \frac{\partial \Psi_a}{\partial t_a} ~. \end{equation}
这个过程也可以简单认为是把式 7 中的 $\hbar$ 替换为数字 1。

   我们再看波函数的归一化公式

\begin{equation} 1 = \int \left\lvert \Psi \right\rvert ^2 \,\mathrm{d}{x} = \beta_\Psi^2 \beta_x \int \left\lvert \Psi_a \right\rvert ^2 \,\mathrm{d}{x_a} ~, \end{equation}
为了使归一化公式的形式不变,必须令
\begin{equation} \beta_\Psi = \beta_x^{-1/2}~. \end{equation}
同理,对 $N$ 维波函数有 $\beta_\Psi = \beta_x^{-N/2}$。

例 3 动量算符

   国际单位中 $x$ 方向动量的本征方程为

\begin{equation} - \mathrm{i} \hbar \frac{\mathrm{d}}{\mathrm{d}{x}} \psi(x) = p \psi(x)~, \end{equation}
解得本征态(已归一化)为平面波
\begin{equation} \psi(x) = \frac{1}{\sqrt{2\pi\hbar}} \exp\left( \mathrm{i} p x/\hbar\right) ~. \end{equation}
例 2 同理,在原子单位中式 12 变为
\begin{equation} - \mathrm{i} \frac{\mathrm{d}}{\mathrm{d}{x_a}} \psi_a(x_a) = k_a \psi_a(x_a)~, \end{equation}
其中 $k_a$ 就是原子单位中的动量,它和波数 $k_a = 2\pi/\lambda_a$ 数值相同。从这个意义上,原子单位中,动量等于波数,能量等于频率。

   动量本征态式 13 变为

\begin{equation} \psi_a(x) = \frac{1}{\sqrt{2\pi}} \mathrm{e} ^{ \mathrm{i} kx}~. \end{equation}
注意这就是傅里叶变换中的基底。所以在原子单位中,把波函数投影的动量本征态上就是做傅里叶变换

   由以上约束条件,上面出现的所有转换常数中只剩下两个自由的,例如只要确定 $\beta_x$ 和 $\beta_m$,剩下的 $\beta$ 也就确定了。

3. Hartree 原子单位

   原子单位最常见的版本是 Hartree 原子单位,即定义 $\beta_m$ 等于电子的质量,$\beta_x$ 等于玻尔半径(式 4 ),再由式 2 式 11 确定 $\beta_E, \beta_t, \beta_\Psi$,如表 1 所示4。注意许多常数都与氢原子的玻尔模型(原子核不动)的基态(表中简称基态)有关。表中还定义了一些其他的物理量的转换常数,它们的定义可以使以下无单位公式成立(以后我们在不至于混淆的情况下省略角标 $a$)

表1:原子单位转换常数表(参考 “物理学常数”)
物理量 $\beta$ 描述 数值(国际单位)
质量 $m$ $m_e$ 电子质量 $9.10938215 \times 10^{-31} \,\mathrm{kg} $
长度 $x$ $a_0 = \dfrac{4\pi \epsilon_0 \hbar ^2}{m_e e^2}$ 玻尔半径 $5.2917721067 \times 10^{-11} \,\mathrm{m} $
速度 $v$ $\dfrac{\hbar}{m_e a_0}$ 基态电子速度 $2.1876912633 \times 10^{6} \,\mathrm{m/s} $
时间 $t$ $m_e a_0^2/\hbar$ 长度除以速度 $2.418884326 \times 10^{-17} \,\mathrm{s} $
角频率 $\omega$ $\dfrac{\hbar}{m_e a_0^2}$ 时间的倒数 $4.134137335 \times 10^{16} \,\mathrm{s^{-1}} $
能量 $E$ $\dfrac{\hbar^2}{m_e a_0^2} = \dfrac{e^2}{4\pi \epsilon_0 a_0}$ 基态电子势能大小 $4.3597446499 \times 10^{-18} \,\mathrm{J} $
角动量 $L$ $m_e v_0 a_0 = \hbar$ 长度乘以动量 $1.0545718176 \times 10^{-34} \,\mathrm{kg\cdot m^2/s} $
电荷 $q$ $e$ 或 $q_e$ 电子电荷 $1.6021766208 \times 10^{-19} \,\mathrm{C} $
电场强度 $\mathcal{E}$ $\dfrac{e}{4\pi \epsilon_0 a_0^2}$ 基态轨道电场强度 $5.1422067070 \times 10^{11} \,\mathrm{N/C} $
磁感应强度 $B$ $\dfrac{\hbar}{ea_0^2}$ $2.350517567 \times 10^{5} \,\mathrm{T} $
电势 $V$ $\dfrac{e}{4\pi\epsilon_0 a_0}$ 基态轨道电势 $27.211386019 \,\mathrm{J/C} $

   常见公式

\begin{align} &\omega = \frac{2\pi}{T}~,\\ &x = v t~,\\ & \boldsymbol{\mathbf{L}} = m \boldsymbol{\mathbf{r}} \boldsymbol\times \boldsymbol{\mathbf{v}} \qquad \text{(角动量)}~,\\ &\mathcal{E} = \frac{q}{r^2} \qquad \text{(点电荷电场)}~,\\ &\rho_{E} = \frac{1}{8\pi} ( \left\lvert \boldsymbol{\mathbf{\mathcal{E}}} \right\rvert ^2 + \left\lvert c \boldsymbol{\mathbf{B}} \right\rvert ^2) \qquad \text{(电磁场能量密度)}~,\\ & \boldsymbol{\mathbf{F}} = q \boldsymbol{\mathbf{v}} \boldsymbol\times \boldsymbol{\mathbf{B}} \qquad \text{(洛伦兹力)}~,\\ &U = \frac{q}{r} = \mathcal{E} x \qquad \text{(点电荷电势)}~,\\ & V = qU = -q\mathcal{E} x \qquad \text{(匀强电场电势能)}~,\\ & \boldsymbol{\mathbf{\mathcal{E}}} = - \boldsymbol\nabla \varphi - \frac{\partial}{\partial{t}} \boldsymbol{\mathbf{A}} \qquad \text{(标势和矢势)}~,\\ & \boldsymbol{\mathbf{B}} = \boldsymbol{\nabla}\boldsymbol{\times} \boldsymbol{\mathbf{A}} \qquad \text{(标势和矢势)}~,\\ & \boldsymbol{\mathbf{s}} = \boldsymbol{\mathbf{\mathcal{E}}} \boldsymbol\times \boldsymbol{\mathbf{B}} /(4\pi\alpha^2) \qquad \text{(坡印廷矢量)}~,\\ & -\frac{1}{2m} \frac{\partial^{2}{\Psi}}{\partial{x}^{2}} + V\Psi= \mathrm{i} \frac{\partial \Psi}{\partial t} \qquad \text{(薛定谔方程)}~. \end{align}
注意当考察对象为电子时,式中 $m = 1$,可省略。

例 4 匀强电场中电子的薛定谔方程

   令式 27 中 $m = 1$,$q = -1$,再将式 23 代入,得

\begin{equation} -\frac12 \frac{\partial^{2}{\Psi}}{\partial{x}^{2}} + \mathcal{E} x \Psi= \mathrm{i} \frac{\partial \Psi}{\partial t} ~, \end{equation}
其中电场 $\mathcal{E}$ 可以是时间的函数。

习题 1 氢原子的基态能量

   计算原子单位下玻尔模型中氢原子基态的能量,假设原子核固定不动(答案:$-1/2$)。

电磁学相关

   把精细结构常数记为 $\alpha$。国际单位中有 $a_0 = \hbar/(\alpha c m_e)$。原子单位中光速时精细结构常数的倒数

\begin{equation} c = \frac{1}{\alpha} \approx 137.036~. \end{equation}
另外 $\epsilon_0 = 1/(4\pi)$,$\mu_0 = 4\pi\alpha^2$,$\mu_0\epsilon_0 = \alpha^2$。对平面电磁波 $E = cB$。

   注意原子单位制和经典电磁学中的高斯单位制是两种不兼容的单位制,不能混用。高斯单位制中 $\hbar\ne 1$。

4. 其他原子单位

   当问题涉及一主要角频率 $\omega_0$ 的时候(例如研究原子在单频激光中的变化),可选择 $\beta_E = \hbar\omega_0$ 做能量单位。同样令 $\beta_m$ 等于电子质量,$\beta_q$ 等于元电荷,由式 2

\begin{equation} \beta_x = \sqrt{\frac{\hbar}{m_e\omega_0}}~, \qquad \beta_t = \frac{1}{\omega_0}~. \end{equation}
为了使式 23 成立,得电场的转换常数为
\begin{equation} \beta_\mathcal{E} = \frac{\hbar\omega_0}{e \beta_x}~. \end{equation}
一种常见的情况是平面电磁波中的电场用国际单位表示为 $\mathcal{E}(t) = \mathcal{E}_0 \cos\left(\omega_0 t\right) $,而这种原子单位下为
\begin{equation} \mathcal{E}(t) = \mathcal{E}_0\cos t~, \end{equation}
注意右边不含 $\omega$,形式更简洁。

   另一个常见的例子是简谐振子,若使用其振动的固有频率来定义 $\beta_t$,其薛定谔方程为

\begin{equation} -\frac12 \frac{\partial^{2}{\Psi}}{\partial{x}^{2}} + \frac12 x^2 \Psi= \mathrm{i} \frac{\partial \Psi}{\partial t} ~, \end{equation}
能级为
\begin{equation} E_n = \frac12 + n \qquad (n = 0, 1, 2\dots)~, \end{equation}
归一化的基态波函数为
\begin{equation} \psi_0(x) = \pi^{-1/4} \mathrm{e} ^{-x^2/2}~. \end{equation}

5. 转换为国际单位制

   使用 “无量纲的物理公式” 中的方法,容易把原子单位的物理量和公式换回到国际单位,前提是要明确每个物理量的量纲。

例 5 薛定谔方程

   要将薛定谔方程式 33 转换为国际单位的公式。即先把所有无单位的物理量替换成有单位的物理量除以对应的 $\beta$ 常数,得(两边已同乘 $\beta_\Psi$)

\begin{equation} -\beta_x^2\frac12 \frac{\partial^{2}{\Psi}}{\partial{x}^{2}} + \frac{1}{\beta_x^2}\frac12 x^2 \Psi= \beta_t \mathrm{i} \frac{\partial \Psi}{\partial t} ~. \end{equation}
式 30 代入,两边乘以 $\omega\hbar$ 得国际单位下的简谐振子薛定谔方程
\begin{equation} -\frac{\hbar^2}{2m} \frac{\partial^{2}{\Psi}}{\partial{x}^{2}} + \frac12 m\omega^2 x^2 \Psi= \mathrm{i} \hbar \frac{\partial \Psi}{\partial t} ~. \end{equation}
类似地,也可以将式 34 变为
\begin{equation} E = \left(\frac12 + n \right) \beta_E = \left(\frac12 + n \right) \omega\hbar \qquad (n = 0, 1, 2\dots)~. \end{equation}


1. ^ 参考 Wikipedia 相关条目。以及 [1] 的附录。
2. ^ 例如在扇形面积公式 $S = \theta R^2/2$ 中的 $\theta$ 可以看作具有单位 “Rad”,但面积的单位却只是 “$ \,\mathrm{m^2} $” 而无需记为 “$ \,\mathrm{Rad\cdot m^2} $”。
3. ^ 根据偏微分的定义,常数可以移到偏微分算符外,如 $ \partial^{2}/\partial {(\beta_x x_a)}^{2} = (1/\beta_x^2) \partial^{2}/\partial {x_a}^{2} $
4. ^ 为了区别能量与电场,以下用 $E$ 表示能量,用 $\mathcal{E}$ 表示电场。


[1] ^ Bransden, Physics of Atoms and Molecules, 2ed

                     

© 小时科技 保留一切权利