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

24小时热门版块排行榜    

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

dxcharlary

专家顾问 (职业作家)

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

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

» 猜你喜欢
1楼 2010-03-23 16:26:41
已阅   回复此楼   关注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 个回答

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的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
小木虫(金币+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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿哈工大，初试329，求环境科学与工程调剂！ +10 余未辛 2026-04-06 10/500 2026-04-07 15:44 by 上岸快快
[考研] 化工学硕 285求调剂 +16 Wisjxn 2026-04-07 16/800 2026-04-07 15:06 by 一只xx浩
[考研] 085404 293求调剂 +8 勇远库爱314 2026-04-06 9/450 2026-04-07 13:05 by flydream1314
[考研] 297求调剂 +13 GENJIOW 2026-04-07 14/700 2026-04-07 12:25 by 1018329917
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +10 Naiko 2026-04-04 10/500 2026-04-05 09:40 by sam3303
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 292分，材料与化工，申请调剂 +22 程晴之 2026-04-01 26/1300 2026-04-04 22:03 by hemengdong
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 材料295 +13 小英11 2026-04-03 14/700 2026-04-04 09:02 by 来看流星雨10
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 数二英二348求调剂 +4 hxdzj1 2026-04-03 5/250 2026-04-03 21:25 by zhq0425
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
[考博] 申博求助 +3 Reee1Llll 2026-04-01 3/150 2026-04-02 22:29 by 这是一个无聊的�
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 一志愿北京科技，085601总分305求调剂 +9 半生瓜！ 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[考研] 08生物与医药专硕初试346找调剂 +6 dianeeee 2026-04-01 7/350 2026-04-02 08:23 by guoweigw
[考研] 环境工程调剂 +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
信息提示
关闭
请填处理意见
关闭