24小时热门版块排行榜    

查看: 1837  |  回复: 23

改王之王

金虫 (正式写手)

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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点坐标。
思想重于技巧,内涵重于表象
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是在哪一步中止的吧~
思想重于技巧,内涵重于表象
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不能并行计算……
思想重于技巧,内涵重于表象
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有没有生成
思想重于技巧,内涵重于表象
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的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
信息提示
请填处理意见