First principles molecular dynamics simulations have been performed for the liquid-vapor interfaces of liquid Li, Mg, Al and Si. We analize the oscillatory ionic and valence electronic density profiles obtained, their wavelengths and the mechanisms behind their relative phase-shift.PACS numbers: 61.25. Mv, 64.70.Fx, 71.15.Pd X-ray reflectivity measurements on the surface of liquid metals and alloys, along with other techniques like diffuse scattering or grazing incidence diffraction, have shown the existence of layering in the ionic density profile. 1,2,3,4,5,6,7,8 Its origin is yet not clear and several reasons have been mooted. Rice and coworkers 9,10,11 have pointed that the reason for surface layering in metals is the interconnection between the ionic and electronic densities and that the abrupt decay of the electron density at the surface induces an effective wall against which the ions, behaving like hard-spheres, stack. Recently, it has been suggested that surface layering is a rather universal phenomenon, although in most cases it is frustrated by solidification; 12,13 therefore it only appears in systems whose melting temperature is very low compared with the critical temperature.The reflectivity experiments probe the total electronic density profile. Therefrom, the ionic density profile is derived by a superposition approximation where the the total electron density is taken as the sum of atomic-like electron densities around each of the nuclei in the system. Whereas this approach may be vindicated for the tightly bound core electrons, the case of valence electrons is more subtle. Early calculations of Lang and Kohn for semiinfinite step surfaces 14 showed that the valence electron density does not decay so abruptly, but displays some spill-out outside the surface; moreover inside the bulk its exhibits the so called Friedel oscillations.Computer simulations of liquid surfaces can evaluate directly the ionic density profile. But only a fully ab initio method can deliver an electronic density selfconsistent with the discrete nature of the ions. Orbital based ab initio simulations are scarce (just for Si 15 and Na 16 ) due to the huge computational demands they pose. Orbital free ab initio simulations are less demanding, although still expensive, and recent orbital free ab initio molecular dynamics (OFAIMD) calculations with 2000 and 3000 particles have studied the surface properties of liquid Li, Na, Na 0.3 K 0.7 and Li 0.4 Na 0.6 . 17,18 These studies showed that the superposition approximation produces a valence electronic density profile very similar to the fully selfconsistent one, except for the width of the interface due to the spill-out.This communication reports results for the liquidvapor interface of liquid Li, Mg, Al and Si near their respective melting points. The calculations were per- formed by the OFAIMD method where the forces acting on the nuclei are computed from electronic structure calculations, based on the density functional theory (DFT), which are performed as the ionic trajecto...