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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 自己使用QHA时计算退出
52/1返回列表
查看: 1945  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

改王之王

金虫 (正式写手)

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

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

» 猜你喜欢
1楼 2011-09-03 19:13:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★
贺仪(金币+3): 多谢指导 2011-09-07 20:47:25
改王之王(金币+16): 大感谢………问题终于解决了…………呵呵 2011-09-09 18:12:14
引用回帖:
17楼: Originally posted by 改王之王 at 2011-09-06 20:05:41:
去了EOF，还是出现上面的错，后来我自己又写了遍计算AlAs的文件，就没有问题…………
又重新把ZnO的计算文件删了重写还是错，我写一下步骤：
1.从tetrahedra文件夹中找到相应的ttrinp_hcp文件，该文件名为ttri ...

感觉你做的很不仔细，关键文件是否存在根本没提及。
我说一下详细步骤：
首先是替换ttrinp中的数据，然后生成K点，把这些K点粘贴到matdyn.in的后面（记得把那个多余的EOF去掉），然后matdyn.x生成frequency，同时生成了matdyn.modes。
之后，保证替换后的ttrinp还在当前文件夹中。创建温度输入文件和德拜温度输入文件，然后创建phdos1.in文件。运行Partial_phonon_DOS.x 然后做好phdos.in文件，运行phonon_dos.x 以上所有文件都在同文件夹中，并且文件名和大小写不要轻易更换。最好完全保持与例子中的相同。
CODE:
$ phonon_dos.x natoms==           4
irec====         168
  0.75000000000000000     
Ti1 Ti2
It seems you have imaginary frequences.                     
Hopefully you know what you are doing
nstep====         805
  114.5737  114.5743  147.1030
  122.1805  134.7498  138.2352
  125.8343  131.7529  138.9782
  128.8438  128.8446  139.4173
  125.2058  128.6042  153.5996
  130.8430  131.7825  148.8230
  132.6401  133.3911  146.9807
  132.0160  134.5900  146.1146
  133.6986  137.3024  143.4115
  134.1994  140.3732  140.3739
  118.5830  121.3872  162.3758
E_min=   0.0000000000000000         E_max=   581.27710000000002     
nstep====         805
****** input tetrahedra for BZ-integration ******
   6   3  12
  0.000000  0.000000  0.000000
  0.000000  0.577350  0.000000
-0.333330  0.577350  0.000000
  0.000000  0.000000  0.806500
  0.000000  0.577350  0.806500
-0.333330  0.577350  0.806500
****** input tetrahedra for BZ-integration ******
********************** generate_tetra *********************
  NT0=           3  NTETMX=        1728
    0.0000   0.0000  -0.3333  -0.3333
    0.0000   0.5774   0.5774   0.5774
    0.0000   0.0000   0.0000   0.8065
    1.0000   1.0000   1.0000   1.0000
  volume of tetrahedron =  0.02587
    0.0000   0.0000   0.0000  -0.3333
    0.0000   0.0000   0.5774   0.5774
    0.0000   0.8065   0.8065   0.8065
    1.0000   1.0000   1.0000   1.0000
  volume of tetrahedron =  0.02587
    0.0000   0.0000  -0.3333   0.0000
    0.0000   0.5774   0.5774   0.5774
    0.0000   0.8065   0.8065   0.0000
    1.0000   1.0000   1.0000   1.0000
  volume of tetrahedron =  0.02587
  total volume of BZ is =   0.0776  omg48=  25.7716
****************** end of generate_tetra ******************
before integration:  E_min=   0.0000000000000000         E_max=   581.27710000000002     
Freq,Tot_DOS,DOS==   0.000000       0.000000       0.000000       0.000000       0.000000       0.000000       0.000000   
Freq,Tot_DOS,DOS==  0.7500000      0.2083669E-06  0.8334674E-06  0.2084666E-06  0.4792305E-07  0.1302991E-06  0.3024444E-07
Freq,Tot_DOS,DOS==   1.500000      0.1666935E-05  0.3333870E-05  0.8342646E-06  0.2019754E-06  0.4941277E-06  0.1381614E-06
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
18楼2011-09-07 13:58:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

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的回帖
查看全部 24 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭