Evolution of neutron stars with toroidal magnetic fields: Axisymmetric simulation in full general relativity

Abstract

We study the stability of neutron stars with toroidal magnetic fields by magnetohydrodynamic simulation in full general relativity under assumption of axial symmetry. Nonrotating and rigidly rotating neutron stars are prepared for a variety of magnetic field configuration. For modeling the neutron stars, the polytropic equation of state with the adiabatic index $\Gamma =2$ is used for simplicity. It is found that nonrotating neutron stars are dynamically unstable for the case that toroidal magnetic field strength varies $\propto {\varpi }^{2k-1}$ with $k\ge 2$ (here $\varpi$ is the cylindrical radius), whereas for $k=1$ the neutron stars are stable. After the onset of the instability, unstable modes grow approximately in the Alfvén time scale and, as a result, a convective motion is excited to change the magnetic field profile until a new state, which is stable against axisymmetric perturbation, is reached. We also find that rotation plays a role in stabilization, although the instability still sets in in the Alfvén time scale when the ratio of magnetic energy to rotational kinetic energy is larger than a critical value $\sim 0.2$. Implication for the evolution of magnetized protoneutron stars is discussed.

accepted to PRD

References in corpus (6)

• Evolution of magnetized, differentially rotating neutron stars: Simulations in full general relativity
• Magnetorotational collapse of massive stellar cores to neutron stars: Simulations in full general relativity
• The Proto-neutron Star Phase of the Collapsar Model and the Route to Long-soft Gamma-ray Bursts and Hypernovae
• Magnetically-driven explosions of rapidly-rotating white dwarfs following Accretion-Induced Collapse
• General relativistic simulations of pasive-magneto-rotational core collapse with microphysics
• Origin of neutron star magnetic fields