24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jmlv

金虫 (著名写手)

应助: 50 (小学生)
金币: 4640.2
散金: 20
红花: 6
帖子: 1357
在线: 129小时
虫号: 83527
注册: 2005-07-30
性别: GG
专业: 催化化学

[求助] Ni团簇的磁矩为什么总是波动？已有1人参与

我创建了一个气态游离(free standing)、79个原子的截角八面体Ni团簇，对其进行自旋极化计算是，出现咄咄怪事：

第一步，对从体相截出的Ni79团簇进行完全弛豫，收敛门槛设为EDIFF=1E-5，得到磁矩83.10，总能量-121.724 eV。

第二步，对第一步得到的优化构型进行完全弛豫，收敛门槛仍设为EDIFF=1E-5，按理说应该得到同样的磁矩和总能量，但事与愿违：磁矩82.64，总能量-137.479 eV。

总能量差别如此之大，令人无法信任也无法使用这种结果呀。有达人碰到过类似的问题吗？求解！谢谢。

附上我的INCAR参数：
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ENCUT  = 400.00 eV
PREC = Normal
LREAL  = Auto

EDIFF  = 1E-5
NSW = 1000
IBRION = 2
ISIF = 2

ISMEAR = 2
SIGMA  = 0.2

ISPIN  = 2             # spin polarized

# DFT+U settings
LDAU = .TRUE.    # switch on the +U corrections
LDAUTYPE = 2       # Dudarev's approach
LDAUL  = 2       # -1 no U; 1 for p electrons; 2 d and 3 f.
LDAUU  = 6.0    # U values
LDAUJ  = 0.0    # do not specify them for Dudarev's approach

# output controls
LCHARG = .FALSE.
LWAVE  = .FALSE.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

回复此楼

» 猜你喜欢

三无产品还有机会吗已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有7人回复
压汞仪和BET测气凝胶孔隙率已经有4人回复
博士申请都是内定的吗？已经有14人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复
投稿精细化工已经有6人回复

1楼2015-11-25 09:54:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jmlv

金虫 (著名写手)

应助: 50 (小学生)
金币: 4640.2
散金: 20
红花: 6
帖子: 1357
在线: 129小时
虫号: 83527
注册: 2005-07-30
性别: GG
专业: 催化化学

引用回帖:

5楼: Originally posted by WanderingHeart at 2015-11-26 11:06:26
如果收敛条件是原子所受的最小的力的话，好理解，只要它小于一个阈值就可以认为收敛到稳定结构（当然这也仅仅是local minimal）；那么再算的时候程序会不会在最开始判断这个力已经满足收敛条件，然后就不开始计算呢 ...

我的经验，使用能量作为收敛的判据，和使用力作为收敛的判据，并无优劣之分。

能量作为收敛判据，是说如果最后两个或若干个离子步的能量差低于用户设定值，则弛豫停止。而这两个离子步的差异，来自于两者根据原子间力而移动原子位置后的结构不同。换言之，能量差越来越小，也同时意味着内部应力越来越小。这两者是同步的。