请问用G4方法计算电离能，分子中含碘原子，基组该如何选择？高斯文件咋输入呢？

flybird52888 · 发表于 Post on 2023-5-16 09:25:31

本帖最后由 flybird52888 于 2023-5-17 09:18 编辑

如题，请问用G4方法计算绝热电离能，分子中含碘原子（体系较小，就一个六元环），基组该如何选择？高斯输入文件咋写呢？

flybird52888 · 发表于 Post on 2023-5-17 13:33:58

wxsxjj 发表于 2023-5-17 11:41
直接CCSD(T)就完事了

谢谢！

flybird52888 · 发表于 Post on 2023-5-17 13:33:00

sobereva 发表于 2023-5-17 10:51
以标准做法计算两个态的U(0)求差就完了
B3LYP/def2-TZVP做opt freq得到结构和ZPE，CCSD(T)/def2-QZVP（ ...

好的，谢谢！

wxsxjj · 发表于 Post on 2023-5-17 11:41:55

flybird52888 发表于 2023-5-17 09:13
那还有别的方法可以算G4吗？求告知！谢谢！

直接CCSD(T)就完事了

sobereva · 发表于 Post on 2023-5-17 10:51:20

flybird52888 发表于 2023-5-17 09:10
需要算绝热电离能，那用啥方法比较好呢？您有啥推荐的方法吗？我这个分子就一个六元环，体系比较小

以标准做法计算两个态的U(0)求差就完了
B3LYP/def2-TZVP做opt freq得到结构和ZPE，CCSD(T)/def2-QZVP（算不动就def2-TZVPP）算电子能量，二者相加得到U(0)。对阳离子和中性状态都这么算，然后求差

flybird52888 · 发表于 Post on 2023-5-17 09:13:42

wxsxjj 发表于 2023-5-16 12:45
如果你是用Gaussian计算G4，不能设定基组。

那还有别的方法可以算G4吗？求告知！谢谢！

flybird52888 · 发表于 Post on 2023-5-17 09:10:13

本帖最后由 flybird52888 于 2023-5-17 09:16 编辑

sobereva 发表于 2023-5-16 22:07
根本没必要用G4
并且说清楚是垂直电离还是绝热电离能

需要算绝热电离能，那用啥方法比较好呢？您有啥推荐的方法吗？我这个分子就一个六元环，体系比较小

sobereva · 发表于 Post on 2023-5-16 22:07:44

根本没必要用G4
并且说清楚是垂直电离还是绝热电离能

mfdsrax2 · 发表于 Post on 2023-5-16 15:13:51

flybird52888 发表于 2023-5-16 10:46
哦哦，好的，谢谢回复！是一个方法！现在问题是我怎么输入基组呢？在输入文件第三行输入g4/genecp, 输入 ...

详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

wxsxjj · 发表于 Post on 2023-5-16 12:45:23

flybird52888 发表于 2023-5-16 10:46
哦哦，好的，谢谢回复！是一个方法！现在问题是我怎么输入基组呢？在输入文件第三行输入g4/genecp, 输入 ...

如果你是用Gaussian计算G4，不能设定基组。

flybird52888 · 发表于 Post on 2023-5-16 10:46:03

fatpig 发表于 2023-5-16 10:19
找G4的原始文献，一步步拆解了手动分步算。
找几个小分子试一下，保证分步后每一步能量加起来的总能量一定 ...

哦哦，好的，谢谢回复！是一个方法！现在问题是我怎么输入基组呢？在输入文件第三行输入g4/genecp, 输入原子坐标以后I用SDD基组，那轻元素原子的基组咋输入呢？

fatpig · 发表于 Post on 2023-5-16 10:19:42

找G4的原始文献，一步步拆解了手动分步算。
找几个小分子试一下，保证分步后每一步能量加起来的总能量一定和原G4最后给出的能量一致。

此时换成你要算的分子，在每一步中把 I 的基组换成相对应的级别，分步算，加起来...

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

[Gaussian/gview] 请问用G4方法计算电离能，分子中含碘原子，基组该如何选择？高斯文件咋输入呢？