In this contribution, the thermodynamic analysis of α- and β-pinene epoxide isomerization over Fe and Cu supported on MCM-41 is presented using computational chemistry and group contribution methods (GCMs). Some physical–chemical data (Tc, Pc, vc, Zc, ω, Tb, Tfus) and thermodynamic (S°298.15, Cp,298.15°, Cv,298.15°, ΔHf,298.15°, ΔGf,298.15°, ΔHvb°, ΔHfus, CpL) properties obtained by different GCMs are reported for several monoterpenes and monoterpenoids, which significantly contribute to the knowledge of the properties of these compounds. Density functional theory (DFT), PBE-D3/6-311G(d,p), was employed for determining the Gibbs free energy and the heat of reaction associated with the transformation of monoterpene epoxides into aldehydes, ketones, and related oxygenated compounds in the presence of different solvents and at several temperatures. The calculations were compared with available data reported and the experimental results of the catalytic reactions. The transformation of α- and β-pinene epoxides into aldehydes appears to be more spontaneous and favorable than their transformations into alcohols in a wide range of temperatures. These results are in agreement with the experiments over Fe/MCM-41 and Cu/MCM-41, where α-pinene epoxide isomerization yields campholenic aldehyde (50–80% selectivity) as the main product. The 1.7Fe/MCM-41 material was more active in all solvents than 1.3Cu/MCM-41 for both α- and β-pinene epoxide isomerization. However, perillyl alcohol (20–70% selectivity) was the most favored for the isomerization reaction, except when ethyl acetate was the solvent. Enthalpy and Gibbs free energy of the studied reactions estimated by both GCMs and DFT calculations did not show large differences for most of the reactions at evaluated temperatures.