24小时热门版块排行榜    

查看: 3588  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

赵奇一

新虫 (小有名气)

[求助] 第一次运用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编译器)
回复此楼
天使可以飞起来是因为他们把自己看得很轻
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

引用回帖:
10楼: Originally posted by 赵奇一 at 2015-06-14 14:17:37
那这个计算过程对吗?
第一步计算的内容也是scf.第二步计算nscf。第三步就是phonon计算。第四步编辑phonon.第五步画图...

不一定需要计算nscf。编辑Phonon不知道是什么意思。
11楼2015-06-14 17:30:41
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见