| 查看: 2765 | 回复: 15 | |||
huazhorg铁杆木虫 (正式写手)
|
[交流]
【求助】ph.x计算振动频率的问题已有5人参与
|
| 我想计算一个吸附在硅表面上的水分子的振动频率。所有的原子都计算的话，我试过，太慢了，由于我只想计算水分子的振动频率，就把表面硅原子全部固定。方法是在ph.x的输入文件里加上这样一项：nat_todo=3，并给出水分子三个原子的原子系数。能跑起来，但是算完后，在out文件里却没有任何振动频率的结果。反复试了几次都是如此，还以为是bug,换用最新的QE4.2版本算也是同样的结果。问题出在哪里？请高人指点迷津... |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
各位大侠，我可以优化结构，但只要一计算频率就出现下面的问题？
已经有4人回复- 【交流】振动模式与红外光谱、拉曼光谱的原理与计算问题的讨论
已经有14人回复
- 振动模式，声子支与对称性分类--期待热心的虫友
已经有11人回复
- 【求助】请教一下关于Gaussian计算振动频率的问题
已经有7人回复
- 【求助】gaussian频率计算问题
已经有9人回复
- 【求助】MS5.0中Dmol3计算频率问题
已经有21人回复
- 【求助】如何分析ph.x计算振动频率的结果
已经有4人回复
- 【求助】用MS中的castep计算声子振动模频率，具体怎么算
已经有10人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2010-05-31 12:52:42
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
ph.x 的输入是这样的： vibrations of H2O &inputph nat_todo=3, tr2_ph=4.0d-17, prefix='H2O', outdir='～/tmp', amass(1)=28.d0, amass(2)=16.d0, amass(3)=1.d0, fildyn='H2O.dyn', asr=.true., zue=.true., epsil=.true., trans=.true., / 0.0 0.0 0.0 25 26 27 最下面的三个数，即是水的三个原子在超胞中的系号.（即1-24号原子都是硅原子) 其中输出文件最后的部分内容为： ………………………………………… iter # 10 total cpu time : 54288.4 secs av.it.: 21.0 thresh= 0.256E-07 alpha_mix = 0.700 |ddv_scf|^2 = 0.362E-13 iter # 11 total cpu time : 54509.6 secs av.it.: 21.0 thresh= 0.190E-07 alpha_mix = 0.700 |ddv_scf|^2 = 0.429E-14 End of self-consistent calculation Convergence has been achieved Stopping because representation 1 is not done PHONON : 15h 8m CPU time, 15h11m wall time INITIALIZATION: phq_setup : 1.27s CPU phq_init : 6207.15s CPU phq_init : 6207.15s CPU ……………………………………………………………… 可以看出，循环结束后，按照常规，本应该给出声子结果的，却没有，而是直接结束了，没有给出任何声子的数据。 （补充：asr=.true., zue=.true., epsil=.true., trans=.true. 这几项可以不要，但结果也是一样的：没有声子。） [ Last edited by huazhorg on 2010-5-31 at 14:40 ] |
3楼2010-05-31 14:38:00
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★
huazhorg(金币+2): 2010-05-31 18:40:35
qasd(金币+2):谢谢 2010-05-31 20:42:24
huazhorg(金币+2): 2010-05-31 18:40:35
qasd(金币+2):谢谢 2010-05-31 20:42:24
|
为计算出频率的原因，可能是你的tr2_ph的收敛精度设置太高，而未收敛，适当降低它。 在做ph.x之前进行的pw.x中，使用多个点。比如下面的： K_POINTS {automatic} 6 6 1 1 1 1 将ph.x的输入文件中下面的： vibrations of H2O &inputph tr2_ph=4.0d-14, prefix='H2O', outdir='～/tmp', amass(1)=28.d0, amass(2)=16.d0, amass(3)=1.d0, asr=.true., epsil=.false., trans=.true., lraman=.false., elop=.false., fildyn='H2O.dyn', fildrho='H2O.drho', nat_todo=3, / 0.0 0.0 0.0 25 26 27 上面的调整是针对4.1.2版本的。试试看。 [ Last edited by valenhou001 on 2010-5-31 at 17:00 ] |
4楼2010-05-31 16:57:58
等待的心
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 819.7
- 散金: 620
- 沙发: 1
- 帖子: 372
- 在线: 91.8小时
- 虫号: 947336
- 注册: 2010-01-23
- 性别: GG
- 专业: 物理无机化学
5楼2010-05-31 18:40:17
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2010-05-31 18:44:55
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2010-05-31 18:48:03
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2010-05-31 19:34:52
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★
huazhorg(金币+2): 2010-05-31 22:09:04
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 10:48:04
huazhorg(金币+2): 2010-05-31 22:09:04
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 10:48:04
|
可以的话，把你pw.x和ph.x的输入文件发到valenhou@gmail.com.我再帮你看看。 |
9楼2010-05-31 21:33:43
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4993.8
- 散金: 5
- 红花: 2
- 帖子: 381
- 在线: 363.2小时
- 虫号: 629899
- 注册: 2008-10-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2010-05-31 22:24:50