AbstractMolecular dynamics method was applied to simulate the physical properties of liquid metals: lead, bismuth and a binary alloy-- lead-bismuth eutectic (LBE). The embedded atom method (EAM), an empirical model rooted in density-functional theory, was used to represent the many-body interaction within the liquid metals. The atomic-scale interactions, structure and thermal physical properties of lead, bismuth and LBE were obtained through the simulation, and then compared to the available experimental results. The theoretical results of the physical properties calculated through the MD simulations are in good agreements with the available experimental data.