计算化学公社

标题: 请问用G4方法计算电离能，分子中含碘原子，基组该如何选择？高斯文件咋输入呢？ [打印本页]

作者
Author: flybird52888 时间: 2023-5-16 09:25
标题: 请问用G4方法计算电离能，分子中含碘原子，基组该如何选择？高斯文件咋输入呢？
本帖最后由 flybird52888 于 2023-5-17 09:18 编辑

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

作者
Author: fatpig 时间: 2023-5-16 10:19
找G4的原始文献，一步步拆解了手动分步算。
找几个小分子试一下，保证分步后每一步能量加起来的总能量一定和原G4最后给出的能量一致。

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

作者
Author: flybird52888 时间: 2023-5-16 10:46

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

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

作者
Author: wxsxjj 时间: 2023-5-16 12:45

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

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

作者
Author: mfdsrax2 时间: 2023-5-16 15:13

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

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

作者
Author: sobereva 时间: 2023-5-16 22:07
根本没必要用G4
并且说清楚是垂直电离还是绝热电离能

作者
Author: flybird52888 时间: 2023-5-17 09:10
本帖最后由 flybird52888 于 2023-5-17 09:16 编辑

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

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

作者
Author: flybird52888 时间: 2023-5-17 09:13

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

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

作者
Author: sobereva 时间: 2023-5-17 10:51

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

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

作者
Author: wxsxjj 时间: 2023-5-17 11:41

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

直接CCSD(T)就完事了

作者
Author: flybird52888 时间: 2023-5-17 13:33

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

好的，谢谢！

作者
Author: flybird52888 时间: 2023-5-17 13:33

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

谢谢！

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