“…An accurate description of these phenomena is required and is becoming a primary concern for industrial applications. The heating rate can be computed with more sophisticated hydrodynamic models which are able to describe off-equilibrium regimes (Muscato et al, 1998;Muscato and Di Stefano, 2008, 2011a, c, d, 2012Di Stefano and Muscato, 2012) but we shall not pursue this subject in the sequel. Another possibility for computing a detailed picture of energy dissipation, is to adapt the well known Monte Carlo (MC) simulation method, which has been originally developed for studying hot electron effects when the lattice is considered a thermal reservoir (Muscato, 2000;Majorana et al, 2004;Muscato and Wagner, 2005;Muscato et al, 2010Muscato et al, , 2011).…”