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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 声子谱 虚频 求助
51/1返回列表
查看: 4839  |  回复: 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的回帖

leila3409

捐助贵宾 (正式写手)
引用回帖:
6楼: Originally posted by xhzha at 2013-12-22 13:50:38
EnergyCutoff 的值好大， 有验证需要那么大么，不浪费时间么？没有看到你力的收敛标准， 可以在优化的时候将力的收敛标准设高一些。 如若再出现虚频， 可以考虑将你的虚频的对应的坐标（对应在Eigenvector里)加到原 ...

嗯，很可能是晶胞不够大，我昨天也想到了。我用原胞在各方向扩大5倍，但因为基失不正交，所以由周期性引起的晶胞间相互作用可能还是比较大，我已经将扩展胞正交，现在有250个原子了。
energycutoff 会对别的性质有影响，这里看来对声子影响不大，我接下来还用300就好了。
另外你说的按照波矢移动原子位置，我不懂，可否详细解释一下？谢谢！

扩展胞坐标
0 0 0
0.5 0.5 0.5
虚频位置大概0.2 (2 pi /a)
你说的移动应该如何操作？谢谢！

[ 发自小木虫客户端 ]
赞一下 回复此楼
7楼2013-12-23 08:11:20
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 一志愿武理材料工程348求调剂 +3 ￣^￣゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
关闭
请填处理意见
关闭