“…Models of surface temperature distribution with purely poloidal magnetic fields (Page 1995;Geppert et al 2004) do predict non-uniform surface temperature distributions, , but such inhomogeneities are not strong enough to produce such small X-ray emitting regions surrounded by large cold regions detectable in the optical band as observed. However, inclusion of a toroidal component of the magnetic field, confined to the neutron star crust, has a dramatic effect (Pérez-Azorín et al 2006;Geppert et al 2006; see also Pons 2007;Page 2007): this field component inhibits heat from the stellar core to flow to the surface through most of the crust, except for small domains surrounding the magnetic axis, and results in highly non-uniform surface temperature distributions producing good fits to the observed thermal spectra, from the optical up to the X-ray band. These models of small hot regions, detected in the Xray band, surrounded by large cold regions, detected in the optical band, which allow to reproduce the entire observed thermal spectrum and results in large radii for the star are in contradiction with the proposed strange star interpretation of RX J1856.5-3754, which was based on the small radius detected in the X-ray band.…”