<div dir="ltr"><br>Dear all,<br><br>Over the past decades, significant effort has been devoted to<br>understanding the evolution of isolated neutron stars during their<br>lifetimes, which presents challenges from both physical and numerical<br>perspectives. The temperature and magnetic field of neutron stars<br>evolve together in a complex way. Extremely strong magnetic fields<br>influence how heat flows through the star’s solid crust, channeling<br>thermal energy along preferred directions and insulating some regions<br>while exposing others. At the same time, these magnetic fields are not<br>static: they can change their geometry due to the Hall drift,  and as<br>they gradually decay, they release energy that can partially reheat<br>the star. Understanding this coupled magnetothermal evolution is<br>therefore essential to correctly interpret neutron-star observations,<br>and especially their apparent diversity.<br><br>Several magneto-thermal codes in the past had to choose to either<br>employ detailed microphysics in simplified (axially symmetric)<br>geometries or addressing fully three-dimensional configurations at the<br>expense of using more approximate physical assumptions. In the past<br>few years, we tried to overcome this by developing MATINS (MAgneto<br>Thermal evolutIon of Isolated NSs), funded by the European Research<br>Council (ERC 'MAGNESIA' No. 817661): a three-dimensional numerical<br>framework designed to model the magnetothermal evolution of the crust<br>of isolated neutron stars. MATINS computes the neutron star structure<br>by solving the Tolman–Oppenheimer–Volkoff equations for a range of<br>nuclear matter equations of state, enabling systematic studies across<br>different stellar masses. The code solves the induction equation in<br>the crust, including both Ohmic dissipation and Hall drift, to<br>describe magnetic-field evolution. This is coupled to a fully 3D<br>thermal cooling model based on state-of-the-art microphysical inputs:<br>conductivities, heat capacities, neutrino processes, with the options<br>of different superfluid and superconductive gap models. To solve the<br>governing equations in three-dimensional spherical geometry, MATINS<br>employs a finite-volume scheme on a cubed-sphere grid. Each radial<br>shell is covered by six smoothly connected patches, analogous to<br>inflating a cube into a sphere, thereby avoiding the coordinate<br>singularities at the axis of standard spherical coordinates. This new<br>framework stems directly from the 2D magneto-thermal code developed in<br>the past 20 years at the University of Alicante, adapted to a 3D grid<br>and with additional features including parallelization with openMP.<br>The code can be used to model a variety of observable properties of<br>isolated neutron stars, including X-ray thermal emission, surface<br>magnetic fields, temperature maps, and rotational evolution (spin<br>period and its derivative).<br><br>We are happy to announce the public release of the MATINS code that<br>can be now accessed at:<br><a href="https://ice-csic-astroexotic.github.io/code/matins">https://ice-csic-astroexotic.github.io/code/matins</a> .<br><br>Details on how to install and run the code are also available at the<br>above website, together with the reference articles explaining more<br>technical details. For any queries you can write at:<br><a href="mailto:matins@ice.csis.es">matins@ice.csis.es</a> . Apologies if you receive this message more than<br>once.<br><br>Best wishes,<br><br>Stefano Ascenzi<br>Clara Dehman<br>Davide De Grandis<br>José A. Pons<br>Nanda Rea<br>Daniele Viganò<br><br><br><br>______________________________________________________<br><br>Prof. Nanda Rea  | CSIC Research Professor<br>Instituto de Ciencias del Espacio (ICE, CSIC)<br>Institut d'Estudis Espacials de Catalunya (IEEC)<br>Email: <a href="mailto:rea@ice.csic.es">rea@ice.csic.es</a>; Tel: +34937379788 (ext: 430423)<br>_______________________________________________________</div>