24小时热门版块排行榜    

查看: 2010  |  回复: 8

youqin5912

金虫 (正式写手)

[求助] 关于finite displacement声子谱计算的问题

本人用MS6.0计算类似TiC结构声子热容,但是用linear response老是出错,改用finite displacement后计算成功,但是finite displacement极其耗费内存,一般都是满载内存运行!
想咨询一下finite displacement比较linear response有何缺陷么?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwangan

金虫 (初入文坛)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢指教,新年快乐 2013-01-01 09:01:05
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 15:14:06
这个和VASP中IBRION参数的选择是一样的。当选用linear response(IBRION=7 or 8) 时,计算力常数矩阵的方法是密度泛函扰动理论(DFPT);而选用finite displacement(IBRION= 5 or 6)时,力常数矩阵计算则采用的是finite difference方法(即笛卡尔坐标下将每个离子移动一个微小位移,算出离子受力后得到Hessian矩阵)
Tomorrowisanotherday
2楼2012-12-30 15:59:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)


sunyang1988: 金币+1, 专家考核, 谢谢交流 2013-01-03 12:11:09
DFPT比finite difference要准确,但更耗内存
3楼2012-12-31 05:39:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youqin5912

金虫 (正式写手)

引用回帖:
3楼: Originally posted by jghe at 2012-12-31 05:39:55
DFPT比finite difference要准确,但更耗内存

不对吧,明明FD更加耗内存啊
4楼2012-12-31 11:47:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)

那可能是Castep的问题,finite difference只是移动一下位置计算原子受力而已怎么可能比linear response更耗内存呢?建议你试试VASP吧
5楼2013-01-01 05:53:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
youqin5912: 金币+5, ★★★很有帮助 2013-01-03 10:10:16
sunyang1988: 金币+1, 专家考核, 谢谢指点 2013-01-03 12:11:39
引用回帖:
5楼: Originally posted by jghe at 2013-01-01 05:53:19
那可能是Castep的问题,finite difference只是移动一下位置计算原子受力而已怎么可能比linear response更耗内存呢?建议你试试VASP吧

因为finite displacement是在超胞中移动原子位置计算原子受力的,而不是在原始晶胞中。
至于linear response,计算过程复杂了些,但是并不需要多耗内存,耗时确实是比较恐怖,尤其当模型提及比较大时。不过linear response的结果的确也好些,尤其对于有LO-TO splitting的体系。
思想重于技巧,内涵重于表象
6楼2013-01-03 00:17:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youqin5912

金虫 (正式写手)

引用回帖:
6楼: Originally posted by souledge at 2013-01-03 00:17:00
因为finite displacement是在超胞中移动原子位置计算原子受力的,而不是在原始晶胞中。
至于linear response,计算过程复杂了些,但是并不需要多耗内存,耗时确实是比较恐怖,尤其当模型提及比较大时。不过linear ...

什么是LO-TO splitting的体系?过渡金属例如TiC是么?TiC用那种方式比较好呢?
7楼2013-01-03 10:11:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwangan

金虫 (初入文坛)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 15:15:02
LO-TO splitting一般是指在强离子型的晶体中,由于库伦力(F~r-1)相比范德瓦尔斯力(F~r-6)是长程作用力,所以这会破坏原有的谐振子模型。解决方法就是引入一个非解析项,这就造成了光学声子在Gamma点的横波和纵波的分离(LO-TO splitting)
Tomorrowisanotherday
8楼2013-01-06 09:34:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoaosky

木虫 (小有名气)

同赞

[ 发自手机版 http://muchong.com/3g ]
9楼2013-01-17 04:21:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 youqin5912 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 01:04 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:03 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见