Here we report on experimental studies of femtosecond laser induced surface metal alloying. We demonstrate that layers of different metals can be mixed in a certain range of laser pulse energies. Numeric simulations demonstrate that the sub-surface melting and mixing is advantaged through the difference in the electron-phonon coupling constants of the metals in the multi-layer system. Dependence of the depth of the mixed layer on the number of laser pulses per unit area is studied. Numeric simulations illustrate physical picture of the laser alloying process. (C) 2011 Elsevier B.V. All rights reserved.