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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】scf的计算出的压强与晶向的选取有关吗？
51/1返回列表
查看: 700  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +4 兮兮夜夜 2026-04-09 4/200 2026-04-09 08:02 by may_新宇
[考研] 0703化学 +29 妮妮ninicgb 2026-04-04 32/1600 2026-04-08 22:50 by 猪会飞
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 专硕310求调剂 +7 捞捞我…. 2026-04-04 8/400 2026-04-08 20:24 by yutian743
[考研] 生物学308分求调剂（一志愿华东师大）做过分子实验 +6 相信必会光芒万� 2026-04-07 7/350 2026-04-08 16:49 by tjzhao
[考研] 一志愿郑州大学085600求调剂 +21 吃的不少 2026-04-05 24/1200 2026-04-08 16:47 by sunhuadong
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 312 化工或制药调剂 +9 小小墨123 2026-04-02 10/500 2026-04-08 10:30 by screening
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 0855求调剂材料 +11 红桃灼灼 2026-04-04 12/600 2026-04-06 10:26 by 蓝云思雨
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 278求调剂 +3 依旧！ 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
[考研] 一志愿北京科技大学材料工程085601，求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 320求调剂 +3 农业工程与信息� 2026-04-03 3/150 2026-04-03 11:40 by 土木硕士招生
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 346求调剂 +5 郑诚乐 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
信息提示
关闭
请填处理意见
关闭