The following examples illustrate cooling curves computed with the BSK24 equation of state over one million year for a neutron star of 1.6 M⊙, under varying surface magnetic field strengths ranging from 1012 G to 5×1015 G. The models assume a simple magnetic field configuration consisting of a poloidal dipole plus a toroidal quadrupole component, and include two types of magnetized envelope compositions: heavy-element and light-element envelopes.
The cooling curves and evolutionary track data are available for download here: