Charge order at magnetite Fe3O4(001): surface and Verwey phase transitions

At ambient conditions, the Fe3O4(0 0 1) surface shows a (root 2 x root 2) R45 degrees reconstruction that has been proposed as the surface analog of the bulk phase below the Verwey transition temperature, T-V. The reconstruction disappears at a high temperature, T-S, through a second order transition. We calculate the temperature evolution of the surface electronic structure based on a reduced bulk unit cell of P2/m symmetry that contains the main features of the bulk charge distribution. We demonstrate that the insulating surface gap arises from the large demand of charge of the surface O, at difference with that of the bulk. Furthermore, it is coupled to a significant restructuration that inhibits the formation of trimerons at the surface. An alternative bipolaronic charge distribution emerges below T-S, introducing a competition between surface and bulk charge orders below T-V.