| 查看: 2627 | 回复: 16 | ||
| 【悬赏金币】回答本帖问题，作者13s154019将赠送您 5 个金币 | ||
13s154019新虫 (初入文坛)
|
[求助]
超晶格结构优化时INCAR参数设置问题已有2人参与
|
|
|
学第一性原理照着文献学习计算呢，不过在结构优化时就碰到了问题。 文献上的BaTiO3/SrTiO3平衡晶格参数a=3.95, c=7.9不过我优化的结果是a=3.99, b=8。挺愁人的 下面是文献，还有我的INCAR设置，各位大神看看我是哪点出现问题了？？  1111.png  2222222.png  333.png |
回复此楼» 猜你喜欢
院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度
已经有10人回复
-
哎
已经有3人回复
-
25申博记录贴
已经有4人回复
-
紫外光谱
已经有3人回复
-
事业单位还是大学好？
已经有16人回复
-
考研英语带学
已经有5人回复
-
298求调剂
已经有8人回复
-
家乡二本高校/沿海传统私企，如何抉择？
已经有13人回复
-
下雨了
已经有19人回复
-
小孩7岁，上一年级，不自信，学习较差
已经有30人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 二维体系加电场的ISIF取2还是3？？？
已经有0人回复
- 富锂材料的一些零零碎碎
已经有118人回复
- 表面体系的结构优化
已经有12人回复
- 新手求教pw算声子譜提示缺少.dyn1文件
已经有9人回复
- CASTEP 咱们板块最近不少讨论CASTEP模块的，我开辟一帖子，供大伙讨论
已经有4人回复
- 2010年度中俄（NSFC－RFBR）合作项目批准通知
已经有0人回复
- 【原创】第一原理计算的一些心得
已经有160人回复
- 【转帖】Materials-Studio 论坛问答全集（精选众多论坛讨论贴）
已经有122人回复
- 【分享】经典计算书籍《冶金和材料计算物理化学》（乔芝郁）
已经有49人回复
- 【转帖】VASP 初学者必读
已经有280人回复
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:29:45
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:29:45
|
另外说实话你的结果和人家的挺接近的，只是两个原子沿着z方向平移了一点 发自小木虫IOS客户端 |
5楼2015-10-08 17:49:44
LJ云水禅心
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 3008.7
- 散金: 50
- 红花: 6
- 帖子: 412
- 在线: 393.9小时
- 虫号: 2341134
- 注册: 2013-03-12
- 性别: MM
- 专业: 半导体材料
10楼2015-10-09 09:10:37
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:29:51
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:29:51
|
晶格常数相差1%以内是正常的。 如果要优化晶格常数需要ISIF=3。 另外POTIM设的太小了 发自小木虫IOS客户端 |
2楼2015-10-07 23:06:05
13s154019
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7
- 帖子: 31
- 在线: 56.5小时
- 虫号: 3155035
- 注册: 2014-04-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2015-10-08 11:41:50
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:29:57
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:29:57
|
首先你的ISIF=2只会优化原子位置，不会优化晶格常数，如果想优化晶格常数那么就需要ISIF=3。 自己算的和别人算不一样很正常，看你想要误差在百分之多少。GGA的晶格常数误差一般在1%。 有些因素会影响你的结果，比如pseudopotential的选择，K点密度，cutoff大小，还有范德华力修正 发自小木虫IOS客户端 |
4楼2015-10-08 17:48:05
13s154019
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7
- 帖子: 31
- 在线: 56.5小时
- 虫号: 3155035
- 注册: 2014-04-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
差之毫厘，谬之千里啊。按照我这个优化结果接着计算声子谱，冻结声子，再得到新的相，差别就大了。声子谱根本不一样。 其中k点，还有cutoff都好选。赝势PAW和交换互联函数GGA下的PBE不太清楚。你再看看这个计算过程文件nohup.out。这个怎么还显示是LDA呢？？？不会是这点错了吧 running on 1 nodes distr: one band on 1 nodes, 1 groups vasp.5.2.12 26May11 complex POSCAR found : 5 types and 10 ions LDA part: xc-table for Pade appr. of Perdew POSCAR, INCAR and KPOINTS ok, starting setup FFT: planning ...( 4 ) WAVECAR not read entering main loop N E dE d eps ncg rms rms(c) |
6楼2015-10-08 18:05:43
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
13s154019: 金币+10, ★有帮助 2015-10-08 18:37:08
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:30:04
13s154019: 金币+10, ★有帮助 2015-10-08 18:37:08
liliangfang: 金币+1, 谢谢交流 2015-10-10 10:30:04
|
PBE就是这么显示的没有错。你的声子谱有虚频吗。别光看原子坐标，可能会有平移的，你的晶格常数还有键长和别人差多少 发自小木虫IOS客户端 |
7楼2015-10-08 18:26:51
13s154019
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7
- 帖子: 31
- 在线: 56.5小时
- 虫号: 3155035
- 注册: 2014-04-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2015-10-08 22:46:56
13s154019
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7
- 帖子: 31
- 在线: 56.5小时
- 虫号: 3155035
- 注册: 2014-04-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-10-08 22:53:53
10