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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于缩放系数的优化
51/1返回列表
查看: 1862  |  回复: 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的回帖

匿名

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

泽润东方

金虫 (正式写手)
引用回帖:
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的回帖

匿名

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

对马拉二胡

至尊木虫 (职业作家)
我个人觉得其实楼主算的差不多了，因为每个人采用的参数都不同，最后得出的基本性质符合实际就行，文献上的也不一定是对的。
赞一下 回复此楼
相互学习，共同进步。
5楼2014-11-21 19:11:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 泽润东方 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[硕博家园] 深圳大学硕士招生（2026秋，传感器方向，仅录取第一志愿） +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭