24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 1668  |  回复: 6
本帖产生 1 个 1ST强帖 ,点击这里进行查看

发光的橘子

铜虫 (小有名气)

[求助] 新手求助~关于VASP位置优化的问题

刚接触VASP,还有很多不懂的~
我试图计算锐钛矿的LDA+U计算,优化之后看CONTCAR文件
   1.00000000000000     
     3.9426097708347903    0.0000000000000000    0.0000000000000000
     0.0000000000000000    3.9426097708347903    0.0000000000000000
     0.0000000000000000    0.0000000000000000    9.7285115413474781
   4   8
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000   T   T   T
  0.0000000000000000  0.5000000000000000  0.2500000000000000   T   T   T
  0.5000000000000000  0.5000000000000000  0.5000000000000000   T   T   T
  0.5000000000000000  0.0000000000000000  0.7500000000000000   T   T   T
  0.0000000000000000  0.5000000000000000  0.0419691282295330   T   T   T
  0.0000000000000000  0.0000000000000000  0.2080308717704670   T   T   T
  0.5000000000000000  0.5000000000000000  0.2919691282295330   T   T   T
  0.0000000000000000  0.5000000000000000  0.4580308717704670   T   T   T
  0.5000000000000000  0.0000000000000000  0.5419691282295330   T   T   T
  0.5000000000000000  0.5000000000000000  0.7080308717704670   T   T   T
  0.0000000000000000  0.0000000000000000  0.7919691282295330   T   T   T
  0.5000000000000000  0.0000000000000000  0.9580308717704670   T   T   T
      我觉得有些奇怪,优化之后的所有原子的前两个坐标相对于对应的晶格的大小来说都没变(优化之前也是0.5,0.25,0.75之类的),尤其是前四个原子,三个坐标都没变,只有后面8个原子的最后一个坐标改变了。
   另外,我用的是LDA+U,它的晶格优化的结果与别人计算的结果相差很大,体积相差10%了,在不用LDA+U的情况下,也有上面的原子的相对位置几乎没动的现象。

求助~~~~~~~

[ Last edited by 发光的橘子 on 2011-5-4 at 21:40 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

还记得那个凉凉的深秋
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaods101@163

金虫 (小有名气)

【答案】应助回帖


youzhizhe(金币+1): 谢谢提示。 2011-05-05 10:20:41
这得看你的INCAR设置,让大家看这个,好像看不出什么来吧,呵呵
2楼2011-05-05 09:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

发光的橘子

铜虫 (小有名气)

我的INCAR如下


  NPAR = 6
  LPLANE = .TRUE.
  LSCALU = .FALSE.
  NSIM = 4

  ISTART = 0
  PREC = High
  GGA = PE

  NSW = 200
  IBRION = 2
  ISIF = 3

  LORBIT = 11
  ENCUT = 550.00 eV
  EDIFF = 1.0E-5
  EDIFFG = -0.01
  NELM = 60
  NELMIN = 2

  ISMEAR = -5
  SIGMA = 0.08

  LMAXMIX-tag = 4
  LDAU = .TRUE.
  LDAUL  = 2 -1
  LDAUU = 6.0 0.0
  LDAUJ = 0.0 0.0
U取6.0是随便试试的,还试过其他好多数值,基本上是U越大晶格体积偏离越大
还记得那个凉凉的深秋
3楼2011-05-05 12:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

【答案】应助回帖


发光的橘子(金币+1): 总算弄清楚了,谢谢了 2011-05-05 22:49:49
youzhizhe(金币+1, 1ST强帖+1): 谢谢交流。 2011-05-05 22:58:01
原子坐标不动是很正常的事情,如果一个原子是在高对称位置上,那么弛豫的时候位置一般都不会变化,正如你前面4个原子都是高对称位置。
其次是,LDA+U如果不是计算磁学性质的话,一般晶格优化是不需要的,影响也不应该大到10%
4楼2011-05-05 14:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)

youzhizhe: 谢谢提示,呵呵。 2011-05-05 22:58:26
楼上说的很详细了
心之所向无惧无悔
5楼2011-05-05 14:10:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

发光的橘子

铜虫 (小有名气)

引用回帖:
Originally posted by bingmou at 2011-05-05 14:00:03:
原子坐标不动是很正常的事情,如果一个原子是在高对称位置上,那么弛豫的时候位置一般都不会变化,正如你前面4个原子都是高对称位置。
其次是,LDA+U如果不是计算磁学性质的话,一般晶格优化是不需要的,影响也不 ...

也就是说VASP在优化晶格的时候是会考虑到对称性的了?

看起来你是对的,我在原来的POSCAR每个坐标后面都加上一个很小的误差,那么就不会出现这种情况了,LDA+U计算的体积也几乎没有什么变化了。
说明我之前的计算中,VASP以牺牲计算准确性来换取对称性的保证。

THANKS~
还记得那个凉凉的深秋
6楼2011-05-05 22:49:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

【答案】应助回帖


youzhizhe(金币+1): 谢谢交流。 2011-05-27 10:38:28
引用回帖:
Originally posted by 发光的橘子 at 2011-05-05 22:49:16:
也就是说VASP在优化晶格的时候是会考虑到对称性的了?

看起来你是对的,我在原来的POSCAR每个坐标后面都加上一个很小的误差,那么就不会出现这种情况了,LDA+U计算的体积也几乎没有什么变化了。
说明我之前 ...

这刻通过调节ISYM这个参数的来控制的,查看手册
7楼2011-05-06 16:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 发光的橘子 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑州大学085600求调剂 +21 吃的不少 2026-04-05 24/1200 2026-04-08 16:47 by sunhuadong
[考研] 324求调剂 +16 想上学求调 2026-04-03 16/800 2026-04-08 16:07 by luoyongfeng
[考研] 调剂 +7 月@163.com 2026-04-08 7/350 2026-04-08 14:32 by 418490947
[考研] 0703化学调剂325分 +13 15771691647 2026-04-04 14/700 2026-04-08 10:28 by screening
[考研] 293分求调剂,外语为俄语 +5 加一一九 2026-04-07 7/350 2026-04-07 23:04 by JourneyLucky
[考研] 一志愿郑州大学材料与化工085600,求调剂 +34 吃的不少 2026-04-02 34/1700 2026-04-07 20:01 by lrll?l
[考研] 材料调剂 +12 一样YWY 2026-04-02 14/700 2026-04-07 12:52 by Sammy2
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 0703化学 +9 goldtt 2026-04-02 11/550 2026-04-06 10:35 by 无际的草原
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 302求调剂一志愿华中师范大学 +8 小江小江江江 2026-04-02 8/400 2026-04-04 19:50 by 蓝云思雨
[考研] 11408 一志愿西电,277分求调剂 +4 zhouzhen654 2026-04-03 4/200 2026-04-04 18:10 by 猪会飞
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 11408,284分,二战真诚求调剂 +4 12.27 2026-04-02 4/200 2026-04-03 14:14 by dxiaoxin
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
信息提示
请填处理意见