The two-temperature model has been widely used to predict the electron and phonon temperature distributions in ultra-short laser processing of metals. This paper extends the use of existing two-temperature model to high electron temperatures by using full-run quantum treatments, including the electron heat capacity, electron relaxation time and electron conductivity . A dual-hyperbolic two temperature model is solved numerically using an accurate and stabile forward-difference scheme, in which one kind of artificial viscosities and adaptive time steps method are introduced. This model predicts the damage thresholds more accurately compared with published experimental results.© (2007) COPYRIGHT SPIE--The International Society for Optical Engineering. Downloading of the abstract is permitted for personal use only.