24小时热门版块排行榜    

查看: 1169  |  回复: 4

ljyang320

木虫 (正式写手)

小书虫爱吃糖

[求助] 请教一个castep几何优化的问题

向各位做计算的虫子们请教一个问题:
在“Structural and electrolyte properties of Li4P2S6” which published in Solid State Ionics (Solid State Ionics 284 (2016) 61–70)这篇文章中,作者利用density functional theory the projector augmented wave (PAW) formalism,LDA近似。计算了Li4P2S6的集合结构,得出的结果比实验值低1.2%。但是我在重复此项工作的时候利用CASTEP,GGA-PBE/PBEsol,LDA,采用截断能700eV,k8x8x7,计算结果a~6.25 A,c~10A a大约高出4%,c值高出44%。请大家帮忙分析一下此种情况是何缘故?
     PS:作者两篇已经发表的文献:
1.Structural and electrolyte properties of Li4P2S6 (Solid State Ionics 284 (2016) 61–70);
2.Larry E. Rush“First principles investigation of the structural and electrochemical
properties of Na4P2S6 and Li4P2S6(Solid State Ionics 286 (2016) 45–50)
回复此楼

» 猜你喜欢

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

但使岁月静好,现世安稳!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

【答案】应助回帖

★ ★ ★ ★ ★
ljyang320(蛋蛋小童鞋代发): 金币+5, 谢谢参与,小木虫鼓励积极回答问题者 2017-02-13 10:05:00
蛋蛋小童鞋: 应助指数+1 2017-02-13 10:05:22
1、LDA也有好多种,CASTEP采用的是CA-PZ一种,文献采用的是PW92,虽然原则上不同的LDA区别不大;
2、可能是FFT mesh不够,原则上使用超软赝势的话FFT mesh至少是fine以上结果才比较可靠;
3、可能是赝势和PAW的差异,不排除部分的赝势可移植性比较差,你可以换一下模守恒(Norm-Conserving)赝势试试
不一定挂在论坛,计算问题问题欢迎留言。
2楼2017-02-12 18:12:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljyang320

木虫 (正式写手)

小书虫爱吃糖

引用回帖:
2楼: Originally posted by 卡开发发 at 2017-02-12 18:12:08
1、LDA也有好多种,CASTEP采用的是CA-PZ一种,文献采用的是PW92,虽然原则上不同的LDA区别不大;
2、可能是FFT mesh不够,原则上使用超软赝势的话FFT mesh至少是fine以上结果才比较可靠;
3、可能是赝势和PAW的差 ...

多谢!先试探计算一下。
但使岁月静好,现世安稳!
3楼2017-02-12 22:17:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljyang320

木虫 (正式写手)

小书虫爱吃糖

引用回帖:
2楼: Originally posted by 卡开发发 at 2017-02-12 18:12:08
1、LDA也有好多种,CASTEP采用的是CA-PZ一种,文献采用的是PW92,虽然原则上不同的LDA区别不大;
2、可能是FFT mesh不够,原则上使用超软赝势的话FFT mesh至少是fine以上结果才比较可靠;
3、可能是赝势和PAW的差 ...

还是不行计算结果跟上面差不多。
但使岁月静好,现世安稳!
4楼2017-02-13 09:37:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
4楼: Originally posted by ljyang320 at 2017-02-13 09:37:55
还是不行计算结果跟上面差不多。...

模守恒赝势也验证了?如果是这样,可能得用QE去重复一下他的计算了,ATOMPAW的作者参与了该文献。
不一定挂在论坛,计算问题问题欢迎留言。
5楼2017-02-13 10:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ljyang320 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:58 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:14 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:59 by winsaint
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 3/150 2026-08-29 23:14 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见