24小时热门版块排行榜    

查看: 3599  |  回复: 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的回帖

赵奇一

新虫 (小有名气)

引用回帖:
7楼: Originally posted by KalaShayminS at 2015-06-13 23:58:56
对的。...

谢谢您一直耐心帮我,我对这个一窍不通,看手册总怕理解不到位,所以想从照猫画虎开始理解。
我想再请教一下计算步骤的事。对应研究对象确定之后,这位朋友是计算能带的步骤是这样的:1.scf.2nscf.3.band计算4.编辑bands.5.画图
想请问您一下,是不是所有的计算前两个步骤都是一样的?比方说我打算计算声子谱。第一步计算的内容也是scf.第二步计算nscf。第三步就是phonon计算。第四步编辑phonon.第五步画图
我这么对应过来不知道正确不正确?
还有,这位朋友在前三步采用的功能都是pw.x 我想知道一下是不是计算方面都采用pw.x或者说计算的内容和功能有对应关系?
天使可以飞起来是因为他们把自己看得很轻
8楼2015-06-14 08:53:31
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见