OpenMX Ver. 3.9 では、全エネルギーを各基底関数の寄与に分解することが可能です。 これは、基底関数への射影によって行われます [68]。 全エネルギーの軌道分解は、以下のキーワードを加えるだけで実施可能です。
Energy.Decomposition on # on|off, default=off
例として、メタン分子の全エネルギーの軌道分解を以下に示します。
% mpirun -np 5 openmx Methane_ED.dat > met_ed.std &
ここで使用している「Methane_ED.dat」は「work」ディレクトリに収められています。
計算が終わると「met_ed.out」が得られ、以下の様な情報が書き出されます。
*************************************************************************
*************************************************************************
Decomposed energies in Hartree unit
Utot = Utot(up) + Utot(dn)
= Ukin(up) + Ukin(dn) + Uv(up) + Uv(dn)
+ Ucon(up)+ Ucon(dn) + Ucore+UH0 + Uvdw
Uele = Ukin(up) + Ukin(dn) + Uv(up) + Uv(dn)
Ucon arizes from a constant potential added in the formalism
up: up spin state, dn: down spin state
*************************************************************************
*************************************************************************
Total energy (Hartree) = -8.216132481346387
Decomposed.energies.(Hartree).with.respect.to.atom
Utot Utot(up) Utot(dn) Ukin(up) Ukin(dn) Uv(up) ....
1 C -6.132295762765 -3.066148 -3.066148 2.076016 2.076016 -2.957459 ....
2 H -0.520959186503 -0.260480 -0.260480 0.300675 0.300675 -0.499086 ....
3 H -0.520959174111 -0.260480 -0.260480 0.300675 0.300675 -0.499086 ....
4 H -0.520959173764 -0.260480 -0.260480 0.300675 0.300675 -0.499086 ....
5 H -0.520959184204 -0.260480 -0.260480 0.300675 0.300675 -0.499086 ....
Decomposed.energies.(Hartree).with.respect.to.atomic.orbital
1 C Utot Utot(up) Utot(dn) Ukin(up) Ukin(dn) Uv(up)
multiple
none -4.483770 -2.241885 -2.241885 0.000000 0.000000 0.000000 ....
s 0 -0.675699 -0.337849 -0.337849 0.203145 0.203145 -0.556473 ....
s 1 0.003690 0.001845 0.001845 0.036240 0.036240 -0.034310 ....
px 0 -0.325884 -0.162942 -0.162942 0.496144 0.496144 -0.673031 ....
py 0 -0.325912 -0.162956 -0.162956 0.496166 0.496166 -0.673068 ....
pz 0 -0.325884 -0.162942 -0.162942 0.496144 0.496144 -0.673031 ....
px 1 0.005096 0.002548 0.002548 0.107318 0.107318 -0.104552 ....
...
..
.
全エネルギーが各原子の各軌道からの寄与に分解されているのが分かります。
不純物原子や空孔による局所的な摂動が隣接原子の安定性/不安定性にどのように影響するかを
議論する際に、本機能は有用です。
また、本機能は幅広い現象に対する物理的な描像を解析するのに役立つものと期待されます。
しかしながら、この機能の公開はまだ試験的なものです。
我々はこの機能を近い将来にさらに発展させることを考えています。