The geometry optimization requires calculations of energy derivatives with respect to nuclear coordinates of molecules. The finite difference method is one approach for energy derivative calculations.