Development of the temperature-dependent interatomic potential for molecular dynamics simulation of metal irradiated with an ultrashort pulse laser

Abstract

Laser ablation is often explained by a two-temperature model (TTM) with different electron and lattice temperatures. To realize a classical molecular dynamics simulation of the TTM, we propose an extension of the embedded atom method to construct an interatomic potential that is dependent on the electron temperature. This method is applied to copper, and its validity is demonstrated by comparison of several physical properties, such as the energy–volume curve, phonon dispersion, electronic heat capacity, ablation threshold, and mean square displacement of atoms, with those of finite-temperature density functional theory.