24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2728  |  回复: 18

fzuvivi

金虫 (小有名气)

[求助] 结构超胞后替位掺杂,如何加快超胞的弛豫速度?在线等回复~

各位,我试图对原胞扩2×2×2的胞后,进行替代原子2%的微量掺杂,如图选定原子位置。众所周知,这样掺杂后,超胞的对称性降低,从而导致掺杂后的胞在弛豫时收敛十分缓慢,甚至不收敛。望各位有好的方法指导!谢谢了~
我的INCAR中已经增加了号称对弛豫有利的maxing参数。如:
Amix = 0.2
Bmix = 0.0001
Amix_MAG = 0.8
Bmix_MAG = 0.0001

PS:体系具有金属磁性。另外,扩胞之前,原胞已充分弛豫(EDIFF=1E-5;EDIFFG=-1E-4);扩胞后再掺杂,设置弛豫精度为:EDIFF=1E-4;EDIFFG=-1E-3.



[ Last edited by fzuvivi on 2013-1-4 at 09:34 ]
回复此楼

» 猜你喜欢

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

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 15:40:05
有可能是你的参数设置的不合理    体相的计算收敛一般都是很快的
你说的”弛豫收敛十分缓慢“ 是每个离子步耗时很大 还是跑了上百步了还是不收敛? 如果每个离子步耗时很大 那应该是你k点设置太大了。如果是后者,这可能是参数设置不合理   另外,mixing方法不是万能的,有时候使用它反而不好收敛
2楼2013-01-04 12:52:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)

引用回帖:
2楼: Originally posted by emilyoyang at 2013-01-04 12:52:01
有可能是你的参数设置的不合理    体相的计算收敛一般都是很快的
你说的”弛豫收敛十分缓慢“ 是每个离子步耗时很大 还是跑了上百步了还是不收敛? 如果每个离子步耗时很大 那应该是你k点设置太大了。如果是后者, ...

谢谢回复~其实在掺杂之前,收敛性是很好的。基本上20步左右一定能收敛的。可是微量替位掺杂后,胞的对称性变低了。然后收敛很慢,就是需要80步都不一定能收敛。不知道您收的参数不合理指的是哪一块?
K点的话,我还是用掺杂之前的,因为怕减低K点的话,不利于做微量掺杂。不知道是不是跟我设置的弛豫方式有关系(IBRION=1;POTIM=0.2 );这两天也在尝试做不加mixing的弛豫,等结果中。。。
3楼2013-01-05 10:45:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cclw

新虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 15:40:27
引用回帖:
3楼: Originally posted by fzuvivi at 2013-01-05 10:45:28
谢谢回复~其实在掺杂之前,收敛性是很好的。基本上20步左右一定能收敛的。可是微量替位掺杂后,胞的对称性变低了。然后收敛很慢,就是需要80步都不一定能收敛。不知道您收的参数不合理指的是哪一块?
K点的话,我 ...

修改mixing参数一般是针对电子迭代不收敛的情况,如果一个离子步里电子迭代很多步不收敛,表示电子密度混合的方式可能不合理,可以ALGO=38或fast,如果不能搞定再,修改mixing参数;
如果是跑了很多个离子步,体系很难弛豫到基态构型,很可能是算法不够好,可以调IBRION和POTIM,观察每一步的能量,如果能量下降的趋势很慢,POTIM可以设大点;如果有明显的震荡,改小POTIM。
另外仅做结构优化的话,楼主的收敛精度太细了。初步的弛豫,EDIFFG=-0.1就够了。
4楼2013-01-07 08:10:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)

引用回帖:
4楼: Originally posted by cclw at 2013-01-07 08:10:19
修改mixing参数一般是针对电子迭代不收敛的情况,如果一个离子步里电子迭代很多步不收敛,表示电子密度混合的方式可能不合理,可以ALGO=38或fast,如果不能搞定再,修改mixing参数;
如果是跑了很多个离子步,体系 ...

非常感谢cclw~看来唯有的办法就只有是弛豫,弛豫再弛豫。不能一步求登天啊~呵呵
5楼2013-01-07 08:38:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)

饼干队长

本帖仅楼主可见
6楼2013-01-07 10:22:50
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

fzuvivi

金虫 (小有名气)

引用回帖:
6楼: Originally posted by 旭旭宝宝 at 2013-01-07 10:22:50
请问,什么情况下建立超包?  超包在计算后 处理方法 和原包有何不同?  我 新手 讨论下

我的情况是要做替位掺杂,所以需要将原胞扩大掺杂。超胞前后的处理方法基本大同小异。
7楼2013-01-07 15:57:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 21:44:08
引用回帖:
3楼: Originally posted by fzuvivi at 2013-01-05 10:45:28
谢谢回复~其实在掺杂之前,收敛性是很好的。基本上20步左右一定能收敛的。可是微量替位掺杂后,胞的对称性变低了。然后收敛很慢,就是需要80步都不一定能收敛。不知道您收的参数不合理指的是哪一块?
K点的话,我 ...

你这种情况下 采用IBRION=2 效果要比IBRION=1好点 你可以试试
8楼2013-01-07 16:21:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖


franch: 金币+1, 谢谢回帖交流,, 2013-01-07 21:44:20
引用回帖:
3楼: Originally posted by fzuvivi at 2013-01-05 10:45:28
谢谢回复~其实在掺杂之前,收敛性是很好的。基本上20步左右一定能收敛的。可是微量替位掺杂后,胞的对称性变低了。然后收敛很慢,就是需要80步都不一定能收敛。不知道您收的参数不合理指的是哪一块?
K点的话,我 ...

另外, k点稍微减低点是可以的  计算能量是再提高到正在k点
9楼2013-01-07 16:22:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzuvivi

金虫 (小有名气)

引用回帖:
9楼: Originally posted by emilyoyang at 2013-01-07 16:22:20
另外, k点稍微减低点是可以的  计算能量是再提高到正在k点...

谢谢回复~
我弛豫时,试过IBRION=2,算20步离子步就自动停掉。。
K点减低后,优化得到的结果,你在算能量的时候再提高K点,看OUTCAR时,你会发现,体系的stress force还是比较大的。所以,我不太认同你的想法。。。
10楼2013-01-07 17:12:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fzuvivi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕311分求调剂 +17 梁富贵险中求 2026-04-04 18/900 2026-04-08 22:53 by may_新宇
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 311求调剂 +5 surte 2026-04-08 11/550 2026-04-08 22:10 by 土木硕士招生
[考研] 求调剂 一志愿西南交通大学085701环境工程 282分 +15 多多爱吃汉堡 2026-04-04 16/800 2026-04-08 11:39 by i_cooler
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 一志愿南科大生物学297分,求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 一志愿北京化工085600 310分求调剂 +20 0856材料与化工3 2026-04-04 22/1100 2026-04-07 15:14 by 上岸快快
[考研] 一志愿北京2,材料与化工308求调剂 +16 熊二想上岸 2026-04-04 17/850 2026-04-07 14:58 by BruceLiu320
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 生物学求调剂 +5 15064154688 2026-04-03 5/250 2026-04-06 11:56 by lijunpoly
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[考研] 301求调剂 +18 骆驼男人 2026-04-02 18/900 2026-04-04 20:33 by 蓝云思雨
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 一志愿沪985,326分求调剂 +3 刘墨墨 2026-04-03 3/150 2026-04-04 11:16 by 悲伤的芋头
[考研] 26调剂 086003 +6 失活的细胞 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 生物学硕341求调剂 +4 你笑起来像云朵 2026-04-03 4/200 2026-04-03 10:32 by macy2011
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
信息提示
请填处理意见