应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp NEB 如何选取原胞大小 ？
111/2返回列表12下一页
查看: 3468  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fangyongxinxi

新虫 (正式写手)

[求助] vasp NEB 如何选取原胞大小 ？

如题，想做一个体系两种结构 A B 之间的过渡态情况。A B 两个体系最小原胞是6个原子，自己做了测试：
将最小原胞设置成过渡态搜索的初始和终态，得到一种过渡态。
然后，将A B 都扩展成18个原子的体系，得到另一个过渡态。
问题：这两种过渡态得到的趋势和过渡态结构都不一样，为什么？自己的想法，想是不是体系初始结构的大小，会严重影响过渡态的计算 ？对于过渡态计算，体系应该选多大合适些 ？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-06-15 10:45:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胡同学

金虫 (正式写手)
关注
回复此楼
人生苦短，抓紧科研，抓紧找对象。
2楼2011-06-15 12:10:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖

★
fangyongxinxi(金币+2): 2011-06-15 13:54:43
gzqdyouxia(金币+1): 谢谢交流 2011-06-17 09:29:36
主要是因为你6个原子和18个原子的原子对应方式不同，这样可能导致了不同的相变路径。 自然可能得到两个过渡态。
赞一下(2人) 回复此楼
真理是一点点接近的！
3楼2011-06-15 13:22:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)
引用回帖:
Originally posted by gavinliu7390 at 2011-06-15 13:22:12:
主要是因为你6个原子和18个原子的原子对应方式不同，这样可能导致了不同的相变路径。 自然可能得到两个过渡态。

对应方式不同，
你指的是初始和终态，两个构型原子的对应方式不一样吗 ？
即：小体系中 a1-b1,a2-b ,a3-b3,a4-b4,a5-b5,a6-b6;
但在大体系中：a1-b2, a2-b1,(其余部分对应正确)..........

我仔细看了两个体系原子的对应方式，都一样。
所以上面说的那种情况不存在。

另外，大家做NEB选取原胞一般多大 ？
赞一下 回复此楼
4楼2011-06-15 13:34:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★
fangyongxinxi(金币+5): 2011-06-15 13:51:32
gzqdyouxia(金币+2): 谢谢指教 2011-06-17 09:29:52
你确定你两个系统的初始路径相同吗？如果相同的话，说明18个原子的体系相变过程中有对称性破坏，中间态不能用6个原子来表达，那么应该是18个原子的原子势垒要小于6个原子的。应该以18个原子的为准。不过这种情况在NEB体系里面很少见。我很怀疑你是否效验了初始路径，用的image过少或者最后没收敛。
另外提醒你一点，你用的是vtst的code吗？那个code是算分子体系，溶液或者是表面吸附体系的。不能算晶体相变的。除非你的初态和末态的元胞一模一样。他要求元胞必须是个不变的框框。见下面的链接。not yet availble for public use。不过这个yet已经好久了。
http://theory.cm.utexas.edu/vtsttools/neb/
赞一下(1人) 回复此楼
5楼2011-06-15 13:38:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

fangyongxinxi(金币+1): 2011-06-15 13:54:58
引用回帖:
Originally posted by fangyongxinxi at 2011-06-15 13:34:11:
对应方式不同，
你指的是初始和终态，两个构型原子的对应方式不一样吗 ？
即：小体系中 a1-b1,a2-b ,a3-b3,a4-b4,a5-b5,a6-b6;
但在大体系中：a1-b2, a2-b1,(其余部分对应正确)..........

我仔细看了两个 ...

大体系和小体系有一个原子错位的话，路径不一样就很正常了
赞一下 回复此楼
6楼2011-06-15 13:41:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)
引用回帖:
Originally posted by gavinliu7390 at 2011-06-15 13:22:12:
主要是因为你6个原子和18个原子的原子对应方式不同，这样可能导致了不同的相变路径。 自然可能得到两个过渡态。

可能你说的是对的，我仔细看了类似体系大原胞初始和终态的POSCAR文件，发现里面有这样的情况：
初始：
Direct
  0.0000000000000000  0.7500000000000000  0.2500000000000000
  0.5000000000000000  0.7500000000000000  0.2500000000000000
0.0000000000000000  0.2500000000000000  0.7500000000000000
  0.5000000000000000  0.2500000000000000  0.7500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000
终态：
Direct
  0.0427919664457328  0.6726506146310324  0.2500000000000000
  0.5427919664457328  0.6726506146310324  0.2500000000000000
  0.4572080335542673  0.3273493853689677  0.7500000000000000
  0.9572080335542672  0.3273493853689677  0.7500000000000000
  0.1790221628519589  0.2141590376332636  0.2500000000000000
  0.6790221628519589  0.2141590376332636  0.2500000000000000
  0.3209778371480411  0.7858409623667364  0.7500000000000000
  0.8209778371480411  0.7858409623667364  0.7500000000000000

是不是这样的对应方式不行 ？
赞一下 回复此楼
7楼2011-06-15 13:44:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)
小体系的 初始和终态的POSCAR如下：（大体系的POSCAR,为在小体系的构型上在 a 方向扩大2倍，同时将相应原子坐标做变换（不知道里面有错没有））
初始
Direct
0    0.75    0.25
0    0.25    0.75
0.5  0.25    0.25
0.5  0.75    0.75
终态
Direct
0.086    0.673    0.25
0.914    0.327    0.75
0.358  0.214    0.25
0.642  0.786    0.75

大体系（扩大2倍后）
初始：
Direct
  0.0000000000000000  0.7500000000000000  0.2500000000000000
  0.5000000000000000  0.7500000000000000  0.2500000000000000
0.0000000000000000  0.2500000000000000  0.7500000000000000
  0.5000000000000000  0.2500000000000000  0.7500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000
终态：

Direct
  0.0427919664457328  0.6726506146310324  0.2500000000000000
  0.5427919664457328  0.6726506146310324  0.2500000000000000
  0.4572080335542673  0.3273493853689677  0.7500000000000000
  0.9572080335542672  0.3273493853689677  0.7500000000000000
  0.1790221628519589  0.2141590376332636  0.2500000000000000
  0.6790221628519589  0.2141590376332636  0.2500000000000000
  0.3209778371480411  0.7858409623667364  0.7500000000000000
  0.8209778371480411  0.7858409623667364  0.7500000000000000


so， 高手看看，问题出在哪里。
非常感谢。
赞一下 回复此楼
8楼2011-06-15 13:49:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangyongxinxi

新虫 (正式写手)
引用回帖:
Originally posted by dxcharlary at 2011-06-15 13:38:50:
你确定你两个系统的初始路径相同吗？如果相同的话，说明18个原子的体系相变过程中有对称性破坏，中间态不能用6个原子来表达，那么应该是18个原子的原子势垒要小于6个原子的。应该以18个原子的为准。不过这种情况在 ...

谢谢，麻烦再看看后面我的回复。
另外，自己用的是 nebmaker.perl 来生成的NEB中间点。
赞一下 回复此楼
9楼2011-06-15 13:52:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-06-16 09:46:59
fangyongxinxi(金币+2): 2011-06-16 20:31:37
小体系里面终态第二个原子应该写0.914-1，大体系的23行显然反了，而且同小体系0.9那个要减1。
初始路径非常重要，一定要仔细核对。另外可以用nebmovie.pl生成路径，然后用xcrysden观察。
另外我很好奇，你做的是个什么东西？方便透露一起讨论吗？
赞一下(1人) 回复此楼
10楼2011-06-15 20:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fangyongxinxi 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否申诉？ +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:37 by hZiFeudZyoGR
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:35 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:53 by hZiFeudZyoGR
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 17:36 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:14 by hZiFeudZyoGR
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭