mrsf(nstate=70)/bhhlyp/6-31g* energy(S0) tdhf(ixcore=2)
geom="""
C      0.000000     0.000000    -0.505994
N      0.000000     0.000000     0.657328
H      0.000000     0.000000    -1.565330
"""