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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 新手求助~关于VASP位置优化的问题
71/1返回列表
查看: 1632  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
还记得那个凉凉的深秋
1楼 2011-05-04 21:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaods101@163

金虫 (小有名气)

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢提示。 2011-05-05 10:20:41
这得看你的INCAR设置，让大家看这个，好像看不出什么来吧，呵呵
赞一下(1人) 回复此楼
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%
赞一下(1人) 回复此楼
4楼2011-05-05 14:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
youzhizhe: 谢谢提示，呵呵。 2011-05-05 22:58:26
楼上说的很详细了
赞一下(1人) 回复此楼
心之所向无惧无悔
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这个参数的来控制的，查看手册
赞一下(1人) 回复此楼
7楼2011-05-06 16:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 发光的橘子 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[考研] 085601专硕，总分342求调剂，地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0703化学调剂 ，六级已过，有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 288求调剂，一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 307求调剂 +3 冷笙123 2026-03-17 3/150 2026-03-18 09:55 by macy2011
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 334求调剂 +3 志存高远意在机� 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
关闭
请填处理意见
关闭