计算化学公社

标题: 求助计算一维周期结构的开壳层单重态 [打印本页]

作者
Author: diffcwmol 时间: 2023-6-2 06:01
标题: 求助计算一维周期结构的开壳层单重态
各位同仁：

请问我在计算一个反铁磁耦合体系，关键词为：
#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
##################

作者
Author: zjxitcc 时间: 2023-6-2 10:05
我猜可以算，但不支持stable=opt

作者
Author: diffcwmol 时间: 2023-6-2 10:20

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

谢谢哥们

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

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