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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】scf的计算出的压强与晶向的选取有关吗？
51/1返回列表
查看: 782  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dxcharlary

专家顾问 (职业作家)

[交流] 【求助】scf的计算出的压强与晶向的选取有关吗？

一个基态的结构，我把它用material studio成最简元胞之后，带到scf里面去算怎么会有60KBar左右的压强呢？原来的结构没有。
还有为什么vc-relax算出来的压强和scf算出来的压强会有10KBar左右的差距？
回复此楼

» 猜你喜欢
1楼 2010-03-23 16:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+2):谢谢解答！ 2010-03-24 11:38
&cell
cell_dynamics = 'bfgs'
press=94.34
press_conv_thr=0.1

你为什么给出94.34Kbar的压力？你这个是在加压条件下优化的。
赞一下(2人) 回复此楼
4楼2010-03-23 20:38:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

goldenfisher

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+2):谢谢解答！ 2010-03-24 11:37
有可能是参数选择的问题，我算出来几乎都一样。
在构造原胞时一定要注意边界原子，在pwscf中这类原子的处理相当重要。
最好把输入文件贴出来看看。
赞一下(2人) 回复此楼
2楼2010-03-23 17:31:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)
引用回帖:
Originally posted by goldenfisher at 2010-03-23 17:31:27:
有可能是参数选择的问题，我算出来几乎都一样。
在构造原胞时一定要注意边界原子，在pwscf中这类原子的处理相当重要。
最好把输入文件贴出来看看。

貌似我的权限不能上传附件，只能把输入文件贴出来了。先贴那个relax问题的。
ZnO.relax.in:
&control
    calculation='vc-relax'
    restart_mode='from_scratch',
    tstress = .true.
    tprnfor = .true.
    prefix='ZnO',
    pseudo_dir = '/home/dx/pseudo_pwscf/',
    outdir='/home/dx/Image-29/temp'
    nstep = 200
/
&system   
    nosym=.true.,
    ibrav= 0,
    A = 0.529177249,
    nat= 4,
    ntyp= 2,
    ecutwfc = 60.0
/
&electrons
    conv_thr =  1.0d-10
    mixing_beta = 0.7
/
&ions
   ion_dynamics = 'bfgs'
/
&cell
cell_dynamics = 'bfgs'
press=94.34
press_conv_thr=0.1
/
ATOMIC_SPECIES
Zn  1    Zn.pz-van_ak.UPF
O   1     O.pz-van_ak.UPF
ATOMIC_POSITIONS {crystal}
    O    -0.007730710828  -0.008164084083  -0.000128479944
    O     0.531057323447   0.513284399685   0.504978366043
   Zn     0.534941222596   0.531208870365   0.000599504368
   Zn    -0.011852429318  -0.025943383010   0.504255886213
K_POINTS { automatic }
  6 6 4 0 0 0
CELL_PARAMETERS { cubic }
   5.635398814464  -0.028104864643   0.053770344840
   0.731974786133   5.474870714882   0.051201281737
  -0.075971213739  -0.067189850828   7.888057990345

优化后的结果是   
Final enthalpy =    -317.7756731328 Ry
     new unit-cell volume =    243.27897 a.u.^3 (    36.05022 Ang^3 )

CELL_PARAMETERS (alat)
   5.546435726  -0.379100058   0.050116256
   0.384917387   5.547182443   0.049721204
  -0.075313599  -0.065158857   7.868525935

ATOMIC_POSITIONS (crystal)
O        0.011795021   0.002466959   0.002513451
O        0.511423363   0.502713555   0.502339320
Zn       0.511545204   0.502481456   0.002526922
Zn       0.011651818   0.002723832   0.502325584



ZnO.scf.in

&control
    calculation='scf'
    restart_mode='from_scratch',
    tstress = .true.
    tprnfor = .true.
    prefix='ZnO',
    pseudo_dir = '/home/dx/pseudo_pwscf/',
    outdir='/home/dx/ZnO_test/vc-NEB/Image-60/temp'
/
&system   
    nosym=.true.,
    ibrav= 0,
    A = 0.529177249,
    nat= 4,
    ntyp= 2,
    ecutwfc = 60.0
/
&electrons
    conv_thr =  1.0d-6
    mixing_beta = 0.7
/
ATOMIC_SPECIES
Zn  1    Zn.pz-van_ak.UPF
O   1     O.pz-van_ak.UPF
ATOMIC_POSITIONS {crystal}
    O     0.011795021000   0.002466959000   0.002513451000
    O     0.511423363000   0.502713555000   0.502339320000
   Zn     0.511545204000   0.502481456000   0.002526922000
   Zn     0.011651818000   0.002723832000   0.502325584000
K_POINTS { automatic }
  6 6 4 0 0 0
CELL_PARAMETERS { cubic }
   5.546435726000  -0.379100058000   0.050116256000
   0.384917387000   5.547182443000   0.049721204000
  -0.075313599000  -0.065158857000   7.868525935000


结果计算的压强是
    entering subroutine stress ...

          total   stress  (Ry/bohr**3)                   (kbar)     P=  114.74
   0.00078423  -0.00000628  -0.00000012        115.36     -0.92     -0.02
  -0.00000628   0.00078345  -0.00000006         -0.92    115.25     -0.01
  -0.00000012  -0.00000006   0.00077232         -0.02     -0.01    113.61
赞一下 回复此楼
3楼2010-03-23 18:33:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ice_rain(金币+2):xiexie~ 2010-03-24 22:02
relax中用的是 conv_thr =  1.0d-10
scf中是       conv_thr =  1.0d-6

这个会造成一定差别。另外，你确实如“老虎大王 ”说的，目标压强是94.34 Kbar，所以你scf算出来的114.74kbar的压强一点也不奇怪。

另外，vc-relax中，你大多数都用了默认的参数，这个会减慢收敛速度。
赞一下(2人) 回复此楼
5楼2010-03-24 21:10:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 35/1750 2026-08-12 14:38 by sdfapple719
[基金申请] 长年满屏的广告，版主太不责任了。基金也等的急 +3 gltch 2026-08-12 3/150 2026-08-12 13:25 by dcqxinyang
[基金申请] filecode +11 documentary 2026-08-10 12/600 2026-08-12 12:42 by kudofaye
[基金申请] 不应该看fileCode +4 且听虎啸 2026-08-12 4/200 2026-08-12 12:40 by 行无人
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +5 Tide man 2026-08-11 5/250 2026-08-12 09:24 by kangban
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
关闭
请填处理意见
关闭