I am calculating OER and HER reaction pathway for photocatalytic acitivity of a 2D material with a supercell size of 2x2x1. Hence, I need to find free energy correction for several intermediate states. Using Dmol3, a geometry optimization followed by a frequency calculation gives the free energy values relatively quickly. In comparison, to obtain the same results using castep, geometry optimization followed by phonon calculation takes insanely longer time. I want to know which method is followed as standard i.e. gives more reliable values.