计算化学公社

标题: 利用optDFTw在高斯中对开壳层计算时不收敛情况（已解决） [打印本页]

作者
Author: lyc460360210 时间: 2018-4-26 22:15
标题: 利用optDFTw在高斯中对开壳层计算时不收敛情况（已解决）
情况说明：             关键词： lc-wpbe  6-31g(d)  IOp(3/107=0413900000,3/108=0413900000)

Charge = -1 Multiplicity = 2

            使用sob社长的optDFTw工具对一个70多个原子的分子优化w参数，在第12圈的时候由于不收敛停止（l502 Convergence failure）,参考博文http://sobereva.com/61 ，目前已尝试了以下方法：

1. SCF(novaracc,noincfock) int=acc2e=12

2. int=ultrafine
3. SCF=noDIIS

            不想降低收敛限，同一体系的另外3个分子都顺利跑完得到w参数，只有这个很奇怪。

以下是optDFTw运行日志的末尾,在第11圈的时候 J^2 突然变大很多，不清楚原因。
Iteration:    9
Running: g09 < N.gjf > N.out
Running: g09 < N-1.gjf > N-1.out
Running: g09 < N+1.gjf > N+1.out
w: 0.416951 J:       0.0219834800 J^2:       0.0002426689

Iteration: 10
Running: g09 < N.gjf > N.out
Running: g09 < N-1.gjf > N-1.out
Running: g09 < N+1.gjf > N+1.out
w: 0.413334 J:       0.0217516700 J^2:       0.0002375672

Iteration: 11
Running: g09 < N.gjf > N.out
Running: g09 < N-1.gjf > N-1.out
Running: g09 < N+1.gjf > N+1.out
w: 0.412480 J:       0.0410186800 J^2:       0.0010016469

Iteration: 12
Running: g09 < N.gjf > N.out
Running: g09 < N-1.gjf > N-1.out
Running: g09 < N+1.gjf > N+1.out
Error: segmentation violation
rax 0000000000000000, rbx ffffffffffffffff, rcx 0000148b9f4a2767
rdx 000000000001591c, rsp 00007ffe34d94878, rbp 00007ffe34d94df0
rsi 000000000000000b, rdi 000000000001591c, r8  0000000000000020
r9  0000000000000401, r10 000000000000022d, r11 0000000000000206
r12 00007ffe34d94e38, r13 0000000000000000, r14 0000000000000000
r15 00000000000003e6
  /lib/x86_64-linux-gnu/libpthread.so.0(+0x11390) [0x148b9fb51390]
  /lib/x86_64-linux-gnu/libc.so.6(kill+0x7) [0x148b9f4a2767]
  /home/g16/l502.exe() [0x4afeb9]
  /home/g16/l502.exe() [0x4d1e0a]
  /home/g16/l502.exe() [0x54e80c]
  /home/g16/l502.exe() [0x52d7bf]
  /home/g16/l502.exe() [0x4529db]
  /home/g16/l502.exe() [0x450c2b]
  /home/g16/l502.exe() [0x450b34]
  /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0) [0x148b9f48d830]
  /home/g16/l502.exe() [0x450a29]
Aborted (core dumped)
Error in finding Alpha virt. eigenvalues

         呃呃呃，刚要写求助的时候，发现刚提交加入SCF(novaracc,noincfock) int=acc2e=12 后顺利收敛，额，嗯，社长大法好！

作者
Author: lyc460360210 时间: 2018-4-27 11:05
后来发现当w值过大时SCF总会不收敛，然后通过调整w的初猜值为0.18，w顺利收敛。

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