24小时热门版块排行榜

>论坛更新日志 (2574)
>虫友互识 (138)
>导师招生 (70)
>文献求助 (56)
>休闲灌水 (47)
>硕博家园 (40)
>考研 (26)
>考博 (22)
>博后之家 (18)
>基金申请 (17)
>论文投稿 (17)
>论文道贺祈福 (16)
>找工作 (15)
>教师之家 (14)
>公派出国 (14)
>金属 (4)

返回列表

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

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

» 收录本帖的淘帖专辑推荐

学习

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

一定压力下如何计算弹性常数的问题，请各位高手同志们帮忙解决！谢谢你们啦！已经有5人回复
请问VASP计算，怎么样得到形成焓？已经有12人回复
如何计算物质某个压强下的焓值已经有9人回复

凌云

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
专业: 原子和分子物理

引用回帖:

7楼: Originally posted by weiyongkai at 2013-05-13 11:30:54
我试了其他的赝势，但老提示错误，比如用H.blyp-van_ak.UPF和Al.blyp-n-van_ak.UPF提示：
task# 6
from read-pseudo-nl: error #6
Reading pseudofile (qfunc)
难道从官方网站下下来的势函数不能直接用么？...

肯定可以直接用的。你的设置可能有些问题。另外，带ak字样的是从CPMD中转换过来的，有些小问题。