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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于缩放系数的优化
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1893  |  回复: 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的回帖

匿名

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

泽润东方

金虫 (正式写手)
引用回帖:
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楼2014-11-21 19:11:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物与医药273求调剂 +8 荔题南墙 2026-04-05 8/400 2026-04-06 22:40 by chenzhimin
[考研] 277求调剂 +3 倪建设 2026-04-06 3/150 2026-04-06 22:39 by qlm5820
[考研] 材料调剂 +14 小刘同学吖吖 2026-04-06 15/750 2026-04-06 22:37 by qlm5820
[考研] 求调剂 +5 熊二想上岸 2026-04-06 5/250 2026-04-06 22:27 by chenzhimin
[考研] 生物学调剂 可调剂到生物与医药 +3 李政莹 2026-04-06 3/150 2026-04-06 19:02 by macy2011
[考研] 308求调剂 +13 倘若起风了呢 2026-04-05 13/650 2026-04-06 14:20 by 蒋皓禹
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 312求调剂 +3 Say Never 2026-04-04 3/150 2026-04-05 11:19 by guoweigw
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 0705理学294求调剂 +3 成果成果cg5 2026-04-03 3/150 2026-04-03 14:04 by simons1972
[考研] 一志愿安徽大学0817化学工程与技术，求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 0856初试324分求调剂 +6 想上学求调 2026-04-01 6/300 2026-04-02 11:42 by 星空星月
[考研] 271求调剂 +15 勒布朗@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 材料调剂 +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
信息提示
关闭
请填处理意见
关闭