加载中…[转载]如何在VASP计算中指定自旋多重度
(2014-05-23 21:28:49)
标签：
转载 |
原文地址：如何在VASP计算中指定自旋多重度作者：原子修补匠
在用gaussian时,如果分子的基态不是单重态，可以通过在title
section里指定spin来计算其他组态如三重态的能量。
比如以下的计算指定氧原子为单重态
--------------------------
#
hf/6-31G(d)
Oxygen atom
0 1
O 0.0000 0.0000 0.0000
--------------------------
而以下的输入则指定氧原子为三重态
--------------------------
#
hf/6-31G(d)
Oxygen atom
0 3
O 0.0000 0.0000 0.0000
--------------------------
在VASP中，控制自旋的参数是INCAR中的ISPIN,
ISPIN=1表示进行非自旋极化的计算，而ISPIN=2表示进行自旋极化的计算。为了进一步得到我们所需要的组态，可以通过NUPDOWN指定自旋向上和自旋向下的电子数差，此方法与直接指定自旋多重度是一样的。
VASP里，C的参考态是S2P2，即C的基态，ISPIN=2应该是这个态。
值得注意的是，在进行固定自旋多重态的计算中，如果初始的电荷密度的自旋极矩和INCAR中设置的NUPDOWN不同的话，收敛会很慢。如果是从读入的初始波函数开始计算，则没有类似问题。
以下是vasp手册上的原话：
There is
a word of caution required: If NUPDOWN is set in the INCAR file the
initial moment for the charge density should be the same. Otherwise
convergence can slow down. When starting from atomic charge density
(ICHARG=2), VASP will try to do this automatically by setting
MAGMOM to NUPDOWN/NIONS. The user can of course overwrite this
default by specifying a different MAGMOM (which should still result
in the correct total moment). If one starts from the wavefunctions,
the initial moment will be always correct, because VASP will
``push'' the required number of electrons from the down to the up
component. If starting from a chargedensity supplied in the CHGCAR
file (ICHARG=1), the initial moment is usually
incorrect!
喜欢
0
赠金笔