Of all the metals, gold is the most well-known and widely used material in scientific research, industrial production and, more recently, in biomedicine problems. In the temperature range 300 ≤ T ≤ 2000 K, including the region of the meltingcrystallization phase transition, the results of modeling the phonon thermal conductivity of gold are presented. Phonon thermal conductivity plays an important role in modeling the mechanisms of interaction of pulsed laser radiation with gold in the framework of the two-temperature continuum model. In the region of the phase transition, overheating-undercooling of the solid phase occurs, the substance changes its structure. These phenomena are associated with changes in the phonon subsystem of gold, therefore, for mathematical modeling of heatingcooling, it is necessary to know the characteristic of heat transfer as the thermal conductivity of the phonon subsystem of gold. Obtaining the temperature dependence of phonon thermal conductivity in such a wide temperature range from experiment is problematic. In this work, phonon thermal conductivity was obtained by the direct non-equilibrium method in the framework of molecular dynamics modeling using the EAM potential.