关于激发态的优化求助

lastzealot

Sob老师您好，在一次C6H5N双自由基的激发态优化的计算中每次都出错我不知道是为什么，麻烦老师帮忙看一下

（这个结构是Ph-N:，C2v点群（A1，B1，B2，A2），基态是b1轨道和b2轨道各一个单电子，即波函数对称性为B1*B2=3A2。我想计算的电子组态是1A2，即在不可约表示B2和A2的轨道中分别有一个单电子且自旋相反）

谢谢老师

输入文件
***, c6h5n
geomtyp=xyz
geometry={
C          0.0000000000       -1.2102338407       -1.0503576403
C          0.0000000000       -1.2251057123        0.3285536652
C          0.0000000000        0.0000000000        1.0566247904
C          0.0000000000        1.2251057123        0.3285536652
C          0.0000000000        1.2102338407       -1.0503576403
C          0.0000000000        0.0000000000       -1.7471550171
H          0.0000000000       -2.1449242583       -1.5957303914
H          0.0000000000       -2.1538392973        0.8815671517
H          0.0000000000        2.1538392973        0.8815671517
H          0.0000000000        2.1449242583       -1.5957303914
H          0.0000000000        0.0000000000       -2.8287071576
N          0.0000000000        0.0000000000        2.3840718537}
basis=cc-pvtz
{hf
wf, 48, 4, 0
print, orbital=5}
{multi
wf, 48,4,0
}
optg


输出文件末尾
PROGRAM * RHF-SCF (CLOSED SHELL)       Authors: W. Meyer, H.-J. Werner


NUMBER OF ELECTRONS:      24+   24-    SPACE SYMMETRY=4    SPIN SYMMETRY: Singlet
CONVERGENCE THRESHOLDS:    1.00E-06 (Density)    1.00E-07 (Energy)
MAX. NUMBER OF ITERATIONS:       60
INTERPOLATION TYPE:            DIIS
INTERPOLATION STEPS:              2 (START)      1 (STEP)
LEVEL SHIFTS:                  0.00 (CLOSED)  0.00 (OPEN)


SYYMETRY= 4  MS2= 0
? Error
? Open-shell singlet not implemented
? The problem occurs in scfpro

ymygca

双自由基是open-shell system，
报错信息告诉你了Open-shell singlet not implemented

lastzealot

ymygca 发表于 2019-10-23 01:55
双自由基是open-shell system，
报错信息告诉你了Open-shell singlet not implemented

我看到了，但是我想知道对于开壳层的单线态应该如何设置才能不出问题