Let's say there is a rxn : A+B -> C+D
i) Enthalpy of reaction/gibbs free energy of a reaction can be calculated using the thermochemistry output from Gaussian.
deltaG = sum(Gc +Gd) - sum(Ga +Gb).
ii) On the other hand, if we use IRC. We can obtain the reaction mechanism (A+B) --> TS --> (C+D)
-----------------------------------------------------------
The results I obtained from
i) showed that the reaction is thermodynamically unfavourable. (deltaG > 0)
ii) showed that (C+D) is at a lower energy than (A+B) which means the reaction is thermodynamically favourable.
-----------------------------------------------------------
I understand that by using i) to determine the reaction feasibility, we assume that A and B are located far enough to have the interaction neglected. Same as C and D. Will this incur some errors here?
As far as I am concerned, the distance of A and B we set in Gaussian has an effect on the output energy.
So in this case, IRC is a better representation of the reaction?
Another question, for example I want to know whether the synthesis of the new chemical is feasible or not. In this case is it appropriate to use the deltaG = sum(Gc +Gd) - sum(Ga +Gb) to get a quick answer to the feasibility of the reaction using Gaussian.
Thank you.