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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于缩放系数的优化
51/1返回列表
查看: 1867  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

泽润东方

金虫 (正式写手)

[求助] 关于缩放系数的优化 已有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. 出现计算误差的原因有哪些 ？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
希望就在不远处，成功就在路尽头。
1楼 2014-10-21 16:36:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

泽润东方

金虫 (正式写手)
引用回帖:
2楼: Originally posted by zyl雨田 at 2014-10-23 21:12:30
计算优化得到的结果是否可以 一般都是跟实验的参数进行比较，越接近肯定越好。你优化得到的C参数 跟实验值跟别人已经计算得到的值差距还挺大，这样看来至少说明你没有优化好结构。
优化结构，不也就那几个参数了？ ...

最近按照侯老师的建议做了一下。取在实验值的左右取10来个体积值，进行ISIF=4的计算，得到体积---总能的10来对数据，然后状态方程拟合，得到一个能量最低所对应 的体积值（即平衡态下的体积值），然后以此体积值作为基础，计算出一个初始的a, b, c，再对此体积下结构进行ISIF=4计算。
但是计算出来的结果与实验值相差还是很大。但是跟2013年JPCC上面发表的一篇文章的结果比较接近。我想这多半不是我的方法的问题。
赞一下 回复此楼
希望就在不远处，成功就在路尽头。
3楼2014-11-18 13:41:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

匿名

用户注销 (职业作家)
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2014-10-26 19:13:30
泽润东方: 金币+5, ★有帮助 2014-11-18 13:38:20
本帖仅楼主可见
赞一下
2楼2014-10-23 21:12:30
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

匿名

用户注销 (职业作家)
本帖仅楼主可见
赞一下
4楼2014-11-19 10:40:28
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

对马拉二胡

至尊木虫 (职业作家)
我个人觉得其实楼主算的差不多了，因为每个人采用的参数都不同，最后得出的基本性质符合实际就行，文献上的也不一定是对的。
赞一下 回复此楼
相互学习，共同进步。
5楼2014-11-21 19:11:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +3 小贾同学123 2026-03-15 5/250 2026-03-16 20:39 by zhq0425
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-16 4/200 2026-03-16 20:13 by Wangjingyue
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[硕博家园] 深圳大学硕士招生（2026秋，传感器方向，仅录取第一志愿） +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
关闭
请填处理意见
关闭