24小时热门版块排行榜    

查看: 2923  |  回复: 12

xhzha

铁杆木虫 (正式写手)

[求助] 用siesta画声子谱的问题

大家好,我用siesta中的vibra画声子谱,但画出来的声子谱和文献中的声子谱差距很大,文献中声子谱线条很少,而我画出来的声子谱线条很多,看起来跟能带似的。我以为是我K点取得太多,就减少了K点的数目,但画出来的还是和文献有很大差距。下面是我画出来的图以及文献图,还有我的输入文件,请大家帮我看看是什么原因,谢谢了。
文献图

我用比较少k点的图



[ Last edited by xhzha on 2011-11-2 at 17:06 ]
回复此楼

» 收录本帖的淘帖专辑推荐

siesta问题

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


xhzha(金币+2): 2011-11-03 09:04:46
zzy870720z(金币+1): 谢谢提示 2011-11-05 12:23:24
现在能看出的是你的k点取错了
嘟啦啦嘟啦啦嘟
2楼2011-11-02 21:10:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by 贺仪 at 2011-11-02 21:10:38:
现在能看出的是你的k点取错了

谢谢,我一时粗心大意将算管状结构的k点拿过来没有修改,我将k点改为10,10,1重新计算,但我还是想问下,你用siesta算声子谱能得到像文献那样比较少线条的声子谱么,因为我之前算管状结构的声子谱也是很多线条构成的。而频率范围差不多。
3楼2011-11-03 09:04:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

★ ★
zzy870720z(金币+2): 谢谢指教 2011-11-05 12:23:39
线条多少是和原子数目有关的 siesta是用超胞法计算声子的 线条数目应该是明显多于线性响应方法的
4楼2011-11-03 16:28:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
4楼: Originally posted by identation at 2011-11-03 16:28:58:
线条多少是和原子数目有关的 siesta是用超胞法计算声子的 线条数目应该是明显多于线性响应方法的

哦,你是说这是软件的问题,非常感谢,我看了很多文献以及固体物理书上的振动谱线条都很少,所以一直很迷惑。非常感谢! 最近在开始学quantum espresso 看到用DFPT的方法的例子中只有算单个点的声子谱(如Gamma,X),可以算特定方向的声子谱么?希望能给更多指导和建议。谢谢。
5楼2011-11-03 19:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
4楼: Originally posted by identation at 2011-11-03 16:28:58:
线条多少是和原子数目有关的 siesta是用超胞法计算声子的 线条数目应该是明显多于线性响应方法的

抱歉,请点应助回帖,我这样给不了金币·~
6楼2011-11-03 19:22:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

【答案】应助回帖


贺仪(金币+1): 多谢指教! 2011-11-03 21:51:05
xhzha(金币+5): 2011-11-04 12:23:24
引用回帖:
5楼: Originally posted by xhzha at 2011-11-03 19:19:29:
哦,你是说这是软件的问题,非常感谢,我看了很多文献以及固体物理书上的振动谱线条都很少,所以一直很迷惑。非常感谢! 最近在开始学quantum espresso 看到用DFPT的方法的例子中只有算单个点的声子谱(如Gamma, ...

PWSCF都能算,可以计算声子色散曲线 PWSCF的作者就是线性响应理论的发明人 可以试试
7楼2011-11-03 21:37:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

★ ★
xhzha(金币+3): 2011-11-04 12:29:27
zzy870720z(金币+2): 鼓励交流 2011-11-05 12:23:59
引用回帖:
1楼: Originally posted by xhzha at 2011-11-02 17:01:46:
大家好,我用siesta中的vibra画声子谱,但画出来的声子谱和文献中的声子谱差距很大,文献中声子谱线条很少,而我画出来的声子谱线条很多,看起来跟能带似的。我以为是我K点取得太多,就减少了K点的数目,但画出来 ...

你取一个单胞算算,看看是不是声子谱少不少啊。虽然SIESTA计算中采用超胞法,但是最后的谱线条数依赖与你体系中的原子的数目(以及你选择的基组的大小,当然,你所选定的频率范围内的普线条数基本不依赖基组大小,但是位置可能会发生变化)。当你选择的原子数目多的时候,也就说不仅仅包含一个最下原胞的时候,谱线会有折叠的。这个是我的理解。
弘德明志博学笃行
8楼2011-11-04 07:57:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
7楼: Originally posted by identation at 2011-11-03 21:37:04:
PWSCF都能算,可以计算声子色散曲线 PWSCF的作者就是线性响应理论的发明人 可以试试

哦,谢谢,我再试试。
9楼2011-11-04 12:22:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)

引用回帖:
8楼: Originally posted by zhangguangping at 2011-11-04 07:57:24:
你取一个单胞算算,看看是不是声子谱少不少啊。虽然SIESTA计算中采用超胞法,但是最后的谱线条数依赖与你体系中的原子的数目(以及你选择的基组的大小,当然,你所选定的频率范围内的普线条数基本不依赖基组大小, ...

哦,我再试试,之前有算过单胞的声子谱,但总是会出现虚频,用SIESTA.即使力的收敛精度很高。看文献有说frozen phonon(直接法)的方法的原子数最好不要低于100,而且我用超单胞也能消除虚频。所以一直没想过算单个单胞的声子谱了。
10楼2011-11-04 12:28:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xhzha 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 24/1200 2026-08-16 10:16 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
请填处理意见