计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3606|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[波函数分析求助] 高斯计算三重态的自旋不对,如何解决

[复制链接 Copy URL]

17

帖子

0

威望

187

eV
积分
204

Level 3 能力者

如图,高斯计算后的log文件,查找spin后,Cu的自旋为一正一负,我计算的是三重态,按道理来说应该是两个接近0.5的正数,请问各位大佬们,有什么解决办法

4289

帖子

4

威望

9544

eV
积分
13913

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2021-11-5 10:00:11 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2021-11-5 10:12 编辑

(1)scfcon=5算出来结果没有实际意义,不能用于发表,作为粗结果、后续再恢复正常精度计算是可以的,但直接拿来发表的话是严重的垃圾数据。不是说程序不报错、结果就能用,这是错误的观点。
(2)复杂结构应该先算单点,检验波函数稳定性(stable=opt),检查自旋布居与期望的是否有偏差。然后再读取波函数做结构优化。如果波函数不稳定,布居分析是没有意义的。

由于你已经做了一次结构优化,不妨就在此结构上(记得恢复正常收敛精度)算单点、检验波函数稳定性,看看结果如何。
自动做多参考态计算的程序MOKIT

17

帖子

0

威望

187

eV
积分
204

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-11-5 10:11:40 | 只看该作者 Only view this author
zjxitcc 发表于 2021-11-5 10:00
(1)scfcon=5算出来结果没有实际意义,不能用于发表,作为粗结果、后续再恢复正常精度计算是可以的,但直 ...

老师,请问若自旋与期望值有偏差,以什么办法解决呢

4289

帖子

4

威望

9544

eV
积分
13913

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2021-11-5 10:13:37 | 只看该作者 Only view this author
LYFE 发表于 2021-11-5 10:11
老师,请问若自旋与期望值有偏差,以什么办法解决呢

可能性:
(1)自己的化学直觉有误,猜想有误;
(2)自己的化学直觉无误,SCF存在多个解,甚至多个稳定的解,需要找出其他的解。除了确保波函数稳定外,还可以借助《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82 找出其他的SCF解
自动做多参考态计算的程序MOKIT

17

帖子

0

威望

187

eV
积分
204

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-11-5 10:16:36 | 只看该作者 Only view this author
zjxitcc 发表于 2021-11-5 10:13
可能性:
(1)自己的化学直觉有误,猜想有误;
(2)自己的化学直觉无误,SCF存在多个解,甚至多个稳 ...

老师,还有一个问题,高斯计算中scfcon是代表精度嘛,scfcon等于的数值越高越好嘛

4289

帖子

4

威望

9544

eV
积分
13913

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2021-11-5 10:21:42 | 只看该作者 Only view this author
LYFE 发表于 2021-11-5 10:16
老师,还有一个问题,高斯计算中scfcon是代表精度嘛,scfcon等于的数值越高越好嘛

高斯官网都有http://gaussian.com/scf,相当于scf(conver=N),N越大越精确,太大没必要,徒增耗时,但不要比默认值粗糙。
自动做多参考态计算的程序MOKIT

17

帖子

0

威望

187

eV
积分
204

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-11-5 10:22:46 | 只看该作者 Only view this author
zjxitcc 发表于 2021-11-5 10:21
高斯官网都有http://gaussian.com/scf,相当于scf(conver=N),N越大越精确,太大没必要,徒增耗时,但不 ...

好的,感谢老师!!!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

8#
发表于 Post on 2021-11-5 16:54:06 | 只看该作者 Only view this author
你都不给出体系具体结构,这个问题根本没法回答,别人都不知道两个Cu的化学环境是否相同、位置是否是对称的
但凡用gen的时候,具体用的什么基组都必须说清楚,下文明确强调了。仔细看下文再提问
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620http://bbs.keinsci.com/thread-25787-1-1.html

不懂的关键词一律不要随便加,scf=conver=x(等同于scfcon=x)怎么用好好看下文
解决SCF不收敛问题的方法
http://sobereva.com/61

写iop(1/8=6)完全莫名其妙。有maxstep现成的关键词不用用这么抽象的IOp做什么,而且不加notrust的话根本起不到实际作用。maxstep怎么用、什么时候才该用仔细看
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

绝对不要从别人那里瞎抄关键词。我担心你在入门阶段看了不少乱七八糟的东西已被误导,建议看
谈谈学量子化学如何正确地入门
http://sobereva.com/355http://bbs.keinsci.com/thread-4447-1-1.html

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

17

帖子

0

威望

187

eV
积分
204

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-11-7 23:18:28 | 只看该作者 Only view this author
sobereva 发表于 2021-11-5 16:54
你都不给出体系具体结构,这个问题根本没法回答,别人都不知道两个Cu的化学环境是否相同、位置是否是对称的 ...

好的,我是个新手,才学习,谢谢老师的教导!!!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 06:59 , Processed in 0.172524 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list