版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

suspend2014

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 202.5
红花: 2
帖子: 72
在线: 28.5小时
虫号: 2684265
注册: 2013-09-24
专业: 理论和计算化学

[求助] 请问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

回复此楼

» 猜你喜欢

之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有12人回复
博士读完未来一定会好吗已经有27人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复

1楼2015-07-22 16:15:11

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

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

仔細想一下其�?我們兩個的作法都是對的

對於你第二種結構表示法
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
是有問題的

我認為 3.937 應該要乘進去 lattice vector，變成
La2O3
1.000000
  3.409542014699336  -1.968500000000000 0.000000000000000
  0.000000000000000 3.937000000000000 0.000000000000000
  0.000000000000000 0.000000000000000 6.150000000000002
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

wmy8802217 認為應原子位置應該要除以 3.937，變成
La2O3
3.937000
  0.866025403784439  -0.500000000000000 0.000000000000000
  0.000000000000000 1.000000000000000 0.000000000000000
  0.000000000000000 0.000000000000000 1.562103124206249
2 3
Cartesian
0.288675134594813 0.500000000000000 0.385995681991364
0.577350269189626 0.000000000000000 1.176107442214885
0.000000000000000 0.000000000000000 0.000000000000000
0.288675134594813 0.500000000000000 1.007868935737871
0.577350269189626 0.000000000000000 0.554234188468377

�?際上這兩個結構是同一個結構

只是我講的其中有個地方是不正確的，也就是 "Cartesian 會忽略第二行" 這件事，?K不正確
事�?上，第二行還是會被考慮進?� (詳情請看手冊第四十頁)

簡單說一下
假設 lattice constant (POSCAR第二行，也可叫 scale factor) = a
      lattice vector (POSCAR 第三四五行) = V
   atom position (Direct or Cartesian 以下的所有數字) = x
原子真�?位置 R
在 Direct 情況下
R=a * V * x
在 Cartesian 情況下
R = a * x

赞一下

回复此楼

25楼2015-07-23 21:02:08

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

查看全部 28 个回答

wmy8802217

木虫 (正式写手)

应助: 89 (初中生)
金币: 2877
红花: 12
帖子: 350
在线: 461.9小时
虫号: 2935682
注册: 2014-01-16
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

感谢参与，应助指数 +1

下面一个的原子位子全部要除以3.937000才是对的吧？

赞一下

回复此楼

凝聚态理论与计算物理

2楼2015-07-22 19:40:44

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

suspend2014

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 202.5
红花: 2
帖子: 72
在线: 28.5小时
虫号: 2684265
注册: 2013-09-24
专业: 理论和计算化学

引用回帖:

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

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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 的兩個結構之 CONTCAR
就知道~ 這兩種 POSCAR 產生出?淼� CONTCAR 是不同的
主要原因是在於你在POSCAR內用了 Cartesian 這個東西
在VASP計算遇到 POSCAR 內的 Cartesian
VASP 會根據 POSCAR 內給的 lattice vectors ?矶x座標
切記，VASP "只" 會根據 POSCAR 內給的 lattice vectors ?矶x座標
它?K不會去理會你給的 lattice constant (第二行的數字 3.937000)
換句話說就是，你若是想在 POSCAR 內用 Cartesian
那第二行 lattice constant 一定要是 1.0000000
DONE

赞一下

回复此楼

4楼2015-07-22 20:55:48

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

查看全部 28 个回答