1 Introduction Geometry optimization is an essential technique for studying the chemical properties of molecules. Conventional geometry optimization uses energy gradients obtained from quantum mechanical (QM) calculations 1-3. Geometry optimization using energy gradients is applicable not only to the global minimum (GM) but also to local minima (LMs), transition states (TSs), and conical intersections (CIs).