We calculate the mass-radius relationship of quark stars with the magnetized densitydependent quark mass model in this work, considering two magnetic field geometries: a statistically isotropic, tangled field and a force-free configuration. In both cases, magnetic field production decreases in the case of maximum quark star mass. Furthermore, a tangled, isotropic magnetic field has a relatively smaller impact on the mass and radius, compared to the force-free configuration, which implies that the geometry of the interior magnetic field is at least as important as the field strength itself when the influence of the strong magnetic field on the mass and radius is assessed.
Rotochemical heating originates in the deviation from beta equilibrium due to spin-down compression, which is closely related to the dipole magnetic field. We numerically calculate the deviation from chemical equilibrium and thermal evolution of neutron stars with decaying magnetic fields. We find that the power-law long term decay of the magnetic field slightly affects the deviation from chemical equilibrium and surface temperature. However, the magnetic decay leads to older neutron stars that could have a different surface temperature with the same magnetic field strength. That is, older neutron stars with a low magnetic field(108G) could have a lower temperature even with rotochemical heating in operation, which probably explains the lack of other observations on older millisecond pulsars with higher surface temperature,except millisecond pulsar J0437–4715.