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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
31522/32返回列表首页上一页202122232425下一页末页
查看: 15954  |  回复: 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下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 电子信息270求调剂 +12 terminal469 2026-04-07 12/600 2026-04-11 11:19 by zhq0425
[考研] 085400 328分 求调剂 +9 喂你一个大橙子 2026-04-09 13/650 2026-04-11 11:15 by zhq0425
[考研] 326求调剂 +5 Shansyn 2026-04-10 5/250 2026-04-10 22:23 by 猪会飞
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 0856专硕求调剂 希望是a区院校 +21 好好休息好不好 2026-04-09 24/1200 2026-04-10 16:58 by luoyongfeng
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 招收有机化学、化工，药学，食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +11 emokidd 2026-04-08 12/600 2026-04-09 09:24 by wyf236
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 327求调剂 +12 Xxjc1107. 2026-04-06 12/600 2026-04-08 16:46 by luoyongfeng
[考研] 071000生物学，一志愿深圳大学296分，求调剂 +12 TIckLw 2026-04-06 13/650 2026-04-07 20:34 by lijunpoly
[考研] 生物工程求调剂 +13 喜欢还是不甘心 2026-04-05 13/650 2026-04-07 16:55 by Ecowxq666！
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 327考研调剂推荐 +6 呜呜呜呜呢 2026-04-06 6/300 2026-04-06 21:39 by 啵啵啵0119
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭