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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问VASP计算时，为什么我的POSCAR稍微改了下，结果能差那么多
111/2返回列表12下一页
查看: 1189  |  回复: 27
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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
回复此楼

» 猜你喜欢
1楼2015-07-22 16:15:11
已阅   回复此楼   关注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的回帖

suspend2014

新虫 (小有名气)
引用回帖:
4楼: Originally posted by jpchou at 2015-07-22 20:55:48
La2O3  
3.937000
  0.866025403784439  -0.500000000000000   0.000000000000000
  0.000000000000000   1.000000000000000   0.000000000000000
  0.000000000000000   0.000000000000000   1.56210312420624 ...

高手啊，膜拜~~~
那请问，如果我想用 3.937000，是不是一定要用DIRECT呢
赞一下 回复此楼
5楼2015-07-22 21:06:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suspend2014

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

嗯，谢谢，我去试试看，因为想用那个murn.f来算晶格常数，3.937比较容易改
赞一下 回复此楼
7楼2015-07-22 21:29:49
 已阅   回复此楼   关注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的回帖

suspend2014

新虫 (小有名气)
引用回帖:
10楼: Originally posted by wmy8802217 at 2015-07-23 08:47:11
嗯，第二行如果不是1.0000，下面的原子坐标必须所有的都除以这个数字，Cartesian也一样。...

好像说是Cartesian表示的话第二行只能是1
赞一下 回复此楼
11楼2015-07-23 10:43:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suspend2014

新虫 (小有名气)
引用回帖:
9楼: Originally posted by jpchou at 2015-07-23 03:44:10
比如我想算LA2O3，最终得到的TOTEN，是不是要减去2*LA的参考组态的能量与3*O的参考组态的能量才是我要求的LA2O3的能量呢？

> 我想.....
   你這裡指的 "LA2O3的能量" 是想要問 "LA2O3 的 for ...

不是，就想问的是LA2O3本身的能量，应该用MS算的有好千，VASP却只有40多，我想是不是要减去POTCAR里的EATOM
赞一下 回复此楼
12楼2015-07-23 10:52:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suspend2014

新虫 (小有名气)
引用回帖:
14楼: Originally posted by wmy8802217 at 2015-07-23 10:55:30
可以去看看VASP手册，里边都没这个说法。
看手册40页：

Cubic BN
3.57
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1 1
Selective dynamics
Cartesian
0.00 0.00 0.00 T T F
0.25 0.25 0.25 F F F

很明显 ...

可是我的POSCAR已经全部除过了啊，换成Direct就没错
赞一下 回复此楼
15楼2015-07-23 13:03:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suspend2014

新虫 (小有名气)
引用回帖:
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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suspend2014

新虫 (小有名气)
引用回帖:
27楼: Originally posted by wmy8802217 at 2015-07-23 22:36:53
总的来说一句话，第二行的数字对下面所有的数字都有影响。在计算过程中，都得乘以这个数。

嗯嗯，所以说cif2cell这个脚本没有考虑到这一点，谢谢
赞一下 回复此楼
28楼2015-07-24 15:07:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 suspend2014 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭