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

wzkchem5 · 发表于 Post on 2023-5-23 15:28:38

Alicean 发表于 2023-5-23 08:22
各位老师好，我目前在对二甲醚CH3OCH3的初始氢提取反应进行过渡态搜索，但是无论是使用TS还是QST3方式都会 ...

先确定这个反应是不是无垒反应。某些热力学足够有利的氢提取反应根本没有过渡态，因此不管用什么方法也优化不出来过渡态

sxh3231796 · 发表于 Post on 2023-5-23 15:39:34

sobereva 发表于 2023-5-23 15:17
CAS里加上lanczos

谢谢sob老师，我刚试了重算了另一个casscf(10,10)的结果。还是没有找到相应的组态信息，能麻烦您再给看一下吗？

Alicean · 发表于 Post on 2023-5-23 16:14:35

wzkchem5 发表于 2023-5-23 15:28
先确定这个反应是不是无垒反应。某些热力学足够有利的氢提取反应根本没有过渡态，因此不管用什么方法也优 ...

感谢老师的解答，那么我如何确定它是否为无垒反应呢？是有系统的方法，还是需要我查阅文献来确定呢？

hbi.trans · 发表于 Post on 2023-5-23 16:17:48

博主您好，我使用Gaussian16 Linux版本优化了我的结果（茂铁和碘）。也产生了chk文件，但是每次打开的时候总是出现下图的问题，请问怎么解决。附件附上我的输入文件和输出文件。
输出文件我在楼下给出期待各位大佬的帮助

hbi.trans · 发表于 Post on 2023-5-23 16:18:17

本帖最后由 hbi.trans 于 2023-5-23 16:23 编辑

我无法上传输出文件，这个是输出文件的最后内容，不知道是否对解决问题有帮助。

zjxitcc · 发表于 Post on 2023-5-23 16:23:49

Linux下计算产生的chk文件，要在Linux下用formchk小程序将其转化为fch文件，然后到Windows下才能用GaussView打开fch文件。

sobereva · 发表于 Post on 2023-5-23 16:55:45

hbi.trans 发表于 2023-5-23 16:17
博主您好，我使用Gaussian16 Linux版本优化了我的结果（茂铁和碘）。也产生了chk文件，但是每次打开的时候 ...

GaussView打开chk文件时会自动调用本机的formchk将之转换成fch再打开
Windows机子上的formchk没法转Linux算的chk文件，此问题下文都说了
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379（http://bbs.keinsci.com/thread-6020-1-1.html）

所以哪算的就在哪转成fch

sobereva · 发表于 Post on 2023-5-23 16:56:32

sxh3231796 发表于 2023-5-23 15:39
谢谢sob老师，我刚试了重算了另一个casscf(10,10)的结果。还是没有找到相应的组态信息，能麻烦您再给看一 ...

我之前明确说了加上IOp(4/43=3)

kimariyb · 发表于 Post on 2023-5-23 17:29:32

hbi.trans 发表于 2023-5-23 16:18
我无法上传输出文件，这个是输出文件的最后内容，不知道是否对解决问题有帮助。

这是正常结束任务，chk文件是二进制文件，应该用 formchk 命令转化为 fchk 文件，在 gview 里打开

wzkchem5 · 发表于 Post on 2023-5-23 17:34:52

Alicean 发表于 2023-5-23 09:14
感谢老师的解答，那么我如何确定它是否为无垒反应呢？是有系统的方法，还是需要我查阅文献来确定呢？

做柔性扫描。如果从反应物可以连续扫到产物（或反之），两边键长扫得足够远，扫描步长足够小，扫描曲线平滑，且中间没有极大值点，那就可以认为是无垒反应

Alicean · 发表于 Post on 2023-5-23 18:32:00

wzkchem5 发表于 2023-5-23 17:34
做柔性扫描。如果从反应物可以连续扫到产物（或反之），两边键长扫得足够远，扫描步长足够小，扫描曲线平 ...

好的，感谢老师的解答

KCM@SHU · 发表于 Post on 2023-5-23 18:40:24

各位老师我想问一下我在做离子的结合能的时候发现，离子在低配位的时候结合能为正值，在高配位的时候为负值，我想问一下我该怎么判断是我计算参数设置问题导致它低配为为正还是它本身在低配位就是正呢？

sxh3231796 · 发表于 Post on 2023-5-23 18:46:46

本帖最后由 sxh3231796 于 2023-5-23 18:58 编辑

sobereva 发表于 2023-5-23 16:56
我之前明确说了加上IOp(4/43=3)

不好意思老师，传错了文件，是下面这个文件：

wzkchem5 · 发表于 Post on 2023-5-23 19:47:30

KCM@SHU 发表于 2023-5-23 11:40
各位老师我想问一下我在做离子的结合能的时候发现，离子在低配位的时候结合能为正值，在高配位的时候为负值 ...

你这个问题等价于“怎么知道我做的计算是对的”，这个没法回答，因为我们没办法把你所有可能犯的错误列举出来让你一一检查。。。一定要问的话，你必须把你的详细操作流程列出来，然后我们才能告诉你是否做对了

sobereva · 发表于 Post on 2023-5-23 23:41:37

sxh3231796 发表于 2023-5-23 18:46
不好意思老师，传错了文件，是下面这个文件：

诸如Configuration 1651 Symmetry 1 11ab00a1b0这就是组态函数信息，1代表相应轨道双占据，a和b对应单占据alpha和beta电子

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

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

高斯16优化问题