| 查看: 3104 | 回复: 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 |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有12人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教各位大侠，优化石墨烯的晶格常数时，IBRION与ISIF是如何设置的呢？
已经有4人回复
- vasp 优化大概的模型 isif=?合适啊
已经有18人回复
- VASP结构优化ISIF=3，为什么晶格常数不变
已经有12人回复
- 在vasp计算里，有没有好的只优化原子位置和晶胞体积的方法？
已经有8人回复
- 求助成功：VASP进行体积弛豫时用ISIF＝7，如何看弛豫后的结果
已经有7人回复
- 【求助】vasp优化时选择了isif=3 怎么限制晶格的一个维度呢？
已经有16人回复
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
|
现在呢？ 第一张图就是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
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
|
是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
10楼2014-01-16 22:05:48
40