24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1732  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 急需调剂 +5 绝不放弃22 2026-04-15 5/250 2026-04-15 23:11 by Tide man
[考研] 22408 312求调剂 +20 门路摸摸 2026-04-14 21/1050 2026-04-15 20:39 by cuisz
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 调剂 +12 月@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 材料085601调剂 +32 何润采123 2026-04-10 34/1700 2026-04-14 08:47 by 木木mumu~
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
信息提示
请填处理意见