24小时热门版块排行榜    

查看: 4705  |  回复: 21

sunyang1988

至尊木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
bingmou(金币+1):感谢分享 2010-12-08 19:05:09
引用回帖:
Originally posted by lgxyz at 2010-12-08 14:50:44:


弹簧力是沿着两个image连线的切线方向,是为了保证每个image之间的联系;而这里讲到的力是构型在势能面上,是为了让这个构型优化为局域极小;这二者同时作用才能找到一级鞍点。我是这样理解的,不对之处,请多 ...

按道理说应该是这样的,但VASP自带的EB计算如果取力的收敛标准不知道它指的力是什么力,能量标准时几个镜像点能量之和,如果受力也是几个镜像点之和弹簧力应该能够消掉。但我算过几次无一例外都是在采用力的收敛标准时结果出现了严重错误。用CLNEB是可以的
具体为什么会这样还没搞懂
心之所向无惧无悔
11楼2010-12-08 15:44:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

嗯,我说的就是CI-NEB
12楼2010-12-08 19:05:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wqg1806

木虫 (著名写手)

科研农民工

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-12-08 21:22:55
我以前也遇到过,我觉得是路径设计的不合理。中间可能还要经过一些过渡态
一路走到黑
13楼2010-12-08 21:07:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

引用回帖:
Originally posted by wqg1806 at 2010-12-08 21:07:27:
我以前也遇到过,我觉得是路径设计的不合理。中间可能还要经过一些过渡态

多谢!我隐约也觉得是这样的!
14楼2010-12-08 21:20:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

04nylxb

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by hakuna at 2010-12-08 13:39:49:


不对吧,可以用能量标准的。CI-NEB收敛后一般情况下能量高的image就是过渡态了,个别情况会是二阶鞍点,消掉多出来的虚频就可以了,不理解你为什么那样做....

请教,如何消掉多出来的虚频?非常感谢
集中精力发文章
15楼2010-12-09 08:57:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
sunyang1988(金币+1):谢谢指点 2010-12-09 11:21:46
怎么调我就不知道了,但NEB确实是以image受的力为收敛标准的。NEB的力是包含两部分:切线方向的弹性和垂直于切线方向的真实力。
16楼2010-12-09 11:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-12-10 00:10:27
引用回帖:
Originally posted by sunyang1988 at 2010-12-08 14:00:45:

请教,如何消掉多出来的虚频?非常感谢

这没什么标准方法,经验,直接调坐标、部分优化......

[ Last edited by hakuna on 2010-12-10 at 08:16 ]
17楼2010-12-09 23:47:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿克飞顿

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
引用回帖:
7楼: Originally posted by hakuna at 2010-12-08 13:39:49
不对吧,可以用能量标准的。CI-NEB收敛后一般情况下能量高的image就是过渡态了,个别情况会是二阶鞍点,消掉多出来的虚频就可以了,不理解你为什么那样做.......

您好,请问您一个问题。
我用CI-NEB找过渡态,分成三步计算的,第一步EDIFFG设为-0.5,经历二三十步达到力的收敛判据后,第二部设为-0.1.一般要经历100多步才能达到力收敛,第三部设为-0.05一般要两三百才能收敛。期间我用nebconverge.pl查看力的收敛趋势,发现中间多次出现跳跃,即本来已经收敛到-0.07了,结果突然跳跃到-0.3然后又慢慢收敛下来。
请问我这样得到的过渡态可信吗???
另外加入我插入了四个构型,1 2 4都达到力判据了,3还差-0.01左右就能达到连判据,而2才是想要的过渡态,这是我能不能就停止计算,用2作为我的过渡态或者把2的构型再拿出去结构优化完作为过渡态??
18楼2015-09-10 10:54:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
18楼: Originally posted by 阿克飞顿 at 2015-09-10 10:54:01
您好,请问您一个问题。
我用CI-NEB找过渡态,分成三步计算的,第一步EDIFFG设为-0.5,经历二三十步达到力的收敛判据后,第二部设为-0.1.一般要经历100多步才能达到力收敛,第三部设为-0.05一般要两三百才能收敛。 ...

1.TS搜索过程中,力或能量出现跳跃应该都是正常现象,搜索TS本质上也是优化,而目前的优化都是与梯度相关的局域搜索方法,所以出现跳跃很正常.......
2.建议还是等NEB收敛后再决定

» 本帖已获得的红花(最新10朵)

19楼2015-09-10 11:24:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿克飞顿

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
引用回帖:
19楼: Originally posted by hakuna at 2015-09-10 11:24:44
1.TS搜索过程中,力或能量出现跳跃应该都是正常现象,搜索TS本质上也是优化,而目前的优化都是与梯度相关的局域搜索方法,所以出现跳跃很正常.......
2.建议还是等NEB收敛后再决定...

那这种要经过两三百步才能达到力判据收敛的状况正常吗?为什么有的人说超过五六十就说明没必要继续算了呢?它们说的是电子步吧!
这是我的INCAR:
SYSTEM = TS1
ENCUT = 400
EDIFF  = 1E-4
EDIFFG = -0.05
GGA=PE
ISTART = 0 ; ICHARG = 2
ISMEAR = 0;SIGMA=0.2
ISIF=2
LREAL = AUTO
NSW=1000
IBRION=3
POTIM=0.0
NELM = 500
IALGO=48
LDIAG=.TRUE.
ALGO = Fast
LWAVE = .TRUE.
LCHARG = .TRUE.
NPAR = 5;NSIM = 4
LPLANE = .TRUE.
LSCALU = .TRUE.
IOPT=3
ICHAIN=0
MAXMOVE=0.2
SPRING=-5
IMAGES=3
LCLIMB = .TRUE.

KPOINT
我用的2×2×1
体系有80多个原子共七层,下面四层固定

发自小木虫Android客户端
20楼2015-09-10 11:49:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgxyz 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +8 调剂的考研学生 2026-03-09 8/400 2026-03-15 22:14 by Winj1e
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 337一志愿华南理工材料求调剂(有希望2吗?) +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 求调剂! +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 数二英二309分请求调剂 +3 dtdxzxx 2026-03-09 4/200 2026-03-09 19:56 by yuningshan
信息提示
请填处理意见