| 查看: 2545 | 回复: 11 | ||
孙远慧木虫 (小有名气)
|
[求助]
VASP优化几何时，ISIF=3时，如何只改变晶格，而不改变位矢方向.已有2人参与
|
|
各位大神，我在ms建了一个立方晶系，之后用ISIF=2优化几何，发现OUTCAR里面的FORCE on cell下面的YY的值并不是负值，而是正的11点多（图一）。之后试了ISIF=3优化几何发现YY降到负的了。但是伴随产生的另一个新问题就是CONTCAR的坐标发生变化，如图（图二、三）所示。有没有人知道设置什么参数，可以只改变a,b,c的长度，而不改变各个位矢的夹角，就是原来是零的位置优化完之后还是零，但是可以改变晶胞的三个边长。希望各位大神指教！！ 1.jpg  2.jpg  3.jpg |
回复此楼» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求教各位大侠，优化石墨烯的晶格常数时，IBRION与ISIF是如何设置的呢？
已经有4人回复- vasp 优化大概的模型 isif=?合适啊
已经有18人回复
- VASP结构优化ISIF=3，为什么晶格常数不变
已经有12人回复
- 在vasp计算里，有没有好的只优化原子位置和晶胞体积的方法？
已经有8人回复
- 求助成功：VASP进行体积弛豫时用ISIF＝7，如何看弛豫后的结果
已经有7人回复
- 【求助】vasp优化时选择了isif=3 怎么限制晶格的一个维度呢？
已经有16人回复
孙远慧
木虫 (小有名气)
- 应助: 20 (小学生)
- 金币: 1450.3
- 散金: 282
- 红花: 5
- 帖子: 229
- 在线: 224.3小时
- 虫号: 1375586
- 注册: 2011-08-22
- 性别: GG
- 专业: 理论和计算化学
9楼2014-01-10 16:30:02
kjsun
木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 5557.5
- 散金: 100
- 红花: 5
- 帖子: 434
- 在线: 469.8小时
- 虫号: 285070
- 注册: 2006-10-14
- 性别: GG
- 专业: 催化化学
★ ★
sunyang1988: 金币+2, 谢谢交流 2014-01-17 13:52:49
sunyang1988: 金币+2, 谢谢交流 2014-01-17 13:52:49
|
vasp 书上的例子稍微变通一下就能满足你的要求。 对比vasp的例子，仔细看下面的变化 rm WAVECAR for i in 3.7 3.8 3.9 4.0 4.1 do cat >POSCAR <<! fcc: 1 $i 0.0 0.0 0.0 3.7 0.0 0.0 0.0 3.7 1 cartesian 0 0 0 ! echo "a= $i" ; vamp E=`tail -1 OSZICAR` ; echo $i $E >>SUMMARY.fcc done cat SUMMARY.fcc |
11楼2014-01-17 09:16:31
2楼2014-01-09 11:24:36
孙远慧
木虫 (小有名气)
- 应助: 20 (小学生)
- 金币: 1450.3
- 散金: 282
- 红花: 5
- 帖子: 229
- 在线: 224.3小时
- 虫号: 1375586
- 注册: 2011-08-22
- 性别: GG
- 专业: 理论和计算化学
|
现在呢？ 第一张图就是FORCE on cell 里面YY竖行和 in total横行对应的值是11点多，第二张到第三张的是 —————————————————————————— 20.000000000 0.000000000 0.000000000 0.000000000 45.000000000 0.000000000 0.000000000 0.000000000 7.571900000 —————————————————————————— 变成图三 —————————————————————————— 19.1459010205 -0.007219492 0.0378856652 -0.0141628126 45.608278082 0.0187786055 0.01463468126 0.031684302 7.5643364137 —————————————————————————— 我想让图二中20 45 7.5719这三个数据变化，而原来是零的不变化，该添加什么参数啊？  2014-01-09_11-59-19.jpg |
3楼2014-01-09 12:09:31
孙远慧
木虫 (小有名气)
- 应助: 20 (小学生)
- 金币: 1450.3
- 散金: 282
- 红花: 5
- 帖子: 229
- 在线: 224.3小时
- 虫号: 1375586
- 注册: 2011-08-22
- 性别: GG
- 专业: 理论和计算化学
4楼2014-01-09 12:11:48
kjsun
木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 5557.5
- 散金: 100
- 红花: 5
- 帖子: 434
- 在线: 469.8小时
- 虫号: 285070
- 注册: 2006-10-14
- 性别: GG
- 专业: 催化化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
孙远慧: 金币+5, ★★★很有帮助 2014-01-10 09:28:20
sunyang1988: 金币+2, 谢谢交流 2014-01-10 11:17:44
感谢参与，应助指数 +1
孙远慧: 金币+5, ★★★很有帮助 2014-01-10 09:28:20
sunyang1988: 金币+2, 谢谢交流 2014-01-10 11:17:44
|
如果你看过vasp的源程序，就知道如果不改源程序的话，是不可能只改变边长，不改变角度的。原因是ISIF=3相对于ISIF=0，是将晶格参数的9个变量作为和原子坐标一样的参数并列，然后优化。也就是他们和原子坐标等同，这样他就有额外的9个独立的自由度。你要的是限制某些自由度（角度），vasp里面是没有类似的限制的。 解决办法， 1 （推荐）. 如果想只改变边长不改变角度，最好的办法应该是写个小脚本，脚本的写法看vasp手册最后。 2. （不推荐），改变源程序，找到你的用的算法那部分就是IBRION=1or 2or3 对应的脚本部分，加入你的限制条件，重新编译。 |
5楼2014-01-09 12:26:26
孙远慧
木虫 (小有名气)
- 应助: 20 (小学生)
- 金币: 1450.3
- 散金: 282
- 红花: 5
- 帖子: 229
- 在线: 224.3小时
- 虫号: 1375586
- 注册: 2011-08-22
- 性别: GG
- 专业: 理论和计算化学
|
是110页那个吗？脚本能控制在计算过程中只改变20 45 7.5719这三个值的大小，而不改变其余六个的值吗？不需要成比例变化，我只要保证还是原来的立方晶系就行。能不能麻烦你给我发一个满足要求的脚本，或者告诉我在书上多少页？syh5852767@163.com，谢了。。 |
6楼2014-01-10 10:40:25
kjsun
木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 5557.5
- 散金: 100
- 红花: 5
- 帖子: 434
- 在线: 469.8小时
- 虫号: 285070
- 注册: 2006-10-14
- 性别: GG
- 专业: 催化化学
★
sunyang1988: 金币+1, 谢谢交流 2014-01-17 13:52:36
sunyang1988: 金币+1, 谢谢交流 2014-01-17 13:52:36
|
现在有网络，我好久不看书了，不知道多少页。 脚本参照这个例子 http://cms.mpi.univie.ac.at/vasp/guide/node174.html |
7楼2014-01-10 12:09:15
8楼2014-01-10 13:04:40
孙远慧
木虫 (小有名气)
- 应助: 20 (小学生)
- 金币: 1450.3
- 散金: 282
- 红花: 5
- 帖子: 229
- 在线: 224.3小时
- 虫号: 1375586
- 注册: 2011-08-22
- 性别: GG
- 专业: 理论和计算化学
10楼2014-01-16 22:05:48