计算化学公社

 找回密码 Forget password
 注册 Register

使用Gaussian做镧系金属配合物的量子化学计算

查看数: 50908 | 评论数: 44 | 收藏 Add to favorites 34
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-1-13 05:00

正文摘要:

使用Gaussian做镧系金属配合物的量子化学计算Using Gaussian to perform quantum chemistry calculations for lanthanide metal complexes 文/Sobereva@北京科音  2021-Jan-13 0 前言 过渡金属配 ...

回复 Reply

yjb 发表于 Post on 2023-5-28 21:30:16
sobereva 发表于 2023-5-28 01:59
可行
这总比你只用小核赝势算的电子能量来讨论强

谢谢,卢老师。
yjb 发表于 Post on 2023-5-28 21:23:14
wzkchem5 发表于 2023-5-28 16:15
有一种特殊情况下我感觉可能需要手动校正一下电子熵。假如你的体系不仅有稀土,还有其他未成对电子,而且 ...

谢谢,王老师
wzkchem5 发表于 Post on 2023-5-28 16:18:06
护卫天使 发表于 2022-10-20 03:51
老师好,我想再请教一下,1、计算水中结合能是不是计算吉布斯自由能的差值,焓值差要计算吗?我的想法是 ...

1. 先明确你要算结合能还是结合自由能。结合能可以视具体用途用零点能校正后的能量、内能或者焓,但不能用自由能;结合自由能只能用自由能。
2. 不是。焓和内能差RT,焓和Gibbs自由能差TS
wzkchem5 发表于 Post on 2023-5-28 16:15:22
本帖最后由 wzkchem5 于 2023-5-28 18:31 编辑
yjb 发表于 2023-5-27 15:30
老师,再请教一个问题,由于初始结构的优化和单点能计算的赝势用的不一样,还有就是由于大小核基组导致sp ...

有一种特殊情况下我感觉可能需要手动校正一下电子熵。假如你的体系不仅有稀土,还有其他未成对电子,而且未成对电子数目发生了改变,那么大核赝势的电子熵的变化会不对。比如你的体系有一个Gd(III),然后有一根共价键均裂:[Gd-A-B] -> [Gd-A·] + B·,并且假设Gd(III)和A·是反铁磁耦合。那么按大核赝势算,电子熵变等于Rln4;按小核赝势算,电子熵变等于Rln(7/4)。所以严格的处理应该需要手动计算大核赝势因为忽略f电子对电子熵引入的误差。不过我不知道有没有文献是真的这样做的。

如果反应前后只有一个稀土、反应不产生新的单电子、也没有单电子消失,那么容易发现大核赝势的电子熵虽然是错的,但是电子熵变是对的,此时不需要校正。
sobereva 发表于 Post on 2023-5-28 01:59:59
yjb 发表于 2023-5-27 11:31
老师,如果用上面的方法在求反应能垒的时候,得到一个粗略的结果,大核赝势结构优化下的热力学能量加上小核 ...

可行
这总比你只用小核赝势算的电子能量来讨论强
yjb 发表于 Post on 2023-5-27 22:30:31
wzkchem5 发表于 2023-5-27 15:02
那叫热力学校正量,不叫热力学能量
如果已经算了热力学校正量,尽量加热力学校正量,哪怕是很粗略级别下 ...

老师,再请教一个问题,由于初始结构的优化和单点能计算的赝势用的不一样,还有就是由于大小核基组导致spin也不同,这样会影响热力学矫正量和单点能相加之后的结果吗?
yjb 发表于 Post on 2023-5-27 21:51:29
wzkchem5 发表于 2023-5-27 15:02
那叫热力学校正量,不叫热力学能量
如果已经算了热力学校正量,尽量加热力学校正量,哪怕是很粗略级别下 ...

好的 谢谢老师
wzkchem5 发表于 Post on 2023-5-27 15:02:07
yjb 发表于 2023-5-27 04:31
老师,如果用上面的方法在求反应能垒的时候,得到一个粗略的结果,大核赝势结构优化下的热力学能量加上小核 ...

那叫热力学校正量,不叫热力学能量
如果已经算了热力学校正量,尽量加热力学校正量,哪怕是很粗略级别下计算得到的热力学校正量,也比不加好
yjb 发表于 Post on 2023-5-27 11:31:59
老师,如果用上面的方法在求反应能垒的时候,得到一个粗略的结果,大核赝势结构优化下的热力学能量加上小核赝势下的单点能这样可行吗?还是直接用小核赝势计算出来的单点能。
护卫天使 发表于 Post on 2022-10-20 10:51:57
护卫天使 发表于 2022-9-24 18:30
老师好,对于12-冠-4-五水合钬(III)和12-冠-4-五水合镧(III)两个配合物。想简单比较一下配合物稳定性,我的 ...

老师好,我想再请教一下,1、计算水中结合能是不是计算吉布斯自由能的差值,焓值差要计算吗?我的想法是每个结构带上隐式溶剂算单点,再用您开发的Shermo程序输出。
2、还有就是吉布斯自由能和焓值的热化学量是不是就差个RT?
sobereva 发表于 Post on 2022-10-8 18:30:28
护卫天使 发表于 2022-10-7 10:49
老师好,以您给的例子为例。在CDA和ETS-NOCV分析中,不想考虑水分子对配体或者对稀土的相互作用,能否在 ...

我不建议去掉水。因为配位水对Ho的电子结构状态影响极大,没了水可能更难以靠CDA/ETS-NOCV反映Ho和12-冠-4在实际配合物中的相互作用
护卫天使 发表于 Post on 2022-10-7 10:49:36
sobereva 发表于 2022-8-18 17:00
如果你对不同配体间的相互作用不感兴趣,可以把所有配体一起当个片段用Multiwfn做CDA和ETS-NOCV分析。

老师好,以您给的例子为例。在CDA和ETS-NOCV分析中,不想考虑水分子对配体或者对稀土的相互作用,能否在已经优化到极小值的结构中,去掉水分子,提取12-冠-4-Ho的结构,然后作单点能计算,此时把12-冠-4作为一个片段,Ho作为一个片段,再去讨论?
sobereva 发表于 Post on 2022-9-26 05:42:55
护卫天使 发表于 2022-9-25 10:35
感谢老师的回答。假如我探究两个片段分别是稀土和全部的配体,那么对于ETS-NOCV和CDA分析,我疑惑的是: ...

ETS-NOCV能直接给出轨道相互作用能,这是配体结合能的重要组成部分
护卫天使 发表于 Post on 2022-9-25 10:35:25
sobereva 发表于 2022-9-25 02:11
严格的判断指标是水中的结合自由能
去分析相互作用特征和电荷转移,推荐用Multiwfn做ETS-NOCV(http://s ...

感谢老师的回答。假如我探究两个片段分别是稀土和全部的配体,那么对于ETS-NOCV和CDA分析,我疑惑的是:轨道相互作用能和电荷转移的分析,最终能否根据相互作用能的大小和转移电荷的大小,反映两者稳定性的差异?我写的文章是实验加理论,理论部分就是想验证其稳定性。

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

GMT+8, 2024-11-26 00:52 , Processed in 0.476961 second(s), 26 queries , Gzip On.

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