量子三体问题： 理论和数学

量子力学三体问题

• 基于此篇文章: 库伦量子三体问题，前半部分是数学推导，后半部分是数值计算模拟。理论部分来源于项武义教授的研究成果

• 对于库伦作用的三体问题，薛定谔方程是个9D的偏微分方程，考虑到两个系统守恒量：动量，利用Jacobi坐标可以去除质心坐标运动，9D->6D；角动量，可以利用($L^2,Lz$)算符降维,6D->3D！基于转动不变量的耦合偏微分方程组具有内在的球对称性，利用球坐标，球对称性变得明显，进而根据角度和径向算符的特征函数，耦合偏微分方程组被进一步化为矩阵特征值问题

• 最终化为一个无限维的矩阵特征值求解问题，矩阵结构超级稀疏，相关矩阵可以表示成Kronecker直积$A\otimes B$形式，矩阵维度由三个截断整数参数(m,n,p)控制，(m,n,p)分别代表展开基函数的下标指数

推导过程：公式

• 基本公式。记 { M i , Z i , i = 1 , 2 , 3 } \{M_i,Z_i, i=1,2,3\} {Mi​,Zi​,i=1,2,3} 为三个粒子的质量和电荷量，无量纲的三体带电粒子的薛定谔方程为: $$-\frac{1}{2M}\Delta \Psi +V\Psi =E\Psi$$ 其中$\Delta ={\sum }_{i=1}^3\frac{{\Delta }_i}{m_i}$, 其中$m_i=\frac{M_i}{M}$ 是粒子的质量比，库伦势能$$V=\frac{Z_2{Z}_3}{|{\mathbf{r}}_2-{\mathbf{r}}_3|}+\frac{Z_3{Z}_1}{|{\mathbf{r}}_3-{\mathbf{r}}_1|}+\frac{Z_1{Z}_2}{|{\mathbf{r}}_1-{\mathbf{r}}_2|}$$整个系统由6个参数 { m i , Z i , i = 1 , 2 , 3 } \{m_i,Z_i, i=1,2,3\} {mi​,Zi​,i=1,2,3}完全确定。
• 移除质心坐标。采用雅可比坐标
  $$\{\begin{array}{l}\mathbf{R}=m_1{\mathbf{r}}_1+m_2{\mathbf{r}}_2+m_3{\mathbf{r}}_3\\ \mathbf{x}=\sqrt{\frac{m_1}{1-m_1}}(m_2({\mathbf{r}}_1-{\mathbf{r}}_2)+m_3({\mathbf{r}}_1-{\mathbf{r}}_3))\\ \mathbf{y}=\sqrt{\frac{m_2{m}_3}{1-m_1}}({\mathbf{r}}_2-{\mathbf{r}}_3)\end{array}$$
  此时动能算子$$\Delta ={\Delta }_{\mathbf{R}}+{\Delta }_{\mathbf{x}}+{\Delta }_{\mathbf{y}}$$在质心坐标系下$$\Delta ={\Delta }_{\mathbf{x}}+{\Delta }_{\mathbf{y}}$$角动量算子$$\mathbf{L}={\mathbf{L}}_{\mathbf{x}}+{\mathbf{L}}_{\mathbf{y}}=-i(\mathbf{x}\times {▽}_{\mathbf{x}}+\mathbf{y}\times {▽}_{\mathbf{y}})$$
• 转动不变量。对两个向量$\mathbf{x}$和$\mathbf{y}$,可以定义三个转动不变量$$f_1=\mathbf{x}\cdot \mathbf{x},f_2=\mathbf{y}\cdot \mathbf{y},f_3=\mathbf{x}\cdot \mathbf{y}$$进一步采用Jacobi球坐标
  $$\begin{gathered} \rho =\sqrt{f_1+f_2}, \\ \cos \theta =2(f_1{f}_2-f_3^2{)}^{1/2}(f_1+f_2{)}^{-1}, \\ \tan \beta =2f_3(f_2-f_1{)}^{-1},\sin \beta =2f_3[(f_2-f_1{)}^2+4f_3^2{]}^{-1/2} \\ (0\le \rho <\infty ,0\le \theta \le \pi /2,-\pi \le \beta \le \pi ) \end{gathered}$$记$\xi =\sin \theta$,反变换关系为
  $$f_1=\frac{{\rho }^2}{2}(1-\xi \cos \beta ),f_2=\frac{{\rho }^2}{2}(1+\xi \cos \beta ),f_3=\frac{{\rho }^2}{2}\xi \sin \beta$$
• 粒子间距离表示.
  $$\{\begin{array}{l}|{\mathbf{r}}_{23}{|}^2=k_1{\rho }^2(1+\xi \cos \beta ),\\ |{\mathbf{r}}_{31}{|}^2=k_2{\rho }^2(1+\xi \cos (\beta +2{\beta }_2)),\\ |{\mathbf{r}}_{12}{|}^2=k_3{\rho }^2(1+\xi \cos (\beta +2{\beta }_3)),\end{array}$$具体推导和参数定义可在论文中找到
• 波函数表示
  
  具体细节请参考文章。该波函数是角动量算子${\mathbf{L}}^2,{\mathbf{L}}_3$,和Parity算子的本征函数，依次对应特征值$l(l+1),l$ 和$(-1{)}^{l+\lambda }$. 带入薛定谔方程，利用角动量基函数的正交性，得到一组转动不变量分量的耦合方程
  
  至此，PDE的维度已经降为3维. 利用Jacobi球坐标$(\rho ,\xi ,\beta )$, 耦合方程进一步化为
  
  其中
  
  
  算符$A^l$的本征函数是
  
  对应本征量
  
  具体推导细节参考文章。在算子$A^l$的本征函数空间，不变量分量可以进一步展开
  
  算子$A^l,B_1,B_2$都已经被对应矩阵表示，分量只依赖于坐标$\rho$。其中库伦势能函数
  
  此处
  
  $U(\xi ,\beta )$ 有三个奇点：$(1,\pi )$,$(1,\pi -2{\beta }_2)$,和$(1,\pi -2{\beta }_3)$。这三处奇点完全刻画了库伦三体问题的本质，会显著影响数值计算的收敛性(需要足够大的本征函数空间，才能更好的描述这些cusp尖点)。在质心坐标系下${\rho }^2=m_1|{\mathbf{r}}_1{|}^2+m_2|{\mathbf{r}}_2{|}^2+m_3|{\mathbf{r}}_3{|}^2$,刻画了系统的惯性质量性质。此外库伦势中，$\rho$和$(\xi ,\beta )$是可以分离的，为进一步化简提供可能。
• 几个特殊情况
  1. $l=0$
     
  2. $l=1$, 偶Parity
     
  3. $l=1$, 奇Parity
     
     借助标记
     
     耦合方程可写成更简洁形式
     
• $l\ge 2$情况，借助矩阵标记，耦合方程可以写成
  
  符号代表的具体含义可在论文中找到。至此我们可以把关注重心聚焦于涉及$\rho$变量的微分算符，接下来会找一组对应的$\rho$本征函数,耦合方程会在该本征函数集合上进一步展开，最终化为一个无限维的求解矩阵特征值问题

计算过程：数值设计

计算源码

应用

中性的氢原子会和一个电子结合成束缚态(负能量)系统，有两个束缚态，其中基态能量：

唯一偶宇称束缚态能量:

该状态中相当于两个P态的电子组合，电子之间的作用很显著，变分方法需要的基函数数目会比基态大很多，对能量估计可以比较准确，但对真实波函数的估计精度会下降很多，尤其考虑质子有限质量情况下

---

附注

Helium原子

一个氦原子核和两个电子组成的三体系统。假定氦原子核质量无穷大。两个电子坐标分别记为${\mathbf{r}}_1$和${\mathbf{r}}_2$, 雅可比作坐标取为:
$$\begin{gathered} \mathbf{x}=\frac{{\mathbf{r}}_1-{\mathbf{r}}_2}{\sqrt{2}} \\ \mathbf{y}=\frac{{\mathbf{r}}_1+{\mathbf{r}}_2}{\sqrt{2}} \end{gathered}$$