This work presents a Lagrangian meshless method for modeling transient heat transfer and thermo-capillary effects at the interface of two-phase incompressible fluids. To simulate the thermo-capillary effect, we add the Marangoni force to the continuum surface force (CSF) model, by considering the forces caused by the surface tension gradient which acts tangentially to the interface position. In the current study, we extend our previously proposed model which uses a fully explicit incompressible smoothed particle hydrodynamics (EISPH) approach along with corrected SPH and VKF kernel function, to obtain stable and accurate results in non-isothermal single and multi-phase problems.The model is validated for several test cases, including transient heat conduction, natural convection heat transfer, and thermo-capillary droplet migration, against conventional numerical methods. The validated model is used to study the effects of several dimensionless parameters on thermo-capillary droplet migration. The results show that the proposed EISPH method can model accurately complex heat transfer problems such as thermo-capillary induced motion.