简单量子化学问题答疑专帖（旧版，已关帖）

wzkchem5 · 发表于 Post on 2020-10-25 19:12:09

kyuu 发表于 2020-10-24 20:48
真能强词夺理无理取闹

报错确实和你的不一样。
Cycle 127  Pass 1  IDiag  1:
E= -150.280377184120    Delta-E=       0.000712579022 Rises=F Damp=F
DIIS: error= 1.10D-02 at cycle 127 NSaved= 6.
NSaved= 6 IEnMin= 1 EnMin= -150.281115491133    IErMin= 3 ErrMin= 7.23D-03
ErrMax= 1.10D-02  0.00D+00 EMaxC= 1.00D-01 BMatC= 2.46D-03 BMatP= 1.05D-03
IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00
Rare condition: small coef for last iteration:  0.000D+00
Coeff-En: 0.484D-02 0.000D+00 0.000D+00 0.000D+00 0.507D+00 0.488D+00
Coeff:    0.484D-02 0.000D+00 0.000D+00 0.000D+00 0.507D+00 0.488D+00
Gap=    0.000 Goal= None Shift= 0.000
RMSDP=2.74D-02 MaxDP=1.61D-01 DE= 7.13D-04 OVMax= 2.60D-03

Cycle 128  Pass 1  IDiag  1:
E= -150.281110174776    Delta-E=    -0.000732990655 Rises=F Damp=F
DIIS: error= 7.26D-03 at cycle 128 NSaved= 7.
NSaved= 7 IEnMin= 1 EnMin= -150.281115491133    IErMin= 3 ErrMin= 7.23D-03
ErrMax= 7.26D-03  0.00D+00 EMaxC= 1.00D-01 BMatC= 1.06D-03 BMatP= 1.05D-03
IDIUse=3 WtCom= 9.29D-01 WtEn= 7.15D-02
Large coefficients: NSaved=  7 BigCof= 22.28 CofMax= 10.00 Det=-1.42D-08
Inversion failed.  Reducing to  6 matrices.
Large coefficients: NSaved=  6 BigCof= 17.70 CofMax= 10.00 Det=-1.43D-08
Inversion failed.  Reducing to  5 matrices.
EnCoef did 100 forward-backward iterations
Coeff-Com:  0.190D+01-0.277D+01 0.532D+00-0.111D+01 0.244D+01
Coeff-En: 0.372D+00 0.852D-01 0.000D+00 0.126D+00 0.417D+00
Coeff:    0.179D+01-0.257D+01 0.494D+00-0.102D+01 0.230D+01
Gap=    0.010 Goal= None Shift= 0.000
RMSDP=2.25D-04 MaxDP=1.96D-03 DE=-7.33D-04 OVMax= 2.94D-03

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RM06L) =  -150.281110175    A.U. after  129 cycles
         NFock=128  Conv=0.23D-03    -V/T= 2.0019
KE= 1.499902970266D+02 PE=-4.121198954422D+02 EE= 8.362570365842D+01
Convergence failure -- run terminated.
Error termination via Lnk1e in /work/g16/l502.exe at Sun Oct 25 19:06:48 2020.
Job cpu time:    0 days  0 hours  0 minutes 14.2 seconds.
Elapsed time:    0 days  0 hours  0 minutes 14.8 seconds.
File lengths (MBytes):  RWF=    17 Int=    0 D2E=    0 Chk=    2 Scr=    2

我是l502报错，你是l501报错。
我的输入文件是
#p opt 6-311g(d,p) m06l  td scf=(NoVarAcc,noincfock) int=superfine

O2

0 1
O 0. 0. 0.
O 0. 0. 1.2

高斯版本是G16 revision A.03。

所以现在算是清楚了，你要是真的连输出文件都懒得贴，起码得告诉我输入文件的坐标，以及你用的高斯的版本。我现在已经证明给你看了，当我既不知道你的分子坐标，也不知道你的高斯版本的时候，我报的错和你报的错是不同的，所以我不可能解决你的问题，其他人也不可能。

sophieding · 发表于 Post on 2020-10-25 21:52:38

本帖最后由 sophieding 于 2020-10-25 21:57 编辑

各位老师好，我在学习计算重组能，也看了sob老师的帖子http://sobereva.com/330，有几个地方没搞明白，望得到解答。
1.是否也可用于自由基的计算？
2.如污染物为A，其电子转移的过程如下：A-+B•—A•+B-，想计算重组能，用于能垒图的绘制。请问是否是分别优化（A-）和（A•）的结构，，然后算E1(A-在带负电荷状态的能量)和E2（A-自由基状态下的能量）；E3(A•在带负电荷状态的能量)和E4（A•自由基状态下的能量）？最后根据贴吧里的公式进行加减
3.不知上述的E1/E2/E3/E4的能量是否描述的对？

snljty · 发表于 Post on 2020-10-25 23:10:51

请问如何在orca中测试双杂化泛函的KS部分波函数的稳定性？直接加

%scf
HFTyp UHF
STABPerform true
end

复制代码

会有报错。谢谢！

kyuu · 发表于 Post on 2020-10-26 02:13:26

本帖最后由 kyuu 于 2020-10-26 03:34 编辑

wzkchem5 发表于 2020-10-25 19:12
报错确实和你的不一样。
Cycle 127 Pass 1 IDiag 1:
E= -150.280377184120 Delta-E= 0 ...

证明什么了，scf都没能收敛，连个正版的g16都没有这不是无理取闹是什么

sobereva · 发表于 Post on 2020-10-26 05:38:28

snljty 发表于 2020-10-25 23:10
请问如何在orca中测试双杂化泛函的KS部分波函数的稳定性？直接加
会有报错。谢谢！

我这里倒是没有报错

h2o.out (24.88 KB, 下载次数 Times of downloads: 2)

sobereva · 发表于 Post on 2020-10-26 05:43:10

sophieding 发表于 2020-10-25 21:52
各位老师好，我在学习计算重组能，也看了sob老师的帖子http://sobereva.com/330，有几个地方没搞明白，望得 ...

可以

我不知道贴吧指的是什么玩意

对

sophieding · 发表于 Post on 2020-10-26 09:03:23

sobereva 发表于 2020-10-26 05:43
可以

我不知道贴吧指的是什么玩意

好的，谢谢老师

snljty · 发表于 Post on 2020-10-26 11:22:32

sobereva 发表于 2020-10-26 05:38
我这里倒是没有报错

谢谢卢老师。我测的是一个对称性破缺的计算，用Gaussian得到了SCF级别的稳定对称性破缺波函数，传给ORCA做初猜然后算的双杂化。去掉波函数稳定性测试部分的关键词后计算正常。报错文件贴上来了，不知道您（和其他各位老师）是否有时间看一下。多谢！
主要报错信息是
Electron density on the grid ... found on disk
Building xc-kernel on the grid ... [file orca_tools/Tool-Numint/qcfxc.cpp, line 15455, Process 0]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 15455, Process 1]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

C2H6.zip (35.83 KB, 下载次数 Times of downloads: 1)

axuanyan · 发表于 Post on 2020-10-26 12:07:07

请教一下，
在隐式溶剂模型下分别计算反应物、产物和过渡态的能量并通过校正获得不同温度下的自由能计算得到的自由能垒，和计算气态时各个结构的能量再加上溶解自由能得到的自由能垒相比，两者的数据都是合理的吗？

wzkchem5 · 发表于 Post on 2020-10-26 20:22:16

kyuu 发表于 2020-10-26 02:13
证明什么了，scf都没能收敛，连个正版的g16都没有这不是无理取闹是什么

我只能说，想跟人吵架，还能显得自己占一点理，办法有的是。但你现在这么做不行。基本的吵架技巧你还是得学一学，比如你自己scf也没收敛，你去攻击别人拿你的输入文件也没收敛，就显得可笑了，等于在攻击你自己。
你要是因为生活不顺心，想跟人吵一架发泄一下，那我可以继续奉陪，我不会往心里去。要能疏解你心里的压力，那也是我功德一件。

wzkchem5 · 发表于 Post on 2020-10-26 20:25:42

axuanyan 发表于 2020-10-26 12:07
请教一下，
在隐式溶剂模型下分别计算反应物、产物和过渡态的能量并通过校正获得不同温度下的自由能计算得 ...

当分子结构受溶剂影响不明显时，后者比较合理；当分子结构受溶剂影响明显时，前者比较合理
原因是一般的溶剂模型是用分子的气相结构拟合的，但拟合时选用的分子都是气相结构和溶剂里的结构定性相符的分子

kyuu · 发表于 Post on 2020-10-26 22:46:42

本帖最后由 kyuu 于 2020-10-26 22:50 编辑

wzkchem5 发表于 2020-10-26 20:22
我只能说，想跟人吵架，还能显得自己占一点理，办法有的是。但你现在这么做不行。基本的吵架技巧你还是得 ...

你这是心理疾病，得治 btw 你太恶心了没人陪你玩你好自为之

axuanyan · 发表于 Post on 2020-10-27 02:36:48

wzkchem5 发表于 2020-10-26 20:25
当分子结构受溶剂影响不明显时，后者比较合理；当分子结构受溶剂影响明显时，前者比较合理
原因是一般的 ...

感谢，但这里我还有一些疑惑，
除了从实验数据推断溶剂对分子结构和构象影响，从理论上怎么进行判断呢？

sobereva · 发表于 Post on 2020-10-27 07:15:51

axuanyan 发表于 2020-10-27 02:36
感谢，但这里我还有一些疑惑，
除了从实验数据推断溶剂对分子结构和构象影响，从理论上怎么进行判断呢？

溶剂和真空下分别优化进行对比。
更严格的是用molclus（http://www.keinsci.com/research/molclus.html）在真空和溶剂环境下分别做构象搜索，看看最稳定的构象的结构、构象分布如何

sobereva · 发表于 Post on 2020-10-27 07:26:03

snljty 发表于 2020-10-26 11:22
谢谢卢老师。我测的是一个对称性破缺的计算，用Gaussian得到了SCF级别的稳定对称性破缺波函数，传给ORCA ...

这是程序目前的限制，对于PWPB95等极个别的双杂化做不了稳定性测试

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

[综合交流] 简单量子化学问题答疑专帖（旧版，已关帖）

评分 Rate

浏览过的版块