Rodrigo Negreiros, Jeffrey Peterson, Veronica Dexheimer

We study the cooling evolution of neutron stars with strong poloidal magnetic fields (with strength not far from observed values) using the full general relativity 2-dimensional \textit{Astreus} code, which solves consistently Einstein's and Maxwell's equations. We find that central magnetic fields with strengths $3-4\times10^{17}$ G, corresponding to surface magnetic fields $7-8\times10^{16}$, can significantly modify the cooling behavior of neutron stars, leading to stars with similar masses but different magnetic fields to exhibit different thermal evolution. We show a non-linear increase in the thermal relaxation time with increasing magnetic fields and that this behavior is associated with the reduction of the Direct Urca process in stars with strong magnetic fields. This is a novel result in which we can observe the magnetic field influence on the thermal evolution of stars, even if it is not strong enough to affect the Fermi distribution of particles.