24小时热门版块排行榜    

查看: 1847  |  回复: 23
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

改王之王

金虫 (正式写手)

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

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

改王之王

金虫 (正式写手)

好的谢谢指点…………我的是纤锌矿hcp结构,应该支持吧…………呵呵
5楼2011-09-04 20:46:21
已阅   回复此楼   关注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的回帖

改王之王

金虫 (正式写手)

引用回帖:
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的回帖

改王之王

金虫 (正式写手)

上传到这里了…………麻烦帮忙看一下…………自洽文件里有原子种类和坐标…………
http://u.115.com/file/dn1p8hdm
12楼2011-09-05 21:07:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

上传完了…………附上了我写的qha文件,和scf、ph文件…………
14楼2011-09-05 22:14:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

引用回帖:
15楼: Originally posted by souledge at 2011-09-06 07:48:40:
嫩滴matdyn.in文件写错鸟,最后多了个EOF~所以就都悲剧鸟~

原来是这样…………大感谢呀…………我再去试试…………呵呵
16楼2011-09-06 08:16:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

引用回帖:
15楼: Originally posted by souledge at 2011-09-06 07:48:40:
嫩滴matdyn.in文件写错鸟,最后多了个EOF~所以就都悲剧鸟~

去了EOF,还是出现上面的错,后来我自己又写了遍计算AlAs的文件,就没有问题…………
又重新把ZnO的计算文件删了重写还是错,我写一下步骤:
1.从tetrahedra文件夹中找到相应的ttrinp_hcp文件,该文件名为ttrinp,把里面的'X'字母用fc文件里的数值代替
2..../QHA/bin/tetra.x 命令来产生k点,并把产生的k点复制到matdyn.in文件夹中,其中matdyn.in的文件我是照着example写的。
3.执行.../bin/matdyn.xmatdyn.out命令,重新计算频率
4.执行../QHA/bin/Partial_phonon_DOS.x < phdos1.in,,
执行../QHA/bin/phonon_dos.x 然后就提示上面帖子出现的错误,不知道是哪有错,我把例子里的frequency拷到ZnO的文件夹里计算../QHA/bin/phonon_dos.x 另外我的ZnO文件夹是放在QHA的example文件夹下进行计算的……
17楼2011-09-06 20:05:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
信息提示
请填处理意见