计算化学公社

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

[新手求助] 含Fe的离子团簇自旋多重度确定求助

[复制链接 Copy URL]

31

帖子

0

威望

352

eV
积分
383

Level 3 能力者

本帖最后由 wuyeyangguang 于 2021-9-14 19:21 编辑

请教一下各位老师,
1. 该离子团簇总电荷为-1,自旋多重度是否为1,自己在b3lyp/6-311+g**级别下试算过自旋多重度为1,3,5,7,9,11的单点,只有为1的时候收敛了,其他都是SCF不收敛,所以没法比较能量。顺便问下SCF不收敛是因为b3lyp根本不适合算激发态吗
2.另外想请教以下几种情况的自旋多重度又分别是多少?单个[FeCl4]-带一个负电荷情况;[FeCl4]-与咪唑阳离子组合成的中性分子;两个咪唑阳离子与一个[FeCl4]-组成的带一个正电荷的离子团;

202109141655091361..jpg (25.93 KB, 下载次数 Times of downloads: 27)

离子团

离子团

4289

帖子

4

威望

9554

eV
积分
13923

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2021-9-14 16:53:38 | 只看该作者 Only view this author
你这算的啥[FeCl4]-,[Fe4Cl]-
自动做多参考态计算的程序MOKIT

150

帖子

0

威望

3925

eV
积分
4075

Level 6 (一方通行)

3#
发表于 Post on 2021-9-14 16:53:43 | 只看该作者 Only view this author
你画的这个是Fe4Cl......

31

帖子

0

威望

352

eV
积分
383

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2021-9-14 17:05:09 | 只看该作者 Only view this author
本帖最后由 wuyeyangguang 于 2021-9-14 19:24 编辑
zjxitcc 发表于 2021-9-14 16:53
你这算的啥[FeCl4]-,[Fe4Cl]-

不好意思画错图了, 已修正。

31

帖子

0

威望

352

eV
积分
383

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-9-14 17:05:39 | 只看该作者 Only view this author
RAL 发表于 2021-9-14 16:53
你画的这个是Fe4Cl......

画错了,不好意思。

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

6#
发表于 Post on 2021-9-14 19:32:32 | 只看该作者 Only view this author
这种含有不止一个自旋中心的体系,应该用broken symmetry来算,除了最高自旋的态可以直接设自旋多重度以外,较低自旋的态都应该要么用片段波函数,要么用guess=mix等方法产生初猜,不能直接设个自旋多重度就跑。一定要直接设个自旋多重度就跑的话,也得做stable=opt。
不收敛就解决不收敛的问题,按照http://sobereva.com/61里的方法,没有任何遗漏地试一遍,仍然不收敛,才叫不收敛。不能一发现默认的SCF设置不收敛就放弃
这种体系应当用配体场理论预测自旋多重度,不然我们告诉你你也是死记硬背,换个体系你又不会判断了
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

31

帖子

0

威望

352

eV
积分
383

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-9-14 21:07:02 | 只看该作者 Only view this author
wzkchem5 发表于 2021-9-14 19:32
这种含有不止一个自旋中心的体系,应该用broken symmetry来算,除了最高自旋的态可以直接设自旋多重度以外 ...

谢谢了,回头仔细学下配体场理论。对于当前这个体系,带一个负电荷的阴离子团似乎1,3,11都有可能,1可能概率大点?单独的一个阴离子貌似是6;单独的一个阳离子肯定是1,阴阳离子组合的中性分子的话是2吗。不知道对不对,再次感谢。

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

8#
发表于 Post on 2021-9-14 21:45:08 | 只看该作者 Only view this author
wuyeyangguang 发表于 2021-9-14 14:07
谢谢了,回头仔细学下配体场理论。对于当前这个体系,带一个负电荷的阴离子团似乎1,3,11都有可能,1可 ...

建议再仔细学习一下自旋多重度的定义。你这里如果阴离子的自旋多重度是6,阳离子是1,那么一个阴离子和一个阳离子的非共价复合物必定是6。只要了解自旋多重度的定义就会知道为什么。
至于包含两个阴离子的体系,取决于两个阴离子之间的J值是正的还是负的,1、11都有可能,应当分别计算,其中1用broken symmetry算(而不是直接设自旋多重度为1),取能量较低的(严格一些的做法是取自由能较低的)
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125189

管理员

公社社长

9#
发表于 Post on 2021-9-15 01:47:55 | 只看该作者 Only view this author
注意看下文了解对称破缺计算是怎么回事
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

SCF不收敛和当前泛函适不适合计算没有必然关系
北京科音自然科学研究中心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

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

10#
发表于 Post on 2021-9-15 23:14:34 | 只看该作者 Only view this author
wzkchem5 发表于 2021-9-14 21:45
建议再仔细学习一下自旋多重度的定义。你这里如果阴离子的自旋多重度是6,阳离子是1,那么一个阴离子和一 ...

如果两边离子分别5alpha5beta的话光用u+guess=mix能收敛到正确的态嘛,感觉要设置片段初猜比较稳一些吧

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

11#
发表于 Post on 2021-9-15 23:31:55 | 只看该作者 Only view this author
喵星大佬 发表于 2021-9-15 16:14
如果两边离子分别5alpha5beta的话光用u+guess=mix能收敛到正确的态嘛,感觉要设置片段初猜比较稳一些 ...

对,高斯里用片段初猜比较稳妥,orca里做flipspin就行了
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

31

帖子

0

威望

352

eV
积分
383

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2021-9-21 20:10:41 | 只看该作者 Only view this author
谢谢社长及各位大佬。在学习了社长帖子之后,做了如下计算:
首先在b3lyp/6-311+g**, em=gd3下对闭壳层单重态进行了几何优化,得到稳定构型并获得该计算级别下的单点能。之后针对优化后的构型,在同样计算级别下指定U并采用guess=mix nosymm stable=opt关键词计算了对称破缺单重态能量,以及自旋多重度为11的能量,得到的能量顺序为:
-6554.48293069(闭壳层单重态)>6554.57287364 (对称破缺单重态) > -6554.64174730 (自旋多重度为11的开壳层)
虽然b3lyp算能量有点次了,但能量相差还是比较大的,是不是也可以说明该体系的基态其实就是自旋多重度为11?
谢谢

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

13#
发表于 Post on 2021-9-21 20:55:53 | 只看该作者 Only view this author
wuyeyangguang 发表于 2021-9-21 13:10
谢谢社长及各位大佬。在学习了社长帖子之后,做了如下计算:
首先在b3lyp/6-311+g**, em=gd3下对闭壳层单 ...

guess=mix主要适用于单电子数等于2的情形,对于单电子数大于2的体系,guess=mix成功率比较低,有比较大概率收敛到错误的态。这种体系算完以后不能只看能量,必须先看自旋布居数,确定自旋布居是对的,再比较能量。如果自旋布居不对,需要换其他方法找对称破缺单重态
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

31

帖子

0

威望

352

eV
积分
383

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2021-9-22 11:29:49 | 只看该作者 Only view this author
wzkchem5 发表于 2021-9-21 20:55
guess=mix主要适用于单电子数等于2的情形,对于单电子数大于2的体系,guess=mix成功率比较低,有比较大概 ...

谢谢大佬指点。根据提示做了自旋布居分析,由于计算中加了弥散函数,因此用了Becke方法。闭壳层单重态下体系总自旋布居数是-0.00000110, 对称破缺态是-0.00000024,自旋多重度11时自旋布居数是9.99999840。是不是可以下结论是11了。

51

帖子

0

威望

835

eV
积分
886

Level 4 (黑子)

15#
发表于 Post on 2021-9-22 11:53:09 | 只看该作者 Only view this author
总布居数是什么鬼……要看的是原子布居或者片段布居吧

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 09:38 , Processed in 0.177352 second(s), 23 queries , Gzip On.

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