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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【其他】在线答疑:VASP的问题
31522/32返回列表首页上一页202122232425下一页末页
查看: 15824  |  回复: 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下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] (出分啦!!!)26 考研调剂 20w+调剂信息火速更新中！【捡漏】qq 群：1087823172 +3 db58 2026-02-28 79/3950 2026-02-28 17:50 by nbf68
[考研] 311求调剂 +6 南迦720 2026-02-28 6/300 2026-02-28 17:17 by 刘茉莉
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考博] 博士推荐 +4 花儿笑? 2026-02-21 5/250 2026-02-28 16:55 by mumin1990
[考研] 0856材料求调剂 +8 hyf hyf hyf 2026-02-28 9/450 2026-02-28 16:27 by etapple
[考研] 285求调剂 +3 满头大汗的学生 2026-02-28 3/150 2026-02-28 16:22 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考博] 26申博 +3 想申博！ 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 材料类求调剂 +3 wana_kiko 2026-02-28 3/150 2026-02-28 15:03 by lature00
[考研] 295求调剂 +4 19171856320 2026-02-28 4/200 2026-02-28 13:39 by ms629
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[基金申请] 什么是人一生最重要的？ +10 瞬息宇宙 2026-02-21 10/500 2026-02-27 08:46 by tfang
[基金申请] 面上可以超过30页吧？ +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[教师之家] 版面费该交吗 +14 苹果在哪里 2026-02-22 17/850 2026-02-26 11:55 by AGI智创机器人
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭