Molecular dynamics simulations of the liquid-vapor interfaces in simple sp-bonded liquid metals have been performed using first-principles methods. Results are presented for liquid Li, Na, K, Rb, Cs, Mg, Ba, Al, Tl, and Si at thermodynamic conditions near their respective triple points, for samples of 2000 particles in a slab geometry. The longitudinal ionic density profiles exhibit a pronounced stratification extending several atomic diameters into the bulk, which is a feature already experimentally observed in liquid K, Ga, In, Sn, and Hg. The wavelength of the ionic oscillations shows a good scaling with the radii of the associated Wigner-Seitz spheres. The structural rearrangements at the interface are analyzed in terms of the transverse pair correlation function, the coordination number, and the bond-angle distribution between nearest neighbors. The valence electronic density profile also shows (weaker) oscillations whose phase, with respect to those of the ionic profile, changes from opposite phase in the alkalis to almost in-phase for Si.

• Received 30 November 2005

DOI:https://doi.org/10.1103/PhysRevB.74.014207