版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前主题已经存档。

dxcharlary

专家顾问 (职业作家)

专家经验: +83
1ST强帖: 2
应助: 251 (大学生)
金币: 10398.6
散金: 155
红花: 65
帖子: 3137
在线: 757.1小时
虫号: 936718
注册: 2010-01-04
专业: 凝聚态物性I:结构、力学和
管辖: 第一性原理

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

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

回复此楼

» 猜你喜欢

1楼2010-03-23 16:26:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1551
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gavinliu7390(金币+2):谢谢解答！ 2010-03-24 11:37

有可能是参数选择的问题，我算出来几乎都一样。
在构造原胞时一定要注意边界原子，在pwscf中这类原子的处理相当重要。
最好把输入文件贴出来看看。

赞一下(2人)

回复此楼

2楼2010-03-23 17:31:27

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

专家经验: +83
1ST强帖: 2
应助: 251 (大学生)
金币: 10398.6
散金: 155
红花: 65
帖子: 3137
在线: 757.1小时
虫号: 936718
注册: 2010-01-04
专业: 凝聚态物性I:结构、力学和
管辖: 第一性原理

引用回帖:

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

老虎大王

木虫 (著名写手)

应助: 26 (小学生)
贵宾: 0.17
金币: 4774.1
散金: 8
红花: 42
帖子: 1361
在线: 215.2小时
虫号: 659094
注册: 2008-11-21
专业: 金属结构材料

★ ★ ★
小木虫(金币+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

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1551
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

★ ★ ★
小木虫(金币+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 的主题更新

返回列表