利用optDFTw在高斯中对开壳层计算时不收敛情况（已解决）

lyc460360210

情况说明：               关键词： 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 后顺利收敛，额，嗯，社长大法好！
         

lyc460360210

后来发现当w值过大时SCF总会不收敛，然后通过调整w的初猜值为0.18，w顺利收敛。