24小时热门版块排行榜    

查看: 15584  |  回复: 314

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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见