24小时热门版块排行榜    

查看: 1268  |  回复: 27
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

suspend2014

新虫 (小有名气)

[求助] 请问VASP计算时,为什么我的POSCAR稍微改了下,结果能差那么多 已有2人参与

我的POSCAR是用不同软件从CIF转的格式,一个是
La2O3
1.0
        3.9370000362         0.0000000000         0.0000000000
       -1.9685000181         3.4095420461         0.0000000000
        0.0000000000         0.0000000000         6.1500000954
    2    3
Cartesian
     0.000000000         2.273028098         1.519664997
     1.968499901         1.136513948         4.630335099
     0.000000000         0.000000000         0.000000000
     0.000000000         2.273028098         3.967980148
     1.968499901         1.136513948         2.182019947
另一种是   
La2O3  
3.937000
  0.866025403784439  -0.500000000000000   0.000000000000000
  0.000000000000000   1.000000000000000   0.000000000000000
  0.000000000000000   0.000000000000000   1.562103124206249
   2   3
Cartesian
  1.136514004899778   1.968500000000000   1.519665000000000
  2.273028009799556   0.000000000000000   4.630335000000001
  0.000000000000000   0.000000000000000   0.000000000000000
  1.136514004899778   1.968500000000000   3.967980000000000
  2.273028009799556   0.000000000000000   2.182020000000000
INCAR是用脚本自动生成的rm WAVECAR
for i in 300 350 400 450 500
do
cat >INCAR<<!
SYSTEM=LA2O3
ENCUT=$i
ISTART=0;ICHARG=2
ISMEAR=0;sigma=0.1
NELM=90
NELMIN=8
ISPIN=2
EDIFF=10E-5
LREAL=.FALSE.
PREC=Accurate
!
echo "ENCUT =$i eV"; time vasp
E=$(grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}')
echo $i $E >>comment
done
可为什么结果却是以下两种呢,上面的POSCAR是从同一个CIF中转换而来的,不过一个是用VESTA,另一个是用CIF2CELL,而且我看也没什么问题
300 -13.822602
350 -13.591244
400 -13.714225
450 -13.539265
500 -13.536550

300 -42.346665
350 -42.039226
400 -41.951873
450 -41.911078
500 -41.899045
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suspend2014

新虫 (小有名气)

引用回帖:
6楼: Originally posted by jpchou at 2015-07-22 21:21:14
嗯~
如果堅持第二行一定要是 3.937
那一定要用 DIRECT 才行

大神,我再问个小白问题,比如我想算LA2O3,最终得到的TOTEN,是不是要减去2*LA的参考组态的能量与3*O的参考组态的能量才是我要求的LA2O3的能量呢?
8楼2015-07-22 21:34:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 28 个回答

wmy8802217

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
下面一个的原子位子全部要除以3.937000才是对的吧?
凝聚态理论与计算物理
2楼2015-07-22 19:40:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suspend2014

新虫 (小有名气)

引用回帖:
2楼: Originally posted by wmy8802217 at 2015-07-22 19:40:44
下面一个的原子位子全部要除以3.937000才是对的吧?

好像是全部除的3.937吧,只是基矢的选择略有差别吧
3楼2015-07-22 20:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
suspend2014: 金币+10, ★★★★★最佳答案 2015-07-22 21:07:25
La2O3  
3.937000
  0.866025403784439  -0.500000000000000   0.000000000000000
  0.000000000000000   1.000000000000000   0.000000000000000
  0.000000000000000   0.000000000000000   1.562103124206249
   2   3
Cartesian
  1.136514004899778   1.968500000000000   1.519665000000000
  2.273028009799556   0.000000000000000   4.630335000000001
  0.000000000000000   0.000000000000000   0.000000000000000
  1.136514004899778   1.968500000000000   3.967980000000000
  2.273028009799556   0.000000000000000   2.182020000000000
?@?? POSCAR ??????}??


???????ε????? 300 eV cutoff energy ?????Y??? CONTCAR
?????~ ?@??N POSCAR ?a???????? CONTCAR ??????
???????????????POSCAR?????? Cartesian ?@???|??
??VASP??????? POSCAR ??? Cartesian
VASP ?????? POSCAR ??o?? lattice vectors ?????x????
?????VASP "?" ?????? POSCAR ??o?? lattice vectors ?????x????
???K??????????o?? lattice constant (????е???? 3.937000)
?Q????f??????????????? POSCAR ???? Cartesian
?????? lattice constant ?????? 1.0000000
DONE
4楼2015-07-22 20:55:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +6 嘉年新程 2026-03-15 6/300 2026-03-21 17:07 by Dream007008
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿武理材料工程348求调剂 +3  ̄^ ̄゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
信息提示
请填处理意见