The melting mechanism in binary Pd_0.25Ni_0.75 nanoparticles: molecular dynamics simulations

Abstract

The melting mechanism for Pd_0.25Ni_0.75 alloy nanoparticles (NPs) was investigated using molecular dynamics (MD) simulations with quantum Sutton-Chen many-body potentials. NPs of six different sizes ranging from 682 to 22,242 atoms were studied to observe the effect of size on the melting point. The melting temperatures of the NPs were estimated by following the changes in both the thermodynamic and structural quantities such as the total energy, heat capacity and Lindemann index. We also used a thermodynamics model to better estimate the melting point and to check the accuracy of MD simulations. We observed that the melting points of the NPs decreased as their sizes decreased. Although the MD simulations for the bulk system yielded higher melting temperatures because of the lack of a seed for the liquid phase, the melting temperatures determined for both the bulk material and the NPs are in good agreement with those predicted from the thermodynamics model. The melting mechanism proceeds in two steps: firstly, a liquid-like shell is formed in the outer regions of the NP with increasing temperature. The thickness of the liquid-like shell increases with increasing temperature until the shell reaches a critical thickness. Then, the entire Pd–Ni NP including core-related solid-like regions melts at once.

Keywords:

• Melting mechanism
• classical molecular dynamics simulations
• bond-order parameter
• Lindemann index
• nanoparticle
• self-diffusion coefficient

Disclosure statement

No potential conflict of interest was reported by the authors.

Acknowledgments

The authors would like to thank M. Colakogullari (Trakya University) for sharing his experiences on the molecular dynamics and useful discussions.