应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » QE计算在某一压强下的晶体参数要比相同条件下的VASP小的很多，是怎么回事呢？
101/1返回列表
查看: 3063  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

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

学习

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
凌云
1楼 2013-05-10 11:28:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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赝势会低估晶格常数。
赞一下(1人) 回复此楼
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的回帖

weiyongkai

木虫 (职业作家)
引用回帖:
3楼: Originally posted by huazhorg at 2013-05-12 10:36:53
LDA赝势会低估晶格常数。

是的，可是从结果看来低估的有些大，需要转换势函数么？可以用abinit势函数转化过来用么？
赞一下 回复此楼
凌云
5楼2013-05-13 09:16:29
已阅   回复此楼   关注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的回帖

weiyongkai

木虫 (职业作家)
引用回帖:
6楼: Originally posted by goldenfisher at 2013-05-13 10:18:42
80%是赝势的问题。
QE里面Al和H的赝势那么多，选择也非常多、

我试了其他的赝势，但老提示错误，比如用H.blyp-van_ak.UPF和Al.blyp-n-van_ak.UPF提示：
task#  6
from  read-pseudo-nl: error #6
Reading pseudofile (qfunc)
难道从官方网站下下来的势函数不能直接用么？
赞一下 回复此楼
凌云
7楼2013-05-13 11:30:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
weiyongkai: 金币+10, ★★★很有帮助, 谢谢帮助！ 2013-05-14 08:52:51
引用回帖:
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)
难道从官方网站下下来的势函数不能直接用么？...

你试试PBE赝势，比如Al.pbe-n-van.UPF和H.pbe-van_ak.UPF
赞一下 回复此楼
8楼2013-05-14 08:36:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
引用回帖:
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中转换过来的，有些小问题。
赞一下 回复此楼
9楼2013-05-14 16:07:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyongkai

木虫 (职业作家)
★
ben_ladeng: 金币+1, 欢迎分享您的经验 2013-05-15 07:37:29
引用回帖:
9楼: Originally posted by goldenfisher at 2013-05-14 16:07:54
肯定可以直接用的。你的设置可能有些问题。另外，带ak字样的是从CPMD中转换过来的，有些小问题。...

我也怀疑是我的设置的问题，另外我觉得势函数的选择应该要根据所要计算的要求进行选择，谢谢
赞一下 回复此楼
凌云
10楼2013-05-14 16:54:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 weiyongkai 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-11 08:12 by 1Nb0IVlmV6v0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 07:08 by Z7qfCsdIiD1T
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 07:03 by Z7qfCsdIiD1T
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-11 06:57 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 by Z7qfCsdIiD1T
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 04:12 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:54 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:32 by Z7qfCsdIiD1T
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:12 by Z7qfCsdIiD1T
[考博] 申请环境，微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭