Abstract Magnetars, neutron stars thought to be with ultrastrong magnetic fields of 10 14–15 G, are observed much hotter than ordinary pulsars ∼10 12 and additional heating sources required. One possibility is by the ambipolar diffusion in stellar core. This scenario examined calculating models using relativistic thermal evolutionary code without making isothermal approximation. The results sho...