请问VASP计算时，为什么我的POSCAR稍微改了下，结果能差那么多 - 第一性原理 - Vasp&MedeA

凝聚态理论与计算物理

21楼2015-07-23 16:33:18

jpchou

专家顾问 (著名写手)

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

引用回帖:

20楼: Originally posted by wmy8802217 at 2015-07-23 16:29:25
你可以算算，拿手算，这种体系刚好两种坐标是等价的。或者你用VESTA之类的看看，体系完全一样。
0.25*【（0 0.5 0.5）+（0.5 0 0.5）+（0.5 0.5 0）】=（0.25 0.25 0.25）...

你的意思是說

Cubic BN
3.57
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
B N
1 1
Selective Dynamics
Direct
0.00 0.00 0.00  T T T
0.25 0.25 0.25  T T T

這樣的結構 "剛好" 與下面這個等價

Cubic BN
3.57
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
B N
1 1
Selective Dynamics
Cartesian
0.00 0.00 0.00  T T T
0.25 0.25 0.25  T T T

你是認真的嬤? 你真的有試過嬤?
確定那一行 Cartesian or Direct 在這系統剛好不會影響結果嬤?

我把這兩種?G進 vesta 畫
怎麼看都是不同結構...........................

22楼2015-07-23 17:55:30

wmy8802217

木虫 (正式写手)

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

引用回帖:

22楼: Originally posted by jpchou at 2015-07-23 17:55:30
你的意思是說

Cubic BN
3.57
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
B N
1 1
Selective Dynamics
Direct
0.00 0.00 0.00 T T T
0.25 0.25 0.25 T T T

這樣的結構 "剛好" 與下面這個等價 ...

那我只能说，我怎么看都是一样的结果。如果你真的懂这些数字到底什么意思，就知道我刚才写的0.25*【（0 0.5 0.5）+（0.5 0 0.5）+（0.5 0.5 0）】=（0.25 0.25 0.25）是什么意思

凝聚态理论与计算物理

23楼2015-07-23 18:35:12

wmy8802217

木虫 (正式写手)

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

引用回帖:

你也可以试试用VESTA把两个结构都转换一下成相同的结构，你看看结果一不一样

凝聚态理论与计算物理

24楼2015-07-23 18:43:09

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

suspend2014

新虫 (小有名气)

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

引用回帖:

25楼: Originally posted by jpchou at 2015-07-23 21:02:08
仔細想一下其�?我們兩個的作法都是對的

對於你第二種結構表示法
La2O3
3.937000
  0.866025403784439  -0.500000000000000 0.000000000000000
  0.000000000000000 1.000000000000000 0.00000000000 ...

哦，原来如此，彻底明白了，谢谢两位大神的讨论

26楼2015-07-23 21:09:02