求助锕系金属体系在隐式溶剂模型下的计算问题

yuanliu18

大家好，我想计算锕系金属，比如说ThCl4的结构优化和溶解自由能计算，但是出现一个问题，我在用Tpssh（用这个方法是根据这篇sobereva老师贴子：http://sobereva.com/272）和stuttgart小核赝势，在真空状态和PCM（ethanol）的环境下，都很快的优化计算完成，但是我根据sobereva老师这两篇贴子（http://sobereva.com/327，http://sobereva.com/613）中，将得到的结构在SMD和SM18的方法下作单点能计算，发现scf不收敛，然后我就根据这个贴子(http://sobereva.com/61)去解决,我并没有尝试所以的方法，但我尝试了用PCM的波函数作为初猜，改变fermi能级，增加scf次数，发现都没有用，最后用xqc的方法得到收敛，但我不确定这个能量是不是可以用。所以我又查看了这篇关于隐式模型溶剂的贴子(http://sobereva.com/42) 和SMD的原文，发现PCM和SMD的区别在于对deltaG(nonpolar)的定义，SMD对于这部分是对实验的拟合，然后我查看它的训练集里面没有锕系元素。所以我想是不是这个原因导致我的SMD计算不收敛。如果是这个原因，那再用SMD计算我的这个体系是否还合理？而且我做了个简单谷歌文献搜索，发现最近的关于锕系的隐式计算用的都是cpcm或者pcm。基于此，是不是用PCM的方法就可以了，或者是不是还有别的方法？


最后，谢谢大家的解惑！！！

wzkchem5

ThCl4放在乙醇里，乙醇会和Th配位吧？如果是的话，配位的乙醇必须显式描述。第一配位层以外的溶剂用PCM、CPCM还是SMD倒是次要的。
SMD的训练集里没有锕系元素，可能会导致结果不如非锕系体系准，但不会直接导致不收敛。此外PCM、CPCM更是没有用锕系元素来训练，所以这也不算SMD相对于其他溶剂模型的缺点

hebrewsnabla

wzkchem5 发表于 2024-1-20 17:11
ThCl4放在乙醇里，乙醇会和Th配位吧？如果是的话，配位的乙醇必须显式描述。第一配位层以外的溶剂用PCM、CP ...

PCM和SMD相比，基本没有和溶质有关的经验参数，所以只有SMD需要拟合这些参数。SMD的这个拟合问题挺大的，不只是没有拟合锕系，连一个金属都没有……或者说并不是作者的问题，只是后人扩大了SMD的使用范围。

当然，第一配位层用显式的话这就不重要了。

sobereva

(IEF)PCM和SMD的区别不仅在于后者额外明确定义了非极性部分，还在于溶质孔洞的构造，SMD自己有一套原子半径用来构造孔洞。而非极性部分并不会影响SCF收敛
用IEFPCM也可以要求Gaussian计算非极性部分的贡献，只不过对于SMD参数化时考虑的那些体系（都不含过渡金属），肯定没有SMD好。而对于镧系金属配合物，则不一定

wzkchem5

hebrewsnabla 发表于 2024-1-20 15:11
PCM和SMD相比，基本没有和溶质有关的经验参数，所以只有SMD需要拟合这些参数。SMD的这个拟合问题挺大的， ...

没错，所以SMD算金属，总体上来说肯定不如SMD算非金属靠谱。但是SMD算金属和PCM算金属哪个好就不好说了，我个人还是倾向于用SMD，起码当金属被配体包得足够严实时SMD肯定是好的

yuanliu18

谢谢各位老师的帮忙，我看能否找到实验数据去计算一下哪一种方法更好。