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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 第一次运用QE计算貌似出现问题，求解答
51/1返回列表
查看: 4268  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

赵奇一

新虫 (小有名气)

[求助] 第一次运用QE计算貌似出现问题，求解答 已有3人参与

昨天刚把QE安装，完成，想模仿论坛里一个朋友的计算过程进行一次运算，求解决。首先建立一个纯文本文档，内容直接将这位朋友编辑好的粘贴过来，内容如下：
&CONTROL
                 calculation = 'scf' ,
                      outdir = './' ,
                  pseudo_dir = './' ,
                      prefix = 'silicon' ,
/
&SYSTEM
                       ibrav = 2,
                   celldm(1) = 10.2,
                         nat = 2,
                        ntyp = 1,
                     ecutwfc = 12.0 ,
                        nbnd = 8,
            exxdiv_treatment = 'gygi-baldereschi' ,
/
&ELECTRONS
/
ATOMIC_SPECIES
   Si   28.08600  Si.vbc.UPF
ATOMIC_POSITIONS
   Si      0.000000000    0.000000000    0.000000000
   Si      0.250000000    0.250000000    0.250000000
K_POINTS
2
   0.250000000    0.250000000    0.750000000      3.000000000   
   0.250000000    0.250000000    0.250000000      1.000000000


执行的命令为：$espresso-5.3.2/bin/pw.x<si.scf.in>si.scf.out


输入终端显示的结果为：
zhaoqiyi@zhaoqiyi-HP-G42-Notebook-PC:~$ espresso-5.1.2/bin/pw.x<si.scf.in>si.scf.out
STOP 1

同时产生的一个文档si.scf.out，这个应该是输出文件吧，里边显示的内容为：
Program PWSCF v.5.1.2 starts on 13Jun2015 at 19:54:28

     This program is part of the open-source Quantum ESPRESSO suite
     for quantum simulation of materials; please cite
         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
          URL http://www.quantum-espresso.org",
     in publications or presentations arising from this work. More details at
     http://www.quantum-espresso.org/quote

     Serial version
     Waiting for input...
     Reading input from standard input
     Message from routine read_cards :
     DEPRECATED: no units specified in ATOMIC_POSITIONS card
     Message from routine read_cards :
     ATOMIC_POSITIONS: units set to alat

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine readpp (2):
     file ./Si.vbc.UPF not found
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...

请问这是我软件没安装好？还是说我用的编译器不合适（因为论坛里有朋友说我不应该用gfortran而应该用intel编译器）
回复此楼

» 猜你喜欢
天使可以飞起来是因为他们把自己看得很轻
1楼 2015-06-13 20:16:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵奇一

新虫 (小有名气)
引用回帖:
9楼: Originally posted by KalaShayminS at 2015-06-14 10:42:22
不对，band是电子计算，只需要pw.x。声子计算需要ph.x，而且输入文件也完全不一样。这种时候需要固体物理基础知识外加慢慢看文档来做。...

那这个计算过程对吗？
第一步计算的内容也是scf.第二步计算nscf。第三步就是phonon计算。第四步编辑phonon.第五步画图
赞一下 回复此楼
天使可以飞起来是因为他们把自己看得很轻
10楼2015-06-14 14:17:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
赵奇一: 金币+20, ★★★很有帮助 2015-06-14 08:29:47
需要在pseudo_dir指定的目录里面（这里是当前目录）放入赝势文件（例如Si.vbc.UPF）
自带的pseudo文件夹或者官网上都有赝势下载
赞一下 回复此楼
2楼2015-06-13 20:18:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuguiqiu

禁虫 (文坛精英)
感谢参与，应助指数 +1
本帖内容被屏蔽
4楼2015-06-13 21:16:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 赵奇一 at 2015-06-13 21:14:47
不好意思，我的过程都是按一个朋友的帖子照猫画虎操作的。弱弱的问一句，在5.1.2的文件夹里有pseud文件夹，并且找到一个Si.pz-vbc.UPF文件（这个应该就是您说的Si.vbc.UPF吧...），该怎么样将这个文件放在pseudo_d ...

输入文件里写了pseudo_dir = './'   就是说是当前目录
输入文件里的UPF文件和放到目录里的文件名要对应。

直接复制UPF文件到当前目录就好了。
要用Si.pz-vbc.UPF的话就把输入文件里面的Si.vbc.UPF改成Si.pz-vbc.UPF
赞一下 回复此楼
5楼2015-06-13 21:22:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:30 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:16 by cITaGg3p5edr
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:06 by cITaGg3p5edr
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-18 01:40 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 9/450 2026-09-17 23:05 by cuPZTDXS3VOd
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +7 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 22:53 by cuPZTDXS3VOd
[基金申请] 国社科系统bug了，是不是要放榜了？ +5 kynobel 2026-09-16 6/300 2026-09-17 21:37 by tosogo
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 5BDX0d0WFp7t 2026-09-15 4/200 2026-09-17 16:29 by 0ngsZVSCTS7z
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-17 14:17 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 40ms4Wlo4umh 2026-09-14 6/300 2026-09-17 13:53 by ObVzqQQh1rrE
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 13:39 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-17 13:15 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 12:29 by ObVzqQQh1rrE
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
信息提示
关闭
请填处理意见
关闭