计算化学公社

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

[新手求助] 一篇文章计算casscf时进行定域化,其系数和电子组态的疑问

[复制链接 Copy URL]

10

帖子

0

威望

134

eV
积分
144

Level 2 能力者

    各位老师好,最近看到一篇文章10.1002/anie.202423347,我对其中的figure3产生了疑问,发邮件给通讯作者一周多也没理我,所以来论坛问问各位老师。
    在论坛查找了关于casscf方法的帖子,能比较浅显地知道左半部分的意思,但是对于中间的等式和右边“localized orbitals”不太理解。按照我理解的意思,右边的图应该是:第一步用uPBE-D3/def2-TZVP guess=mix得到alpha和beta两组电子的信息,第二步是用uPBE-D3/def2-TZVP guess=(read,naturalorbitals,save,only)得到轨道系数,第三步确定了活性空间后用CASSCF(2,2)/def2-TZVP计算得到相应的特征值和系数。
    右边的图我是这样理解的:在上述过程第一步时,就进行定域化,得到相应的轨道图像。但是我不理解的是下面的ψ1、ψ2、ψ3和c1、c2、c3在右半部分代表什么意思,并且是怎么计算出来的呢?
    希望各位老师批评指正。


1.pdf (1.75 MB, 下载次数 Times of downloads: 5)




1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

2#
发表于 Post on 2025-2-10 15:57:50 | 只看该作者 Only view this author
上面已经写了c_i和ψ_i是什么意思了,令i=1,2,3就知道ψ1、ψ2、ψ3和c1、c2、c3是什么意思了
定域化的方式不止一种,作者应当在正文或SI的某处写了具体是怎么做的,如果确定作者没写,那么是作者的责任,只能问作者
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

10

帖子

0

威望

134

eV
积分
144

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-2-10 21:06:17 | 只看该作者 Only view this author
wzkchem5 发表于 2025-2-10 15:57
上面已经写了c_i和ψ_i是什么意思了,令i=1,2,3就知道ψ1、ψ2、ψ3和c1、c2、c3是什么意思了
定域化的方 ...

谢谢老师解答,作者在文中和SI都没提及定域化的方式,我应该再发邮件询问。
关于右图,可能我的问题描述的不太清楚,我想问老师的是:1、右边的localized orbitals和图中的两条式子有什么关系呢;2、如果相应的c_i是用casscf计算出来的,那么这在gaussian上用什么关键词实现呢?

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

4#
发表于 Post on 2025-2-11 00:32:02 | 只看该作者 Only view this author
jim_lee_ 发表于 2025-2-10 21:06
谢谢老师解答,作者在文中和SI都没提及定域化的方式,我应该再发邮件询问。
关于右图,可能我的问题描述 ...

定域化轨道可以从分子轨道做酉变换(属于线性变换)得到,下文做了介绍。如果没说定域化方法,大概率是很常用的Pipek-Mezey定域化
Multiwfn的轨道定域化功能的使用以及与NBO、AdNDP分析的对比
http://sobereva.com/380http://bbs.keinsci.com/thread-6053-1-1.html

任何程序做CASSCF计算都可以给出组态系数。Gaussian也是如此。参考北京科音高级量子化学培训班(http://www.keinsci.com/KAQC)讲Gaussian做CASSCF计算的幻灯片的例子:

北京科音自然科学研究中心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

10

帖子

0

威望

134

eV
积分
144

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-2-11 12:11:10 | 只看该作者 Only view this author
sobereva 发表于 2025-2-11 00:32
定域化轨道可以从分子轨道做酉变换(属于线性变换)得到,下文做了介绍。如果没说定域化方法,大概率是很 ...

谢谢老师解答。经过站内查找我明白做CASSCF计算后如何在输出文件中找到相应的组态系数和贡献值。但不明白的是如何将定域化的轨道(如右图的ψ1)和另外两条轨道(右图的ψ2,ψ3)放在一起去做CASSCF计算呢。

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

6#
发表于 Post on 2025-2-11 17:18:49 | 只看该作者 Only view this author
jim_lee_ 发表于 2025-2-11 05:11
谢谢老师解答。经过站内查找我明白做CASSCF计算后如何在输出文件中找到相应的组态系数和贡献值。但不明白 ...

图中画的ψ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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

7#
发表于 Post on 2025-2-12 05:21:42 | 只看该作者 Only view this author
jim_lee_ 发表于 2025-2-11 12:11
谢谢老师解答。经过站内查找我明白做CASSCF计算后如何在输出文件中找到相应的组态系数和贡献值。但不明白 ...

ψ1代表1号CSF(configurational state function,组态函数)的展开系数
如果不了解组态函数,看Szabo的Modern quantum chemistry相关部分
北京科音自然科学研究中心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

10

帖子

0

威望

134

eV
积分
144

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2025-2-12 09:49:50 | 只看该作者 Only view this author
wzkchem5 发表于 2025-2-11 17:18
图中画的ψ1是基于定域化轨道构建的组态,不是定域化轨道本身

好的谢谢老师!我去认真学习一下!

10

帖子

0

威望

134

eV
积分
144

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2025-2-12 09:50:11 | 只看该作者 Only view this author
sobereva 发表于 2025-2-12 05:21
ψ1代表1号CSF(configurational state function,组态函数)的展开系数
如果不了解组态函数,看Szabo的 ...

好的谢谢老师!我认真研读一下!

本版积分规则 Credits rule

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

GMT+8, 2025-8-17 02:13 , Processed in 0.733137 second(s), 23 queries , Gzip On.

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