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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
31522/32返回列表首页上一页202122232425下一页末页
查看: 15745  |  回复: 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下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 3/150 2026-02-08 03:05 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 3/150 2026-02-08 02:52 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:27 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:25 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:12 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 02:07 by vs90ilomwc
[教师之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 02:05 by vs90ilomwc
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 01:46 by vs90ilomwc
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 01:45 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 01:32 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 01:12 by vs90ilomwc
[教师之家] 有院领导为了换新车，用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人） +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[硕博家园] 博士延得我，科研能力直往上蹿 +7 偏振片 2026-02-02 7/350 2026-02-04 17:36 by 陈氏帝国
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
关闭
请填处理意见
关闭