Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.
Magnetic fields are observed in massive Ap/Bp stars and are presumably present in the radiative zone of solar-like stars. To date, there is no clear understanding of the dynamics of the magnetic field in stably stratified layers. A purely toroidal magnetic field configuration is known to be unstable, developing mainly non-axisymmetric modes. Rotation and a poloidal field component may lead to stabilization. Here we perform global MHD simulations with the EULAG-MHD code to explore the evolution of a toroidal magnetic field located in a layer whose Brunt-Väisälä frequency resembles the lower solar tachocline. Our numerical experiments allow us to explore the initial unstable phase as well as the long-term evolution of such field. During the first Alfven cycles, we observe the development of the Tayler instability with the prominent longitudinal wavenumber, m = 1. Rotation decreases the growth rate of the instability and eventually suppresses it. However, after a stable phase, energy surges lead to the development of higher-order modes even for fast rotation. These modes extract energy from the initial toroidal field. Nevertheless, our results show that sufficiently fast rotation leads to a lower saturation energy of the unstable modes, resulting in a magnetic topology with only a small fraction of poloidal field, which remains steady for several hundreds of Alfven traveltimes. The system then becomes turbulent and the field is prone to turbulent diffusion. The final toroidal–poloidal configuration of the magnetic field may represent an important aspect of the field generation and evolution in stably stratified layers.
more » « less- NSF-PAR ID:
- 10401254
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 521
- Issue:
- 1
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 1415-1428
- Size(s):
- p. 1415-1428
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
Abstract We provide exact analytical solutions for the magnetic field produced by prescribed current distributions located inside a toroidal filament of finite thickness. The solutions are expressed in terms of toroidal functions, which are modifications of the Legendre functions. In application to the MHD equilibrium of a twisted toroidal current loop in the solar corona, the Grad–Shafranov equation is decomposed into an analytic solution describing an equilibrium configuration against the pinch-effect from its own current and an approximate solution for an external strapping field to balance the hoop force. Our solutions can be employed in numerical simulations of coronal mass ejections (CMEs). When superimposed on the background solar coronal magnetic field, the excess magnetic energy of the twisted current loop configuration can be made unstable by applying flux cancellation to reduce the strapping field. Such loss of stability accompanied by the formation of an expanding flux rope is typical for the Titov & Démoulin eruptive event generator. The main new features of the proposed model are as follows: the filament is filled with finite
β plasma with finite mass and energy, the model describes an equilibrium solution that will spontaneously erupt due to magnetic reconnection of the strapping magnetic field arcade, and there are analytic expressions connecting the model parameters to the asymptotic velocity and total mass of the resulting CME, providing a way to connect the simulated CME properties to multipoint coronograph observations. -
more » « less
ABSTRACT Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean-field dynamo theories used to model the growth of large-scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coefficients with successively higher order spatial derivatives of the mean magnetic field. One ingredient of standard models is the kinematic coefficient of the zeroth-order term, namely the averaged kinetic pseudo-tensor $\boldsymbol \alpha$, that converts toroidal to poloidal fields. Here we demonstrate an efficient way to calculate this quantity for rotating stratified turbulence, whereby the pre-averaged quantity is calculated for the motion of a single plume, and the average is then taken over an ensemble of plumes of different orientations. We calculate the plume dynamics in the most convenient frame, before transforming back to the lab frame and averaging. Our concise configuration space calculation gives essentially identical results to previous lengthier approaches. The present application exemplifies what is a broadly applicable method.
-
more » « less
ABSTRACT Convection in massive main-sequence stars generates large-scale magnetic fields in their cores that persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components that are expected to stabilize each other. Previous analytic and numerical calculations did not find any evidence for instability of the Prendergast field over short time-scales. In this paper, we present numerical simulations which show a long time-scale, linear instability of this magnetic field. We find the instability to be robust to changes in boundary conditions and it is not stabilized by strong stable stratification. The instability is a resistive instability, and the growth rate has a power-law dependence on the resistivity, in which the growth rate decreases as the resistivity decreases. We estimate the growth rate of the instability in stars by extrapolating this power law to stellar values of the resistivity. The instability is sufficiently rapid to destabilize the magnetic field on time-scales shorter than the stellar evolution time-scale, indicating that the Prendergast field is not a good model to use in studies of magnetic fields in stars.
-
more » « less
Abstract The proton–alpha drift instability is a possible mechanism of the alpha-particle deceleration and the resulting proton heating in the solar wind. We present hybrid numerical simulations of this instability with particle-in-cell ions and a quasi-neutralizing electron fluid for typical conditions at 1 au. For the parameters used in this paper, we find that fast magnetosonic unstable modes propagate only in the direction opposite to the alpha-particle drift and do not produce the perpendicular proton heating necessary to accelerate the solar wind. Alfvén modes propagate in both directions and heat the protons perpendicularly to the mean magnetic field. Despite being driven by the alpha temperature anisotropy, the Alfvén instability also extracts the energy from the bulk motion of the alpha particles. In the solar wind, the instabilities operate in a turbulent ambient medium. We show that the turbulence suppresses the Alfvén instability but the perpendicular proton heating persists. Unlike a static nonuniform background, the turbulence does not invert the sense of the proton heating associated with the fast magnetosonic instability and it remains preferentially parallel.
