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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 自己使用QHA时计算退出
241/3返回列表123下一页
查看: 2156  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +6 yuleib84 2026-08-25 7/350 2026-08-25 14:40 by 六两废铜
[基金申请] 什么时候开奖？ +12 CrisMessi 2026-08-18 13/650 2026-08-25 14:40 by ZJTJZ
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 今天基金会出结果吗？20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭