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

lym563557413

请教优化的问题，谢谢！优化时死在了l303,如下：
Ernie: Thresh=  0.10000D-02 Tol=  0.10000D-05 Strict=F.
There are  3200 symmetry adapted cartesian basis functions of A   symmetry.
There are  3200 symmetry adapted basis functions of A   symmetry.
  3200 basis functions,  5600 primitive gaussians,  3200 cartesian basis functions
   600 alpha electrons      600 beta electrons
       nuclear repulsion energy     47609.3102299344 Hartrees.
IExCor=    0 DFT=F Ex=HF Corr=None ExCW=0 ScaHFX=  1.000000
ScaDFX=  1.000000  1.000000  1.000000  1.000000 ScalE2=  1.000000  1.000000
IRadAn=      0 IRanWt=     -1 IRanGd=            0 ICorTp=0 IEmpDi=  4
NAtoms=  320 NActive=  320 NUniq=  320 SFac= 1.00D+00 NAtFMM=   60 NAOKFM=T Big=T
Integral buffers will be    131072 words long.
Raffenetti 1 integral format.
Two-electron integral symmetry is turned on.
Leave Link  301 at Sat Jul 25 15:54:06 2020, MaxMem=    33554432 cpu:         0.9
(Enter /share/apps/gaussian/g09d01/g09/l302.exe)
NPDir=0 NMtPBC=     1 NCelOv=     1 NCel=       1 NClECP=     1 NCelD=      1
         NCelK=      1 NCelE2=     1 NClLst=     1 CellRange=     0.0.
One-electron integrals computed using PRISM.
PrsmSu:  requested number of processors reduced to:   1 ShMem   1 Linda.
PrsmSu:  requested number of processors reduced to:   1 ShMem   1 Linda.
GSVD:  LWork=   -12585890 too small for GESVD, short by    53568290 words or    53568290 for optimal perf.
NBasis=  3200 RedAO= T EigKep=  3.28D-04  NBF=  3200
NBsUse=  3200 1.00D-06 EigRej= -1.00D+00 NBFU=  3200
Leave Link  302 at Sat Jul 25 15:55:33 2020, MaxMem=    33554432 cpu:      2024.9
(Enter /share/apps/gaussian/g09d01/g09/l303.exe)
DipDrv:  MaxL=1.
Out-of-memory error in routine ShPair-CalcS2-2 (IEnd=      18917376 MxCore=      13056829)
Use %mem=38MW to provide the minimum amount of memory required to complete this step.
Error termination via Lnk1e in /share/apps/gaussian/g09d01/g09/l303.exe at Sat Jul 25 15:55:34 2020.
Job cpu time:       0 days  0 hours 34 minutes  3.1 seconds.
File lengths (MBytes):  RWF=    555 Int=      0 D2E=      0 Chk=      1 Scr=      1

下面是输入文件：
%Nproc=24
mem=24GB
%chk=1.chk
#p hf/6-31g(d,p) 5d opt

pdb

0 1
C    -1              -3.0090  14.0650  -2.3780
O    -1              -3.4230  12.9310  -2.5350
C    -1              -3.9580  15.2150  -2.0310
H     0              -4.9740  14.9210  -2.2450
H     0              -3.7020  16.0830  -2.6200
H     0              -3.8650  15.4540  -0.9810
N    -1              -1.7410  14.3450  -2.4980
C    -1              -0.7780  13.2630  -2.8340
C    -1              -0.6420  12.3160  -1.6570
O    -1              -0.3680  11.1440  -1.8190
C    -1               0.5450  13.9680  -3.1320
C    -1               1.6490  12.9250  -3.3180
O    -1               2.2680  12.5630  -2.3320
O    -1               1.8560  12.5060  -4.4460
H     0              -1.4230  15.2580  -2.3680
H     0              -1.1170  12.7230  -3.6920
H     0               0.4450  14.5530  -4.0340
H     0               0.8010  14.6170  -2.3080
N    -1              -0.8700  12.7960  -0.4740
C    -1              -0.7920  11.8920   0.6880
C    -1              -1.8360  10.7830   0.4710
O    -1              -1.7250   9.6980   0.9860
C    -1              -1.2060  12.7390   1.8990
O    -1              -1.7250   9.6980   0.9860
C    -1              -1.2060  12.7390   1.8990
C    -1              -2.6360  12.4290   2.1080
C    -1              -3.6390  13.0000   1.4380
C    -1              -3.2120  11.3630   2.9030
N    -1              -4.8200  12.3900   1.8060
C    -1              -4.6020  11.3680   2.7140
C    -1              -2.6660  10.4200   3.7810
C    -1              -5.4240  10.4510   3.3610
C    -1              -3.4880   9.4970   4.4430
C    -1              -4.8680   9.5100   4.2290
H     0              -1.1250  13.7330  -0.3570
H     0               0.1970  11.4910   0.8130
H     0              -0.6260  12.4580   2.7680
H     0              -1.0820  13.7890   1.6850
H     0              -3.5300  13.8060   0.7310
H     0              -5.7050  12.6340   1.4770
H     0              -1.6090  10.4020   3.9420
H     0              -6.4730  10.4490   3.1730
H     0              -3.0520   8.7710   5.1130
H     0              -5.5020   8.8000   4.7380
N    -1              -2.8870  11.1010  -0.2500
C    -1              -3.9690  10.1100  -0.4750
C    -1              -3.4660   9.0710  -1.4400
O    -1              -3.4920   7.8900  -1.1580
C    -1              -5.1450  10.8850  -1.0430
C    -1              -6.3070  10.7760  -0.0600
C    -1              -7.2010  12.0120  -0.1800
C    -1              -7.1260   9.5230  -0.3780
H     0              -2.9710  12.0010  -0.6210
H     0              -4.2480   9.6540   0.4550
H     0              -2.9710  12.0010  -0.6210
H     0              -4.2480   9.6540   0.4550
H     0              -4.8680  11.9210  -1.1700
H     0              -5.4340  10.4610  -1.9910
H     0              -5.9100  10.7040   0.9500
H     0              -6.6650  12.7960  -0.6950
H     0              -7.4770  12.3540   0.8070
H     0              -8.0920  11.7600  -0.7360
H     0              -7.1470   9.3690  -1.4470
H     0              -8.1340   9.6500  -0.0130
H     0              -6.6750   8.6670   0.1020
N    -1              -2.9200   9.4960  -2.5480
C    -1              -2.3230   8.4930  -3.4610
C    -1              -1.3900   7.6710  -2.5770
O    -1              -1.1400   6.5020  -2.7940
C    -1              -1.5430   9.2930  -4.5070
C    -1              -0.6080   8.3560  -5.2720
C    -1              -0.1050   9.0560  -6.5370
C    -1              -0.6330   8.3210  -7.7710
N    -1               0.5210   8.2690  -8.7110
H     0              -2.8440  10.4580  -2.7300
H     0              -3.0820   7.8800  -3.9170
H     0              -2.2350   9.7530  -5.1960
H     0              -0.9610  10.0570  -4.0150
H     0               0.2340   8.0960  -4.6450
H     0              -1.1420   7.4590  -5.5480
H     0              -0.4570  10.0770  -6.5460
H     0               0.9740   9.0460  -6.5490
H     0              -0.9510   7.3220  -7.5060
H     0              -1.4470   8.8720  -8.2170
H     0               0.8240   9.2360  -8.9430
H     0              -1.4470   8.8720  -8.2170
H     0               0.8240   9.2360  -8.9430
H     0               0.2370   7.7750  -9.5820
H     0               1.3090   7.7600  -8.2640
N    -1              -0.9170   8.3180  -1.5360
C    -1              -0.0380   7.6620  -0.5460
C    -1              -0.8420   6.6640   0.3140
O    -1              -0.4010   5.5530   0.5370
C    -1               0.4900   8.8150   0.3020
H     0              -1.1650   9.2560  -1.3970
H     0               0.7800   7.1730  -1.0360
H     0               0.1010   8.7310   1.3060
H     0               0.1670   9.7520  -0.1300
H     0               1.5670   8.7840   0.3260
N    -1              -2.0190   7.0260   0.8030
C    -1              -2.7930   6.0580   1.6270
C    -1              -3.3130   4.9220   0.7490
O    -1              -3.5890   3.8380   1.2210
C    -1              -3.9520   6.8650   2.2110
C    -1              -4.7860   5.9720   3.0970
C    -1              -4.2150   5.3780   4.2290
C    -1              -6.1320   5.7380   2.7880
C    -1              -4.9870   4.5490   5.0510
C    -1              -6.9050   4.9090   3.6100
C    -1              -6.3330   4.3150   4.7420
H     0              -2.3950   7.9200   0.6320
H     0              -2.1790   5.6700   2.4180
H     0              -3.5620   7.6870   2.7910
H     0              -4.5640   7.2480   1.4080
H     0              -3.1770   5.5580   4.4670
H     0              -6.5730   6.1960   1.9150
H     0              -3.1770   5.5580   4.4670
H     0              -6.5730   6.1960   1.9150
H     0              -4.5460   4.0910   5.9240
H     0              -7.9430   4.7290   3.3710
H     0              -6.9300   3.6760   5.3760
N    -1              -3.4470   5.1580  -0.5260
C    -1              -3.9450   4.0880  -1.4310
C    -1              -2.9120   2.9620  -1.5190
O    -1              -3.1630   1.8400  -1.1260
C    -1              -4.1210   4.7650  -2.7870
C    -1              -5.2300   4.0830  -3.5520
C    -1              -4.9520   2.9490  -4.3250
C    -1              -6.5360   4.5840  -3.4870
C    -1              -5.9820   2.3170  -5.0340
C    -1              -7.5640   3.9520  -4.1960
C    -1              -7.2870   2.8180  -4.9690
O    -1              -8.3010   2.1950  -5.6680
H     0              -3.2190   6.0360  -0.8890
H     0              -4.8880   3.7130  -1.0840
H     0              -4.3740   5.8040  -2.6370
H     0              -3.2020   4.6930  -3.3460
H     0              -3.9460   2.5630  -4.3750
H     0              -6.7500   5.4590  -2.8910
H     0              -5.7680   1.4420  -5.6310
H     0              -8.5710   4.3380  -4.1460
H     0              -9.0710   2.7670  -5.6460
N    -1              -1.7470   3.2550  -2.0270
C    -1              -0.7030   2.2000  -2.1330
C    -1              -0.3690   1.6570  -0.7440
O    -1              -0.0360   0.4990  -0.5840
C    -1               0.5110   2.9000  -2.7470
O    -1              -0.0360   0.4990  -0.5840
C    -1               0.5110   2.9000  -2.7470
C    -1               1.7060   1.9460  -2.7420
O    -1               2.3010   1.7780  -1.6900
O    -1               2.0080   1.3990  -3.7900
H     0              -1.5590   4.1650  -2.3360
H     0              -1.0380   1.4080  -2.7750
H     0               0.2840   3.1890  -3.7630
H     0               0.7520   3.7790  -2.1680
N    -1              -0.4650   2.4800   0.2630
C    -1              -0.1650   2.0050   1.6370
C    -1              -1.2530   1.0340   2.0920
O    -1              -1.0180   0.1420   2.8850
C    -1              -0.1700   3.2650   2.5050
C    -1              -0.0890   2.8690   3.9810
C    -1               0.8350   3.8380   4.7200
C    -1               1.7300   3.0560   5.6840
N    -1               2.4380   4.0970   6.4810
H     0              -0.7410   3.4040   0.1160
H     0               0.8020   1.5360   1.6690
H     0               0.6800   3.8800   2.2490
H     0              -1.0810   3.8170   2.3320
H     0              -1.0760   2.9060   4.4180
H     0               0.3040   1.8660   4.0640
H     0               1.4490   4.3660   4.0050
H     0               0.2420   4.5470   5.2770
H     0               1.1290   2.4290   6.3280
H     0               2.4440   2.4620   5.1360
H     0               3.0690   4.6410   5.8590
H     0               2.9980   3.6390   7.2290
H     0               1.7420   4.7380   6.9110
H     0               2.9980   3.6390   7.2290
H     0               1.7420   4.7380   6.9110
N    -1              -2.4470   1.1990   1.5890
C    -1              -3.5490   0.2840   1.9880
C    -1              -3.4050  -1.0500   1.2540
O    -1              -3.6600  -2.1020   1.8030
C    -1              -4.8480   1.0170   1.6110
C    -1              -5.3700   0.5590   0.2440
C    -1              -5.9100   0.7330   2.6750
H     0              -2.6130   1.9190   0.9480
H     0              -3.5190   0.1290   3.0470
H     0              -4.6520   2.0750   1.5810
H     0              -5.8760  -0.3890   0.3510
H     0              -4.5450   0.4490  -0.4410
H     0              -6.0610   1.2930  -0.1430
H     0              -5.6560   1.2540   3.5860
H     0              -5.9520  -0.3290   2.8680
H     0              -6.8720   1.0730   2.3210
N    -1              -2.9870  -1.0120   0.0200
C    -1              -2.8140  -2.2750  -0.7410
C    -1              -1.5840  -3.0250  -0.2260
O    -1              -1.4700  -4.2280  -0.3630
C    -1              -2.6120  -1.8410  -2.1860
H     0              -2.7820  -0.1530  -0.4010
H     0              -3.6920  -2.8820  -0.6650
H     0              -3.0770  -2.5610  -2.8420
H     0              -1.5560  -1.7870  -2.3990
H     0              -3.0630  -0.8730  -2.3340
N    -1              -0.6650  -2.3170   0.3730
C    -1               0.5560  -2.9760   0.9030
C    -1               0.2220  -3.7250   2.1910
C    -1               0.5560  -2.9760   0.9030
C    -1               0.2220  -3.7250   2.1910
O    -1               0.6060  -4.8620   2.3780
C    -1               1.5340  -1.8320   1.1820
C    -1               2.6530  -2.3270   2.1010
C    -1               3.8190  -1.3370   2.0640
O    -1               4.2490  -1.0000   0.9730
O    -1               4.2620  -0.9330   3.1260
H     0              -0.7810  -1.3550   0.4790
H     0               0.9680  -3.6430   0.1690
H     0               1.9590  -1.4870   0.2500
H     0               1.0100  -1.0190   1.6620
H     0               2.2800  -2.4060   3.1120
H     0               2.9940  -3.2940   1.7650
N    -1              -0.4970  -3.0960   3.0790
C    -1              -0.8600  -3.7750   4.3510
C    -1              -1.8140  -4.9320   4.0670
O    -1              -1.7000  -5.9970   4.6380
C    -1              -1.5420  -2.7000   5.2010
C    -1              -2.3290  -3.3650   6.3320
C    -1              -2.8760  -2.2900   7.2740
C    -1              -2.5950  -2.6910   8.7240
N    -1              -3.8770  -2.4500   9.4440
H     0              -0.7980  -2.1800   2.9050
H     0               0.0220  -4.1340   4.8440
H     0              -0.7930  -2.0450   5.6190
H     0              -2.2180  -2.1280   4.5830
H     0              -3.1510  -3.9300   5.9150
H     0              -1.6780  -4.0270   6.8840
H     0              -2.3940  -1.3470   7.0610
H     0              -3.9410  -2.1920   7.1300
H     0              -2.3940  -1.3470   7.0610
H     0              -3.9410  -2.1920   7.1300
H     0              -2.3220  -3.7370   8.7770
H     0              -1.8150  -2.0760   9.1420
H     0              -3.6820  -2.2640  10.4470
H     0              -4.4870  -3.2890   9.3570
H     0              -4.3570  -1.6260   9.0280
N    -1              -2.7480  -4.7370   3.1830
C    -1              -3.6960  -5.8290   2.8560
C    -1              -2.9700  -6.9080   2.0490
O    -1              -3.3500  -8.0670   2.0520
C    -1              -4.7890  -5.1660   2.0210
C    -1              -5.8380  -4.5530   2.9500
C    -1              -5.3790  -3.1640   3.3960
C    -1              -7.1700  -4.4320   2.2040
H     0              -2.8200  -3.8770   2.7300
H     0              -4.1170  -6.2410   3.7520
H     0              -4.3510  -4.3920   1.4090
H     0              -5.2560  -5.9050   1.3920
H     0              -5.9660  -5.1850   3.8160
H     0              -5.5740  -2.4500   2.6090
H     0              -4.3200  -3.1860   3.6080
H     0              -5.9170  -2.8750   4.2860
H     0              -7.9670  -4.2760   2.9140
H     0              -7.3540  -5.3400   1.6470
H     0              -7.1260  -3.5960   1.5220
N    -1              -1.9110  -6.5330   1.3760
C    -1              -1.1410  -7.5240   0.5800
C    -1              -0.4610  -8.4950   1.5330
O    -1              -0.5410  -9.6980   1.3810
C    -1              -0.1100  -6.6990  -0.1920
O    -1              -0.5410  -9.6980   1.3810
C    -1              -0.1100  -6.6990  -0.1920
C    -1               1.0150  -7.6120  -0.6870
C    -1               2.1410  -6.7610  -1.2770
C    -1               2.2920  -7.0730  -2.7670
N    -1               3.7610  -7.1710  -2.9930
H     0              -1.6130  -5.6010   1.4090
H     0              -1.7900  -8.0440  -0.1020
H     0              -0.5880  -6.2250  -1.0380
H     0               0.3040  -5.9420   0.4570
H     0               1.3970  -8.1930   0.1400
H     0               0.6320  -8.2760  -1.4470
H     0               1.9050  -5.7130  -1.1500
H     0               3.0660  -6.9850  -0.7680
H     0               1.8100  -8.0120  -3.0030
H     0               1.8790  -6.2750  -3.3640
H     0               3.9420  -7.5080  -3.9590
H     0               4.1730  -7.8390  -2.3080
H     0               4.1940  -6.2350  -2.8670
N    -1               0.1860  -7.9760   2.5350
C    -1               0.8490  -8.8630   3.5240
C    -1              -0.2310  -9.6580   4.2520
O    -1              -0.2030 -10.8710   4.3030
C    -1               1.5750  -7.9230   4.4870
C    -1               2.0340  -8.7080   5.7180
C    -1               3.1020  -7.9070   6.4660
O    -1               2.7550  -6.8870   7.0400
O    -1               4.2470  -8.3270   6.4530
H     0               0.2160  -7.0030   2.6460
H     0               1.5490  -9.5180   3.0360
H     0               2.4350  -7.4930   3.9920
H     0               1.5490  -9.5180   3.0360
H     0               2.4350  -7.4930   3.9920
H     0               0.9060  -7.1340   4.7950
H     0               1.1900  -8.8810   6.3700
H     0               2.4490  -9.6550   5.4080
N    -1              -1.1940  -8.9680   4.8010
C    -1              -2.3070  -9.6600   5.5160
C    -1              -2.7220 -10.9280   4.7580
O    -1              -2.4480 -12.0310   5.1890
C    -1              -3.4530  -8.6480   5.5300
H     0              -1.1880  -7.9880   4.7350
H     0              -2.0140  -9.9000   6.5260
H     0              -3.1490  -7.7670   6.0750
H     0              -4.3150  -9.0880   6.0090
H     0              -3.7040  -8.3750   4.5160
N    -1              -3.3810 -10.7880   3.6340
C    -1              -3.8020 -11.9920   2.8730
C    -1              -3.0220 -12.0900   1.5590
O    -1              -3.5730 -12.4330   0.5320
C    -1              -5.2930 -11.7950   2.6030
C    -1              -5.5080 -10.4710   1.9150
C    -1              -5.4300 -10.3860   0.5210
C    -1              -5.7810  -9.3250   2.6720
C    -1              -5.6270  -9.1560  -0.1180
C    -1              -5.9780  -8.0960   2.0340
C    -1              -5.9010  -8.0110   0.6400
H     0              -3.5990  -9.8990   3.2930
H     0              -3.6500 -12.8700   3.4660
H     0              -5.6540 -12.5920   1.9710
H     0              -5.8330 -11.8050   3.5390
H     0              -5.2200 -11.2690  -0.0640
H     0              -5.8330 -11.8050   3.5390
H     0              -5.2200 -11.2690  -0.0640
H     0              -5.8410  -9.3910   3.7480
H     0              -5.5680  -9.0900  -1.1940
H     0              -6.1890  -7.2120   2.6180
H     0              -6.0510  -7.0630   0.1480
N    -1              -1.7500 -11.8010   1.5480
H     0              -1.3040 -11.5250   2.3750
H     0              -1.2410 -11.8620   0.7130

447951397

输入文件没有指定内存用量，最少38MW，大约320M。你%men少了个%号。Use %mem=38MW to provide the minimum amount of memory required to complete this step.

哎呀妈呀

%mem

snljty

文件里写的很清楚内存不够。
现在不要用HF这种已经没有任何实际意义的说法

sobereva

lym563557413 发表于 2020-7-25 16:02
请教优化的问题，谢谢！优化时死在了l303,如下：
Ernie: Thresh=  0.10000D-02 Tol=  0.10000D-05 Strict ...

很长的输入文件别直接贴，应当上传附件方式提供。
mem前头加上%，否则不生效
绝对不要用HF，毫无实际意义

sobereva

阿di达思 发表于 2020-7-23 14:52
@sobereva 请问老师：
1，算不同过渡金属（Fe3+和Zn2+）与有机分子的螯合作用，应根据高斯输出文件中的哪 ...

1 要求不高的话通过电子能量变化考察，更严格是通过自由能变化考察
2 看
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471（http://bbs.keinsci.com/thread-12488-1-1.html）

ljk

Inversion failed.  Reducing to  9 matrices.
Coeff-Com: -0.205D+01 0.154D+01-0.363D+00 0.365D+00-0.467D+00 0.285D+00
Coeff-Com:  0.803D+00 0.582D+00 0.300D+00
Coeff:     -0.205D+01 0.154D+01-0.363D+00 0.365D+00-0.467D+00 0.285D+00
Coeff:      0.803D+00 0.582D+00 0.300D+00
Gap=     0.065 Goal=   None    Shift=    0.000
RMSDP=5.62D-06 MaxDP=9.54D-04 DE= 3.20D-10 OVMax= 2.66D-06

Cycle  35  Pass 1  IDiag  1:
RMSU=  1.94D-06    CP:  1.00D+00  3.00D+00 -3.00D+00 -3.00D+00 -3.00D+00
                    CP: -3.00D+00  3.00D+00 -3.00D+00  3.00D+00 -3.00D+00
                    CP: -3.00D+00  3.00D+00  7.48D-01  3.00D+00  3.00D+00
E= -9309.66796365872     Delta-E=        0.000000001011 Rises=F Damp=F
DIIS: error= 5.97D-08 at cycle  35 NSaved=  10.
NSaved=10 IEnMin= 4 EnMin= -9309.66796366272     IErMin=10 ErrMin= 5.97D-08
ErrMax= 5.97D-08 EMaxC= 1.00D-01 BMatC= 7.77D-12 BMatP= 1.34D-11
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.252D+01 0.315D+01-0.157D+01 0.710D-01-0.542D+00 0.468D+00
Coeff-Com:  0.657D+00 0.125D+00 0.444D+00 0.710D+00
Coeff:     -0.252D+01 0.315D+01-0.157D+01 0.710D-01-0.542D+00 0.468D+00
Coeff:      0.657D+00 0.125D+00 0.444D+00 0.710D+00
Gap=     0.065 Goal=   None    Shift=    0.000
RMSDP=1.04D-06 MaxDP=1.51D-04 DE= 1.01D-09 OVMax= 1.77D-06

Cycle  36  Pass 1  IDiag  1:
RMSU=  5.70D-07    CP:  1.00D+00  3.00D+00  5.33D-01 -3.00D+00 -3.00D+00
                    CP: -3.00D+00  3.00D+00 -3.00D+00  3.00D+00 -3.00D+00
                    CP: -3.00D+00  3.00D+00  8.51D-01  3.00D+00  3.00D+00
                    CP:  9.56D-01
E= -9309.66796365954     Delta-E=       -0.000000000826 Rises=F Damp=F
DIIS: error= 3.06D-08 at cycle  36 NSaved=  11.
NSaved=11 IEnMin= 4 EnMin= -9309.66796366272     IErMin=11 ErrMin= 3.06D-08
ErrMax= 3.06D-08 EMaxC= 1.00D-01 BMatC= 2.33D-12 BMatP= 7.77D-12
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.157D+00 0.695D+00-0.310D+00 0.683D-02-0.455D+00 0.363D+00
Coeff-Com: -0.662D-02-0.269D+00 0.219D+00 0.416D+00 0.499D+00
Coeff:     -0.157D+00 0.695D+00-0.310D+00 0.683D-02-0.455D+00 0.363D+00
Coeff:     -0.662D-02-0.269D+00 0.219D+00 0.416D+00 0.499D+00
Gap=     0.065 Goal=   None    Shift=    0.000

sobereva

ljk 发表于 2020-7-26 20:47
Inversion failed.  Reducing to  9 matrices.
Coeff-Com: -0.205D+01 0.154D+01-0.363D+00 0.365D+00-0 ...

Inversion failed. Reducing to 9 matrices不算报错提示

当前明显还没算完，继续等着就完了

xiaobchine

这是输入文件：
%nprocshared=16
%chk=I-anion.chk
%mem=8GB
#p BP86/genecp opt=(noeigen) freq=noraman scrf(smd,solvent=DiChloroEthane)

Title Card Required

-1 1
I                 -3.43558279    0.21472392    0.00000000

I 0
SDD
****

I 0
SDD

--Link1--
%nprocshared=16
%chk=I-anion.chk
%mem=8GB
#p M06L/genecp int=grid=99770 geom=allcheck guess=read scrf(smd,solvent=DiChloroEthane)

I 0
SDD
****

I 0
SDD

这是输出文件最后几段：
Leave Link  703 at Sun Jul 26 22:51:37 2020, MaxMem=  1073741824 cpu:         1.8
(Enter /share/apps/gauss/G09D01/g09/l716.exe)
Dipole        = 0.00000000D+00 0.00000000D+00 0.00000000D+00
***** Axes restored to original set *****
-------------------------------------------------------------------
Center     Atomic                   Forces (Hartrees/Bohr)
Number     Number              X              Y              Z
-------------------------------------------------------------------
      1       53           0.000000000    0.000000000    0.000000000
-------------------------------------------------------------------
Cartesian Forces:  Max     0.000000000 RMS     0.000000000
Leave Link  716 at Sun Jul 26 22:51:37 2020, MaxMem=  1073741824 cpu:         0.7
(Enter /share/apps/gauss/G09D01/g09/l103.exe)

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GEDIIS/GDIIS optimizer.
Search for a local minimum.
Step number   2 out of a maximum of   20
All quantities printed in internal units (Hartrees-Bohrs-Radians)
RMS Force = .19758D-18 SwitMx=.10000D-02 MixMth= 2
Mixed Optimization -- En-DIIS/RFO-DIIS
Swapping is turned off.
Update second derivatives using D2CorN and points    2
The second derivative matrix:
                          X1        Y1        Z1
           X1           0.00000
           Y1           0.00000   0.00000
           Z1           0.00000   0.00000   0.00000
ITU=  0  0
     Eigenvalues ---
En-DIIS/RFO-DIIS IScMMF=        0 using points:     2    1
RFO step:  Lambda=-5.96046448D-09.
RFO step:  ModMin= 1 Lambda=-5.96046448D-09 EMin= 0.00000000D+00
Using NR instead of RFO step for point 1.
RFO step:  ModMin= 1 Lambda=-5.96046448D-09 EMin= 0.00000000D+00
Using NR instead of RFO step for point 2.
NNeg= 0 NP= 2 Switch=  2.50D-03 Rises=F DC=  2.84D-14 SmlDif=  1.00D-05
RMS Error=               NaN NUsed= 2 EDIIS=F
(END)

谢谢各位老师！

dingniu2

请问老师，平时做高斯计算，可以一次写很多个gjf文件，将文件名写到一个bat文件里，例如：
g16 1.gjf
g16 2.gjf
........
然后sh *.bat&就可以批量计算一系列任务。最近用xtb也想这样，但xtb输出结构文件名都相同，都是xtbopt.xyz，不能放在一个文件夹里算，所以我将几十个xyz文件（1.xyz，2.xyz.......）分别放入不同文件夹（1,2......）里，然后写个bat文件为：
xtb ./1/1.xyz --chrg 0 --uhf 0 --gfn 1 --ohess verytight
xtb ./2/2.xyz --chrg 0 --uhf 0 --gfn 1 --ohess verytight
.......
sh运行之后不行，xtb不会进入每个文件夹去算，请问老师我该如何进行批量的xtb计算？谢谢老师！

ljk

我的电脑16GB 6核12线程
INPUT：
%mem=200MW
%nprocs=8
.......
虽然这样设置但CPU还是用到了12核，内存设置最多200MW，多了就报错

ljk

sobereva 发表于 2020-7-26 20:57
Inversion failed. Reducing to 9 matrices不算报错提示

当前明显还没算完，继续等着就完了

谢谢，我算了一天多，退了重算了，6核12线程，16GB的电脑，但是内存最多只能分配200MW，我想尽可能用到更多内存，设多了就报错，怎么办哪

喵星大佬

在WSL里面用linux版，没有办法，win版的gaussian基本就是个废物，WSL里面可以发挥接近90%的性能

ljk

喵星大佬 发表于 2020-7-27 00:44
在WSL里面用linux版，没有办法，win版的gaussian基本就是个废物，WSL里面可以发挥接近90%的性能

装虚拟机再装linux吗？

喵星大佬

一个原子怎么opt，dft方法写freq=noraman没有意义因为本来就不会算。直接写#p freq sdd scrf=(smd,solvent=dichloroethane) bp86
就完了，写了那么多了没用的。