mrsf(nstate=10)/bhhlyp/6-31+g* ekt(ip=false,ea=true)
geom="""
O      0.000000     0.000000     0.444340
O      0.000000     1.083397    -0.222170
O      0.000000    -1.083397    -0.222170
"""