计算化学公社

标题: pwscf做结构优化帮看看问题所在 [打印本页]

作者
Author: 茶味果冻 时间: 2014-11-5 15:40
标题: pwscf做结构优化帮看看问题所在
计算的为半导体吸附团簇模型，in文件如下：

&CONTROL

calculation= "relax",

restart_mode='from_scratch'

prefix= "cds",

pseudo_dir= "/vol6/home/zhang00/pseudo1",

outdir= "./",

&SYSTEM

ibrav=0,celldm(1)=7.552000,nat=52,ntyp=4,

ecutwfc=40,ecutrho=400,

&ELECTRONS

conv_thr= 1.D-6,

mixing_beta= 0.3,

&IONS

ion_dynamics='bfgs',

bfgs_ndim=1

ATOMIC_SPECIES

O 16.00 O.pbe-n-rrkjus_psl.0.1.UPF

S 32.00 S.pbe-n-rrkjus_psl.0.1.UPF

Ti 48.00 Ti.pbe-sp-van_ak.UPF

Cd 112.40 Cd.pbe-n-van.UPF

ATOMIC_POSITIONS crystal

O 0.000000 0.000000 0.000000

O 0.250000 0.250000 0.264805

。。。此处省略

K_POINTS automatic

2 2 1 0 0 0

CELL_PARAMETERS angstrom

7.5520000 0.0000000 0.0000000

0.0000000 7.5520000 0.0000000

0.0000000 0.0000000 17.9113240

输出文件如下：

(, 下载次数 Times of downloads: 56) (, 下载次数 Times of downloads: 45) (, 下载次数 Times of downloads: 48) 是不是参数值的设置有问题，虚拟内存不足？求帮看看，在此谢过了

作者
Author: helpme 时间: 2014-11-5 16:15
ecutwfc=40,ecutrho=400,还有k，都改小一点看能不能算。

作者
Author: 茶味果冻 时间: 2014-11-5 16:56

helpme 发表于 2014-11-5 16:15
ecutwfc=40,ecutrho=400,还有k，都改小一点看能不能算。

k改成1 截断能25  200 ，结果又报错：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from diropn : error #       10
   error opening cds.wfc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   stopping ...

作者
Author: helpme 时间: 2014-11-5 18:49
这个错已经和前一个不一样了嘛。
你能不能先算一个简单的体系，看你的格式对不对。
现在是打开文件出错了。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)