This article presents a novel coupling of three state-of-the-art particle methods, capable of simulating heat transfer in multiphase flows with minimum error. The idea originates from increasing the robustness of the geometric features required to calculate physical quantities, such as the boundary conditions in the heat equation and the surface tension force in Navier-Stokes equations. The workability of this new approach for handling thermal effects in complex moving interfaces is examined with a numerical test case wherein multiple scenarios such as internal flow, bursting at free surface as well as bubble evolutionary motions occur. Performing a benchmark comparison with the reference model implemented in COMSOL Multiphysics®, a markedly improved mass conservation using the proposed particle-based solver is observed. This, in turn, results in a more realistic and accurate approximation of the temperature.