计算化学公社

标题: 如何使用Gaussian计算Pauli repulsion？ [打印本页]

作者
Author: wrtgcn 时间: 2021-7-5 19:38
标题: 如何使用Gaussian计算Pauli repulsion？
最近看了一篇文献（The Pauli Repulsion-Lowering Concept in Catalysis, Acc. Chem. Res. 2021, 54, 1972−1981），作者提出相比LUMO-lowering理论，Pauli repulsion-lowering理论可能更好的解释一些催化过程。具体来说，作者使用activation strain model将电子能量拆分为activation strain modelthe strain (ΔEstrain) that results from the distortion of the individual reactants以及the interaction (ΔEint) between the deformed reactants along the reaction coordinate defined as the projection of the intrinsic reaction coordinate (IRC) on a critical geometry parameter that directly defines the progress of the studied reaction；之后再用canonical energy decomposition analysis (EDA) method 将ΔEint拆分为electrostatic interaction (ΔVelstat), the Pauli repulsion (ΔEPauli) between closed-shell orbitals which is responsible for steric repulsion, and the stabilizing orbital attractions (ΔEoi) that account, among others, for HOMO−LUMO interactions.

想请问大家，这个Pauli repulsion (ΔEPauli) 的概念能否通过Gaussian计算，又应该如何实现呢？谢谢。

作者
Author: sobereva 时间: 2021-7-5 19:41
Gaussian没法算
用GAMESS-US做Morokuma或LMOEDA能量分解得到

作者
Author: wrtgcn 时间: 2021-7-5 20:04

sobereva 发表于 2021-7-5 19:41
Gaussian没法算
用GAMESS-US做Morokuma或LMOEDA能量分解得到

好的，谢谢Sob老师~

作者
Author: lijiayisjtu 时间: 2021-7-5 21:12
distortion interaction analysis分析完如果要能量分解的化也可以用SAPT做找到TS之后沿着IRC 把每个点的结构都进行SAPT能量分解然后拟合曲线也可以得到类似的曲线
F. Matthias Bickelhaupt和Kendall N. Houk用这个发过很多文章，

作者
Author: wrtgcn 时间: 2021-7-5 22:18

lijiayisjtu 发表于 2021-7-5 21:12
distortion interaction analysis分析完如果要能量分解的化也可以用SAPT做找到TS之后沿着IRC 把每个点的 ...

好的，谢谢~
多问一句，这些概念中有比较近似Pauli repulsion的吗？

作者
Author: sobereva 时间: 2021-7-6 07:50
LMOEDA里有repulsion项，SAPT里有exchange项，都和这本质对应
SAPT并不适合片段间有化学键程度作用的情况，不要盲目效仿一些文章

作者
Author: zjxitcc 时间: 2021-7-6 10:22
本帖最后由 zjxitcc 于 2021-7-6 10:43 编辑

GAMESS里的Morokuma-EDA只能用于RHF波函数（限制性DFT也不行），因此如今没有实际用途。
GAMESS里的LMO-EDA倒是可以做点计算，算有机体系没啥问题，自由基、过渡金属很容易碰到SCF不收敛。
XEDA程序（苏培峰修改版的GAMESS，可发邮件supi[at]xmu.edu.cn向其索要）里可以做现如今的GKS-EDA，且支持读入轨道，不会有SCF不收敛问题。如果不知道怎么写GKS-EDA输入文件，可以用MOKIT自动生成（https://gitlab.com/jxzou/mokit），在PDF手册5.3.3章节有例子，其自动调用高斯做计算然后将关键词、坐标、基组和轨道写入GKS-EDA输入文件，这样GKS-EDA计算一般只需若干分钟（小体系只需要十几秒）。

Pauli-repulsion在GKS-EDA里是exchange term + repulsion term，这两项在GKS-EDA算完后有输出，自己加一下即可。

作者
Author: wrtgcn 时间: 2021-7-6 12:08

zjxitcc 发表于 2021-7-6 10:22
GAMESS里的Morokuma-EDA只能用于RHF波函数（限制性DFT也不行），因此如今没有实际用途。
GAMESS里的LMO-ED ...

好的，非常感谢~

作者
Author: naonao5205 时间: 2021-12-2 19:43

lijiayisjtu 发表于 2021-7-5 21:12
distortion interaction analysis分析完如果要能量分解的化也可以用SAPT做找到TS之后沿着IRC 把每个点的 ...

您好，请问可以告诉下这篇文章的出处吗？

作者
Author: LYX-TLiu 时间: 2022-10-21 22:18

naonao5205 发表于 2021-12-2 19:43
您好，请问可以告诉下这篇文章的出处吗？

https://pubs.rsc.org/en/content/articlelanding/2021/SC/D0SC06628A

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