求助 计算一维周期结构的开壳层单重态

diffcwmol

各位同仁：


请问我在计算一个反铁磁耦合体系，关键词为：
#p uhf/6-31G** pbc nosymm guess=mix stable=opt
我要计算体系的自旋极化单重态，体系具有一维周期性，原子个数不超过20个。

无论我怎么改理论方法和基组，都会显示错误，错误信息如下。
如果我删掉周期性参数TV那一行，程序可以正常结束。
我简单测试了两个H原子的一维周期体系，计算开壳层单重态，仍然会出现下面的错误。

请问是不是G16算不了周期性体系的自旋极化单重态？或者需要怎么特殊处理？

######
UHF ground state
Doing stability rather than CIS.
Two-electron integral symmetry not used.
MDV= 20132659200 DFT=T DoStab=T Mixed=F DoRPA=F DoScal=F NonHer=F
Making orbital integer symmetry assigments:
    12 initial guesses have been made.
Convergence on wavefunction:    0.001000000000000
Davidson Disk Diagonalization:  ConvIn= 1.00D-03 SkipCon=T Conv= 1.00D-03.
Max sub-space:  2000 roots to seek:    12 dimension of matrix:        19698
Iteration     1 Dimension    12 NMult     0 NNew     12
CISAX will form    12 AO SS matrices at one time.
NMat=    12 NSing=    12 JSym2X= 0.
FoFJK:  IHMeth= 1 ICntrl=     500 DoSepK=F KAlg= 0 I1Cent=           0 FoldK=F
IRaf= 760000000 NMat=      12 IRICut=       1 DoRegI=T DoRafI=F ISym2E= 0 IDoP0=0 IntGTp=1.
FoFCou: FMM=T IPFlag=           0 FMFlag=      990000 FMFlg1=        1001
         NFxFlg=           0 DoJE=F BraDBF=F KetDBF=F FulRan=T
         wScrn=  0.000000 ICntrl=       500 IOpCl=  0 I1Cent=           0 NGrid=           0
         NMat0=   12 NMatS0=     12 NMatT0=    0 NMatD0=   12 NMtDS0=    0 NMtDT0=    0
Symmetry not used in FoFCou.
NB too large for integers on this machine.
Error termination via Lnk1e in /home/gaussian/g16/l914.exe at Wed May 24 01:11:56 2023.
Job cpu time:       0 days  4 hours 15 minutes 14.7 seconds.
Elapsed time:       0 days  0 hours  7 minutes  7.6 seconds.
File lengths (MBytes):  RWF=   2200 Int=      0 D2E=      0 Chk=     15 Scr=      1
##################

zjxitcc

我猜可以算，但不支持stable=opt

diffcwmol

zjxitcc 发表于 2023-6-2 10:05
我猜可以算，但不支持stable=opt

谢谢哥们 试了一下删掉stable，就没有报错了。