版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4755)
>基金申请 (974)
>休闲灌水 (918)
>招聘信息布告栏 (791)
>导师招生 (514)
>虫友互识 (408)
>论文投稿 (379)
>文献求助 (374)
>考博 (175)
>硕博家园 (121)
>学术会议 (117)
>新药研发 (108)
>博后之家 (91)
>有奖起名 (79)
>论文道贺祈福 (66)
>教师之家 (44)

返回列表

identation

金虫 (著名写手)

[交流] 【求助】声子计算的疑问

phonon for STMO
&inputph
tr2_ph=1.0d-12,
prefix='PHONON',
fildvscf='PHdv',
amass(1)=,
amass(2)=,
amass(3)=,
amass(4)=,
outdir='./',
fildyn='PH.dyn',
elph=.false.,
trans=.true.,
ldisp=.true.,
nq1=1, nq2=1, nq3=1
/
这是PH.X的输入文件
   Mode symmetry, C_4h (4/m)  point group:

   omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_g          R
   omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_g*          R
   omega(  1 -  1) =       -32.6  [cm-1] --> A_u          I
   omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_u          I
   omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_u*          I
   omega(  2 -  2) =       69.6  [cm-1] --> A_g          R
   omega(  2 -  2) =       69.6  [cm-1] --> ?
   omega(  2 -  2) =       69.6  [cm-1] --> 0E_g          R
   omega(  2 -  2) =       69.6  [cm-1] --> ?
   omega(  2 -  2) =       69.6  [cm-1] --> 0E_g*          R
   omega(  2 -  2) =       69.6  [cm-1] --> ?
   omega(  2 -  2) =       69.6  [cm-1] --> 0E_u          I
   omega(  2 -  2) =       69.6  [cm-1] --> ?
   omega(  2 -  2) =       69.6  [cm-1] --> 0E_u*          I
   omega(  3 -  4) =       77.7  [cm-1] --> E_u          I
   omega(  3 -  4) =       77.7  [cm-1] --> E_u*          I
   omega(  5 -  6) =       86.4  [cm-1] --> E_g          R
   omega(  5 -  6) =       86.4  [cm-1] --> E_g*          R
   omega(  7 -  8) =       90.9  [cm-1] --> E_u          I
   omega(  7 -  8) =       90.9  [cm-1] --> E_u*          I
   omega(  9 -  9) =       148.3  [cm-1] --> ?
   omega(  9 -  9) =       148.3  [cm-1] --> 0E_g          R
   omega(  9 -  9) =       148.3  [cm-1] --> ?
   omega(  9 -  9) =       148.3  [cm-1] --> 0E_g*          R
   omega(  9 -  9) =       148.3  [cm-1] --> A_u          I
   omega(  9 -  9) =       148.3  [cm-1] --> ?
   omega(  9 -  9) =       148.3  [cm-1] --> 0E_u          I
   omega(  9 -  9) =       148.3  [cm-1] --> ?
   omega(  9 -  9) =       148.3  [cm-1] --> 0E_u*          I
   omega( 10 - 11) =       150.7  [cm-1] --> E_g          R
   omega( 10 - 11) =       150.7  [cm-1] --> E_g*          R
   omega( 12 - 12) =       150.9  [cm-1] --> B_g          R
   omega( 12 - 12) =       150.9  [cm-1] --> ?
   omega( 12 - 12) =       150.9  [cm-1] --> 0E_g          R
   omega( 12 - 12) =       150.9  [cm-1] --> ?
   omega( 12 - 12) =       150.9  [cm-1] --> 0E_g*          R
   omega( 12 - 12) =       150.9  [cm-1] --> ?
   omega( 12 - 12) =       150.9  [cm-1] --> 0E_u          I
   omega( 12 - 12) =       150.9  [cm-1] --> ?
   omega( 12 - 12) =       150.9  [cm-1] --> 0E_u*          I
   omega( 13 - 13) =       157.3  [cm-1] --> ?
   omega( 13 - 13) =       157.3  [cm-1] --> 0E_g          R
   omega( 13 - 13) =       157.3  [cm-1] --> ?
   omega( 13 - 13) =       157.3  [cm-1] --> 0E_g*          R
   omega( 13 - 13) =       157.3  [cm-1] --> A_u          I
   omega( 13 - 13) =       157.3  [cm-1] --> ?
   omega( 13 - 13) =       157.3  [cm-1] --> 0E_u          I
   omega( 13 - 13) =       157.3  [cm-1] --> ?
   omega( 13 - 13) =       157.3  [cm-1] --> 0E_u*          I
   omega( 14 - 15) =       163.9  [cm-1] --> E_u          I
   omega( 14 - 15) =       163.9  [cm-1] --> E_u*          I
   omega( 16 - 17) =       274.8  [cm-1] --> E_u          I
   omega( 16 - 17) =       274.8  [cm-1] --> E_u*          I
   omega( 18 - 18) =       277.1  [cm-1] --> ?
   omega( 18 - 18) =       277.1  [cm-1] --> 0E_g          R
   omega( 18 - 18) =       277.1  [cm-1] --> ?
   omega( 18 - 18) =       277.1  [cm-1] --> 0E_g*          R
   omega( 18 - 18) =       277.1  [cm-1] --> B_u
   omega( 18 - 18) =       277.1  [cm-1] --> ?
   omega( 18 - 18) =       277.1  [cm-1] --> 0E_u          I
   omega( 18 - 18) =       277.1  [cm-1] --> ?
   omega( 18 - 18) =       277.1  [cm-1] --> 0E_u*          I
   omega( 19 - 20) =       391.8  [cm-1] --> E_u          I
   omega( 19 - 20) =       391.8  [cm-1] --> E_u*          I
   omega( 21 - 21) =       396.1  [cm-1] --> ?
   omega( 21 - 21) =       396.1  [cm-1] --> 0E_g          R
   omega( 21 - 21) =       396.1  [cm-1] --> ?
   omega( 21 - 21) =       396.1  [cm-1] --> 0E_g*          R
   omega( 21 - 21) =       396.1  [cm-1] --> A_u          I
   omega( 21 - 21) =       396.1  [cm-1] --> ?
   omega( 21 - 21) =       396.1  [cm-1] --> 0E_u          I
   omega( 21 - 21) =       396.1  [cm-1] --> ?
   omega( 21 - 21) =       396.1  [cm-1] --> 0E_u*          I
   omega( 22 - 22) =       426.4  [cm-1] --> B_g          R
   omega( 22 - 22) =       426.4  [cm-1] --> ?
   omega( 22 - 22) =       426.4  [cm-1] --> 0E_g          R
   omega( 22 - 22) =       426.4  [cm-1] --> ?
   omega( 22 - 22) =       426.4  [cm-1] --> 0E_g*          R
   omega( 22 - 22) =       426.4  [cm-1] --> ?
   omega( 22 - 22) =       426.4  [cm-1] --> 0E_u          I
   omega( 22 - 22) =       426.4  [cm-1] --> ?
   omega( 22 - 22) =       426.4  [cm-1] --> 0E_u*          I
   omega( 23 - 24) =       435.1  [cm-1] --> E_g          R
   omega( 23 - 24) =       435.1  [cm-1] --> E_g*          R
   omega( 25 - 25) =       473.0  [cm-1] --> A_g          R
   omega( 25 - 25) =       473.0  [cm-1] --> ?
   omega( 25 - 25) =       473.0  [cm-1] --> 0E_g          R
   omega( 25 - 25) =       473.0  [cm-1] --> ?
   omega( 25 - 25) =       473.0  [cm-1] --> 0E_g*          R
   omega( 25 - 25) =       473.0  [cm-1] --> ?
   omega( 25 - 25) =       473.0  [cm-1] --> 0E_u          I
   omega( 25 - 25) =       473.0  [cm-1] --> ?
   omega( 25 - 25) =       473.0  [cm-1] --> 0E_u*          I
   omega( 26 - 26) =       475.0  [cm-1] --> B_g          R
   omega( 26 - 26) =       475.0  [cm-1] --> ?
   omega( 26 - 26) =       475.0  [cm-1] --> 0E_g          R
   omega( 26 - 26) =       475.0  [cm-1] --> ?
   omega( 26 - 26) =       475.0  [cm-1] --> 0E_g*          R
   omega( 26 - 26) =       475.0  [cm-1] --> ?
   omega( 26 - 26) =       475.0  [cm-1] --> 0E_u          I
   omega( 26 - 26) =       475.0  [cm-1] --> ?
   omega( 26 - 26) =       475.0  [cm-1] --> 0E_u*          I
   omega( 27 - 27) =       504.2  [cm-1] --> ?
   omega( 27 - 27) =       504.2  [cm-1] --> 0E_g          R
   omega( 27 - 27) =       504.2  [cm-1] --> ?
   omega( 27 - 27) =       504.2  [cm-1] --> 0E_g*          R
   omega( 27 - 27) =       504.2  [cm-1] --> A_u          I
   omega( 27 - 27) =       504.2  [cm-1] --> ?
   omega( 27 - 27) =       504.2  [cm-1] --> 0E_u          I
   omega( 27 - 27) =       504.2  [cm-1] --> ?
   omega( 27 - 27) =       504.2  [cm-1] --> 0E_u*          I
   omega( 28 - 29) =       520.0  [cm-1] --> E_u          I
   omega( 28 - 29) =       520.0  [cm-1] --> E_u*          I
   omega( 30 - 30) =       772.8  [cm-1] --> A_g          R
   omega( 30 - 30) =       772.8  [cm-1] --> ?
   omega( 30 - 30) =       772.8  [cm-1] --> 0E_g          R
   omega( 30 - 30) =       772.8  [cm-1] --> ?
   omega( 30 - 30) =       772.8  [cm-1] --> 0E_g*          R
   omega( 30 - 30) =       772.8  [cm-1] --> ?
   omega( 30 - 30) =       772.8  [cm-1] --> 0E_u          I
   omega( 30 - 30) =       772.8  [cm-1] --> ?
   omega( 30 - 30) =       772.8  [cm-1] --> 0E_u*          I
这是计算得到的声子（0.0，0.0，0.0）点的

可是DOS的峰位与ph.out文件中的上述峰位不一致呢，比如说dos图中的虚频在-50cm-1左右，而omega( 1 - 1) = -32.6 [cm-1] --> ?
omega( 1 - 1) = -32.6 [cm-1] --> 0E_g R
请问这是怎么回事呢？

回复此楼

» 猜你喜欢

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

1楼 2010-12-07 19:21:45

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

暖秋

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 12.5
帖子: 70
在线: 36小时
虫号: 797547

★
identation(金币+1):谢谢参与

引用回帖:

Originally posted by identation at 2010-12-07 19:21:45:
phonon for STMO
&inputph
tr2_ph=1.0d-12,
prefix='PHONON',
fildvscf='PHdv',
amass(1)=,
amass(2)=,
amass(3)=,
amass(4)=,
outdir='./',
fildyn='PH.dyn',
elph=.false.,
trans=.true.,
ldi ...

楼主好，顺便请教你一下：下面这些的意义你能解释一下吗？omega(1-1).......omega(2-2).....这些代表什么？后面的－32.6和最后一列是什么意思？谢谢你。
omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_g          R
   omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_g*          R
   omega(  1 -  1) =       -32.6  [cm-1] --> A_u          I
   omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_u          I
   omega(  1 -  1) =       -32.6  [cm-1] --> ?
   omega(  1 -  1) =       -32.6  [cm-1] --> 0E_u*          I
   omega(  2 -  2) =       69.6  [cm-1] --> A_g          R
   omega(  2 -  2) =       69.6  [cm-1] --> ?
   omega(  2 -  2) =       69.6  [cm-1] --> 0E_g          R
   omega(  2 -  2) =       69.6  [cm-1] --> ?
   omega(  2 -  2) =       69.6  [cm-1] --> 0E_g*          R
   omega(  2 -  2) =       69.6  [cm-1] --> ?
   omega(  2 -  2) =       69.6  [cm-1] --> 0E_u          I
   omega(  2 -  2) =       69.6  [cm-1] --> ?

赞一下(1人)

回复此楼

2楼2010-12-08 10:11:07

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

暖秋

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 12.5
帖子: 70
在线: 36小时
虫号: 797547

引用回帖:

楼主，按你说的意思理解，可能你算出来的动力学矩阵文件不止一个，应该有＊.dyn1~*.dynn个吧？你说的这个最低是－30多点，其余的可能负得更多？我不懂，说错了别见怪，同时麻烦你看看上一个回复，教教我。谢谢。。。

赞一下

回复此楼

3楼2010-12-08 10:14:42

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

huazhorg

铁杆木虫 (正式写手)

应助: 21 (小学生)
金币: 4953.8
帖子: 381
在线: 363.2小时
虫号: 629899

★
identation(金币+1):谢谢参与

为何会有这么多虚频？按说虚频在0附近出现可以理解为数值计算的原因，可这些值已经负几十了，应该不可忽视。楼主确保结构是稳定合理的么？还是什么别的原因？

赞一下(1人)

回复此楼

4楼2010-12-08 15:12:45

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

相关版块跳转我要订阅楼主 identation 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 面上合作单位盖章 +6	ssyjh 2026-08-27	9/450	2026-09-01 20:01 by huagongfeihu
[基金申请] 为什么国自然不能直接公布 +5	bjdxyxy 2026-08-26	5/250	2026-09-01 17:54 by 小伟大博士
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +13	Kittylucky 2026-08-27	14/700	2026-09-01 11:06 by feng6531
[论文投稿] 小白求助投论文要求的highlights应该如何写 5+3	l1963982152 2026-08-29	4/200	2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6	叶九微 2026-08-26	6/300	2026-08-31 23:54 by yudaoqian88
[基金申请] 能否申诉？ +7	echo8914667 2026-08-30	8/400	2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7	349506619 2026-08-28	7/350	2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中，共勉 +12	eulota 2026-08-26	12/600	2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知？ +5	tannykie 2026-08-29	5/250	2026-08-30 21:05 by purplejack
[考博] 找导师 +6	yuanjiabao 2026-08-29	7/350	2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13	wangmingqi 2026-08-28	19/950	2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11	董八千 2026-08-26	11/550	2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4	winsaint 2026-08-27	9/450	2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5	a089 2026-08-26	6/300	2026-08-27 14:05 by jonewore
[基金申请] 为什么国际(地区)合作与交流项目没有放榜？ 10+3	majunge000 2026-08-26	11/550	2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 国合里面能看到了 +7	一怀馨秋 2026-08-26	7/350	2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18	Ldrop2023 2026-08-26	18/900	2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4	yanglien 2026-08-26	5/250	2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6	wittyboy 2026-08-26	14/700	2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3	paperzjh 2026-08-26	3/150	2026-08-26 10:41 by LemmonTr