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

24小时热门版块排行榜    

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

souledge

专家顾问 (著名写手)

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-09-05 21:31:28
引用回帖:
9楼: Originally posted by 改王之王 at 2011-09-05 18:29:54:
生成了，但是点击后显示提示“无法显示”…………

那都是二进制文件，肯定是无法显示～实在不成，你把fc文件和原子种类及坐标发上来～
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
11楼2011-09-05 18:45:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

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

souledge

专家顾问 (著名写手)

【答案】应助回帖

★
uuv2010(金币+1): 欢迎参加交流 2011-09-05 21:31:50
引用回帖:
12楼: Originally posted by 改王之王 at 2011-09-05 21:07:01:
上传到这里了…………麻烦帮忙看一下…………自洽文件里有原子种类和坐标…………
http://u.115.com/file/dn1p8hdm

文件拥有者暂时没有权限分享文件。
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
13楼2011-09-05 21:19:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
上传完了…………附上了我写的qha文件，和scf、ph文件…………
赞一下 回复此楼
14楼2011-09-05 22:14:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★
uuv2010(金币+1): 欢迎参加交流 2011-09-06 22:43:07
引用回帖:
14楼: Originally posted by 改王之王 at 2011-09-05 22:14:23:
上传完了…………附上了我写的qha文件，和scf、ph文件…………

嫩滴matdyn.in文件写错鸟，最后多了个EOF~所以就都悲剧鸟~
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
15楼2011-09-06 07:48:40
已阅   回复此楼   关注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的回帖

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的回帖

改王之王

金虫 (正式写手)
引用回帖:
18楼: Originally posted by souledge at 2011-09-07 13:58:33:
感觉你做的很不仔细，关键文件是否存在根本没提及。
我说一下详细步骤：
首先是替换ttrinp中的数据，然后生成K点，把这些K点粘贴到matdyn.in的后面（记得把那个多余的EOF去掉），然后matdyn.x生成frequency，同 ...

我就是这样算的呀…………不过我产生partial_DOS文件是9.1M，我只要把例子里的别的材料的frequency沾过来就能往下算，可能是我从matdyn.x那步就算错了？我去了EOF文件还是不行…………把asr改成simple也不行…………太纠结了…………
赞一下 回复此楼
19楼2011-09-07 18:55:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
不会是程序编译错误吧…………我直接双击的compile…………没有显示编译失败呀…………而且例子也能运行…………
赞一下 回复此楼
20楼2011-09-07 19:38:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
242/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭