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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于缩放系数的优化
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1901  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +14 超级伊昂大王 2026-04-06 14/700 2026-04-08 07:03 by 无际的草原
[考研] 材料科学与工程320求调剂，080500 +9 黄瓜味薯片 2026-04-06 9/450 2026-04-08 00:46 by JourneyLucky
[考研] 化学308分求调剂 +14 你好明天你好 2026-04-07 16/800 2026-04-07 23:09 by JourneyLucky
[考研] 0854电子信息319求调剂（接受跨专业调剂） +5 星星不眨眼喽 2026-04-05 6/300 2026-04-07 22:16 by hemengdong
[考研] 085600材料与化工，求调剂 +5 won_qii 2026-04-07 5/250 2026-04-07 17:10 by 啵啵啵0119
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 化学调剂 +18 艾志恒 2026-04-03 19/950 2026-04-07 16:00 by 起飞的比熊1
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 求调剂 +5 小沢 2026-04-03 5/250 2026-04-06 22:45 by 875465
[考研] 材料334求调剂 +19 Eecho# 2026-04-03 19/950 2026-04-06 08:37 by 小小树2024
[考研] 070300化学学硕311分求调剂 +11 梁富贵险中求 2026-04-04 13/650 2026-04-06 07:24 by houyaoxu
[考研] 320分人工智能调剂 +8 振—TZ 2026-04-03 8/400 2026-04-05 22:33 by 范式思维
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 11408 一志愿西电，277分求调剂 +4 zhouzhen654 2026-04-03 4/200 2026-04-04 18:10 by 猪会飞
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 一志愿北交大材料工程，总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[硕博家园] 考研调剂 +5 骆驼男人 2026-04-01 5/250 2026-04-01 14:28 by syjjj0321
信息提示
关闭
请填处理意见
关闭