Growth of the MRI in accretion discs – the influence of radiation transport

In this paper, we investigate the influence of radiative transport on the growth of the magnetorotational instability (MRI) in accretion discs. The analysis is performed by the use of analytical and numerical means. We provide a general dispersion relation together with the corresponding eigenfunctions describing the growth rates of small disturbances on a homogeneous background shear flow. The dispersion relation includes compressibility and radiative effects in the flux-limited diffusion approximation. By introducing an effective speed of sound, all the effects of radiation transport can be subsumed into one single parameter. It can be shown that the growth rates of the vertical modes – which are the fastest growing ones – are reduced by radiative transport. For the case of non-vertical modes, the growth rates may instead be enhanced. We quantify the effects of compressibility and radiative diffusion on the growth rates for the gas-pressure dominated case. The analytical discussion is supplemented by numerical simulations, which are also used for a first investigation of the non-linear stage of the MRI in gas-pressure dominated accretion discs with radiation transport included.     

The magnetosphere of oscillating neutron stars in general relativity

Just as a rotating magnetized neutron star has material pulled away from its surface to populate a magnetosphere, a similar process can occur as a result of neutron-star pulsations rather than rotation. This is of interest in connection with the overall study of neutron star oscillation modes but with a particular focus on the situation for magnetars. Following a previous Newtonian analysis of the production of a force-free magnetosphere in this way Timokhin et al., we present here a corresponding general-relativistic analysis. We give a derivation of the general relativistic Maxwell equations for small-amplitude arbitrary oscillations of a non-rotating neutron star with a generic magnetic field and show that these can be solved analytically under the assumption of low current density in the magnetosphere. We apply our formalism to toroidal oscillations of a neutron star with a dipole magnetic field and find that the low current density approximation is valid for at least half of the oscillation modes, similarly to the Newtonian case. Using an improved formula for the determination of the last closed field line, we calculate the energy losses resulting from toroidal stellar oscillations for all of the modes for which the size of the polar cap is small. We find that general relativistic effects lead to shrinking of the size of the polar cap and an increase in the energy density of the outflowing plasma. These effects act in opposite directions but the net result is that the energy loss from the neutron star is significantly smaller than suggested by the Newtonian treatment.     

Alfvén polar oscillations of relativistic stars

We study polar Alfvén oscillations of relativistic stars endowed with a strong global poloidal dipole magnetic field. Here, we focus only on the axisymmetric oscillations which are studied by numerically evolving the two-dimensional perturbation equations. Our study shows that the spectrum of the polar Alfvén oscillations is discrete in contrast to the spectrum of axial Alfvén oscillations which is continuous. We also show that the typical fluid modes, such as the f and p modes, are not significantly affected by the presence of the strong magnetic field.     

Kolmogorov dissipation scales in weakly ionized plasmas

In a weakly ionized plasma, the evolution of the magnetic field is described by a 'generalized Ohm's law' that includes the Hall effect and the ambipolar diffusion terms. These terms introduce additional spatial and time-scales which play a decisive role in the cascading and the dissipation mechanisms in magnetohydrodynamic turbulence. We determine the Kolmogorov dissipation scales for the viscous, the resistive and the ambipolar dissipation mechanisms. The plasma, depending on its properties and the energy injection rate, may preferentially select one of these dissipation scales, thus determining the shortest spatial scale of the supposedly self-similar spectral distribution of the magnetic field. The results are illustrated taking the partially ionized part of the solar atmosphere as an example. Thus, the shortest spatial scale of the supposedly self-similar spectral distribution of the solar magnetic field is determined by any of the four dissipation scales given by the viscosity, the Spitzer resistivity (electron–ion collisions), the resistivity due to electron–neutral collisions and the ambipolar diffusivity. It is found that the ambipolar diffusion dominates for reasonably large energy injection rate. The robustness of the magnetic helicity in the partially ionized solar atmosphere would facilitate the formation of self-organized vortical structures.     

Axisymmetric magnetic fields in stars: relative strengths of poloidal and toroidal components

In this third paper in a series on stable magnetic equilibria in stars, I look at the stability of axisymmetric field configurations and, in particular, the relative strengths of the toroidal and poloidal components. Both toroidal and poloidal fields are unstable on their own, and stability is achieved by adding the two together in some ratio. I use Tayler's stability conditions for toroidal fields and other analytic tools to predict the range of stable ratios and then check these predictions by running numerical simulations. If the energy in the poloidal component as a fraction of the total magnetic energy is written as E_p / E , it is found that the stability condition is a ( E / U ) < E_p / E ≲ 0.8 where E /U is the ratio of magnetic to gravitational energy in the star and a is some dimensionless factor whose value is of order 10 in a main-sequence star and of order 10³ in a neutron star. In other words, whilst the poloidal component cannot be significantly stronger than the toroidal, the toroidal field can be very much stronger than the poloidal–given that in realistic stars we expect E / U < 10⁻⁶. The implications of this result are discussed in various contexts such as the emission of gravitational waves by neutron stars, free precession and a 'hidden' energy source for magnetars.     

The Rotating Sunspot in AR 10930

We study the pattern and behavior of a rotating sunspot in Active Region 10930. The rotational angular speed has been extracted from the apparent motions of the sunspot determined by applying a new optical technique – called non-linear affine velocity estimator (NAVE) – to high-resolution G-band images taken by the Solar Optical Telescope (SOT) onboard the Hinode satellite. The structure and dynamics of coronal loops in this active region have been examined using the images obtained by the X-Ray Telescope (XRT) and the spectral data taken by the Extreme-ultraviolet Imaging Spectrometer (EIS), both also onboard Hinode. Our results are summarized as follows: i) The small sunspot of positive polarity rotated counterclockwise about its center by 540° during the period of five days. ii) Its angular velocity varied with the azimuth angle as well as the radial distance, being affected by the asymmetric shape of the umbra. iii) The angular velocity increased up to 8° h⁻¹ until 13 December as the sunspot grew, and then decreased rapidly down to 3° h⁻¹ on the next day as the sunspot decayed. iv) The coronal loops that connected the two sunspots became sigmoidal in shape. v) The coronal emissions from the regions around the rotating sunspot were blueshifted, which may indicate the expansion of the coronal loops. Our results suggest that the rotation of the sunspot may be closely related to the dynamic development of emerging twisted magnetic fields.     

Magnetized relativistic jets and long-duration GRBs from magnetar spin-down during core-collapse supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spin-down of a newly formed millisecond,   B ∼ 10¹⁵ G  , magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spin-down powers  (∼10⁵¹–10⁵² erg s⁻¹)  , the magnetar wind is superfast at almost all latitudes, while for lower spin-down powers  (∼10⁵⁰ erg s⁻¹)  , the wind is subfast but still super-Alfvénic. In all cases, the rates at which the neutron star loses mass, angular momentum and energy are very similar to the corresponding free wind values (≲30 per cent differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated  (∼5–10°)  relativistic jet out along the rotation axis of the star. Nearly all of the spin-down power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.