Optical properties of real surfaces: Local-field effects at oxidized Si(100)(2x2) computed with an efficient numerical scheme

We show the application of an efficient numerical scheme to obtain the independent-particle dynamic polarizability matrix chi((0))(r,r(\textquoteright),omega), a key quantity in modern ab initio excited-state calculations. The method has been applied to the study of the optical response of a realistic oxidized silicon surface, including the effects of crystal local fields. The latter are shown to substantially increase the surface optical anisotropy in the energy range below the bulk band gap. Our implementation in a large-scale ab initio computational code allows us to make a quantitative study of the CPU time scaling with respect to the system size, and demonstrates the real potential of the method for the study of excited states in large systems.