计算化学公社

标题: gaussian计算离子间相互作用力求助 [打印本页]

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

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

作者
Author: zjxitcc 时间: 2021-12-2 17:39
本帖最后由 zjxitcc 于 2021-12-2 17:40 编辑

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

作者
Author: comestronger 时间: 2021-12-2 20:20

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

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

作者
Author: comestronger 时间: 2021-12-2 20:23

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

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

作者
Author: wzkchem5 时间: 2021-12-2 20:23

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

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

作者
Author: comestronger 时间: 2021-12-2 20:44

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

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

作者
Author: wzkchem5 时间: 2021-12-2 22:49

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

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

作者
Author: comestronger 时间: 2021-12-2 23:01

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

非常感谢您！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)