  We simulate the liquid silicon surface with first-principles molecular
dynamics in a slab geometry. We find that the atom-density profile presents a
pronounced layering, similar to those observed in low-temperature liquid metals
like Ga and Hg. The depth-dependent pair correlation function shows that the
effect originates from directional bonding of Si atoms at the surface, and
propagates into the bulk. The layering has no major effects in the electronic
and dynamical properties of the system, that are very similar to those of bulk
liquid Si. To our knowledge, this is the first study of a liquid surface by
first-principles molecular dynamics.
