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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » Vasp截断能测试的问题
111/2返回列表12下一页
查看: 3458  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kwanqsh

金虫 (小有名气)

[求助] Vasp截断能测试的问题已有1人参与

想请教一下各位，在做vasp截断能和k点的测试时，是用实验晶胞参数值还是优化后的，如果是优化后的，那么优化时的截断能和k点要怎么选取？还有，有没有可能测试Encut得到的结果是总能随着Encut值增大而增大，最后趋于稳定值？求解答啊~~谢谢~~~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2015-03-10 09:34:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
方便起见可以先测试，用测试的结果做结构优化后再对优化后的结果进行测试校验。

有可能，甚至有可能震荡着收敛。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
2楼2015-03-10 09:50:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kwanqsh

金虫 (小有名气)
引用回帖:
2楼: Originally posted by 卡开发发 at 2015-03-10 09:50:00
方便起见可以先测试，用测试的结果做结构优化后再对优化后的结果进行测试校验。

有可能，甚至有可能震荡着收敛。

谢谢您回答，但是不是说总能越低体系越稳定吗，如果encut增加总能也增加，虽然收敛了 ，结果可信吗
赞一下 回复此楼
3楼2015-03-10 10:54:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
3楼: Originally posted by kwanqsh at 2015-03-10 10:54:00
谢谢您回答，但是不是说总能越低体系越稳定吗，如果encut增加总能也增加，虽然收敛了 ，结果可信吗...

那个是基函数尺寸不充分大造成的计算误差，肯定得按照收敛的结果来考虑，这样才是可靠的。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
4楼2015-03-10 11:09:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ncbbd

铜虫 (初入文坛)
引用回帖:
4楼: Originally posted by 卡开发发 at 2015-03-10 11:09:16
那个是基函数尺寸不充分大造成的计算误差，肯定得按照收敛的结果来考虑，这样才是可靠的。...

你好，我也是初学者，现在正在进行Co体系结构优化，当ENCUT=350ev,结合能是-11.10477701ev;当ENCUT=400ev,结合能是-11.09002913ev;当ENCUT=450ev,结合能是-11.09053261ev.可是候老师的手册上说过Co原子ENCUT的选取最合适的为267ev，最大也不过是它的1.5倍，现在所取得截断能几乎已经超过了1.5倍，可是却出现了上述状况，请您指导一二，谢谢
赞一下 回复此楼
你不坚强懦弱给谁看
5楼2015-07-13 17:09:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
5楼: Originally posted by ncbbd at 2015-07-13 17:09:00
你好，我也是初学者，现在正在进行Co体系结构优化，当ENCUT=350ev,结合能是-11.10477701ev;当ENCUT=400ev,结合能是-11.09002913ev;当ENCUT=450ev,结合能是-11.09053261ev.可是候老师的手册上说过Co原子ENCUT的选取 ...

你可以再仔细看看侯老师的手册，别说PAW，就是截断能需要比较低的超软赝势有些体系也得600以上，还有50eV间隔略有点短。

之前我自己测试下来如Co3O4体系截断能可能得550以上。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
6楼2015-07-14 00:24:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ncbbd

铜虫 (初入文坛)
引用回帖:
6楼: Originally posted by 卡开发发 at 2015-07-14 00:24:06
你可以再仔细看看侯老师的手册，别说PAW，就是截断能需要比较低的超软赝势有些体系也得600以上，还有50eV间隔略有点短。

之前我自己测试下来如Co3O4体系截断能可能得550以上。...

非常感谢您的帮助，可是我现在又有疑问，希望您可以指点，我现在正在对体系的截断能进行测试，文献所用的截断能是350，可是我们的结果不能很好地与文献值一致，因此在不断地修改截断能大小，当选取截断能为450时，平均每个原子的能量与文献值接近，可是发现磁矩不一样，文献是2，我的计算结果是3，而且，从450以后的截断能所计算的结果几乎不变，我现在也是一头雾水，不知道该怎么办，请指点迷津
赞一下 回复此楼
你不坚强懦弱给谁看
7楼2015-07-14 20:33:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by ncbbd at 2015-07-14 20:33:05
非常感谢您的帮助，可是我现在又有疑问，希望您可以指点，我现在正在对体系的截断能进行测试，文献所用的截断能是350，可是我们的结果不能很好地与文献值一致，因此在不断地修改截断能大小，当选取截断能为450时， ...

要不然把文献和INCAR、OUTCAR发来看看，这样便于找出问题。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
8楼2015-07-15 00:17:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ncbbd

铜虫 (初入文坛)

【答案】应助回帖

引用回帖:
8楼: Originally posted by 卡开发发 at 2015-07-15 00:17:49
要不然把文献和INCAR、OUTCAR发来看看，这样便于找出问题。...

您好，我现在把我阅读的文献给您传上去，并且将相关INCAR文件也传过去了，第一个OUTCAR文件时单个钴原子的结果文件，而第二个OUTCAR文件时Co2的结果文件，之所以计算单个原子的结合能就是为了与文献上的平均每个原子的键能进行比对，但是结果不尽人意，并且我实在是不知道自己的问题出现在哪里，希望您能指点迷津
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 2008_Con-1Mn_and_Con-1V.pdf
    • 2015-07-15 16:05:43, 678.22 K
  • 附件 2 : INCAR
    • 2015-07-15 16:06:56, 212 bytes
  • 附件 3 : OUTCAR
    • 2015-07-15 16:07:51, 83.42 K
  • 附件 4 : OUTCAR
    • 2015-07-15 16:13:09, 86.85 K
你不坚强懦弱给谁看
9楼2015-07-15 16:13:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
9楼: Originally posted by ncbbd at 2015-07-15 16:13:25
您好，我现在把我阅读的文献给您传上去，并且将相关INCAR文件也传过去了，第一个OUTCAR文件时单个钴原子的结果文件，而第二个OUTCAR文件时Co2的结果文件，之所以计算单个原子的结合能就是为了与文献上的平均每个原 ...

没找出很具体的原因，你要么把POSCAR、INCAR、POTCAR传来，我做些验算看看
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
10楼2015-07-15 22:52:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kwanqsh 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭