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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】scf的计算出的压强与晶向的选取有关吗？
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 695  |  回复: 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):给个红包，谢谢回帖交流
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的回帖

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的回帖
相关版块跳转 我要订阅楼主 dxcharlary 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 324求调剂 +16 想上学求调 2026-04-03 16/800 2026-04-08 16:07 by luoyongfeng
[考研] 304求调剂 +10 素年祭语 2026-04-06 17/850 2026-04-08 09:05 by 蓝云思雨
[考研] 专硕085403，291分，有两篇专利，一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 一志愿北京2，材料与化工308求调剂 +16 熊二想上岸 2026-04-04 17/850 2026-04-07 14:58 by BruceLiu320
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 285求调剂 +8 AZMK 2026-04-04 11/550 2026-04-06 13:56 by BruceLiu320
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 296求调剂 +3 汪！？！ 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 315求调剂 +5 ＆123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[考研] 求调剂 +4 晟功? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[考研] 求调剂 +3 ffyyu 2026-04-02 3/150 2026-04-04 19:03 by 蓝云思雨
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 一志愿北京科技大学材料工程085601，求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 295求调剂 +3 尚偌呀 2026-04-03 4/200 2026-04-03 21:23 by zhq0425
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 一志愿北京交通大学材料工程总分358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by 百灵童888
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 化学工程专硕324分，一志愿中国矿业大学求调剂 +7 耿耿1314 2026-04-01 7/350 2026-04-02 07:40 by 尚水阁主
信息提示
关闭
请填处理意见
关闭