简单量子化学问题答疑专帖

toto · 发表于 Post on 2023-2-17 18:38:31

本帖最后由 toto 于 2023-2-17 18:56 编辑

各位好，请教一下，g16同一输入文件里使用dft做精度高一点的opt+freq，关键词除了opt=tight, 还学要添加什么嘛？多谢

NGC626 · 发表于 Post on 2023-2-17 19:05:37

wzkchem5 发表于 2023-2-14 22:13
这是化学吸附啊，吸附能都到6eV的量级了。
建议用orca的embedded cluster功能（https://pubs.acs.org/do ...

老师您好，请问这个结构用RI-B3LYP-D3(BJ)/def2-TZVP(-f)做几何优化是否合适呢？如果用高斯的话用哪一个理论方法和基组优化合适呢？

wzkchem5 · 发表于 Post on 2023-2-17 20:08:45

NGC626 发表于 2023-2-17 12:05
老师您好，请问这个结构用RI-B3LYP-D3(BJ)/def2-TZVP(-f)做几何优化是否合适呢？如果用高斯的话用哪一个 ...

不太好，一般金属团簇用纯泛函比较合适。考虑到碳原子受self-interaction error影响比较大，对于这个体系也可能需要用杂化泛函，但是金属体系用杂化泛函的话PBE0、TPSSh等一般比B3LYP合适。参见http://sobereva.com/272。保险起见可以查一下过渡金属团簇表面碳原子吸附用什么泛函算比较准。
高斯里的理论级别可以和orca一致。
此外注意波函数稳定性问题

wzkchem5 · 发表于 Post on 2023-2-17 20:10:13

先说清楚你这里的精度指的是什么精度，以及提高精度的目的是什么（消虚频？保持体系对称性/平面性？还是其他目的？）。

sobereva · 发表于 Post on 2023-2-18 06:05:49

toto 发表于 2023-2-17 18:38
各位好，请教一下，g16同一输入文件里使用dft做精度高一点的opt+freq，关键词除了opt=tight, 还学要添加什 ...

精度指代不明
理论方法的精度、基组的精度、SCF收敛的精度、几何优化收敛的精度、积分格点精度等等都叫精度
说不清楚就说明实际目的

sobereva · 发表于 Post on 2023-2-18 06:08:55

hzfish 发表于 2023-2-16 22:43
气相中联二苯斜式结构最稳定，重叠式结构次之，交叉式结构能量高。
而联二BN苯类似物只有两种构象，重叠稳 ...

做二面角旋转的扫描，取一些关键结构，用Multiwfn做ETS-NOCV考察片段间轨道相互作用对比讨论

一个阿豪 · 发表于 Post on 2023-2-18 11:05:37

配置文件见图一，

图二为报错信息

调试多次无果，希望各位老师能够指点

量化新人 · 发表于 Post on 2023-2-18 11:28:15

各位老师好，我在寻找一个过渡态，今日优化出来的结果没有虚频，然而查看输出文件是这样的，请问各位老师我该怎么办？

sobereva · 发表于 Post on 2023-2-18 12:07:40

一个阿豪发表于 2023-2-18 11:05
配置文件见图一，
图二为报错信息

格式不对

不熟悉手写QST2输入文件就用GaussView的calculation setup界面生成，确保格式不至于出错

sobereva · 发表于 Post on 2023-2-18 12:08:41

量化新人发表于 2023-2-18 11:28
各位老师好，我在寻找一个过渡态，今日优化出来的结果没有虚频，然而查看输出文件是这样的，请问各位老师我 ...

过渡态初猜不当，或者其它因素考虑不当、关键词不合适等等。没具体信息没法说

lurensan · 发表于 Post on 2023-2-18 15:39:41

本帖最后由 lurensan 于 2023-2-18 15:41 编辑

各位老师好，请问TS1-R-up-Re3与前一步的反应能垒是2.85-(-6.42)还是2.85-(-9.41)呢，我的个人理解是2.85-(-6.42)。谢谢各位老师
不好意思老师，前面网卡了

吴双 · 发表于 Post on 2023-2-18 16:23:26

老师您好，我想请教一下。我想用高斯计算两种分子的HOMO能级，比较这两种分子的失电子能力。分子1电荷和自选多重度是 -3 4。分子2是在分子1的基础上加了一个卤素离子，电荷和自旋多重度是-4 4。在优化中对于分子1加了3个H+使其保持电中性。分子2加了4个H+使其保持电中性。最后对比HOMO能级的时候发现有α-HOMO,和β-HOMO这块应该怎么选择呢？

糖糖糖豆9988 · 发表于 Post on 2023-2-18 16:49:12

请教各位老师，使用“#p UCCSD(T,T1diag) cc-pVTZ”计算单点能，out文件在附件中，这是意外中断了吗？

wzkchem5 · 发表于 Post on 2023-2-18 16:59:20

吴双发表于 2023-2-18 09:23
老师您好，我想请教一下。我想用高斯计算两种分子的HOMO能级，比较这两种分子的失电子能力。分子1电荷和自 ...

不要用HOMO能级比较失电子能力。参见http://sobereva.com/543

wzkchem5 · 发表于 Post on 2023-2-18 17:00:46

lurensan 发表于 2023-2-18 08:39
各位老师好，请问TS1-R-up-Re3与前一步的反应能垒是2.85-(-6.42)还是2.85-(-9.41)呢，我的个人理解是2.85 ...

当且仅当R-up-Re3和R-down-Re4是快平衡时，是2.85-(-9.41)，否则是2.85-(-6.42)

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

[综合交流] 简单量子化学问题答疑专帖

关于gaussian16结构优化和频率分析

求助；使用gaussian16进行qst2过渡态搜索，已纠正原子序号，依然报错

过渡态无虚频，怎么办