| 测试过，这就是能量最低的 |
| 有对铁-EDTA的多重度测试过吗，能量最低的是多重度为6吗？ |
小七丫 发表于 2021-10-21 09:41 不要说"L9999的错误“，导致L9999报错的问题有一大堆，不能排除是因为另外的原因报错的。 IEFPCM严格好于CPCM，CPCM在那些不支持IEFPCM的程序里可以用，但是既然高斯支持IEFPCM，你用CPCM容易被审稿人质疑，因为两者计算量没多大区别，没有任何理由用一个差一些的模型，哪怕差得未必很多。 试一下以opt=loose收敛的结构为基础，继续优化，优化时用比较小的步长 |
| 嗯嗯，我去试试，谢谢老师 |
小七丫 发表于 2021-10-20 08:51 你加的是什么溶剂模型？如果加的是SMD，换成IEFPCM试试（所有物种都要换，即使那些SMD能收敛的也是如此，同样是为了误差抵消的考虑），SMD有数值不稳定性问题，结构优化比较难收敛。单点能可以仍然用SMD。此外注意区分SCF震荡和结构优化震荡，只有对于前者，“读取没带溶剂模型的波函数作为初猜”才有用，对于后者没用；而opt=(MaxStep=5, NoTrustUpdate, GDIIS)只对后者有用，对前者没用。 “除了Fe(H2O)6是用的气相优化的结构算溶解自由能，其余都是溶剂模型优化的结构算的，但是和实验值还是差很多”差得不多才怪。如果理论级别不统一，根本就不要尝试去和实验比，哪怕是因为不收敛，也应该只问收敛问题，而不要因此就去拿理论级别不统一的结果去减。理论级别任何细微的不统一，都会导致能量相减的时候有极大的误差， |
小七丫 发表于 2021-10-17 11:13 （1）统一了泛函和基组，不代表你统一了校正因子，校正因子是否统一取决于你在Shermo里是否对于所有物种都输入了同一个校正因子。假如你是纯有机体系，用的是PBE0/6-311G*，没有Fe原子，那么应当始终统一用PBE0/6-311G*的校正因子，结果比不用校正因子靠谱（但改善可能很小，看不出来）。但因为你现在有Fe原子，Fe原子用的是SDD，所以用PBE0/6-311G*的校正因子并不完全严格，不保证比不用校正因子靠谱，但也不保证比不用校正因子更差，所以此时有两种选项，一个是对所有物种都用PBE0/6-311G*的校正因子，一个是对所有物种都不用校正因子，两个做法都对，因为两个做法都有可能是较好的做法，不确定哪个更好。但混着用一定是错的，因为无法实现误差抵消。 （2）基频校正因子是用来算振动频率的，比如你算红外、拉曼，需要用那个校正因子。算热力学数据要用针对热力学数据的校正因子。 （3）必须要么对所有物种都加溶剂模型，要么都对所有物种都不加。溶剂模型也是理论级别的一部分。我不明白你为什么总是不懂误差抵消的概念，因为你算自由能变化是用非常大的数减去非常大的数，所以必须所有物种都用同样的理论级别才能实现误差抵消，宁可对某些物种用对这个物种来说很不适合的理论级别，也不能让两个物种的理论级别不一致。建议打开你的所有高斯输入文件，一个字一个字对比，有任何差别都统一，除非你明确知道这个差别是可以接受的（比如我之前说的，对于没有Fe的体系可以不写SDD），否则一律统一，一个关键词乃至一个字母都不要差。 （4）M052X/6-31G*只能用来算溶解自由能，严格按照sob老师的方法来，不要自己引申。不能用这个级别优化结构、算频率，更不能算单点能 |
|
本帖最后由 小七丫 于 2021-10-17 18:23 编辑 谢谢老师的解答，我气相中的结构优化现在都统一成了PBE0(PBE1PBE)/6-311G*，对Fe用的SDD,这应该算是校正因子统一吧，您的意思是可以不用校正因子是吗，还有就是关于基准（Fundamentals）频率校正因子,它是用在哪里的呀，我没在Shermo的ini文件里看到。关于溶剂里的情况，现在只有EDTA4价负离子加了，我把其它的也加上试试，我直接在M052X/6-31G*加溶剂的情况下做结构优化和频率计算，得到的能量就等于单点能是吗 |
小七丫 发表于 2021-10-17 07:27 校正因子也是理论级别的一部分，所以也是宁可用同一个差的校正因子，也不能用不同的好的校正因子，所以你可以以有的物种用了SDD为由，统一不加校正因子，也可以以只有一个原子用了SDD、占总体系比例很小为由，统一用PBE0/6-311+G*的校正因子，但是不能只对没有Fe的分子用、对有Fe的分子不用。此外，在校正因子完全统一的前提下，校正因子对结果的影响不是太大，估计最多也就能差出1kcal/mol。 你的结构优化是统一成了哪个理论级别，加没加溶剂？EDTA四价负离子在气相里很可能不稳定，所以优化时必须加溶剂，又根据理论级别必须统一的要求，所有分子优化时都必须加溶剂 |
小七丫 发表于 2021-10-15 01:56 对，优化级别和单点能计算级别必须严格一致。唯一允许不一致的地方就是，如果包含铁的物种你对铁用了SDD，那么不包含铁的物种不用写SDD。不能分别对不同物种用不同的理论级别，宁可用同一个差的理论级别，也不能对每个物种分别用不同的好的理论级别 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-24 13:08 , Processed in 0.196024 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites