Composition-dependent melting behaviour of NaxK55–x core–shell nanoalloys

Molecular dynamics simulations at constant temperature are performed to investigate melting-like transition in Na₁₃K₄₂, Na₁₉K₃₆ and Na₂₆K₂₉ nanoalloys using a second-moment-approximation tight-binding analytic potential to calculate the forces on the constituent atoms. A weighted histogram analysis method is employed to remove non-ergodicity issues due to the complex potential energy surface of these nanoalloys. The heat capacity shows three distinctive steps in melting for Na₁₃K₄₂, while Na₂₆K₂₉ and Na₁₉K₃₆ have two-step and one-step melting transition, respectively. The steepest descent method is used to quench the configurations in a given interval during the simulation and also study the isomerisation processes occurring at different temperatures. Analysing the configuration energies of quenched structures for the entire nanoalloy and the core atoms separately gives more details about the melting mechanism. The Lindemann parameter is also calculated at several temperatures during the simulation which shows a gradual increase for Na₁₃K₄₂ and Na₂₆K₂₉ while a sharp change is observed for Na₁₉K₃₆. These findings are in agreement with the multi-step nature of the phase transition in Na₁₃K₄₂ and Na₂₆K₂₉ and one-step melting of the Na₁₉K₃₆ magic composition.