24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2562  |  回复: 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的回帖

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的回帖
查看全部 10 个回答

huxiao.gl

木虫 (小有名气)

【答案】应助回帖

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

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

good luck
2楼2014-07-19 02:48:21
已阅   回复此楼   关注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的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

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

原则上来说,截断能和k越充分结果是越好的,但是由于DFT的特性,于是不见得是这样。对于某些结构,PBEsol也许比PBE好得多。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-07-19 09:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料科学与工程调剂 +5 深V宿舍吧 2026-03-30 5/250 2026-03-30 13:46 by 唐沐儿
[考研] 332求调剂 +12 032500 2026-03-25 12/600 2026-03-30 12:16 by 学员8dgXkO
[考研] 材料工程专硕求调剂 +6 hyl3153942 2026-03-29 6/300 2026-03-30 12:04 by 1018329917
[考研] 375求调剂 +6 雨夏整夜 2026-03-29 6/300 2026-03-30 10:21 by herarysara
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历 +8 cccchenso 2026-03-29 8/400 2026-03-29 23:53 by 我是小康
[考研] 294分080500材料科学与工程求调剂 +8 柳溪边 2026-03-26 8/400 2026-03-29 20:42 by 唐沐儿
[考研] 求收留 +5 1943443204 2026-03-28 5/250 2026-03-29 12:11 by 无际的草原
[考研] 085600,专业课化工原理,321分求调剂 +5 大馋小子 2026-03-28 5/250 2026-03-29 08:56 by qingfeng258
[考研] 生物学学硕,一志愿湖南大学,初试成绩338 +6 YYYYYNNNNN 2026-03-26 7/350 2026-03-28 20:52 by 唐沐儿
[考研] 071000生物学求调剂,初试成绩343 +7 小小甜面团 2026-03-25 7/350 2026-03-28 20:25 by 唐沐儿
[考研] 322求调剂 +5 旧吢 2026-03-24 5/250 2026-03-28 13:26 by Iveryant
[考研] 328求调剂 +7 嗯滴的基本都 2026-03-27 7/350 2026-03-28 04:19 by fmesaito
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 求调剂323材料与化工 +7 1124361 2026-03-24 7/350 2026-03-27 10:22 by wangjy2002
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
信息提示
请填处理意见