For all topics about errors and inefficiencies with the Psi4 geometry optimizers: if you can’t get a geometry to converge, if a geometry optimization takes suspiciously bad steps, or geometry-specific errors. Errors with electronic structure computations during an optimization belong elsewhere.
How can I calculate potential energy surface with help of PSI4?
Please suggest me any direct method or by python script method.