计算化学公社

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

[新手求助] 双自由基寻找到了稳定的对称破缺波函数,但没有自旋密度

[复制链接 Copy URL]

15

帖子

0

威望

186

eV
积分
201

Level 3 能力者

本帖最后由 天月将白 于 2023-12-7 17:18 编辑

老师您好,请问我实验测到了中性有机分子的EPR信号,计算开壳单重态双自由基的时候找到了稳定的对称破缺波函数,但没有自旋密度
命令行:# UB3LYP/6-31G(d) guess=mix stable   0,1
log文件显示:
The wavefunction is stable under the perturbations considered.
说明是寻找到了稳定的波函数但是


其中<S**2>为0,
且用Multiwfn 自旋密度和布局都是0,是说明是找到了闭壳单重态了吗,但是实验捕捉到了液体和固体的EPR信号




6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2023-12-7 17:43:31 | 只看该作者 Only view this author
尝试wB97XD和M06-2X结合def2-TZVP
北京科音自然科学研究中心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

4289

帖子

4

威望

9544

eV
积分
13913

Level 6 (一方通行)

MOKIT开发者

3#
发表于 Post on 2023-12-7 20:01:52 | 只看该作者 Only view this author
你这波函数是稳定的,但并不是对称破缺波函数。guess=mix只是给一个初猜,不对收敛结果做任何保证。
自动做多参考态计算的程序MOKIT

15

帖子

0

威望

186

eV
积分
201

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-12-7 23:03:02 | 只看该作者 Only view this author
zjxitcc 发表于 2023-12-7 20:01
你这波函数是稳定的,但并不是对称破缺波函数。guess=mix只是给一个初猜,不对收敛结果做任何保证。

谢谢老师,请问怎么判断是不是对称破缺函数呀,是看S**0吗

15

帖子

0

威望

186

eV
积分
201

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-12-7 23:03:51 | 只看该作者 Only view this author
sobereva 发表于 2023-12-7 17:43
尝试wB97XD和M06-2X结合def2-TZVP

谢谢sob老师,试过了没有变化,请问会不会是基态三重态呀,因为EPR确实有数据

92

帖子

0

威望

825

eV
积分
917

Level 4 (黑子)

6#
发表于 Post on 2023-12-8 00:34:50 | 只看该作者 Only view this author
天月将白 发表于 2023-12-7 23:03
谢谢sob老师,试过了没有变化,请问会不会是基态三重态呀,因为EPR确实有数据

算一下看看能量是不是更低

4289

帖子

4

威望

9544

eV
积分
13913

Level 6 (一方通行)

MOKIT开发者

7#
发表于 Post on 2023-12-8 15:20:10 | 只看该作者 Only view this author
天月将白 发表于 2023-12-7 23:03
谢谢老师,请问怎么判断是不是对称破缺函数呀,是看S**0吗

不存在“S**0”。是看<S**2>。
自动做多参考态计算的程序MOKIT

42

帖子

0

威望

263

eV
积分
305

Level 3 能力者

8#
发表于 Post on 2023-12-14 10:36:14 | 只看该作者 Only view this author
我是看自旋密度

5

帖子

0

威望

197

eV
积分
202

Level 3 能力者

9#
发表于 Post on 2023-12-14 16:00:55 | 只看该作者 Only view this author
试一试cam-b3lyp/6-311g*

12

帖子

0

威望

241

eV
积分
253

Level 3 能力者

10#
发表于 Post on 2024-7-14 23:46:24 | 只看该作者 Only view this author
请问最后怎么分析的,我也遇到了这个问题,有guess=mix的和没有的计算出来的单重态能量是一模一样的,并且没有自旋密度,实验上有EPR响应

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

11#
发表于 Post on 2024-7-14 23:48:43 | 只看该作者 Only view this author
dhz 发表于 2024-7-14 16:46
请问最后怎么分析的,我也遇到了这个问题,有guess=mix的和没有的计算出来的单重态能量是一模一样的,并且 ...

单重态双自由基本身没有EPR信号,是和它形成平衡的少量三重态双自由基的信号。这就是为什么证明单重态双自由基往往用温度效应,也就是降低温度EPR信号减小(因为三重态平衡浓度减小)。
应当计算三重态双自由基

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
wjc404 + 2

查看全部评分 View all ratings

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
本团队长期招收研究生,有意者可私信联系

4289

帖子

4

威望

9544

eV
积分
13913

Level 6 (一方通行)

MOKIT开发者

12#
发表于 Post on 2024-7-15 15:04:33 | 只看该作者 Only view this author
补上一下,如果你想做EPR计算,那应该用三重态。但若不做EPR,且体系基态是单重态,若你只想知道单重态极小点结构下的未成对电子数,多组态特征y0值,未成对电子密度等等,则可以用MOKIT做自动CASSCF计算,获得这些分析结果,例子见https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap5-9.html
自动做多参考态计算的程序MOKIT

87

帖子

1

威望

1926

eV
积分
2033

Level 5 (御坂)

吃猫咪的鱼

13#
发表于 Post on 2024-7-15 15:19:56 | 只看该作者 Only view this author
本帖最后由 imasen 于 2024-7-15 15:30 编辑

现在文章里说是基态三重态或者自旋极化单重态的,很多EPR谱图给出的都是个单自由基双重态的信号,一点可信度也没有。极少数确实存在三重态实验证据的文章里也都表明存在大量双重态杂质。如果你要考虑基态三重态的,确保你的实验谱图不是个单峰,可以观察到明显的D值。很多人都拿个双重态单自由基的谱图说是三重态,说D过小测不出来。偶极相互作用就算是7个埃距离粗略估计也有150MHz,对应谱图裂分有10mT/100高斯以上了,不会测不出来。很多情况下只是单自由基信号,在此基础上做变温实验是没有意义的。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
xsc6 + 5 好物!

查看全部评分 View all ratings

12

帖子

0

威望

241

eV
积分
253

Level 3 能力者

14#
发表于 Post on 2024-9-26 09:51:01 | 只看该作者 Only view this author
本帖最后由 dhz 于 2024-9-26 09:52 编辑
wzkchem5 发表于 2024-7-14 23:48
单重态双自由基本身没有EPR信号,是和它形成平衡的少量三重态双自由基的信号。这就是为什么证明单重态双 ...

老师您好,很抱歉现在才看到您的回复,但对于我选的这个体系,我还是有问题不太理解,请您指教,希望您有时间可以回复!
我看到一篇文章,其中一个分子具有EPR响应,但是随温度的升高,响应增强,我判断其为开壳层单重态,但是计算的时候发现这个分子没有对称破缺态。

请问老师这种情况我还是要判断它为三重态双自由基还是闭壳层单重态

2a.png (64.04 KB, 下载次数 Times of downloads: 3)

2a.png

4289

帖子

4

威望

9544

eV
积分
13913

Level 6 (一方通行)

MOKIT开发者

15#
发表于 Post on 2024-9-26 10:50:20 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2024-9-26 12:45 编辑
dhz 发表于 2024-9-26 09:51
老师您好,很抱歉现在才看到您的回复,但对于我选的这个体系,我还是有问题不太理解,请您指教,希望您有 ...

(如果是单重态体系)建议先用CASSCF算个单点看看双自由基特征y0值,别只拿UKS结果下定论。
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 05:43 , Processed in 0.176382 second(s), 24 queries , Gzip On.

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