24小时热门版块排行榜    

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

木虫 (小有名气)

【答案】应助回帖

个人觉得,误差2%~4%不算太大。DFT本来就有误差。
你将截断能增加到600eV,K点增到10*10*4再试试。
对于简单晶体的优化,不需要降低收敛阀值。
如果参数设对了,fine也是很容易收敛的。
good luck
8楼2014-07-19 20:53:54
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +14 junqihahaha 2026-03-26 15/750 2026-03-30 14:12 by 探123
[考研] 南京大学化学调剂 +10 景随风 2026-03-29 15/750 2026-03-30 11:21 by limeifeng
[考研] 一志愿中南大学化学0703总分337求调剂 +6 niko- 2026-03-27 6/300 2026-03-30 10:25 by herarysara
[考研] 375求调剂 +6 雨夏整夜 2026-03-29 6/300 2026-03-30 10:21 by herarysara
[考研] 290求调剂 +3 dfffsar 2026-03-29 3/150 2026-03-29 22:38 by 毛毛毛阿莫2
[考博] 26申博自荐 +6 whh869393 2026-03-24 6/300 2026-03-29 21:05 by nxgogo
[考研] 一志愿北京理工大学本科211材料工程294求调剂 +8 mikasa的围巾 2026-03-28 8/400 2026-03-29 12:48 by 无际的草原
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +7 @taotao 2026-03-29 7/350 2026-03-29 12:03 by longlotian
[考研] 调剂求院校招收 +6 鹤鲸鸽 2026-03-28 6/300 2026-03-29 08:15 by fmesaito
[考研] 调剂 +3 好好读书。 2026-03-28 3/150 2026-03-28 12:04 by 王保杰33
[考研] 求调剂 +4 零八# 2026-03-27 4/200 2026-03-27 18:07 by yu221
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 材料考研求调剂 +3 Dendel 2026-03-23 6/300 2026-03-26 17:51 by fmesaito
[考研] 材料科学与工程 317求调剂 +4 JKSOIID 2026-03-26 4/200 2026-03-26 15:58 by 不吃魚的貓
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
信息提示
请填处理意见