Simulating the long-term evolution of temperature and magnetic fields in neutron stars is a major effort in astrophysics, having significant impact in several topics. A detailed evolutionary model requires, at the same time, the numerical solution of the heat diffusion equation, the use of appropriate numerical methods to control non-linear terms in the induction equation, and the local calculation of realistic microphysics coefficients. Here we present the latest extension of the magneto-thermal 2D code in which we have coupled the crustal evolution to the core evolution, including ambipolar diffusion. It has also gained in modularity, accuracy, and efficiency. We revise the most suitable numerical methods to accurately simulate magnetar-like magnetic fields, reproducing the Hall-driven magnetic discontinuities. From the point of view of computational performance, most of the load falls on the calculation of microphysics coefficients. To a lesser extent, the thermal evolution part is also computationally expensive because it requires large matrix inversions due to the use of an implicit method. We show two representative case studies: (i) a non-trivial multipolar configuration confined to the crust, displaying long-lived small-scale structures and discontinuities; and (ii) a preliminary study of ambipolar diffusion in normal matter. The latter acts on timescales that are too long to have relevant effects on the timescales of interest but sets the stage for future works where superfluid and superconductivity need to be included.

模拟中子星中温度和磁场的长期演变是天体物理学的一项重大工作，在几个主题中产生了重大影响。同时，详细的演化模型需要热扩散方程的数值解，使用适当的数值方法来控制感应方程中的非线性项以及对实际的微物理学系数进行局部计算。在这里，我们介绍了磁热二维代码的最新扩展，其中将地壳演化与核心演化（包括双极性扩散）耦合在一起。它还在模块化，准确性和效率方面得到了提高。我们修改了最合适的数值方法，以准确地模拟类似磁场的磁场，从而再现霍尔驱动的磁间断。从计算性能的角度来看，大部分负载都落在微物理系数的计算上。在较小程度上，热演化部分在计算上也很昂贵，因为由于使用了隐式方法，因此需要大的矩阵求逆。我们展示了两个代表性的案例研究：（i）局限在地壳中的非平凡多极构造，显示了长寿命的小规模结构和不连续性；（ii）对正常物质中双极性扩散的初步研究。后者作用的时间尺度太长，无法对感兴趣的时间尺度产生相关影响，但为将来的工作奠定了基础，在未来的工作中需要包括超流体和超导电性。热演化部分在计算上也很昂贵，因为由于使用了隐式方法，因此需要大的矩阵求逆。我们展示了两个代表性的案例研究：（i）局限在地壳中的非平凡多极构造，显示了长寿命的小规模结构和不连续性；（ii）对正常物质中双极性扩散的初步研究。后者作用的时间尺度太长，无法对感兴趣的时间尺度产生相关影响，但为将来的工作奠定了基础，在未来的工作中需要包括超流体和超导电性。热演化部分在计算上也很昂贵，因为由于使用了隐式方法，因此需要大的矩阵求逆。我们展示了两个代表性的案例研究：（i）局限在地壳中的非平凡多极构造，显示了长期存在的小规模结构和不连续性；（ii）对正常物质中双极性扩散的初步研究。后者作用的时间尺度太长，无法对感兴趣的时间尺度产生相关影响，但为将来的工作奠定了基础，在未来的工作中需要包括超流体和超导电性。显示长期存在的小规模结构和不连续性；（ii）对正常物质中双极性扩散的初步研究。后者作用的时间尺度太长，无法对感兴趣的时间尺度产生相关影响，但为将来的工作奠定了基础，在未来的工作中需要包括超流体和超导电性。显示长期存在的小规模结构和不连续性；（ii）对正常物质中双极性扩散的初步研究。后者作用的时间尺度太长，无法对感兴趣的时间尺度产生相关影响，但为将来的工作奠定了基础，在未来的工作中需要包括超流体和超导电性。