高斯计算1个分子末端基团的原子电荷与2个分子时同一末端基团的电荷相差较大问题

Yaohaohao · 发表于 Post on 2023-7-11 15:38:58

本帖最后由 Yaohaohao 于 2023-7-13 17:24 编辑

请教各位老师，我用高斯将一个分子进行结构优化（文件名为1），体系电荷和自旋多重度设置为-1 3，方法是def2tzvp m06，优化结束后用Multiwfn计算了Hirshfeld charge，将末端基团（原子序号56-75）的charge加和为0.0704。然后将同样2个分子放在一起进行结构优化（文件名为2），体系电荷和自旋多重度设置为-2 1，方法是def2svp m06，溶剂模型用了surface=sas。同样的末端基团（原子序号56-75）的charge和是-0.1155。我想问是因为两个体系设置的自旋多重度以及基组不同导致相同基团的charge相差较大吗。或者说当两个体系所有设置条件相同时，同样的末端基团（原子序号56-75）的charge相差不大才是正常的？请各位老师不吝赐教，感谢。

1.zip (8.12 MB, 下载次数 Times of downloads: 6)
按sob老师回答我将二聚体体系进行优化（因之前报错这次计算没有加em==GD3BJ)，显示不收敛

dimer.gjf (12.66 KB, 下载次数 Times of downloads: 2)

dimer.log (327.41 KB, 下载次数 Times of downloads: 1) ，请老师帮忙看一下如何解决谢谢。

sobereva · 发表于 Post on 2023-7-11 16:09:50

没具体体系信息没法回答，都不知道你算的什么的东西、基团长什么样
而且不要说正数还是负数，必须说具体数值。诸如0.0001和-0.0001根本没本质区别

Yaohaohao · 发表于 Post on 2023-7-11 16:59:23

sobereva 发表于 2023-7-11 16:09
没具体体系信息没法回答，都不知道你算的什么的东西、基团长什么样
而且不要说正数还是负数，必须说具体数 ...

社长好，已经补充了具体信息，请再看一下，谢谢

sobereva · 发表于 Post on 2023-7-12 09:06:44

你先把所有能统一的地方都统一了再比。对于一个没有可比性的数据没法去讨论化学问题

基组必须统一。而且算单体是三重态，为什么算二聚体就成了单重态？
当前你既然直接写了def2-SVP或def2-TZVP，又没写genecp，输入文件后面自定义基组/赝势完全没有任何意义

你这种体系opt freq我更推荐用PBE0-D3(BJ)。基组6-311G*足矣，def2-TZVP太浪费

Yaohaohao · 发表于 Post on 2023-7-12 10:01:11

sobereva 发表于 2023-7-12 09:06
你先把所有能统一的地方都统一了再比。对于一个没有可比性的数据没法去讨论化学问题

基组必须统一。而且 ...

好的社长，我再统一下设置，谢谢

Yaohaohao · 发表于 Post on 2023-7-13 17:26:26

sobereva 发表于 2023-7-12 09:06
你先把所有能统一的地方都统一了再比。对于一个没有可比性的数据没法去讨论化学问题

基组必须统一。而且 ...

老师好，我的二聚体体系计算有问题，已更新了附件，请帮忙看一下，谢谢

sobereva · 发表于 Post on 2023-7-14 11:16:59

关键词里绝对不应当写PBE0，好好看http://bbs.keinsci.com/thread-13660-1-1.html

Yaohaohao · 发表于 Post on 2023-7-14 16:06:06

sobereva 发表于 2023-7-14 11:16
关键词里绝对不应当写PBE0，好好看http://bbs.keinsci.com/thread-13660-1-1.html

学习了，谢谢社长

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 高斯计算1个分子末端基团的原子电荷与2个分子时同一末端基团的电荷相差较大问题