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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
31522/32返回列表首页上一页202122232425下一页末页
查看: 15584  |  回复: 314
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shelay

木虫 (正式写手)
★
小木虫(金币+0.5): 给个红包，谢谢回帖
请问，我优化完TiN结构后，应力为：
Total      -0.00575    -0.00575    -0.00575     0.00000     0.00000     0.00000
  in kB      -0.11958    -0.11958    -0.11958     0.00000     0.00000     0.00000
  external pressure =       -0.12 kB  Pullay stress =        0.00 kB
（INCAR
SYSTEM = TiN
ENCUT = 450
ISTART = 0; ICHARG = 2
ISMEAR = 1;SIGMA = 0.2
NSW = 60; IBRION = 1
ISIF = 3
POTIM = 0.05
EDIFF = 1E-6; EDIFFG = -1E-5
PREC = Accurate）

然后用优化后的结构计算能量的时候应力变为如下：
  Total       0.06807     0.06807     0.06807     0.00000     0.00000     0.00000
  in kB       1.41457     1.41457     1.41457     0.00000     0.00000     0.00000
  external pressure =        1.41 kB  Pullay stress =        0.00 kB
（INCAR
SYSTEM = TiN
ENCUT = 450
ISTART = 0; ICHARG = 2
ISMEAR = 1
PREC = Accurate）
请问是怎么回事阿？我优化的时候不是把应力都优化的很小了吗？为什么直接用优化完的结构计算能量却得到更大的应力，是不是算法有问题阿？还是其它原因阿，清大家赐教！
赞一下 回复此楼
211楼2012-03-14 21:40:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
引用回帖:
19929158楼: Originally posted by shelay at 2012-03-14 21:40:20:
请问，我优化完TiN结构后，应力为：
Total      -0.00575    -0.00575    -0.00575     0.00000     0.00000     0.00000
  in kB      -0.11958    -0.11958    -0.11958     0.00000     0.00000     0.00000 ...

你可以多优化几次。
赞一下 回复此楼
真理是一点点接近的！
212楼2012-03-16 10:25:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+1, 专家考核, 谢谢指教 2012-03-17 21:39:48
引用回帖:
19929158楼: Originally posted by shelay at 2012-03-14 21:40:20:
请问，我优化完TiN结构后，应力为：
Total      -0.00575    -0.00575    -0.00575     0.00000     0.00000     0.00000
  in kB      -0.11958    -0.11958    -0.11958     0.00000     0.00000     0.00000 ...

试过拟合EOS后的结构了么？K点足够么？使用PWscf的结果来看，对于使用primitive cell的TiN，需要K点在15x15x15或者更高的时候才能得到总能变化在1.0E-4Ry级的收敛：
赞一下 回复此楼
思想重于技巧，内涵重于表象
213楼2012-03-17 20:42:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shelay

木虫 (正式写手)
引用回帖:
19929267楼: Originally posted by souledge at 2012-03-17 20:42:13:
试过拟合EOS后的结构了么？K点足够么？使用PWscf的结果来看，对于使用primitive cell的TiN，需要K点在15x15x15或者更高的时候才能得到总能变化在1.0E-4Ry级的收敛：
c5/7e/941486_1331988127_298.png

k点设置为9x9x9
回复此楼
214楼2012-03-17 21:10:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+2, 专家考核, 谢谢指教 2012-03-17 21:39:38
引用回帖:
19929271楼: Originally posted by shelay at 2012-03-17 21:10:12:
k点设置为9x9x9

除此之外，记得有文献提到过，直接用relax的方法优化结构，可能会有Pulay Stress未能正确处理的问题～
赞一下 回复此楼
思想重于技巧，内涵重于表象
215楼2012-03-17 21:23:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shelay

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
19929274楼: Originally posted by souledge at 2012-03-17 21:23:40:
除此之外，记得有文献提到过，直接用relax的方法优化结构，可能会有Pulay Stress未能正确处理的问题～

那怎么优化啊？用脚本优化的话，对于多个变化（键角，键长）同时变化的时候就不好做，难道vasp对于结构优化就不好办吗？但是计算性能的话没有得到正确的优化结构，怎么会有正确的结论呢，现在vasp发表这么多文章啊？
赞一下 回复此楼
216楼2012-03-18 11:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sswfnydx

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
想问一个问题——http://muchong.com/bbs/viewthread.php?tid=4386854&fpage=1 ，望各位大侠回复！谢谢
赞一下 回复此楼
217楼2012-04-15 11:02:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我想算一下SnO2掺杂，
算SnO2体材料的时候是否应该打破对称性ISYM=0？

如何设置X元素的价态呢·？
因为SnO2中+2价的Sn和+4价的Sn的形成能不一样，
那么Sn的+2和+4价如何设置呢·？
赞一下 回复此楼
218楼2012-05-14 09:39:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whw19850730

铁杆木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
在VASP中可以加磁场计算，需要设置下面一些参数，但是我不知道具体怎么设定怎么去算，说明书里没有，希望高手能给予帮助。

ORBITALMAG=     T  #switch on orbital magnetization
LMAGBLOCH =     F  #use Bloch summation to obtain orbital mag.
LCHIMAG   =     F  #perturbation theory with respect to B field
LGAUGE    =     T  #gauge transformation for 0 moment terms
MAGATOM   =     0  #atom on which magnetic dipole is placed
JATOM     =     0  #atom on which magnetic dipole is evaluated
MAGDIPOL  =    0.0000    0.0000    0.0000 #magnetic dipole
AVECCONST =    0.0000    0.0000    0.0000 #constant vector potential (G=0)
LBFCONST  =     T  #constant B field applied
BCONST    =    0.0000    0.0000   10.0000 #applied magnetic field
赞一下 回复此楼
219楼2012-06-02 01:00:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuchangzzw

铜虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我想请问一下，VASP源代码里面，哪些源代码是关于自洽计算的主要代码？我想研究研究这部分代码，来了解自洽的流程。
赞一下 回复此楼
220楼2012-07-13 10:14:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gavinliu7390 的主题更新
31522/32返回列表首页上一页202122232425下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭