din5g 发表于 2021-9-27 17:16 嗯，可以这么理解 |
本帖最后由 din5g 于 2021-9-28 00:19 编辑 wzkchem5 发表于 2021-9-27 22:28 多谢wzkchem5老师提醒，片段处理总是出现Bad data into FinFrg问题，最后我用guess=mix跑完后成功得到了超氧，两个氧原子Mulliken电荷为-0.09和-0.3，可以定性为得到电子，再根据自旋布居数：-0.62和-0.54，可以判断总共得到一个电子左右。这样就获得了超氧阴离子。不知道我的解释对不对 2. 我看H原子上Mulliken电荷为0.34，电荷布居为0.019，这只能说明H原子上没有单电子，不能说明为质子，不过我看全局的单电子分布里，单电子只分布在苄基和O2上，说明催化剂上没有单电子分布，也就说，H原子那里为质子。我的理解对吗？ |
din5g 发表于 2021-9-27 15:01 把O-O键拉长一些，再试试。如果不管怎么尝试都不行，可能意味着超氧负离子只在激发态势能面上存在，你需要做TDDFT计算 |
本帖最后由 din5g 于 2021-9-27 22:22 编辑 wzkchem5 发表于 2021-9-27 21:01 wzkchem5老师的意思是这里得不到一个电子吗？那怎么设置才能让他成为超氧负离子呢另外，突然想到，不知道跑激发态计算可不可以描述甲苯的SP3 C-H键断裂过程以及质子和超氧生成OOH自由基的过程 |
din5g 发表于 2021-9-27 13:49 如果此处氧气是超氧负离子，那就是的。但是这里氧气未必是超氧负离子 |
wzkchem5 发表于 2021-9-27 16:16 好的，多谢wzkchem5老师，我用guess=mix试了一下，找到了破缺态，但我看到两个氧原子（O2）上带相反电荷，相反的自旋单电子，数值很小，小数点第三位。请问wzkchem5老师，正确的分布是不是两个氧原子都带负电荷，并且自旋平行？ |
din5g 发表于 2021-9-27 02:05 1. 对，同一个片段里的单电子必须要么都是alpha，要么都是beta。如果不满足这个条件，结果倒也未必是错误的，但是正确的概率要低很多，因为片段初猜的一个很重要的意义就是把不同自旋划分到不同的片段上 2. 对，自旋布居数、SOMO图、自旋密度图等等 3. 什么叫“所有片段全变成 0 1了”？你从哪里看的？必须用文本编辑器打开高斯输出文件看，不能从gaussview里看 |
本帖最后由 din5g 于 2021-9-27 15:56 编辑 wzkchem5 发表于 2021-9-27 00:29 好的，多谢wzkchem5老师，我这里有两个疑问： 1.“片段初猜计算里，每个片段里的所有单电子自旋必须彼此平行”，我不太理解这句话的意思，因为我在GaussView里设置一个片段自旋时，他只能设置这个片段里的单电子是a还是B，也就是说一个片段里的单电子在设置里不存在一个a和一个B。 2.怎样看单点计算后总体系波函数合不合理呢？是看单电子的分布吗？ 3.另外，我算完单点发现所有片段全变成 0 1了，也就说没找到破缺态，我试了guess=mix后，发现也不行，不知道我的问题在哪 |
|
如果怀疑片段2可能是自旋极化单重态，那就应该把片段2再分成两个片段，一个包含电子，一个包含空穴。片段初猜计算里，每个片段里的所有单电子自旋必须彼此平行。 另外片段初猜建议先做个单点计算，确认收敛的总体系波函数合理后，再另做过渡态搜索，读这个单点任务的波函数。因为你一开始设的自旋多重度，SCF收敛后不一定能保持，所以这样分两步做，容易及时发现问题 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-20 23:40 , Processed in 0.180183 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites