| 查看: 2784 | 回复: 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 |
回复此楼» 猜你喜欢
一志愿北京化工大学070300 学硕336求调剂
已经有8人回复
-
一志愿东华大学化学070300，求调剂
已经有6人回复
-
318求调剂
已经有7人回复
-
求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过
已经有6人回复
-
一志愿武理材料工程348求调剂
已经有5人回复
-
求调剂
已经有7人回复
-
275求调剂
已经有6人回复
-
招08考数学
已经有7人回复
-
303求调剂
已经有6人回复
-
306求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教各位大侠，优化石墨烯的晶格常数时，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
- 专业: 理论和计算化学
4楼2014-01-09 12:11:48
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
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
