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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于缩放系数的优化
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1897  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +11 AZMK 2026-04-05 17/850 2026-04-07 19:05 by 中飞院空管学院�
[考研] 388调剂 +3 玟炡winter 2026-04-01 3/150 2026-04-07 18:33 by ?丶丶丶
[考研] 生物与医药273求调剂 +14 荔题南墙 2026-04-05 15/750 2026-04-07 16:29 by 可口可乐不加冰�
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +9 Mr. Z 2026-04-05 9/450 2026-04-07 15:30 by 尽舜尧1
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 材料工程专业日语生求调剂 +7 111623 2026-04-07 7/350 2026-04-07 14:33 by 诗与自由
[考研] 0702物理学学硕299求调剂 +3 祁柒连 2026-04-06 3/150 2026-04-07 12:48 by 小木虫gyh
[考研] 0703化学 +22 妮妮ninicgb 2026-04-04 25/1250 2026-04-07 10:35 by plum
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 生物学调剂 可调剂到生物与医药 +3 李政莹 2026-04-06 3/150 2026-04-06 19:02 by macy2011
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 0832食品科学与工程学硕282调剂 +6 鱼在水中游a 2026-04-02 9/450 2026-04-05 11:45 by flysky1234
[考研] 考研调剂 +3 mcbbc 2026-04-04 3/150 2026-04-05 10:03 by barlinike
[考研] 材料调剂 +7 dxy调剂 2026-04-04 7/350 2026-04-05 09:15 by 陌秋26
[考研] 一志愿上海大学生物学346 +3 上海大学346调剂 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 初试成绩337找调剂 +3 ??? ?. ? 2026-04-03 3/150 2026-04-03 11:43 by 土木硕士招生
[考研] 0710生物学，325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
信息提示
关闭
请填处理意见
关闭