The total number of geometry cycles required to optimize the same set of molecules Here ASE/Berny and ASE/Sella use redundant internal coordinates, while geomeTRIC uses translation-rotational-internal coordinates, and DL-FIND uses delocalized internal coordinates. Except for ASE/Berny and geomeTRIC, which use a modified version of the original Schlegel update process, as well as ASE/Sella, all methods initialize the Hessian to be the unit matrix. These results reconfirm the two well-known observations that i) internal coordinates outperform the Cartesian coordinate system for molecular geometry optimizations, and ii) good approximations to the initial Hessian are also important to reduce the number of required steps.