Не сходится ССП (фенантрен)

[Vit Nhoc]

У меня никак не получается провести на Гауссиане расчёт фенантрена на PCM B3LYP-D3/AUG-ccpVTZ. Хронически не сходится ССП. Я пробовал и NOSYMM, и SCF(XQC). Можно ввести SCF(SLEAZY), тогда энергия посчитается, но не станет считаться что-то ещё (оптимизация геометрии или ЯМР).

[Гесс]

Если scf сводится хоть как то то прочтите волновую функцию из соответствующего chk файла с помощью Guess=read (имя чекпоинта должно быть такое же как у сошедшегося расчета).
(SLEAZY) я никогда не применял, но как вариант можно попробовать свести в меньшем или менее диффузном базисе.
Насколько я понимаю Sleazy исключительно меняет отсечку, это не должно быть причиной чтоб не шла оптимизация. Там я думаю какая то иная причина. Покажите аутпут или хотя бы последние его строки.

[Vit Nhoc]

Самое страннное, что вроде один и тот же файл то считается нормально, то обрывается.

Вот пример входного файла, для которого я использловал Sleazy:

%CHK=bdp.chk
%NPROCSHARED=3
#p B3LYP/aug-cc-pVTZ EmpiricalDispersion=GD3 SCF=xqc
scrf(pcm,solvent=Chloroform) NOSYMM SCF(Sleazy)

BODIPY

0 1
6 -2.825713000 -0.870742000 0.000000000
6 -3.547811000 0.299863000 0.000000000
6 -1.417345000 -0.858814000 0.000000000
6 -2.871675000 1.528167000 0.000000000
6 -0.725660000 0.385026000 0.000000000
6 -1.494442000 1.565084000 0.000000000
6 -0.676364000 -2.082682000 0.000000000
6 0.677054000 -2.080736000 0.000000000
6 1.416355000 -0.856667000 0.000000000
6 0.726535000 0.385205000 0.000000000
1 -3.335613000 -1.825930000 0.000000000
1 -4.628741000 0.272328000 0.000000000
1 -3.431706000 2.453178000 0.000000000
1 -1.006572000 2.528512000 0.000000000
1 -1.225153000 -3.015614000 0.000000000
1 1.225725000 -3.013893000 0.000000000
6 2.826253000 -0.871544000 0.000000000
6 1.493256000 1.568831000 0.000000000
6 2.870366000 1.528672000 0.000000000
6 3.545986000 0.300852000 0.000000000
1 3.337482000 -1.825562000 0.000000000
1 1.002179000 2.529645000 0.000000000
1 3.430708000 2.454205000 0.000000000
1 4.627085000 0.276278000 0.000000000





--Link1--
%CHK=bdp.chk
%NPROCSHARED=3
#p B3LYP/aug-cc-pVTZ EmpiricalDispersion=GD3 SCF=xqc
OPT(LOOSE) scrf(pcm,solvent=Chloroform)
GUESS(READ) GEOM(ALLCHECK) gfinput




--Link1--
%CHK=bdp.chk
%NPROCSHARED=3
#p B3LYP/aug-cc-pVTZ EmpiricalDispersion=GD3 SCF=xqc
NMR scrf(pcm,solvent=Chloroform)
GUESS(READ) GEOM(ALLCHECK) gfinput

Для этого файла не сходится SCF на второй задаче (оптимизации). Конец файла (я вручную оборвал расчёт):

Cycle 59 Pass 1 IDiag 1:
RMSU= 3.25D-09 CP: 1.00D+00 1.14D+00 8.57D-01 -1.67D-01 5.57D-02
CP: -1.04D-01 1.06D+00 -1.99D-01 1.68D-03 -3.02D-01
CP: -2.70D-01 1.00D-01 -7.00D-03 -4.64D-02 2.81D-04
CP: 7.40D-02 2.47D-01 -6.27D-03
E= -539.745443286354 Delta-E= 0.000000000293 Rises=F Damp=F
DIIS: error= 3.91D-06 at cycle 59 NSaved= 16.
NSaved=16 IEnMin= 5 EnMin= -539.745443326420 IErMin= 8 ErrMin= 1.12D-07
ErrMax= 3.91D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.00D-09 BMatP= 2.42D-12
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Large coefficients: NSaved= 16 BigCof= 0.00 CofMax= 10.00 Det=-6.58D-15
Inversion failed. Reducing to 15 matrices.
Coeff-Com: 0.224D-01 0.233D-01 0.150D+00 0.244D+00 0.332D+00-0.836D+00
Coeff-Com: 0.271D+00-0.282D+01 0.312D+01-0.212D+00 0.832D-01-0.278D+00
Coeff-Com: 0.900D+00-0.742D-02 0.145D-01
Coeff: 0.224D-01 0.233D-01 0.150D+00 0.244D+00 0.332D+00-0.836D+00
Coeff: 0.271D+00-0.282D+01 0.312D+01-0.212D+00 0.832D-01-0.278D+00
Coeff: 0.900D+00-0.742D-02 0.145D-01
Gap= 0.173 Goal= None Shift= 0.000
RMSDP=8.22D-08 MaxDP=4.53D-06 DE= 2.93D-10 OVMax= 2.13D-08

Cycle 60 Pass 1 IDiag 1:
RMSU= 5.96D-10 CP: 1.00D+00 1.10D+00 8.17D-01 -4.97D-02 -2.68D-02
CP: -7.59D-02 1.12D+00 2.23D-01 -1.87D-03 -1.51D-01
CP: -7.53D-02 6.07D-02 -3.04D-03 -4.03D-02 2.29D-04
CP: 2.12D-01 1.03D-01 -3.12D-03 4.22D-01
E= -539.745443286406 Delta-E= -0.000000000053 Rises=F Damp=F
DIIS: error= 3.90D-06 at cycle 60 NSaved= 16.
NSaved=16 IEnMin= 5 EnMin= -539.745443326420 IErMin= 8 ErrMin= 1.12D-07
ErrMax= 3.90D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 9.95D-10 BMatP= 2.42D-12
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Large coefficients: NSaved= 16 BigCof= 0.00 CofMax= 10.00 Det=-2.85D-15
Inversion failed. Reducing to 15 matrices.
Large coefficients: NSaved= 15 BigCof= 10.48 CofMax= 10.00 Det=-2.95D-15
Inversion failed. Reducing to 14 matrices.
Large coefficients: NSaved= 14 BigCof= 10.49 CofMax= 10.00 Det=-2.95D-15
Inversion failed. Reducing to 13 matrices.
Large coefficients: NSaved= 13 BigCof= 10.13 CofMax= 10.00 Det=-2.96D-15
Inversion failed. Reducing to 12 matrices.
Large coefficients: NSaved= 12 BigCof= 10.26 CofMax= 10.00 Det=-2.96D-15
Inversion failed. Reducing to 11 matrices.
Coeff-Com: 0.210D+00 0.232D+01-0.231D+01 0.265D+00-0.674D+01 0.758D+01
Coeff-Com: -0.128D+01-0.446D+00 0.139D+01 0.986D+01-0.985D+01
Coeff: 0.210D+00 0.232D+01-0.231D+01 0.265D+00-0.674D+01 0.758D+01
Coeff: -0.128D+01-0.446D+00 0.139D+01 0.986D+01-0.985D+01
Gap= 0.173 Goal= None Shift= 0.000
RMSDP=9.61D-07 MaxDP=5.04D-05 DE=-5.25D-11 OVMax= 2.36D-07

Cycle 61 Pass 1 IDiag 1:
Restarting incremental Fock formation.
E= -539.745443286105 Delta-E= 0.000000000301 Rises=F Damp=F
DIIS: error= 3.67D-06 at cycle 61 NSaved= 12.
NSaved=12 IEnMin= 2 EnMin= -539.745443326420 IErMin= 5 ErrMin= 1.12D-07
ErrMax= 3.67D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 9.07D-10 BMatP= 2.89D-12
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Large coefficients: NSaved= 12 BigCof= 14.09 CofMax= 10.00 Det=-2.43D-15
Inversion failed. Reducing to 11 matrices.
Coeff-Com: 0.200D+00 0.148D+01-0.212D+01 0.164D+00-0.340D+01 0.422D+01
Coeff-Com: -0.308D+00-0.232D+00 0.985D+00-0.119D+00 0.128D+00
Coeff: 0.200D+00 0.148D+01-0.212D+01 0.164D+00-0.340D+01 0.422D+01
Coeff: -0.308D+00-0.232D+00 0.985D+00-0.119D+00 0.128D+00
Gap= 0.173 Goal= None Shift= 0.000
RMSDP=7.56D-07 MaxDP=4.21D-05 DE= 3.01D-10 OVMax= 2.02D-07

Cycle 62 Pass 1 IDiag 1:
RMSU= 7.56D-07 CP: 1.00D+00
E= -539.745443286513 Delta-E= -0.000000000407 Rises=F Damp=F
DIIS: error= 3.62D-06 at cycle 62 NSaved= 12.
NSaved=12 IEnMin= 2 EnMin= -539.745443326420 IErMin= 5 ErrMin= 1.12D-07
ErrMax= 3.62D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 8.86D-10 BMatP= 2.89D-12
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: 0.212D+00 0.168D+01-0.231D+01 0.121D+00-0.376D+01 0.466D+01
Coeff-Com: -0.373D+00-0.166D+00 0.936D+00-0.870D-01 0.253D+00-0.157D+00
Coeff: 0.212D+00 0.168D+01-0.231D+01 0.121D+00-0.376D+01 0.466D+01
Coeff: -0.373D+00-0.166D+00 0.936D+00-0.870D-01 0.253D+00-0.157D+00
Gap= 0.173 Goal= None Shift= 0.000

[Гесс]

62 итерации это может быть тупо мало
при scf=xqc на 64ой только начинается qc
энергия осцилирует или как?

[Vit Nhoc]

Осциллирует. Минимальная энергия была на 12-м шаге SCF.

[Vit Nhoc]

При SCF(XQC) выдаётся следующее:

CP: -8.52D-03 -3.20D-03 1.73D-02 4.55D-02 -7.12D-02
CP: -2.13D-02 -2.47D-02 -3.77D-02 -3.11D-01
FoF2E skips out because all densities are zero.
ILin=14 X=2.441D-04 Y=-5.397454123382D+02 DE= 2.88D-09 F= -2.18D-10
Restarting incremental Fock formation.
ILin=15 X=1.221D-04 Y=-5.397454123382D+02 DE= 2.90D-09 F= -2.19D-10
RMSU= 2.30D-09 CP: 1.00D+00
ILin=16 X=6.104D-05 Y=-5.397454123382D+02 DE= 2.88D-09 F= -2.19D-10
The polynomial fit failed. Using point 1.
A contracting polynomial of degree 16 produced 0.0000
Search did not lower the energy significantly.
ILin= 1 X=0.000D+00 Y=-5.397454123411D+02 DE= 0.00D+00 F= -1.99D-12
RMSU= 3.65D-07 CP: 1.00D+00 -3.00D+00
ILin= 2 X=1.000D+00 Y=-5.397454123381D+02 DE= 3.06D-09 F= -1.72D-12
RMSU= 2.95D-12 CP: 1.00D+00 -3.00D+00 4.84D-01
ILin= 3 X=5.000D-01 Y=-5.397454123384D+02 DE= 2.71D-09 F= -1.95D-12
RMSU= 6.05D-16 CP: 1.00D+00 -3.00D+00 2.78D-01 8.03D-01
ILin= 4 X=3.000D-01 Y=-5.397454123384D+02 DE= 2.68D-09 F= -2.05D-12
RMSU= 1.14D-17 CP: 1.00D+00 -3.00D+00 7.17D-02 2.65D-01 2.90D+00
FoF2E skips out because all densities are zero.
ILin= 5 X=1.000D-01 Y=-5.397454123387D+02 DE= 2.45D-09 F= -2.14D-12
RMSU= 4.67D-18 CP: 1.00D+00 -3.00D+00 3.30D-02 1.26D-01 2.82D+00
CP: 9.04D-01
FoF2E skips out because all densities are zero.
ILin= 6 X=6.250D-02 Y=-5.397454123386D+02 DE= 2.50D-09 F= -2.16D-12
RMSU= 1.48D-16 CP: 1.00D+00 -3.00D+00 7.55D-04 1.12D-03 3.00D+00
CP: 1.13D+00 5.18D-01
ILin= 7 X=3.125D-02 Y=-5.397454123387D+02 DE= 2.36D-09 F= -2.17D-12
RMSU= 3.76D-18 CP: 1.00D+00 -5.60D-01 3.83D-14 -1.39D-04 5.92D-01
CP: 1.59D-01 3.00D-01 9.16D-05
FoF2E skips out because all densities are zero.
ILin= 8 X=1.563D-02 Y=-5.397454123384D+02 DE= 2.72D-09 F= -2.19D-12
RMSU= 3.54D-18 CP: 1.00D+00 7.20D-01 -8.97D-15 3.20D-05 -1.22D-01
CP: -7.21D-02 3.93D-02 -5.69D-05 1.44D-01
FoF2E skips out because all densities are zero.
ILin= 9 X=7.813D-03 Y=-5.397454123381D+02 DE= 3.00D-09 F= -2.19D-12
RMSU= 3.74D-18 CP: 1.00D+00 1.36D+00 2.08D-14 -7.78D-05 -5.95D-01
CP: -2.25D-01 -1.65D-01 -1.47D-04 2.72D-02 3.25D-01
FoF2E skips out because all densities are zero.
ILin=10 X=3.906D-03 Y=-5.397454123381D+02 DE= 3.01D-09 F= -2.20D-12
RMSU= 4.11D-18 CP: 1.00D+00 1.68D+00 4.92D-14 -1.82D-04 -5.14D-01
CP: -1.59D-01 -3.92D-01 -2.33D-04 -6.48D-02 3.02D-01
CP: 3.43D-01
FoF2E skips out because all densities are zero.
ILin=11 X=1.953D-03 Y=-5.397454123382D+02 DE= 2.95D-09 F= -2.20D-12
RMSU= 3.84D-18 CP: 1.00D+00 1.84D+00 6.63D-14 -2.46D-04 -1.46D+00
CP: -5.28D-01 -1.26D-01 -1.84D-04 3.31D-01 3.00D-01
CP: 2.57D-01 1.49D-02
FoF2E skips out because all densities are zero.
ILin=12 X=9.766D-04 Y=-5.397454123380D+02 DE= 3.12D-09 F= -2.20D-12
RMSU= 3.64D-18 CP: 1.00D+00 1.92D+00 7.61D-14 -2.81D-04 -1.63D+00
CP: -5.85D-01 -1.36D-01 -9.16D-05 3.45D-01 3.74D-01
CP: 2.06D-01 4.83D-02 2.00D-01
FoF2E skips out because all densities are zero.
ILin=13 X=4.883D-04 Y=-5.397454123379D+02 DE= 3.23D-09 F= -2.20D-12
RMSU= 3.46D-18 CP: 1.00D+00 1.96D+00 8.13D-14 -2.99D-04 -1.63D+00
CP: -5.90D-01 -2.25D-01 -1.78D-04 3.66D-01 2.93D-01
CP: 3.19D-01 1.60D-01 3.48D-01 2.41D-01
FoF2E skips out because all densities are zero.
ILin=14 X=2.441D-04 Y=-5.397454123379D+02 DE= 3.25D-09 F= -2.20D-12
RMSU= 2.99D-18 CP: 1.00D+00 1.98D+00 8.38D-14 -3.08D-04 -1.23D+00
CP: -4.67D-01 -2.53D-01 -1.28D-04 1.78D-01 4.68D-01
CP: 4.36D-01 1.47D-01 2.41D-01 4.34D-02 1.85D-03
FoF2E skips out because all densities are zero.
ILin=15 X=1.221D-04 Y=-5.397454123380D+02 DE= 3.06D-09 F= -2.20D-12
RMSU= 3.02D-18 CP: 1.00D+00 1.99D+00 8.49D-14 -3.13D-04 -1.26D+00
CP: -4.61D-01 -2.17D-01 -2.18D-04 1.39D-01 4.69D-01
CP: 3.84D-01 2.90D-01 2.21D-01 1.18D-01 6.61D-02
CP: 5.98D-02
FoF2E skips out because all densities are zero.
ILin=16 X=6.104D-05 Y=-5.397454123379D+02 DE= 3.17D-09 F= -2.20D-12
The polynomial fit failed. Using point 1.
A contracting polynomial of degree 16 produced 0.0000
Search did not lower the energy significantly.
No lower point found -- run aborted.
Error termination via Lnk1e in C:\QuantChem\Gaussian\G09w\l508.exe at Wed May 31 23:21:43 2017.
Job cpu time: 1 days 18 hours 23 minutes 53.0 seconds.
File lengths (MBytes): RWF= 452 Int= 0 D2E= 0 Chk= 193 Scr= 1

[Гесс]

а чтение чекпоинта с этими линками проходит нормально?
Можно попробовать оптимизировать с sleazy критерием и потом применять нормальный на эту геометрию.
Крайне удивительно что qc осцилирует, я такого за собой не помню.

[Yurii]

g16 должен справиться: думаю, что при этом не нужны ни XQC. ни SLEAZY. Хотя использование диффузных функций с растворителем не совсем корректно.

[AlexKsen]

Я бы посоветовал диффузную функцию убрать.

[EvgeniX]

g16 должен справиться: думаю, что при этом не нужны ни XQC. ни SLEAZY. Хотя использование диффузных функций с растворителем не совсем корректно.

А в g18 и подавно сойдётся!

[Гесс]

Я бы посоветовал диффузную функцию убрать.

Всмысле верхнюю, до календарного, или всмысле несколько верхних до maug, или всмысле вообще все, до cc-pVTZ?
Ну сходимость то облегчится но мы незнаем задачи ТС, возможно он хочет в конечном счете таки свести полностью отаугуменченный базис. Но вообще да, это можно использовать как указывалось выше:

как вариант можно попробовать свести в меньшем или менее диффузном базисе.

прочтите волновую функцию из соответствующего chk файла

[Yurii]

.............................................................................
.............................................................................
.............................................................................
#p B3LYP/aug-cc-pVTZ EmpiricalDispersion=GD3 NMR scrf(pcm,solvent=Chlo
roform) GUESS(READ) GEOM(ALLCHECK) gfinput
.............................................................................
.............................................................................
.............................................................................
C,0,-2.825713,-0.870742,0.
C,0,-3.547811,0.299863,0.
C,0,-1.417345,-0.858814,0.
C,0,-2.871675,1.528167,0.
C,0,-0.72566,0.385026,0.
C,0,-1.494442,1.565084,0.
C,0,-0.676364,-2.082682,0.
C,0,0.677054,-2.080736,0.
C,0,1.416355,-0.856667,0.
C,0,0.726535,0.385205,0.
H,0,-3.335613,-1.82593,0.
H,0,-4.628741,0.272328,0.
H,0,-3.431706,2.453178,0.
H,0,-1.006572,2.528512,0.
H,0,-1.225153,-3.015614,0.
H,0,1.225725,-3.013893,0.
C,0,2.826253,-0.871544,0.
C,0,1.493256,1.568831,0.
C,0,2.870366,1.528672,0.
C,0,3.545986,0.300852,0.
H,0,3.337482,-1.825562,0.
H,0,1.002179,2.529645,0.
H,0,3.430708,2.454205,0.
H,0,4.627085,0.276278,0.
.............................................................................
.............................................................................
.............................................................................
SCF GIAO Magnetic shielding tensor (ppm):
1 C Isotropic = 49.7738 Anisotropy = 161.1216
XX= -49.4321 YX= 10.7866 ZX= 0.0000
XY= 16.3060 YY= 41.5388 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 157.2146
Eigenvalues: -51.4064 43.5131 157.2146
2 C Isotropic = 52.8210 Anisotropy = 183.8349
XX= 12.1673 YX= 45.8576 ZX= 0.0000
XY= 44.9803 YY= -29.0819 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 175.3776
Eigenvalues: -58.3398 41.4251 175.3776
3 C Isotropic = 44.8657 Anisotropy = 209.5657
XX= -20.5861 YX= -7.2625 ZX= 0.0000
XY= -9.6536 YY= -29.3930 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 184.5761
Eigenvalues: -34.5252 -15.4539 184.5761
4 C Isotropic = 52.7920 Anisotropy = 182.0186
XX= 11.2375 YX= -44.9610 ZX= 0.0000
XY= -47.2364 YY= -26.9992 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 174.1378
Eigenvalues: -57.7868 42.0251 174.1378
5 C Isotropic = 46.6492 Anisotropy = 197.3458
XX= -12.7448 YX= 13.8903 ZX= 0.0000
XY= 15.7636 YY= -25.5207 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 178.2130
Eigenvalues: -35.2772 -2.9883 178.2130
6 C Isotropic = 56.1448 Anisotropy = 177.8837
XX= -47.8771 YX= -11.9399 ZX= 0.0000
XY= -9.3443 YY= 41.5776 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 174.7340
Eigenvalues: -49.1257 42.8262 174.7340
7 C Isotropic = 51.2169 Anisotropy = 152.1618
XX= -45.5886 YX= -16.0127 ZX= 0.0000
XY= -22.7475 YY= 46.5812 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 152.6582
Eigenvalues: -49.4977 50.4903 152.6582
8 C Isotropic = 51.5388 Anisotropy = 151.7319
XX= -0.9276 YX= -53.4167 ZX= 0.0000
XY= -46.2726 YY= 2.8506 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 152.6934
Eigenvalues: -48.9189 50.8420 152.6934
9 C Isotropic = 45.3011 Anisotropy = 210.6965
XX= -34.4215 YX= -0.1808 ZX= 0.0000
XY= 1.3555 YY= -15.4407 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 185.7654
Eigenvalues: -34.4397 -15.4226 185.7654
10 C Isotropic = 46.8141 Anisotropy = 199.0989
XX= -7.5173 YX= 12.1583 ZX= 0.0000
XY= 10.5043 YY= -31.5871 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 179.5467
Eigenvalues: -36.0821 -3.0223 179.5467
11 H Isotropic = 23.3878 Anisotropy = 7.7564
XX= 23.0209 YX= 0.6265 ZX= 0.0000
XY= 0.2724 YY= 28.5222 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 18.6202
Eigenvalues: 18.6202 22.9844 28.5587
12 H Isotropic = 23.7185 Anisotropy = 5.7447
XX= 26.4427 YX= 1.6767 ZX= 0.0000
XY= 1.6716 YY= 25.0130 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 19.6997
Eigenvalues: 19.6997 23.9075 27.5483
13 H Isotropic = 23.6665 Anisotropy = 6.2006
XX= 26.4011 YX= -1.8706 ZX= 0.0000
XY= -1.9327 YY= 25.2157 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 19.3828
Eigenvalues: 19.3828 23.8165 27.8003
14 H Isotropic = 22.4307 Anisotropy = 10.2901
XX= 22.5534 YX= -1.9526 ZX= 0.0000
XY= -1.0368 YY= 28.9591 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 15.7794
Eigenvalues: 15.7794 22.2218 29.2907
15 H Isotropic = 23.5205 Anisotropy = 7.7229
XX= 23.4291 YX= -0.7019 ZX= 0.0000
XY= -0.5647 YY= 28.5925 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 18.5389
Eigenvalues: 18.5389 23.3526 28.6691
16 H Isotropic = 23.5224 Anisotropy = 7.7428
XX= 26.3160 YX= -2.5773 ZX= 0.0000
XY= -2.7110 YY= 25.7322 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 18.5191
Eigenvalues: 18.5191 23.3639 28.6843
17 C Isotropic = 49.6676 Anisotropy = 161.6071
XX= 24.6924 YX= -35.5290 ZX= 0.0000
XY= -40.6949 YY= -33.0952 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 157.4057
Eigenvalues: -52.0279 43.6251 157.4057
18 C Isotropic = 56.5430 Anisotropy = 179.2704
XX= 2.9509 YX= -44.7882 ZX= 0.0000
XY= -47.0725 YY= -9.3784 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 176.0566
Eigenvalues: -49.5560 43.1285 176.0566
19 C Isotropic = 52.7869 Anisotropy = 181.0847
XX= -57.5376 YX= -0.0874 ZX= 0.0000
XY= 1.5215 YY= 42.3883 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 173.5100
Eigenvalues: -57.5428 42.3935 173.5100
20 C Isotropic = 52.5441 Anisotropy = 182.0860
XX= 27.6726 YX= 33.9473 ZX= 0.0000
XY= 35.5878 YY= -43.9752 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 173.9347
Eigenvalues: -58.0726 41.7700 173.9347
21 H Isotropic = 23.3795 Anisotropy = 7.7709
XX= 27.1697 YX= -2.6275 ZX= 0.0000
XY= -2.2339 YY= 24.3108 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 18.6581
Eigenvalues: 18.6581 22.9204 28.5601
22 H Isotropic = 22.4169 Anisotropy = 10.3012
XX= 25.4878 YX= -3.0831 ZX= 0.0000
XY= -4.0144 YY= 25.9674 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 15.7955
Eigenvalues: 15.7955 22.1708 29.2844
23 H Isotropic = 23.6707 Anisotropy = 6.1820
XX= 23.8404 YX= -0.1528 ZX= 0.0000
XY= -0.1462 YY= 27.7864 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 19.3852
Eigenvalues: 19.3852 23.8347 27.7920
24 H Isotropic = 23.7244 Anisotropy = 5.7369
XX= 27.0497 YX= 1.2184 ZX= 0.0000
XY= 1.2843 YY= 24.4131 ZY= 0.0000
XZ= 0.0000 YZ= 0.0000 ZZ= 19.7105
Eigenvalues: 19.7105 23.9138 27.5490
.............................................................................
.............................................................................
.............................................................................
Normal termination of Gaussian...

[AlexKsen]

Думаю, что можно полностью снять auq, оставив cc-pVTZ.

[Yurii]

А в g18 и подавно сойдётся!

Проблемы у ТС возникли прежде всего из-за наличия растворителя, а g16 в этом вопросе значительно продвинулся по сравнению с g09. Например, ''In Gaussian 16, we use a continuous surface charge formalism that ensures continuity, smoothness and robustness of the reaction field, which also has continuous derivatives with respect to atomic positions and external perturbing fields". Впрочем, почитайте про SCRF в g16 и тогда многие вопросы отпадут сами собой.

[alien308]

Для сходимости имеет смысл поиграть коэффициентом демпфирования (damp ?) и параметрами алгоритма ускорения сходимости. В качестве первого шага можно усилить демпфирование и вообще отключить ускорители сходимости. Не забывайте что при этом скорость сходимости снизится и потребуется больше шагов.
Я не знаю про Гауссиан. Это общие соображения применимые к большинству программ.