#! H2O CISD/6-31G** Optimize Geometry by Energies

memory 250 mb


molecule h2o {
    O
    H 1 1.00
    H 1 1.00 2 103.1
}

set globals {
  basis 6-31G**
  max_force_g_convergence 0.000001
  max_disp_g_convergence  0.0001
}

thisenergy = optimize('cisd', dertype = 0)

