24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2573  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

weiyongkai

木虫 (职业作家)

[求助] 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

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

学习

» 猜你喜欢

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

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

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
weiyongkai: 金币+10, 有帮助 2013-05-13 11:31:03
weiyongkai: 金币+10, ★★★很有帮助, 谢谢提示,正在一个一个尝试:) 2013-05-14 08:54:05
80%是赝势的问题。
QE里面Al和H的赝势那么多,选择也非常多、
6楼2013-05-13 10:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

DFPT

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
weiyongkai: 金币+10, 有帮助 2013-05-13 09:15:05
有没有换算VASP与QE的压强单位的不同?
010-8393-9164
2楼2013-05-10 17:15:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
weiyongkai: 金币+10, 有帮助 2013-05-13 09:16:37
LDA赝势会低估晶格常数。
3楼2013-05-12 10:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyongkai

木虫 (职业作家)

引用回帖:
2楼: Originally posted by DFPT at 2013-05-10 17:15:09
有没有换算VASP与QE的压强单位的不同?

换算了啊,QE中的压强单位为Kbar,设置的时候是按1GPa=10Kbar设置的啊,这也小的太多了啊?
凌云
4楼2013-05-13 09:14:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 考研调剂 +4 watb 2026-03-26 4/200 2026-03-30 10:46 by 唐沐儿
[考研] 一志愿中南大学化学0703总分337求调剂 +6 niko- 2026-03-27 6/300 2026-03-30 10:25 by herarysara
[考研] 0703化学 +13 妮妮ninicgb 2026-03-27 13/650 2026-03-30 09:18 by 王亮_大连医科大
[考研] 化学工程085602 305分求调剂 +23 RichLi_ 2026-03-25 23/1150 2026-03-30 09:01 by Awei099
[考研] 281求调剂 +4 亚克西good 2026-03-26 6/300 2026-03-29 21:06 by 啊李999
[考研] 一志愿华东师范大学有机化学专业,初试351分,复试被刷求调剂! +5 真名有冰 2026-03-29 6/300 2026-03-29 20:53 by 唐沐儿
[考研] 调剂310 +12 温柔的晚安 2026-03-25 13/650 2026-03-29 20:01 by 无际的草原
[考研] 08工科求调剂286 +3 tgs_001 2026-03-28 3/150 2026-03-29 18:29 by 1018329917
[考研] 材料学硕333求调剂 +11 北道巷 2026-03-24 11/550 2026-03-29 08:31 by fmesaito
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 346求调剂 一志愿070303有机化学 +3 萝卜炖青菜 2026-03-28 3/150 2026-03-28 14:11 by 唐沐儿
[考研] 304求调剂 +6 曼殊2266 2026-03-27 6/300 2026-03-28 14:10 by 唐沐儿
[考研] 285求调剂 +4 AZMK 2026-03-27 7/350 2026-03-27 20:59 by AZMK
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 一志愿郑大085600,310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:03 by cmz0325
[考研] 生物学学硕求调剂 +7 小羊睡着了? 2026-03-23 10/500 2026-03-25 02:24 by 清风拂扬。 m
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
信息提示
请填处理意见