24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

weiyongkai

木虫 (职业作家)

应助: 13 (小学生)
金币: 3192
散金: 209
红花: 3
帖子: 3186
在线: 228.1小时
虫号: 1104094
注册: 2010-09-20
专业: 凝聚态物性I:结构、力学和

[求助] QE计算在某一压强下的晶体参数要比相同条件下的VASP小的很多，是怎么回事呢？

输入文件及势函数（官方网站下的）如下：
#! /bin/sh
#
# SCF at dense k-mesh, good enough for electronic DOS
#
cat > al.scf.fit.in <<!
&control
calculation='vc-relax',
restart_mode='from_scratch',
prefix='al',
pseudo_dir = './ps',
outdir='./tmp',
wf_collect =.true.,
etot_conv_thr = 1.0d-6 ,
forc_conv_thr = 1.0d-4 ,
tstress = .true. ,
tprnfor = .true. ,
nstep = 200 ,
/
&system
   ibrav=14,
   a=3.169,b=3.169,c=3.169,
   cosab=0,cosac=0,cosbc=0,
   nat= 8, ntyp= 2,
ecutwfc =100.0,
ecutrho = 800,
occupations='smearing', smearing='methfessel-paxton', degauss=0.05,
la2F=.true.,
/
&electrons
diagonalization='david',
conv_thr =  1.0d-14,
mixing_beta = 0.7,
      mixing_ndim = 10,
diagonalization='david',
/
&IONS
ion_dynamics='bfgs',
ion_positions = 'from_input',
pot_extrapolation = "second_order",
wfc_extrapolation = "second_order",
/
&cell
            cell_dynamics = 'bfgs' ,
            cell_dofree = 'all' ,
            press = 250,
            cell_factor = 5.0,
/
ATOMIC_SPECIES
  H  1.008 H.pz-vbc.UPF
  Al 26.982  Al.pz-vbc.UPF
ATOMIC_POSITIONS {crystal}
  H 0.25 0.00 0.50
  H  -0.25  -0.00 0.50
  H 0.50 0.25  -0.00
  H 0.50  -0.25 0.00
  H 0.00 0.50 0.25
  H 0.00 0.50  -0.25
Al 0.00 0.00 0.00
Al 0.50 0.50 0.50
K_POINTS {automatic}
16 16 16  0 0 0
!
mpirun -np 12 pw.x<al.scf.fit.in > 25.out

B为QE计算的P-V结果，D为VASP计算的P-V结果
Graph1.png

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : ps.rar

2013-05-10 11:28:21, 16.66 K

凌云

1楼 2013-05-10 11:28:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1552
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
weiyongkai: 金币+10, ★有帮助 2013-05-13 11:31:03
weiyongkai: 金币+10, ★★★很有帮助, 谢谢提示，正在一个一个尝试：） 2013-05-14 08:54:05

80%是赝势的问题。
QE里面Al和H的赝势那么多，选择也非常多、