24小时热门版块排行榜    

查看: 1311  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jmlv

金虫 (著名写手)

[求助] 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.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
回复此楼

» 猜你喜欢

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

WanderingHeart

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你的收敛门槛是针对哪个物理量的?
3楼2015-11-25 12:45:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

wangsoft

木虫 (正式写手)

你再对第二步得到的结构迟豫会有啥结果?

发自小木虫Android客户端
2楼2015-11-25 10:03:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

引用回帖:
3楼: Originally posted by WanderingHeart at 2015-11-25 12:45:20
你的收敛门槛是针对哪个物理量的?

EDIFF=1E-5是针对能量的,也可以用EDIFFG来针对力。以我经验,二者对结果影响不大。
4楼2015-11-26 10:08:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WanderingHeart

铁杆木虫 (著名写手)

引用回帖:
4楼: Originally posted by jmlv at 2015-11-26 10:08:29
EDIFF=1E-5是针对能量的,也可以用EDIFFG来针对力。以我经验,二者对结果影响不大。...

如果收敛条件是原子所受的最小的力的话,好理解,只要它小于一个阈值就可以认为收敛到稳定结构(当然这也仅仅是local minimal);那么再算的时候程序会不会在最开始判断这个力已经满足收敛条件,然后就不开始计算呢?

如果收敛条件是能量,我不太理解它是怎么判断结构优化收敛的。
5楼2015-11-26 11:06:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211 初试270分 求调剂 +5 谷雨上岸 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[考研] 299求调剂 +7 某某某某位 2026-03-21 7/350 2026-03-24 15:24 by cuifj
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +5 niko- 2026-03-22 5/250 2026-03-23 22:01 by fuyu_
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap-
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
信息提示
请填处理意见