24小时热门版块排行榜    

查看: 2850  |  回复: 9

jiayoulongnv

新虫 (著名写手)

[求助] Error: Overlap matrix is not positive definite.怎么破?

Dmol3计算氧气分子在金属氧化物表面的结构优化计算中出错信息如下,是什么原因?怎么破?

Error: Overlap matrix is not positive definite.
This most likely means the basis set contains linear dependencies.
Consider increasing the atomic cutoff parameter or changing the basis set.

给出的建议中提到的atomic cutoff parameter是在哪里设置的?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
漫天飘雪: 金币+2, 谢谢交流 2016-08-15 16:17:33
蛋蛋小童鞋: 应助指数+1 2016-08-15 21:32:33
jiayoulongnv: 金币+20, ★★★★★最佳答案, 非常感谢,下次遇到问题再来请教。 2016-08-17 11:25:19
这个就是轨道截断太小导致的。一般DNP级别的基组还不足以产生线性依赖,把Orbital cutoff适当放大就好了。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2016-08-15 14:57:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiayoulongnv

新虫 (著名写手)

嗯嗯,基组DNP不用改,只把cutoff调大就行了吗?精度fine时,cutoff是4.9,迭代11次之后出现上述错误并终止,我在4.6

发自小木虫Android客户端
3楼2016-08-15 15:05:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiayoulongnv

新虫 (著名写手)

我在4.9埃的基础上再调大?

发自小木虫Android客户端
4楼2016-08-15 15:06:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by jiayoulongnv at 2016-08-15 15:05:26
嗯嗯,基组DNP不用改,只把cutoff调大就行了吗?精度fine时,cutoff是4.9,迭代11次之后出现上述错误并终止,我在4.6

一般来说默认是fine应该不大会出现这样的问题,你不妨把计算参数和结构都发来看看。
不一定挂在论坛,计算问题问题欢迎留言。
5楼2016-08-15 22:47:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiayoulongnv

新虫 (著名写手)

引用回帖:
5楼: Originally posted by 卡开发发 at 2016-08-15 22:47:47
一般来说默认是fine应该不大会出现这样的问题,你不妨把计算参数和结构都发来看看。...

我把结构稍微细调了一下,并且把cutoff调为5.1埃,迭代了63次后,于今天早上收敛成功,当然,计算结果我还没有具体分析。我感觉MS的错误提示以及修改建议,还是有非常强的针对性的。谢谢你,在网上看到了你的很多热心回答,下次遇到问题再来请教你。承诺的金币,怎么送给你?
6楼2016-08-17 11:22:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
6楼: Originally posted by jiayoulongnv at 2016-08-17 11:22:08
我把结构稍微细调了一下,并且把cutoff调为5.1埃,迭代了63次后,于今天早上收敛成功,当然,计算结果我还没有具体分析。我感觉MS的错误提示以及修改建议,还是有非常强的针对性的。谢谢你,在网上看到了你的很多热 ...

这个也只是保证了收敛而已,也许只是巧合,我觉得你还是把计算的模型和outmol都发来看下比较稳妥。
不一定挂在论坛,计算问题问题欢迎留言。
7楼2016-08-17 12:54:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiayoulongnv

新虫 (著名写手)

引用回帖:
7楼: Originally posted by 卡开发发 at 2016-08-17 12:54:01
这个也只是保证了收敛而已,也许只是巧合,我觉得你还是把计算的模型和outmol都发来看下比较稳妥。...

好的,你帮我看看 。另外,我想根据outmol文件里的能量信息来计算吸附能,应该是最终总能量和初始总能量做差,对吧?但是问题是,outmol文件中,最初总能量和最终总能量,怎么找到呢?

» 本帖附件资源列表

8楼2016-08-17 17:01:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiayoulongnv

新虫 (著名写手)

引用回帖:
7楼: Originally posted by 卡开发发 at 2016-08-17 12:54:01
这个也只是保证了收敛而已,也许只是巧合,我觉得你还是把计算的模型和outmol都发来看下比较稳妥。...

另外,吸附能是有正负之分的,是吧?对于负的吸附能,怎么比较大小呢?比如,-5eV的吸附能和-2eV的吸附能,谁大谁小呢?
9楼2016-08-17 17:04:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
9楼: Originally posted by jiayoulongnv at 2016-08-17 17:04:39
另外,吸附能是有正负之分的,是吧?对于负的吸附能,怎么比较大小呢?比如,-5eV的吸附能和-2eV的吸附能,谁大谁小呢?...

问题比较多,我列在下面:
1、我不知道你的体系O2你现在的参数优化完是否O2在表面是非解离的状态,如果是的话,体系很有可能是三重态的,或者说磁矩为2,吸附的两个O上应当各自给一个单电子,这个在modify->electonic configuration中设置,我想这个可能是个主要的因素;
2、真空层不够高,DMol3的计算真空层对计算量贡献不大(与平面波完全不同),所以不要太吝啬;
3、体系使用的smearing有点大,一般0.003以下才可能有好的结果,得控制outmol中-kTeS,要足够小;
4、有条件的话Mulitipolar Expansion用Hex吧;
5、吸附能定义为E_ad=E_AB-E_A-EB。能量在total energy里面看,看最后一个。吸附能一般用更正或者更负来形容,严格来说,吸附能没啥确切的意义,得看吸附自由能,仅当熵变贡献不大的时候和吸附自由能才接近,或者,两个吸附体系进行比较的时候熵变贡献接近,仅仅比较能量才有点意义。
不一定挂在论坛,计算问题问题欢迎留言。
10楼2016-08-17 20:34:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiayoulongnv 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:12 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见