Non-equilibrium electron and phonon dynamics in metals under femtosecond laser pulses

Different models for relaxation dynamics of electrons and phonons in a thin metal film heated by femto-pico second laser pulses have been discussed. The traditional two-temperature approach reveals to be inaccurate due to deviations of electrons and phonons from Fermi-Dirac and Bose-Einstein distributions, respectively. Coupled Boltzmann kinetic equations have been adapted for the quantum statistics to study the energy distribution of electrons and phonons in metals. Theoretical details have been discussed and a new solution method has been proposed overcoming numerical problems and improving stability, allowing the study of the dynamics until the complete relaxation. Numerical results have been compared with photoemission spectroscopy experimental data.