gaussian计算离子间相互作用力求助

comestronger · 发表于 Post on 2021-12-2 14:32:18

gaussian计算求助！我想计算一个中性分子吸附H质子后带一个正电荷的结构与金属阳离子之间的作用力，但一直收敛不了，优化过程中，金属阳离子与带正电荷分子之间距离会一直增大，但就是不会收敛！大家有什么好的建议吗？在线求助！

wzkchem5 · 发表于 Post on 2021-12-2 16:43:04

不收敛是正常的。你做的这个计算不是算作用力，而是结构优化，也就是调整结构直至力为0，而因为库伦排斥力是长程的，无论距离增大到多少力都不为零，所以永远不收敛。
还是要明确你要算的是什么，你现在这么算，要么不收敛，只要收敛了作用力一定是零，算不出非零的数

zjxitcc · 发表于 Post on 2021-12-2 17:39:16

本帖最后由 zjxitcc 于 2021-12-2 17:40 编辑

基本概念上有错误，正电荷之间当然是不断排斥远离的，不用算也知道。看sob老师如下博文
《计算化学中的一些常见不良写法和用词》http://sobereva.com/298 第12点
《利用Gaussian或ORCA程序把分子结构拉直的几种方法》http://sobereva.com/594 第3.3 方法3

comestronger · 发表于 Post on 2021-12-2 20:20:46

wzkchem5 发表于 2021-12-2 16:43
不收敛是正常的。你做的这个计算不是算作用力，而是结构优化，也就是调整结构直至力为0，而因为库伦排斥力 ...

谢谢您，想算吡啶分子在酸性溶液与阳离子的结合能，因为吡啶分子在酸性溶液中会结合一个H质子，所以也就是计算带正电荷的H结合后的吡啶分子与阳离子的结合能。

comestronger · 发表于 Post on 2021-12-2 20:23:01

zjxitcc 发表于 2021-12-2 17:39
基本概念上有错误，正电荷之间当然是不断排斥远离的，不用算也知道。看sob老师如下博文
《计算化学中的一 ...

是的，谢谢您，但这种作用力怎么通过能量来比较呢？即计算目的：想算吡啶分子在酸性溶液与阳离子的结合能，因为吡啶分子在酸性溶液中会结合一个H质子，所以也就是计算带正电荷的H结合后的吡啶分子与阳离子的结合能。

wzkchem5 · 发表于 Post on 2021-12-2 20:23:49

comestronger 发表于 2021-12-2 13:20
谢谢您，想算吡啶分子在酸性溶液与阳离子的结合能，因为吡啶分子在酸性溶液中会结合一个H质子，所以也就 ...

如果两者一直互相远离，表示两者根本不会结合，所以结合能没有定义，也可以说结合能等于0。
并不是任意两个分子的结合能都有定义

comestronger · 发表于 Post on 2021-12-2 20:44:52

wzkchem5 发表于 2021-12-2 20:23
如果两者一直互相远离，表示两者根本不会结合，所以结合能没有定义，也可以说结合能等于0。
并不是任意 ...

恩恩，明白了，谢谢您！因为是想通过理论计算来解释实验现象，所以对于阳离子之间的作用力，直接表述为排斥力，也是可以的吧？如果是阴离子和阳离子的结合，这种结合能可以计算出来，作为说明实验结果是合理的了。

wzkchem5 · 发表于 Post on 2021-12-2 22:49:05

comestronger 发表于 2021-12-2 13:44
恩恩，明白了，谢谢您！因为是想通过理论计算来解释实验现象，所以对于阳离子之间的作用力，直接表述为排 ...

可以，就说不管离多远都在排斥，完全没有结合的倾向

comestronger · 发表于 Post on 2021-12-2 23:01:07

wzkchem5 发表于 2021-12-2 22:49
可以，就说不管离多远都在排斥，完全没有结合的倾向

非常感谢您！

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

[新手求助] gaussian计算离子间相互作用力求助