计算化学公社

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

[Gaussian/gview] 隐式溶剂化导致的钯团簇对称性下降应当如何考虑?

[复制链接 Copy URL]

69

帖子

0

威望

1250

eV
积分
1319

Level 4 (黑子)

本帖最后由 妙角不脆 于 2023-4-11 09:43 编辑

大家好,今天看到一篇文章(https://pubs.acs.org/doi/abs/10.1021/acsanm.2c00389)中看到这两个Pd团簇非常奇怪,不明白为什么如此简单的团簇为什么算出来的结果不是最高对称性,尤其是第一个三角形的团簇不是D3H点群而是C2V点群。第一感是心中嘲笑现在计算门槛太低,啥人都敢算,明显对称性不对的东西都敢发表。但出于严谨考虑自己复算了一遍,结果却大大出乎意料。计算软件是Gaussian16,使用的PBE0/Def2-TZVP的计算级别。发现气态下确实是D3H点群,但加了隐式溶剂化后却变成了C2V点群(相关输出文件已经附上)。隐式溶剂化导致的对称性下降令我大呼不解,想请教论坛里朋友这种情况究竟应该如何解决。是应当解释为溶剂模型的噪音还是确实溶剂模型会破坏对称性?还是有什么其他原因?恳请大家不吝赐教!
pd3.log (46.69 KB, 下载次数 Times of downloads: 2)


pd3water.log (54.34 KB, 下载次数 Times of downloads: 1)




127

帖子

0

威望

2765

eV
积分
2892

Level 5 (御坂)

2#
发表于 Post on 2023-4-11 09:51:56 | 只看该作者 Only view this author
这...你觉得真实世界的溶剂里,分子可以个个保持高对称性?

溶剂中的分子绝大部分都是C1,只有极少数有对称性,这才是合理的...

69

帖子

0

威望

1250

eV
积分
1319

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2023-4-11 10:25:46 | 只看该作者 Only view this author
fatpig 发表于 2023-4-11 09:51
这...你觉得真实世界的溶剂里,分子可以个个保持高对称性?

溶剂中的分子绝大部分都是C1,只有极少数有 ...

感谢您的回答。您说的对,计算当然不可能是对真实世界的完全写实描述。但假如给苯加隐式溶剂化就不会导致对称性的下降,一直是D6H。

150

帖子

0

威望

3919

eV
积分
4069

Level 6 (一方通行)

4#
发表于 Post on 2023-4-11 10:51:59 | 只看该作者 Only view this author
这个考虑是隐式溶剂模型增强了Jahn-Teller效应,使分子结构畸变。并不是对称性越高能量越低的

69

帖子

0

威望

1250

eV
积分
1319

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2023-4-11 11:28:16 | 只看该作者 Only view this author
RAL 发表于 2023-4-11 10:51
这个考虑是隐式溶剂模型增强了Jahn-Teller效应,使分子结构畸变。并不是对称性越高能量越低的

感谢您的回复。我学的Jahn-Teller效应都是局限在八面体场中的。之前从未听说过非八面体场中的Jahn-Teller效应,能否请您推荐针对非八面体场中的Jahn-Teller效应一些书籍/文章章节呢?

6万

帖子

99

威望

6万

eV
积分
125137

管理员

公社社长

6#
发表于 Post on 2023-4-11 12:00:45 | 只看该作者 Only view this author
这点偏离D3h差异程度可以忽略不计,照常当D3h就可以。隐式溶剂模型的积分格点分布不严格满足体系对称性,足矣导致这种程度的差异
北京科音自然科学研究中心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

218

帖子

2

威望

2033

eV
积分
2291

Level 5 (御坂)

7#
发表于 Post on 2023-4-11 12:46:13 | 只看该作者 Only view this author
也可以试一试修改溶剂空腔,比如写scrf=read,后面写surface=SAS,有可能能够保持高对称性

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

8#
发表于 Post on 2023-4-11 17:14:55 | 只看该作者 Only view this author
妙角不脆 发表于 2023-4-11 04:28
感谢您的回复。我学的Jahn-Teller效应都是局限在八面体场中的。之前从未听说过非八面体场中的Jahn-Teller ...

I. B. Bersuker, The Jahn-Teller Effect and Vibronic Interactions in Modern Chemistry。这本书可以搜到电子版
此外我觉得这个体系更可能是二阶Jahn-Teller效应。如果是一阶Jahn-Teller效应的话,无法解释为什么真空里完全不畸变。二阶Jahn-Teller效应方面建议看Proc. Nat. Acad. Sci. USA Vol. 72, No. 6, pp. 2104-2106, 1975
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
本团队长期招收研究生,有意者可私信联系

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

9#
发表于 Post on 2023-4-11 20:05:50 | 只看该作者 Only view this author
fatpig 发表于 2023-4-11 02:51
这...你觉得真实世界的溶剂里,分子可以个个保持高对称性?

溶剂中的分子绝大部分都是C1,只有极少数有 ...

溶剂里的分子有对称性的概率,未必比气相里同等大小的分子有对称性的概率小吧?
一般实验上研究的气相分子倾向于比较小(大分子蒸汽压太低,不好研究),而溶剂里的分子倾向于比较大,分子越大有对称性的概率越低。所以搞不好“在溶剂里”和“没有对称性”并非直接的因果关系
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
本团队长期招收研究生,有意者可私信联系

127

帖子

0

威望

2765

eV
积分
2892

Level 5 (御坂)

10#
发表于 Post on 2023-4-12 18:16:10 | 只看该作者 Only view this author
wzkchem5 发表于 2023-4-11 20:05
溶剂里的分子有对称性的概率,未必比气相里同等大小的分子有对称性的概率小吧?
一般实验上研究的气相分 ...

个人认为,自然界里或者实验中的溶剂里,因为各种相互作用的存在,根本就不会有高对称性,能降到C1就一定会降到C1。
哪怕是计算,例如算水团簇,怎么可能各个H2O都是C2v,肯定会降到最低对称性的Cs.
又比如前面楼主提过的苯,单算一个就算加了隐式模型,也是D6h,这很正常。但只要加多几个显式的苯分子,还能保持D6h吗,全部都是C1才是最有可能出现的情况。

而,气相或者真空里,同样的分子,与周围环境其他分子的碰撞概率或者说相互作用,必然是没溶剂中那么强的,更易保持原本的高对称性。

所以,结论就是,在计算中加了溶剂模型,对称性相比真空降低,反而更符合现实世界中的情况。

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

11#
发表于 Post on 2023-4-12 19:07:11 | 只看该作者 Only view this author
fatpig 发表于 2023-4-12 11:16
个人认为,自然界里或者实验中的溶剂里,因为各种相互作用的存在,根本就不会有高对称性,能降到C1就一定 ...

如果是显式溶剂计算,且只取一个snapshot,那溶质分子确实几乎必然是最低对称性。
但是隐式溶剂模型模拟的是显式溶剂计算对溶剂分子做无穷长时间MD取统计平均的结果,而不是显式溶剂计算的一个snapshot。所以显式溶剂计算溶质分子几乎总是最低对称性,不代表隐式溶剂模型也应该优化得到低对称性。
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
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 15:19 , Processed in 0.173303 second(s), 23 queries , Gzip On.

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