24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1728  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

alasq

银虫 (初入文坛)

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

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 5 个回答

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

银虫 (初入文坛)

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854求调剂 +12 门路摸摸 2026-04-15 13/650 2026-04-16 00:14 by Art1977
[考研] 297,工科调剂? +3 河南农业大学-能 2026-04-14 3/150 2026-04-15 22:02 by noqvsozv
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 求助调剂,跨调 +19 X十甫寸Y 2026-04-11 20/1000 2026-04-15 21:18 by cuisz
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +26 慕绝cc 2026-04-09 30/1500 2026-04-14 18:50 by 蔡苏阳
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] 一志愿211 0703化学 346分求调剂 +26 土豆er? 2026-04-09 29/1450 2026-04-13 15:15 by 独醉梦孤城
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
信息提示
请填处理意见