xmcqdpt2/6-31g* cas(active_electrons=2,active_orbitals=2,frozen_core=7) ci(nroot=3,target_spin=singlet) pt2(target_roots="0,1,2",nroot=3,edshft=0.02)
geom="""
C      0.000000     0.000000     0.664311
C      0.000000     0.000000    -0.664311
H      0.000000     0.923944     1.237411
H      0.000000    -0.923944     1.237411
H      0.000000     0.923944    -1.237411
H      0.000000    -0.923944    -1.237411
"""