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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】scf的计算出的压强与晶向的选取有关吗？
51/1返回列表
查看: 675  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +5 z1z2z3879 2026-03-16 7/350 2026-03-17 11:07 by 学员h26Tkc
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[基金申请] 提交后的基金本子，已让学校撤回了，可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
信息提示
关闭
请填处理意见
关闭