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

sobereva · 发表于 Post on 2023-6-19 20:51:16

NGC626 发表于 2023-6-19 15:35
各位老师好，我想用ORCA给一个18碳环和过渡金属吸附的簇模型做optfreq，希望可以尽量节省计算时间，能定性 ...

MN15

sobereva · 发表于 Post on 2023-6-19 20:52:33

hehe321123 发表于 2023-6-19 19:35
各位老师，我已经找到自旋极化单重态，也做了stable，我想做自由基加成的过渡态，需要读自旋极化单重态的波 ...

如果此任务的初始结构与记录了自旋极化单重态波函数文件里的相同，读

jlm · 发表于 Post on 2023-6-19 21:31:22

Gaussina计算时，输入文件的路径，怎么看是不是正确的，如果不正确，怎么修改。

sobereva · 发表于 Post on 2023-6-19 22:12:50

jlm 发表于 2023-6-19 21:31
Gaussina计算时，输入文件的路径，怎么看是不是正确的，如果不正确，怎么修改。

默写Gaussian 100遍

说清楚什么叫你所谓的正确和不正确

如果对Gaussian使用完全没概念，看以下链接
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）
北京科音初级量子化学培训班（http://www.keinsci.com/workshop/KEQC_content.html）

NGC626 · 发表于 Post on 2023-6-20 00:18:47

sobereva 发表于 2023-6-19 20:51
MN15

谢谢老师

NGC626 · 发表于 Post on 2023-6-20 15:04:43

sobereva 发表于 2023-6-19 20:51
MN15

老师好，请问在该级别下计算的结果可否达到发文章的认可度呢？

吴双 · 发表于 Post on 2023-6-20 16:44:42

想跟大家请教一个问题，用高斯做二维柔性扫描势能面的时候什么样的曲面图，才能证明是无垒呢？

sobereva · 发表于 Post on 2023-6-20 18:58:03

NGC626 发表于 2023-6-20 15:04
老师好，请问在该级别下计算的结果可否达到发文章的认可度呢？

我只能告诉你在我来看这是合理的选择，而且也没有明显更好选择。至于实际遇到的审稿人是否comment，那不是我说了算

sobereva · 发表于 Post on 2023-6-20 18:59:42

吴双发表于 2023-6-20 16:44
想跟大家请教一个问题，用高斯做二维柔性扫描势能面的时候什么样的曲面图，才能证明是无垒呢？

没前提的情况下没法回答，什么体系、扫的什么变量都必须非常明确
扫描的变量必须对应反应过程主要改变的变量，如果随着这个变量的变化，能量单调变化，才能说明无能垒

NGC626 · 发表于 Post on 2023-6-20 19:11:10

sobereva 发表于 2023-6-20 18:58
我只能告诉你在我来看这是合理的选择，而且也没有明显更好选择。至于实际遇到的审稿人是否comment，那不 ...

谢谢sob老师

hehe321123 · 发表于 Post on 2023-6-21 10:01:59

各位老师，IRC跑成这样合理吗？

primashen · 发表于 Post on 2023-6-21 11:53:31

本帖最后由 primashen 于 2023-6-21 11:54 编辑

请教各位老师，一醌式结构分子，多种泛函优化闭壳单重态后波函数稳定性测试，结果显示B3LYP和M062x显示波函数稳定，而Cam-b3lyp、wb97xd等存在RHF-UHF不稳定，暂无参考实验数据，请教各位老师如何验证该结构基态为闭壳单重态还是开壳单重态？

wzkchem5 · 发表于 Post on 2023-6-21 13:51:08

primashen 发表于 2023-6-21 04:53
请教各位老师，一醌式结构分子，多种泛函优化闭壳单重态后波函数稳定性测试，结果显示B3LYP和M062x显示波函 ...

无法严格判断，因为闭壳层单重态、开壳层单重态本来就是和泛函有关的，不是体系的客观性质。如果用精确的泛函（准确来说是精确的spin density functional）算，会发现所有单重态都是闭壳层单重态，因为单重态的自旋密度处处为0，所以Kohn-Sham体系的spin density也必须处处为0，所以只可能是闭壳层单重态。如果用full CI算，会发现随着键长拉长，体系的双自由基性质由很小的值逐渐趋近于100%，但不管键长多短，甚至键长在平衡键长附近的时候，双自由基性质都不是0

cokie · 发表于 Post on 2023-6-21 14:00:17

各位老师好，关于Rehm-Weller方程有两个小疑惑想请教大家，

一个是该方程中的εref（也就是在测量氧化还原电势时所用溶剂的介电常数）可以通过guassian计算的方式获得吗（我看到测量溶剂DCM里有加比如0.1M的六氟磷酸叔丁胺盐，文中给出了8.98这样的介电常数但没有给计算方法）？我没有找到一个比较合适的测量方法。
第二个就是我看到一篇文章里用了如此方式的近似（但文中并没有给出可以做此近似的适用范围），请问这个方程在什么情况下可以做这样的近似呢？

sobereva · 发表于 Post on 2023-6-21 14:03:06

hehe321123 发表于 2023-6-21 10:01
各位老师，IRC跑成这样合理吗？

明显不对
注意看下文了解细节和注意事项
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

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

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

关于Gaussinas的输入文件