Role of turbulence and winding in the development of large-scale, strong magnetic fields in long-lived remnants of binary neutron star mergers

We perform a long and accurate large-eddy simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution, and the gradient subgrid-scale model allows us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-Helmholtz instability, they are coherent only over very short spatial scales until t ≳30 ms . Around that time, magnetic winding becomes more efficient leading to a linear growth of the toroidal component and slowly ordering the field to more axisymmetric large scales. The poloidal component only starts to grow at small scales at much later times t ≳90 ms , in a way compatible with the magneto-rotational instability. No strong large-scale poloidal field or jet is produced in the timescales spanned by our simulation, although there is a helicoidal structure gradually developing at late times. We highlight that soon after the merger the magnetic field topology is always strongly dominated by toroidal structures, with a complex distribution in the meridional plane and highly turbulent perturbations. The topology present before the merger is completely erased and forgotten, as long as the initial magnitude of the magnetic field has realistic values. In fact, we advise against the widespread use of unrealistic very strong large-scale initial magnetic fields, since they might affect significantly the postmerger topological outcome, which would be therefore strongly biased by such an artificial choice. Finally, we confirm the universality of the evolved topology, even when starting with very different magnetic fields confined to the outermost layers of each neutron star.