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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 声子谱 虚频 求助
51/1返回列表
查看: 5388  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

leila3409

捐助贵宾 (正式写手)

[求助] 声子谱 虚频 求助 已有4人参与

本人用SIESTA算Nb的声子谱，扩散晶胞为原胞的5*5*5倍，单边长度为3.2985*5=16.5 Ang，共125个原子。EngeryCutoff 选1200 Ry，K点用6*6*6 ，MP展宽用0.3 eV，电子收敛标准为 1.d-5， 原子移动位移为0.04 Bohr。

计算的结果总在GAMMA点附近出现虚频。如下图及与实验值的对比。正确的结果应该是在Gamma附近的声学支出现一个小节，kink，但不应该是虚频。。。

本人所用晶格常数为SIESTA优化后的结果3.2985 Ang,实验值为3.3Ang.

求指教！
声子谱 虚频 求助
计算值


声子谱 虚频 求助-1
实验值
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-20 18:31:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)
★ ★
franch: 金币+2, 谢谢回帖交流 2013-12-23 21:35:04
引用回帖:
7楼: Originally posted by leila3409 at 2013-12-23 08:11:20
嗯，很可能是晶胞不够大，我昨天也想到了。我用原胞在各方向扩大5倍，但因为基失不正交，所以由周期性引起的晶胞间相互作用可能还是比较大，我已经将扩展胞正交，现在有250个原子了。
energycutoff 会对别的性质有 ...

在算声子的输入文件中加上Eigenvectors  True  这个命令，就会得到system.vectors文件，里面会有频率值对应的坐标，找到虚频对应的坐标加到你优化后的原子坐标，然后拿新得到的原子坐标去优化。这样的帖子应该蛮多的。
赞一下 回复此楼
8楼2013-12-23 10:12:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

chenkunkun

铜虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
identation: 金币+2, 谢谢交流 2013-12-21 13:59:47
同求啊 如何消除虚频   请问你这个声子谱是如何算出来的呢  除了siesta 还用其它的辅助软件了吗
赞一下 回复此楼
努力
2楼2013-12-20 22:33:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leila3409

捐助贵宾 (正式写手)
引用回帖:
2楼: Originally posted by chenkunkun at 2013-12-20 22:33:29
同求啊 如何消除虚频   请问你这个声子谱是如何算出来的呢  除了siesta 还用其它的辅助软件了吗

只用SIESTA,没有用其他辅助软件

[ 发自小木虫客户端 ]
赞一下 回复此楼
3楼2013-12-21 06:42:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anionxt

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
identation: 金币+2, 谢谢交流 2013-12-21 13:59:55
Siesta算力常数可靠性比较差,试试Quantum Espresso或者vasp
赞一下 回复此楼
4楼2013-12-21 12:38:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 2rKfVQg2jQww 2026-11-06 15/750 2026-11-09 05:37 by l48hFd8Nwae2
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 05:26 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
关闭
请填处理意见
关闭