RAtom:解决中性原子的非线性Kohn-Sham方程。-开源

RAtom解决了中性原子的非线性Kohn-Sham方程。 实现了基于有限元方法的自适应算法。 微分特征值问题的离散化是通过以Lobatto多项式为基函数的有限元方法完成的。 为了获得原子的总能量，以1E-6 hartree的绝对准确度，使用了高阶高斯求积。 离散化导致广义特征值问题，可通过LAPACK库中的过程解决。 该程序是用C ++实现的，它作为一个线程运行。 下面的论文中描述了实现的算法：1. Z. Romanowski，“原子的Kohn-Sham方程的自适应求解器”，Modeling Simul。 母校科学。 卷17，17，045001（2009）2. Z. Romanowski，“ h自适应，高阶有限元方法在求解径向Schrödinger方程中的应用”，《分子物理学》，第1卷。 107，第1339-1348页（2009）。 进一步的问题：romz@wp.pl