24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

泽润东方

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 1702.8
散金: 406
帖子: 470
在线: 222.7小时
虫号: 684187
注册: 2008-12-29
性别: GG
专业: 基础物理学

[求助] 关于缩放系数的优化已有1人参与

关于VASP计算我些问题想请教各位。
结构弛豫：空间群是Pnma，三个晶格常数互不相等。
实验方面的结构数据如下：
****************************
a = 0.3846 nm, b = 1.1268 nm, c = 1.1321 nm,
α = 90 °, β = 90 °, γ = 90 ° 。
****************************
以实验结果为基础开始优化，采用GGA-PBE泛函。步骤如下：
1.从0.98-1.8之间选取了二十多个数值作为缩放系数进行了一系列的结构弛豫。
2.选择步骤1中能量最低的结果，作为初始结构，首先固定a方向晶格常数，保持b/a和c/a不变（ISIF=4），进行晶格常数优化，但是保持a，b，c之间的夹角不变；其次再调整a方向的常数，重复上面的步骤。反复进行了十多个测试。
3.最后选择步骤2中得到的结构作为最稳定的结构。如下：
****************************
   direct lattice vectors
3.915449903  0.000000000  0.000000000
0.000000000 11.390399806  0.000000000
0.000000000  0.000000000 11.492100194
****************************
4.在步骤3的基础上重新做了结构弛豫，测试了K点、Ecut。最后确定下来结构和Ecut。
之后的静待自洽、能带以及态密度的计算都是步骤的结构为基础进行的。

目前我比较了一下我的结构跟之前发表的理论值与实验值。误差稍微有些大：
   present       EV-GGA（已发表的理论值） Exp1 Exp2    Exp3
a(angstr) 3.87       3.84                               3.85    3.83    3.84
b(angstr) 11.21    11.22                            11.27 11.23 11.22
c(angstr) 11.96    11.30                            11.32 11.29 11.31
Vol(angstr3) 519.8       -                               490.62 -       -

我的问题如下：
a. 在此基础上所做的能带和态密度可信度高么？
b. 如果不可信，应如何从新开始结构参数的优化？
c. 出现计算误差的原因有哪些？

回复此楼