24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2555  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yeahgyz

铜虫 (初入文坛)

[求助] 黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决 已有1人参与

小弟按照文献所述参数对黄铜矿进行了几何结构的优化,但是却一直不能收敛是怎么回事,而且算得的晶胞常数与实验值相差超过了2%,甚至达到了4%或5%,是什么地方出现了问题吗?下面吧具体模拟过程及文献提供给各位大神分析,谢谢!

文献中的设置如下:黄铜矿,理想黄铜矿的晶胞参数为:a=b=5.289 Å,c=10.423 Å,α=β=γ=90°。采用GG-PBE梯度修正函数,在进行几何优化时,选取合适的截断能和k点,收敛标准为:原子位移为0.005 Å,原子间作用力为0.01 eV/nm,原子间的内应力为0.2 GPa,体系总能量的变化不大于5.0×10-4 eV/atom。
截断能和k点选择过程及数据如图一和图二:

我按照245计算,计算到98步进行不下去了,受不了就停止了。
手头上有一个270ev的数据,收敛过程如图三图四,设置参数如图五图六。

另外我算了E点为310,350和400的,k点位5x5x2,很快收敛,但是优化后的晶胞常数误差超过2%,最大的达到4%,应该是不能用的。
请问1.为什么我按照文献设置不能收敛?
       2.为什么提高E点可以收敛,但是晶胞常数误差过大?是不是K点设置问题,还是 其他问题?

谢谢指教!

黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决
360截图20140718223046593.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-1
360截图20140718223114573.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-2
360截图20140718223437784.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-3
360截图20140718223452372.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-4
360截图20140718223536488.jpg


黄铜矿几何结构优化 按照文献参数不能收敛,且晶胞常数误差过大如何解决-5
360截图20140718223610451.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yeahgyz

铜虫 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by huxiao.gl at 2014-07-19 02:48:21
截断能至少350~400eV才合理. 300eV以下算出来的结果基本没有意义。

我看黄铜矿文献上取的245,还有一篇取的是280。如果提高截短能,晶胞常数会不会更趋于实验值,误差更小些?别人为什么截短能200多就能收敛,我的不能呢?

[ 发自小木虫客户端 ]
3楼2014-07-19 09:01:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

huxiao.gl

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
截断能至少350~400eV才合理. 300eV以下算出来的结果基本没有意义。

» 本帖已获得的红花(最新10朵)

good luck
2楼2014-07-19 02:48:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

测试过程建议采用结构单点,一旦确定参数,做结构优化,并在优化结构的单点做参数的检验,这样可以节约时间。最后确定截断能和k不是随便设定一个,而是由测试确定,如果很小的截断能能够满足测试条件,直接使用也无妨(当然我觉得不太可能)。如果测试不能达到标准,或是标准较低,导致的结果就是计算“力”及“应力”会产生偏差,导致结构优化失败。除了这两个参数,空带数目是否充分也很重要。

原则上来说,截断能和k越充分结果是越好的,但是由于DFT的特性,于是不见得是这样。对于某些结构,PBEsol也许比PBE好得多。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-07-19 09:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huxiao.gl

木虫 (小有名气)

引用回帖:
3楼: Originally posted by yeahgyz at 2014-07-19 09:01:42
我看黄铜矿文献上取的245,还有一篇取的是280。如果提高截短能,晶胞常数会不会更趋于实验值,误差更小些?别人为什么截短能200多就能收敛,我的不能呢?
...

245即使与实验值吻合,也是巧合,没有意义的。
不用把时间浪费在这个问题上了。
good luck
5楼2014-07-19 16:42:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历 +8 cccchenso 2026-03-29 8/400 2026-03-29 23:53 by 我是小康
[考研] 327求调剂 +3 小卡不卡. 2026-03-29 3/150 2026-03-29 22:29 by 无际的草原
[考研] 求调剂 +4 QiMing7 2026-03-25 5/250 2026-03-29 21:10 by 唐沐儿
[考研] 各位老师好,我的一志愿为北京科技大学085601材料专硕 +9 Koxui 2026-03-28 9/450 2026-03-29 20:43 by 无际的草原
[考研] 一志愿南昌大学324求调剂 +5 hanamiko 2026-03-29 5/250 2026-03-29 16:46 by 学员8dgXkO
[考研] 一志愿郑州大学,080500学硕,总分317分求调剂 +8 举个栗子oi 2026-03-24 9/450 2026-03-29 13:08 by peike
[考研] 343求调剂 +6 爱羁绊 2026-03-29 6/300 2026-03-29 12:00 by 无际的草原
[考研] 2026年华南师范大学欢迎化学,化工,生物,生医工等专业优秀学子加入! +3 llss0711 2026-03-28 6/300 2026-03-29 10:26 by llss0711
[考研] 一志愿211 初试270分 求调剂 +8 谷雨上岸 2026-03-23 9/450 2026-03-29 01:00 by 我是小康
[考研] 312,生物学求调剂 +3 小译同学abc 2026-03-28 3/150 2026-03-28 15:32 by 落睿可思
[考研] 266分,求材料冶金能源化工等调剂 +7 哇呼哼呼哼 2026-03-27 9/450 2026-03-28 12:22 by zllcz
[考研] 291求调剂 +6 HanBeiNingZC 2026-03-24 6/300 2026-03-28 07:55 by baoball
[考研] 070300化学求调剂 +4 起个名咋这么难 2026-03-27 4/200 2026-03-27 21:39 by 83503孙老师
[论文投稿] Journal of Mechanical Science and Technology +3 Russ_ss 2026-03-25 5/250 2026-03-27 10:49 by 陆小果画大饼
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 317求调剂 +7 蛋黄咸肉粽 2026-03-26 7/350 2026-03-27 02:29 by fmesaito
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
信息提示
请填处理意见