24小时热门版块排行榜    

查看: 2674  |  回复: 8

laplacetalor

银虫 (小有名气)

[求助] vasp52计算弹性常数 已有1人参与

各位好!我想用VASP 5.2计算材料的弹性常数。但算到一般程序终止,提示“Found 2 degree of freedom, adding 6 lattice distortion”。我的结果是面心立方,不知道哪里出了问题。是不是没有写NFREE啊,请指教。谢谢!
POSCAR 如下:
3.07
1 1 0
0 1 1
1 0 1
1 1
D
0.0 0.0 0.0
0.25 0.25 0.25

INCAR 如下:
ISTART=0
ICHARG=2
ISPIN=2
PREC=Accurate
NELM=80
ALGO=V
LREAL=.F.
EDIFF=1E-6
IBRION=6
NSW=100
ISIF=3
ENCUT=355
EDIFFG=-1E-2
ISMEAR=-5
SIGMA=0.1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
2楼2011-12-07 09:03:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖

★ ★
sunyang1988(金币+2): 谢谢热心帮助 2011-12-07 10:51:19
laplacetalor(金币+20): 谢谢 2011-12-07 11:09:37
引用回帖:
1楼: Originally posted by laplacetalor at 2011-12-07 08:39:35:
各位好!我想用VASP 5.2计算材料的弹性常数。但算到一般程序终止,提示“Found 2 degree of freedom, adding 6 lattice distortion”。我的结果是面心立方,不知道哪里出了问题。是不是没有写NFREE啊,请指教。谢 ...

希望对你有帮助

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 弹性常数.rar
  • 2011-12-07 09:26:23, 1.1 M
3楼2011-12-07 09:26:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laplacetalor

银虫 (小有名气)

引用回帖:
3楼: Originally posted by liliangfang at 2011-12-07 09:26:48:
希望对你有帮助

这个我有看到过,不过挺麻烦的。我看说明书上写着VASP52好像可以直接计算弹性常数了,不知道有谁试过没。
4楼2011-12-07 11:10:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

laplacetalor

银虫 (小有名气)

下面是OSZICAR的结果,能量后面的两行不知道什么意思,程序就终止到这里了。
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.956776482372E+02    0.95678E+02   -0.48520E+03  1316   0.827E+02
RMM:   2     0.144197625391E+02   -0.81258E+02   -0.89878E+02  1316   0.201E+02
......
RMM:  20    -0.555847777831E+01   -0.11712E-05   -0.10240E-06  1781   0.889E-03    0.178E-03
RMM:  21    -0.555847864630E+01   -0.86799E-06   -0.66351E-07  1723   0.853E-03
   1 F= -.55584786E+01 E0= -.55584786E+01  d E =0.000000E+00  mag=     0.0000
Found     2 degrees of freedom:
Adding 6 more for cell shape distortion
5楼2011-12-07 11:18:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiangxiao

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by laplacetalor at 2011-12-07 11:18:12
下面是OSZICAR的结果,能量后面的两行不知道什么意思,程序就终止到这里了。
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.956776482372E+02   ...

请问楼主这个问题解决了吗,我也遇到了,能说一下吗,谢谢
6楼2015-08-25 15:26:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

【答案】应助回帖

你用的是原胞,尝试下用晶胞计算。晶胞保持晶体宏观对称性。关于NFREE的问题,自己测试下吧
7楼2015-08-25 22:52:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

【答案】应助回帖

将结构导入MS,BUILD-〉SYMMETRY->FIND SYMMETRY, 找到对称性后,impose symmetry。 得到晶胞结构,导出cif格式,用vesta打开,然后导出vasp格式。
8楼2015-08-25 22:59:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

9楼2024-04-09 08:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 laplacetalor 的主题更新
信息提示
请填处理意见