24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1580  |  回复: 4

alasq

银虫 (初入文坛)

[求助] Pwscf 固定原子(表面体系)算声子谱

计算表面体系的振动频率,因为只关注表面的一些原子的振动频率, 所以算声子的时候将下面几层原子固定了。可是固定原子算出的频率的个数和不固定原子是一样的(即,总原子数的三倍),其中正频率的数目恰好是没固定原子个数的三倍,其余的频率是0,或者是负的。(先只考虑gamma点)。不知道这些正的频率是否就是放开原子的频率呢(相应的振动模式是有其他原子的位移的)。

  没有找到相关的资料,关于固定部分原子算声子的信息很少,所以在此求助,欢迎大家积极讨论指点,非常感谢!
回复此楼

» 猜你喜欢

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

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

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流,, 2013-07-08 22:52:04
QE目前似乎还不支持固定部分原子计算频率。体系不大,就干脆全部计算。对于比较大的体系,比较另类的办法就是直接把你想固定的原子“砍掉”,然后计算剩余原子的频率。
2楼2013-07-06 14:57:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流,, 2013-07-08 22:52:21
alasq: 金币+30, ★★★★★最佳答案, 非常感谢! 2013-07-15 20:50:47
简谐声子是指一个结构内的离子由于原子间力的作用,形成了一个一个的“弹簧”,这些“弹簧”的振动模式。换句话说,就好象要计算带有一定载荷m的给定了虎克系数的弹簧的振动频率,不管有没有“固定”那个载荷,这个“固定”动作根本就没被考虑到计算中,如果说要考虑到,就需要把这个“固定”力作为附加载荷M,最后计算虎克系数为k,载荷为M+m的弹簧振动,得到的结果就根本不再是想要的东西了。
回到LZ的例子中。所谓的固定,是指在relax或者MD的过程,这些离子的质心不变。但是在计算声子过程中,考虑到的是一个模型中原子的所有原子间力相互作用,每个振动模态(共3N个mode)都是所有原子按照一定对称关系受力得到的。所谓的“振动”不过是为了更简单的表述问题,也是为了在frozen phonon方法中实际计算出声子的近似。
最后,至于LZ说的那些实数值的modes,无法说究竟谁贡献多些,但是绝对不是仅仅表面几层“没固定”的离子的。而且,从声子的理论定义来看,那些实数值的modes,到是很可能由中间那些“固定”了的离子贡献的,而纯虚数值的modes,应该是由表面那几层“非固定”离子贡献的。具体的贡献程度,可以通过声子的PDOS进行分析。
思想重于技巧,内涵重于表象
3楼2013-07-07 11:08:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alasq

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by huazhorg at 2013-07-06 14:57:25
QE目前似乎还不支持固定部分原子计算频率。体系不大,就干脆全部计算。对于比较大的体系,比较另类的办法就是直接把你想固定的原子“砍掉”,然后计算剩余原子的频率。

为什么评分评了一个,就没法给您评分和金币了,准备两个回复评一样的分和金币的···非常感谢您的回复。

不过我算的是表面体系,单纯砍掉原子不优化结构的话,体系是不稳定的,应该会有虚频;如果砍掉原子再优化,那底层结构就变了,就不是体相的样子了(表面也不是单纯的切出来优化的,是有重构的表面)。体系还比较大,50多原子。倒是可以考虑搭两个表面对称的结构,层数少一点,我试试看。

谢谢您!
4楼2013-07-15 20:57:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alasq

银虫 (初入文坛)

引用回帖:
3楼: Originally posted by souledge at 2013-07-07 11:08:53
简谐声子是指一个结构内的离子由于原子间力的作用,形成了一个一个的“弹簧”,这些“弹簧”的振动模式。换句话说,就好象要计算带有一定载荷m的给定了虎克系数的弹簧的振动频率,不管有没有“固定”那个载荷,这个 ...

非常感谢您的回帖交流,向您学习了~
固定原子数之后,算出来的频率和不固定是有很大区别的,我知道全放开的话,振动频率大概是取所有原子按3N个方向振动计算的,不清楚固定原子之后具体是怎么计算的。既然程序设计了nat_todo这个参数,应该是有些用的吧~
5楼2013-07-15 21:42:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 alasq 的主题更新
信息提示
请填处理意见