Characteristics of thermionic electron emission during femtosecond laser ablation of gold film are studied numerically. For the rigorous calculation of electron and lattice temperatures, an enhanced two-temperature model with transient thermal and optical properties is developed and it is demonstrated that the model predicts the damage threshold fluences closely matching experimental data. From the calculated electron emission characteristics, quantum efficiency and electron current density are estimated. (C) 2009 Elsevier B.V. All rights reserved.