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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 自己使用QHA时计算退出
241/3返回列表123下一页
查看: 2183  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:01 by mibUvS8DDCwf
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 23:45 by vZfe6xYu34yj
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:21 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:01 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭