薛定谔方程的数值方法

计算科学 量子力学
2021-12-05 03:03:30

我们正在比较各种数值方法的性能,这些方法可用于求解氢原子与强激光脉冲(太强而无法使用微扰方法)相互作用的薛定谔方程。当对径向部分使用离散化方案时,似乎大多数(所有人)将原子放在一个盒子里,只是将半径切掉某个较大的值并求解这些基组。这与将径向变量映射到有限域然后离散化该域(在此过程中,丢弃大部分可用的基组)相比如何?似乎没有人这样做是有原因的吗?

1个回答

贝克等人。1994 年提出了用于原子和分子电子结构计算的径向网格的这种映射。它仍然用于现代电子结构代码,例如FHI-AIMS使用它们,如最近的一篇论文中所述。

即使有这样的映射,同样的问题仍然存在:如果在最外面的网格点之外发生了一些有趣的事情,你就会错过它。然而,这些映射确实具有可以系统地改进网格以包含远处网格点的优点。(这在最近的 FHI-AIMS 论文的第 4.1 节中进行了解释)。