Geometry Optimization: A Comparison of Different Open-Source Geometry Optimizers

17 February 2023, Version 2
This content is a preprint and has not undergone peer review at the time of posting.


Based on a series of energy minimizations with starting structures obtained from the Baker test set of 30 organic molecules, a comparison is made between various open source geometry optimization codes that are interfaced with the open-source QUantum Interaction Computational Kernel (QUICK) program for gradient and energy calculations. The findings demonstrate how the choice of the coordinate system influences the optimization process to reach an equilibrium structure. With fewer steps, internal coordinates outperform Cartesian coordinates while the choice of the initial Hessian and Hessian update method in quasi-Newton approaches made by different optimization algorithms also contributes to the rate of convergence. Overall, DL-FIND, geomeTRIC and ASE, using Hartree--Fock theory with the 6-31G** basis set, needed a comparable number of geometry optimization steps to the approach of Baker using a unit matrix as the initial Hessian to reach the optimized geometry. Due to its performance and the fact that DL-FIND is feature rich we use itas the default optimizer for QUICK.


open-source softwares
geometry optimization

Supplementary materials

Input Files to different Geometry Optimizers
This folder contains input files for different geometry optimizers reported in the paper.


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.