sf(nstate=4)/bhhlyp/6-31g*
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.923944     0.000000    -1.237411
H      0.923944     0.000000    -1.237411
"""