24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1726  |  回复: 7
本帖产生 1 个 1ST强帖 ,点击这里进行查看

从心而觅

新虫 (小有名气)

[求助] 结合使用castep和vasp遇到了问题。

请教大家三个问题,非常感谢!
1:我在vasp对k-point进行了测试,那么这个k-point可以拿到castep中用吗,感觉k点设置和软件没关系。
2:我在计算单个原子能量时,得到H的能量是-13.51ev,与-13.6ev比较接近,说明计算相对可靠。得到C的能量是-148.49ev,而Si的是-103.87ev,这个我无法理解,按说C要大于Si。请问结果合理吗?
3:我想用castep计算声子,而计算声子需要用模守恒贗势,而我优化构型用的是超软,那么如果我用模守恒来重新测试e-cut,进而计算声子,这种做法可行吗?
请求帮助,感激不尽!
回复此楼
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
从心而觅(漫天飘雪代发): 金币+10, 谢谢交流 2016-08-27 16:32:56
蛋蛋小童鞋: 1ST强帖+1, 应助指数+1, 回答详细并且解决了楼主的问题。 2016-08-28 22:02:56
请教大家三个问题,非常感谢!
那我就直接在这里回答一下我自己的看法了
1:我在vasp对k-point进行了测试,那么这个k-point可以拿到castep中用吗,感觉k点设置和软件没关系。
k点的设置理论上似乎和软件没关系,但是其实k点的收敛取值却与每个代码的精度和详细计算有关,所以依然建议对CASTEP的K点进行收敛测试。但是可以从VASP的K点设置作为开端,然后向后(更多)方试几个,如果能证明VASP的K点也能满足CASTEP的收敛那最好,如果不能满足,还要继续找下去
2:我在计算单个原子能量时,得到H的能量是-13.51ev,与-13.6ev比较接近,说明计算相对可靠。得到C的能量是-148.49ev,而Si的是-103.87ev,这个我无法理解,按说C要大于Si。请问结果合理吗?
首先,不清楚LZ说的单个原子能量是什么意思,是建立个诸如1*1*1 nm的框框,中心放个原子进行计算么?如果是的话,这个值没有什么单独拿来看的意义,换句话说,就算你得到的H的能量是-135.1eV也没有关系,这个与赝势的设置、零点参考点的设置有关,各个软件和赝势都有不同。所以我们在实际计算能量时,都要抠掉赝势自带的能量。其实,能量的0点本身在热力学里就没有定义,所以焓值总是定个标准后的相对值(否则谁能告诉我内能究竟有多少呢?)
3:我想用castep计算声子,而计算声子需要用模守恒贗势,而我优化构型用的是超软,那么如果我用模守恒来重新测试e-cut,进而计算声子,这种做法可行吗?
虽然软件是同一个,但是赝势本身的质量还是有差异的,所以不同的赝势给出的优化的结果可能也有小幅的出入。为了免得由于一个更节约时间的优化没有做到位,而导致最后漫长的声子计算出个虚频重新算,还是建议LZ尽量对NC赝势也重新测试e-cut和优化结构(甚至K点),然后计算声子。另外,CASTEP支持线性响应和冻声子两种方法计算声子,前者要求必须NC赝势,后者如果计算资源的内存够大,也是可以用的,不必线性响应慢的,而且不限制赝势类型。
请求帮助,感激不尽!
思想重于技巧,内涵重于表象
3楼2016-08-26 23:40:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

从心而觅

新虫 (小有名气)

再次考虑了一下问题2,觉得结果是可信的。因为最后得到的能量除去了贗势的能量,不知道大家怎么看?
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
2楼2016-08-23 12:47:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)

引用回帖:
3楼: Originally posted by souledge at 2016-08-26 23:40:31
请教大家三个问题,非常感谢!
那我就直接在这里回答一下我自己的看法了
1:我在vasp对k-point进行了测试,那么这个k-point可以拿到castep中用吗,感觉k点设置和软件没关系。
k点的设置理论上似乎和软件没关系, ...

非常详细,太好啦,谢谢老师!
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
4楼2016-08-27 15:43:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)


蛋蛋小童鞋: 金币+1, 谢谢参与,小木虫鼓励积极回帖 2016-08-28 22:04:54
很奇怪!为什么不一直用一个软件用相同的赝势和cutoof ,k点来计算性质而且又可信呢?楼主给人的感觉是,比如你搞对象时三心二意似的,不觉得吗?

[ 发自手机版 http://muchong.com/3g ]
做好自己
5楼2016-08-28 13:03:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)

引用回帖:
5楼: Originally posted by tsglss at 2016-08-28 13:03:13
很奇怪!为什么不一直用一个软件用相同的赝势和cutoof ,k点来计算性质而且又可信呢?楼主给人的感觉是,比如你搞对象时三心二意似的,不觉得吗?

你说的有道理,我这样主要有两个原因。
1:计算条件的限制,比如说我用castep跑不动大的k点,只好用vasp。
2:不同软件在各个方面的优势不一样,比如用casp计算声子比较麻烦,所以为了方便,想结合着使用。
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
6楼2016-08-28 13:13:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shttp125

木虫 (小有名气)

dfasdf

楼主 分享下问题2的单质能量怎么算,我在算熵边上就在这里卡这了,谢谢

发自小木虫Android客户端
dfad
7楼2016-09-01 08:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)

引用回帖:
7楼: Originally posted by shttp125 at 2016-09-01 08:08:56
楼主 分享下问题2的单质能量怎么算,我在算熵边上就在这里卡这了,谢谢

我算的是单个原子,不是单质http://muchong.com/bbs/viewthread.php?tid=10565279
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
8楼2016-09-01 08:53:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 从心而觅 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[考博] 找博士生导师 +3 小代想上岸 2026-05-15 4/200 2026-05-17 09:20 by YuY66
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:55 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:10 by ue3ir18jc3
[考博] 现在不知道怎么办,感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +9 一切都是空工 2026-05-12 20/1000 2026-05-16 19:52 by Equinoxhua
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-16 18:35 by j6b2pdz07o
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[文学芳草园] 裁员滚滚,退居二线 +4 J_wei 2026-05-10 4/200 2026-05-16 10:52 by zh10246
[考博] 2026博士还有哪些学校有名额 +5 小王求读研 2026-05-15 6/300 2026-05-16 10:44 by a4742549
[教师之家] 上海大学实验技术岗位非升即走 +5 嘻嘻哈哈乐呵呵 2026-05-15 5/250 2026-05-16 00:17 by caiyun
[考博] 西南大学考核制博士 +4 lijunjie84 2026-05-11 7/350 2026-05-15 23:20 by 同仁堂教主
信息提示
请填处理意见