| 查看: 3505 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
huazhorg铁杆木虫 (正式写手)
|
[交流]
【求助】ph.x计算振动频率的问题 已有5人参与
|
||
| 我想计算一个吸附在硅表面上的水分子的振动频率。所有的原子都计算的话，我试过，太慢了，由于我只想计算水分子的振动频率，就把表面硅原子全部固定。方法是在ph.x的输入文件里加上这样一项：nat_todo=3，并给出水分子三个原子的原子系数。能跑起来，但是算完后，在out文件里却没有任何振动频率的结果。反复试了几次都是如此，还以为是bug,换用最新的QE4.2版本算也是同样的结果。问题出在哪里？请高人指点迷津... |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位大侠，我可以优化结构，但只要一计算频率就出现下面的问题？
已经有4人回复
- 【交流】振动模式与红外光谱、拉曼光谱的原理与计算问题的讨论
已经有14人回复
- 振动模式，声子支与对称性分类--期待热心的虫友
已经有11人回复
- 【求助】请教一下关于Gaussian计算振动频率的问题
已经有7人回复
- 【求助】gaussian频率计算问题
已经有9人回复
- 【求助】MS5.0中Dmol3计算频率问题
已经有21人回复
- 【求助】如何分析ph.x计算振动频率的结果
已经有4人回复
- 【求助】用MS中的castep计算声子振动模频率，具体怎么算
已经有10人回复
huazhorg
铁杆木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 4953.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 ：电子结构
2楼2010-05-31 12:52:42
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
投票:
