Context <p>The formation of furan, a known carcinogen, in thermally processed foods represents a significant human health risk. Therefore, a clear mechanistic understanding of the transformation process of some common food components to yield furan/furan-derivatives is critical for the development of mitigation efforts or the prevention of furan formation in thermally processed foods. Therefore, in this study, density functional theoretical methods are applied in probing the possibility of furan formation from acetaldehyde, glycolaldehyde and aldotetrose on a step-by-step basis. The thermodynamic parameters associated with the transformation steps in the formation of furan from aldotetrose reveal an overall exergonic process, even at ambient temperature, limited by its second step where ΔG^‡₂₅ = 327.61&#xa0;kJ&#xa0;mol⁻¹. Furan is also generated from the reaction of acetaldehyde with glycolaldehyde via a complex multi-step process, limited by its first step where ΔG^‡₂₅ = 283&#xa0;kJ&#xa0;mol⁻¹. Based on these findings, it is clear that at high temperatures, furan formation is highly probable; hence, kinetic models governing the formation of furan/furan derivatives are proposed.</p> Methods <p>All calculations were carried out at the CAM-RB3LYP /6–311 +  + G(d,p) level of theory, as implemented in the Gaussian-16 software suite. </p>