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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 自己使用QHA时计算退出
241/3返回列表123下一页
查看: 1998  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

改王之王

金虫 (正式写手)

[求助] 自己使用QHA时计算退出

想计算声子的分波态密度，用QHA，编译后例子能正常运行，后来我用自己计算的fc文件，照着例子写了了edit_me文件，但是一点run_me，过一会就马上退出，出来的输出文件都是空白，这是怎么回事？请问有人知道吗？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-03 19:13:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
zzy870720z(金币+2): 谢谢指教 2011-09-04 09:53:24
改王之王(金币+6): 感谢帮助…………呵呵 2011-09-04 11:00:31
你没有先根据给定的路径生成K点吧~查看一下你的frequency文件有没有正常生成，或者matdyn.in文件后有没有非常长的一串K点坐标。
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
2楼2011-09-04 08:06:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
引用回帖:
2楼: Originally posted by souledge at 2011-09-04 08:06:42:
你没有先根据给定的路径生成K点吧~查看一下你的frequency文件有没有正常生成，或者matdyn.in文件后有没有非常长的一串K点坐标。

昨天照脚本逐条看……………………目前用tetra.x能生成k点…………matdyn里有它产生的k点…………不知道为什么不能用脚本运行……现在用终端自己在逐步的输入计算…………看看到底是哪步的毛病…………
一般能正常计算出声子的色散曲线，产生的fc文件应该没问题吧？
赞一下 回复此楼
3楼2011-09-04 11:05:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★
mazuju028(金币+2): 谢谢交流 2011-09-04 13:50:39
引用回帖:
3楼: Originally posted by 改王之王 at 2011-09-04 11:05:11:
昨天照脚本逐条看……………………目前用tetra.x能生成k点…………matdyn里有它产生的k点…………不知道为什么不能用脚本运行……现在用终端自己在逐步的输入计算…………看看到底是哪步的毛病…………
一般能正 ...

fc文件肯定不会有问题的。
不过，也要看你的结构，QHA目前只支持少数几个结构：
FCC，BCC，SC，HCP，ST，BCT，SO这七个吧~如果你的结构不是这其中的一个，那可能也无法得到结果。
如果frequency文件没有问题并且结构也符合要求，而且也有把ttrinp放在计算目录下，那就仔细看看QHA是在哪一步中止的吧~
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
4楼2011-09-04 12:01:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
好的谢谢指点…………我的是纤锌矿hcp结构，应该支持吧…………呵呵
赞一下 回复此楼
5楼2011-09-04 20:46:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★
zzy870720z(金币+1): 鼓励交流 2011-09-05 08:43:29
引用回帖:
5楼: Originally posted by 改王之王 at 2011-09-04 20:46:21:
好的谢谢指点…………我的是纤锌矿hcp结构，应该支持吧…………呵呵

最多限于100个原子的体系~
其实，8个原子的体系，使用QHA时，就够慢的了~需要接近半个小时，因为QHA不能并行计算……
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
6楼2011-09-04 21:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
算到/opt/espresso-4.2.1/QHA/bin/phonon_dos.x < frequency这步出错了，提示如下信息，可是我phdos.in呀？
natoms==           4
irec====         168
  0.750000000000000     
forrtl: severe (24): end-of-file during read, unit 9, file /opt/espresso-4.2.1/QHA/Examples/ZnO/phdos.in
Image              PC                Routine            Line        Source            
phonon_dos.x       000000000048FCCD  Unknown               Unknown  Unknown
phonon_dos.x       000000000048E7D5  Unknown               Unknown  Unknown
phonon_dos.x       0000000000440F40  Unknown               Unknown  Unknown
phonon_dos.x       000000000040778F  Unknown               Unknown  Unknown
phonon_dos.x       0000000000406FC2  Unknown               Unknown  Unknown
phonon_dos.x       00000000004232EB  Unknown               Unknown  Unknown
phonon_dos.x       00000000004009C1  Unknown               Unknown  Unknown
phonon_dos.x       00000000004002EC  Unknown               Unknown  Unknown
phonon_dos.x       000000000049B060  Unknown               Unknown  Unknown
phonon_dos.x       00000000004001B9  Unknown               Unknown  Unknown

Norm of the difference between old and new effective charges:    0.00000000000000000000
Norm of the difference between old and new force-constants:   0.00054944988447139730
     Message from routine matdyn:
     Z* not found in file ZnO.fc, TO-LO splitting at q=0 will be absent!
     A direction for q was not specified:TO-LO splitting will be absent
赞一下 回复此楼
7楼2011-09-04 22:07:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★
zzy870720z(金币+1): 谢谢提示 2011-09-05 17:33:33
引用回帖:
7楼: Originally posted by 改王之王 at 2011-09-04 22:07:48:
算到/opt/espresso-4.2.1/QHA/bin/phonon_dos.x < frequency这步出错了，提示如下信息，可是我phdos.in呀？
natoms==           4
irec====         168
  0.750000000000000     
forrtl: severe (24): ...

看看前一步的eigenv和partile_DOS有没有生成
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
8楼2011-09-05 08:55:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
引用回帖:
8楼: Originally posted by souledge at 2011-09-05 08:55:10:
看看前一步的eigenv和partile_DOS有没有生成

生成了，但是点击后显示提示“无法显示”…………
赞一下 回复此楼
9楼2011-09-05 18:29:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
引用回帖:
6楼: Originally posted by souledge at 2011-09-04 21:19:44:
最多限于100个原子的体系~
其实，8个原子的体系，使用QHA时，就够慢的了~需要接近半个小时，因为QHA不能并行计算……

不过比起计算声子色散串行要8-17个小时好多了…………呵呵
赞一下 回复此楼
10楼2011-09-05 18:31:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
241/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +6 yexuqing 2026-04-19 7/350 2026-04-23 12:32 by yexuqing
[基金申请] 国自然面上和省基金B类撒花 +18 花田半亩～白 2026-04-21 18/900 2026-04-23 11:31 by 12021227
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[论文投稿] 急需审稿人！！！ +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
关闭
请填处理意见
关闭