A concurrent multi-scale modeling method for finite temperature simulation of solids is introduced. The objective is to represent far from equilibrium phenomena using an atomistic model and near equilibrium phenomena using a continuum model, the domain being partitioned in discrete and continuum regions, respectively. An interface sub-domain is defined between the two regions to provide proper coupling between the discrete and continuum models. While in the discrete region the thermal and mechanical processes are intrinsically coupled, in the continuum region they are treated separately. The interface region partitions the energy transferred from the discrete to the continuum into mechanical and thermal components by splitting the phonon spectrum into “low” and “high” frequency ranges. This is achieved by using the generalized Langevin equation as the equation of motion for atoms in the interface region. The threshold frequency is selected such to minimize energy transfer between the mechanical and thermal components. Mechanical coupling is performed by requiring the continuum degrees of freedom (nodes) to follow the averaged motion of the atoms. Thermal coupling is ensured by imposing a flux input to the atomistic region and using a temperature boundary condition for continuum. This makes possible a thermodynamically consistent, bi-directional coupling of the two models.