, where ∆𝐺𝐺 complex corr and ∆𝐺𝐺 𝐴𝐴 corr denote thermal correction to free energy of complex and fragment A evaluated at their respective equilibrium geometries, respectively. It is noticed that Gaussian program employs rigid-rotor harmonic oscillator (RRHO) model to evaluate entropy, while we strongly suggest using Grimme's quasi-RRHO model implemented in our Shermo code instead to calculate the free energy correction when low-frequencies occur, as in this case RRHO may severely overestimate their entropy contribution. If taking all terms mentioned above into account, binding free energy at solvent environment can be represented as