计算化学公社

标题: 高斯view出现 valid structure found代表什么 [打印本页]

作者
Author: 陈亮打陈列 时间: 2021-12-8 19:18
标题: 高斯view出现 valid structure found代表什么
小弟想用qst2算一个SN2取代案例，很小一个体系，但是我用了链接工具确定原子顺序是对的后，提示了该字符。计算也没法继续。想问有大佬遇到过吗，是什么原因，应该怎么解决呀。谢谢

作者
Author: zjxitcc 时间: 2021-12-8 20:00
本帖最后由 zjxitcc 于 2021-12-8 20:01 编辑

截图或复制输出文件末尾信息，不要自行截取几个单词出来。（也不要将成千上万行输出内容直接粘贴进帖子里）
另，最好晒出你的输入文件，看看格式对不对。

作者
Author: 陈亮打陈列 时间: 2021-12-8 21:45

zjxitcc 发表于 2021-12-8 20:00
截图或复制输出文件末尾信息，不要自行截取几个单词出来。（也不要将成千上万行输出内容直接粘贴进帖子里） ...

截图表示如下。
[img]https://s2.loli.net/2021/12/08/bwrWgSjP7hVpKeA.png https://s2.loli.net/2021/12/08/WJ5gsxXDt6mqOFr.png[/img]因为我接触高斯还不够时间，很多操作都是通过高斯view完成的，外加现在高斯正在做一个大任务，没法终止，输入输出文件我明天再贴上来，您看可以吗

作者
Author: sobereva 时间: 2021-12-9 06:53
本来用connection工具自动对两帧原子顺序联配成功后必然显示这个，又不是报错有什么可问的。
跟计算能否继续没有丝毫联系

作者
Author: 陈亮打陈列 时间: 2021-12-9 08:34

sobereva 发表于 2021-12-9 06:53
本来用connection工具自动对两帧原子顺序联配成功后必然显示这个，又不是报错有什么可问的。
跟计算能否继 ...

谢谢卢神。看来是我后面出现的问题

作者
Author: 陈亮打陈列 时间: 2021-12-10 14:35

zjxitcc 发表于 2021-12-8 20:00
截图或复制输出文件末尾信息，不要自行截取几个单词出来。（也不要将成千上万行输出内容直接粘贴进帖子里） ...

你好，我今天上个任务跑完了，就提交了看看，我实验室我现在做的基本就是优化下现有结构，所以找过渡态解耦i书非常差
# opt=(calcfc,qst2) freq b3lyp/6-31+g(d) geom=connectivity
最后几行如下
New curvilinear step not converged.
OIter  1 Iteration    1 Try  1 RMS(Cart)=  0.57740845 RMS(Int)=  0.61887984
OIter  1 Iteration    2 Try  1 RMS(Cart)=  0.19159043 RMS(Int)=  0.64057522
OIter  1 Iteration    3 Try  1 RMS(Cart)=  0.11881261 RMS(Int)=  0.68243972
OIter  1 Iteration    4 Try  1 RMS(Cart)=  0.19103837 RMS(Int)=  0.59903960
Old curvilinear step not converged, using linear step
RedCar/ORedCr failed for GTrans.
Error termination via Lnk1e in D:\gaosi\G09W\l101.exe at Fri Dec 10 14:30:10 2021.
Job cpu time:  0 days  0 hours  0 minutes  0.0 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=    1

问题还是发生了，而且发生非常快，进去就2070错误，能帮忙看一下吗。谢谢

作者
Author: 陈亮打陈列 时间: 2021-12-10 14:39

zjxitcc 发表于 2021-12-8 20:00
截图或复制输出文件末尾信息，不要自行截取几个单词出来。（也不要将成千上万行输出内容直接粘贴进帖子里） ...

%chk=C:\Users\Admin\Desktop\NAOH.chk
# opt=(calcfc,qst2) freq b3lyp/6-31+g(d) geom=connectivity

Title Card Required

0 1
C                -1.91107340 1.12299841 0.85823695
H                -1.55441897 0.11418840 0.85823695
H                -1.55440056 1.62739660 1.73188846
H                -2.98107340 1.12301159 0.85823695
Cl             -1.32439657 1.95266272 -0.57879730
Na             -4.66139032 0.43582779 2.43945796
O                -2.64108603 0.09063202 2.64636480
H                -2.37593929 0.27229157 3.55096251

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5
6 7 1.0
7 8 1.0
8

Title Card Required

0 1
C                -2.18927026 0.85755402 1.47684818
H                -1.48032980 0.19977863 1.01898907
H                -1.69370555 1.75145546 1.79348018
H                -2.95421275 1.10481778 0.77071533
Cl             -0.60176710 0.20948226 3.60137650
Na             -2.36319186 -0.50231328 5.05149026
O                -2.77672319 0.21152836 2.60930124
H                -2.09041964 -0.01031576 3.24284100

1 2 1.0 3 1.0 4 1.0 7 1.0
2
3
4
5 6 1.0
6
7 8 1.0
8
这是输入文件
他最后一行说
New curvilinear step not converged.

作者
Author: zjxitcc 时间: 2021-12-10 14:45
本帖最后由 zjxitcc 于 2021-12-10 14:46 编辑

陈亮打陈列发表于 2021-12-10 14:39
%chk=C:%users\Admin\Desktop\NAOH.chk
# opt=(calcfc,qst2) freq b3lyp/6-31+g(d) geom=connectivity
...

QST2和QST3要求你提供的反应物、产物结构不能相差太大（基本上，一眼看不出如何变化的结构，都可以认为是相差太大），否则无法猜测过渡态，程序就会报错。你的报错就是这个原因。你提供的两个结构相差十万八千里，硬要用QST2的话，需尽可能保证所有原子反应前后的位置不要跑太远，更不能随便点击鼠标在随意位置建立分子。

SN2反应虽简单，但一些SN2过渡态并不好找（概念简单与好算是两码事），要合理构建初始结构。如果你实在不会，可以采用柔性扫描的方式取最高点作为TS初始结构。另外，此反应是不是基元反应需要你自行确认。如果是多步反应，你却当成1步那是很难找到过渡态的。

作者
Author: 陈亮打陈列 时间: 2021-12-10 14:47

zjxitcc 发表于 2021-12-10 14:45
QST2和QST3要求你提供的反应物、产物结构不能相差太大（基本上，一眼看不出如何变化的结构，都可以认为是 ...

我明白了，我现在去修改试试。谢谢您

作者
Author: 陈亮打陈列 时间: 2021-12-10 14:54

zjxitcc 发表于 2021-12-10 14:45
QST2和QST3要求你提供的反应物、产物结构不能相差太大（基本上，一眼看不出如何变化的结构，都可以认为是 ...

谢谢您，已经能够正常跑了

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