Magnetic fields in axisymmetric neutron stars

We derive general equations for axisymmetric Newtonian magnetohydrodynamics and use these as the basis of a code for calculating equilibrium configurations of rotating magnetized neutron stars in a stationary state. We investigate the field configurations that result from our formalism, which include purely poloidal, purely toroidal and mixed fields. For the mixed-field formalism, the toroidal component appears to be bounded at less than 7 per cent. We calculate distortions induced both by magnetic fields and by rotation. From our non-linear work, we are able to look at the realm of validity of perturbative work: we find for our results that perturbative-regime formulae for magnetic distortions agree to within 10 per cent of the non-linear results if the ellipticity is less than 0.15 or the average field strength is less than 10¹⁷ G. We also consider how magnetized equilibrium structures vary for different polytropic indices.     

Equilibrium configurations of strongly magnetized neutron stars with realistic equations of state

We investigate equilibrium sequences of magnetized rotating stars with four kinds of realistic equations of state (EOSs) of SLy, FPS, Shen and LS, employing the Tomimura–Eriguchi scheme to construct the equilibrium configurations. We study the basic physical properties of the sequences in the framework of Newtonian gravity. In addition, we take a new step by taking into account a general relativistic effect to the magnetized rotating configurations. With these computations, we find that the properties of the Newtonian magnetized stars, e.g. structure of magnetic field, highly depends on the EOSs. The toroidal magnetic fields concentrate rather near the surface for Shen and LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. Paying attention to the stiffness of the EOSs, we analyse this tendency in detail. In the general relativistic stars, we find that the difference due to the EOSs becomes small because all the employed EOSs become sufficiently stiff for the large maximum density, typically greater than  10¹⁵ g cm⁻³  . The maximum baryon mass of the magnetized stars with axis ratio   q ∼ 0.7  increases about up to 20 per cent for that of spherical stars. We furthermore compute equilibrium sequences at finite temperature, which should serve as an initial condition for the hydrodynamic study of newly born magnetars. Our results suggest that we may obtain information about the EOSs from the observation of the masses of magnetars.     

On the solution space of differentially rotating neutron stars in general relativity

A highly accurate, multidomain spectral code is used in order to construct sequences of general relativistic, differentially rotating neutron stars in axisymmetry and stationarity. For bodies with a spheroidal topology and a homogeneous or an   N = 1  polytropic equation of state, we investigate the solution space corresponding to broad ranges of degree of differential rotation and stellar densities. In particular, starting from static and spherical configurations, we analyse the changes of the corresponding surface shapes as the rate of rotation is increased. For a sufficiently weak degree of differential rotation, the sequences terminate at a mass-shedding limit, while for moderate and strong rates of differential rotation they exhibit a continuous parametric transition to a regime of toroidal fluid bodies. In this article, we concentrate on the appearance of this transition, analyse in detail its occurrence and show its relevance for the calculation of astrophysical sequences. Moreover, we find that the solution space contains various types of spheroidal configurations, which were not considered in previous work, mainly due to numerical limitations.     

Surface magnetic fields on two accreting T Tauri stars: CV Cha and CR Cha

We have produced brightness and magnetic field maps of the surfaces of CV Cha and CR Cha: two actively accreting G- and K-type T Tauri stars in the Chamaeleon I star-forming cloud with ages of 3–5 Myr. Our magnetic field maps show evidence for strong, complex multipolar fields similar to those obtained for young rapidly rotating main-sequence stars. Brightness maps indicate the presence of dark polar caps and low-latitude spots – these brightness maps are very similar to those obtained for other pre-main-sequence and rapidly rotating main-sequence stars.
Only two other classical T Tauri stars have been studied using similar techniques so far: V2129 Oph and BP Tau. CV Cha and CR Cha show magnetic field patterns that are significantly more complex than those recovered for BP Tau, a fully convective T Tauri star.
We discuss possible reasons for this difference and suggest that the complexity of the stellar magnetic field is related to the convection zone; with more complex fields being found in T Tauri stars with radiative cores (V2129 Oph, CV Cha and CR Cha). However, it is clearly necessary to conduct magnetic field studies of T Tauri star systems, exploring a wide range of stellar parameters in order to establish how they affect magnetic field generation, and thus how these magnetic fields are likely to affect the evolution of T Tauri star systems as they approach the main sequence.     

Magnetic fields of Ap stars as a result of the Tayler instability

Ap star magnetism is often attributed to fossil magnetic fields which have not changed much since the pre‐main‐sequence epoch of the stars. Stable magnetic field configurations are known which could persist probably for the entire mainsequence life of the star, but they may not show the complexity and diversity exhibited by the Ap stars observed. We suggest that the Ap star magnetism is not a result of stable configurations, but is the result of an instability based on strong toroidal magnetic fields buried in the stars. The highly nonaxisymmetric remainders of the instability are reminiscent of the diversity of fields seen on Ap stars. The strengths of these remnant magnetic fields are actually between a few per cent up to considerable fractions of the internal toroidal field; this means field strengths of the order of kGauss being compatible with what is observed. The magnetic fields emerge at the surface rather quickly; rough estimates deliver time‐scales of the order of a few years. Since rotation stabilizes the instability, normal A stars may still host considerable, invisible toroidal magnetic fields (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Can late-type active stars be explained by a dipole magnetic trap?

In this paper we review four different types of X-ray and/or radio observations of active late-type stars. We then consider if a single magnetic source configuration – a toroidal dipole magnetic trap – can possibly explain these various different observations. We conclude that, indeed, dipole magnetic confinement (similar to the magnetic configurations of the Earth's radiation belts and the case of Jupiter and the Io torus) can explain all the diverse observational data. We take this to be very strong observational support for this type of magnetic confinement scheme. We also consider that this magnetic configuration is only likely to be established and maintained in the most active stars.     

Stability of toroidal magnetic fields in the solar tachocline and beneath

Stability of toroidal magnetic field in a stellar radiation zone is considered for the cases of uniform and differential rotation. In the rigidly rotating radiative core shortly below the tachocline, the critical magnetic field for instability is about 600 G. The unstable disturbances for slightly supercritical fields have short radial scales ∼1 Mm. Radial mixing produced by the instability is estimated to conclude that the internal field of the sun can exceed the critical value of 600 G only marginally. Otherwise, the mixing is too strong and not compatible with the observed lithium abundance. Analysis of joint instability of differential rotation and toroidal field leads to the conclusion that axisymmetric models of the laminar solar tachocline are stable to nonaxisymmetric disturbances. The question of whether sun-like stars can posses tachoclines is addressed with positive answer for stars with rotation periods shorter than about two months. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic fields and rotation in white dwarfs and neutron stars

Recent spectropolarimetric observations of Ap and Bp stars with improved sensitivity have suggested that most Ap and Bp stars are magnetic with dipolar fields of at least a few hundred gauss. These new estimates suggest that the range of magnetic fluxes found for the majority of magnetic white dwarfs is similar to that of main-sequence Ap–Bp stars, thus strengthening the empirical evidence for an evolutionary link between magnetism on the main sequence and magnetism in white dwarfs. We draw parallels between the magnetic white dwarfs and the magnetic neutron stars and argue that the observed range of magnetic fields in isolated neutron stars  ( B_p ∼ 10¹¹–10¹⁵ G)  could also be explained if their mainly O-type progenitors have effective dipolar fields in the range of a few gauss to a few kilogauss, assuming approximate magnetic flux conservation with the upper limit being consistent with the recent measurement of a field of   B_p ∼ 1100 G  for θ Orion C.
In the magnetic field–rotation diagram, the magnetic white dwarfs can be divided into three groups of different origin: a significant group of strongly magnetized slow rotators  ( P _rot∼ 50 –100 yr)  that have originated from single-star evolution, a group of strongly magnetized fast rotators  ( P _rot∼ 700 s)  , typified by EUVE J0317–853, that have originated from a merger, and a group of modest rotators ( P _rot∼ hours–days) of mixed origin (single-star and CV-type binary evolution). We propose that the neutron stars may similarly divide into distinct classes at birth , and suggest that the magnetars may be the counterparts of the slowly rotating high-field magnetic white dwarfs.     

Dynamo-generated magnetic fields at the surface of a massive star

Spruit has shown that an astrophysical dynamo can operate in the non-convective material of a differentially rotating star as a result of a particular instability in the magnetic field (the Tayler instability). By assuming that the dynamo operates in a state of marginal instability, Spruit has obtained formulae which predict the equilibrium strengths of azimuthal and radial field components in terms of local physical quantities. Here, we apply Spruit's formulae to our previously published models of rotating massive stars in order to estimate Tayler dynamo field strengths. There are no free parameters in Spruit's formulae. In our models of 10- and  50-M_⊙  stars on the zero-age main sequence, we find internal azimuthal fields of up to 1 MG, and internal radial components of a few kG. Evolved models contain weaker fields. In order to obtain estimates of the field strength at the stellar surface, we examine the conditions under which the Tayler dynamo fields are subject to magnetic buoyancy. We find that conditions for Tayler instability overlap with those for buoyancy at intermediate to high magnetic latitudes. This suggests that fields emerge at the surface of a massive star between magnetic latitudes of about 45° and the poles. We attempt to estimate the strength of the field which emerges at the surface of a massive star. Although these estimates are very rough, we find that the surface field strengths overlap with values which have been reported recently for line-of-sight fields in several O and B stars.     

High-precision magnetic field measurements of Ap and Bp stars

In this paper we describe a new approach for measuring the mean longitudinal magnetic field and net linear polarization of Ap and Bp stars. As was demonstrated by Wade et al., least-squares deconvolution (LSD; Donati et al.) provides a powerful technique for detecting weak Stokes V , Q and U Zeeman signatures in stellar spectral lines. These signatures have the potential to apply strong new constraints to models of stellar magnetic field structure. Here we point out two important uses of LSD Stokes profiles. First, they can provide very precise determinations of the mean longitudinal magnetic field. In particular, this method allows one frequently to obtain 1 σ error bars better than 50 G, and smaller than 20 G in some cases. This method is applicable to both broad- and sharp-lined stars, with both weak and strong magnetic fields, and effectively redefines the quality standard of longitudinal field determinations. Secondly, LSD profiles can in some cases provide a measure of the net linear polarization, a quantity analogous to the broad-band linear polarization recently used to derive detailed magnetic field models for a few stars (e.g. Leroy et al.). In this paper we report new high-precision measurements of the longitudinal fields of 14 magnetic Ap/Bp stars, as well as net linear polarization measurements for four of these stars, derived from LSD profiles.     

Gravitational collapse of polytropic, magnetized, filamentary clouds

For the case in which the gas of a magnetized filamentary cloud obeys a polytropic equation of state, gravitational collapse of the cloud is studied using a simplified model. We concentrate on the radial distribution and restrict ourselves to a purely toroidal magnetic field. If the axial motions and poloidal magnetic fields are sufficiently weak, we could reasonably expect our solutions to be a good approximation. We show that while the filament experiences gravitational condensation and the density at the centre increases, the toroidal flux-to-mass ratio remains constant. A series of spatial profiles of density, velocity and magnetic field for several values of the toroidal flux-to-mass ratio and the polytropic index, is obtained numerically and discussed.     

Relativistic models of magnetars: structure and deformations

We find numerical solutions of the coupled system of Einstein–Maxwell equations with a linear approach, in which the magnetic field acts as a perturbation of a spherical neutron star. In our study, magnetic fields having both poloidal and toroidal components are considered, and higher order multipoles are also included. We evaluate the deformations induced by different field configurations, paying special attention to those for which the star has a prolate shape. We also explore the dependence of the stellar deformation on the particular choice of the equation of state and on the mass of the star. Our results show that, for neutron stars with mass   M = 1.4 M_⊙  and surface magnetic fields of the order of 10¹⁵ G, a quadrupole ellipticity of the order of 10⁻⁶ to 10⁻⁵ should be expected. Low-mass neutron stars are in principle subject to larger deformations (quadrupole ellipticities up to 10⁻³ in the most extreme case). The effect of quadrupolar magnetic fields is comparable to that of dipolar components. A magnetic field permeating the whole star is normally needed to obtain negative quadrupole ellipticities, while fields confined to the crust typically produce positive quadrupole ellipticities.     

Dynamo action in the presence of an imposed magnetic field

Dynamo action within the cores of Ap stars may offer intriguing possibilities for understanding the persistent magnetic fields observed on the surfaces of these stars. Deep within the cores of Ap stars, the coupling of convection with rotation likely yields magnetic dynamo action, generating strong magnetic fields. However, the surface fields of the magnetic Ap stars are generally thought to be of primordial origin. Recent numerical models suggest that a primordial field in the radiative envelope may possess a highly twisted toroidal shape. We have used detailed 3-D simulations to study the interaction of such a twisted magnetic field in the radiative envelope with the core-dynamo operating in the interior of a 2 solar mass A-type star. The resulting dynamo action is much more vigorous than in the absence of such a fossil field, yielding magnetic field strengths (of order 100 kG) much higher than their equipartition values relative to the convective velocities. We examine the generation of these fields, as well as the growth of large-scale magnetic structure that results from imposing a fossil magnetic field. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the r-mode spectrum of relativistic stars: the inclusion of the radiation reaction

We consider both mode calculations and time-evolutions of axial r modes for relativistic uniformly rotating non-barotropic neutron stars, using the slow-rotation formalism, in which rotational corrections are considered up to linear order in the angular velocity Ω. We study various stellar models, such as uniform density models, polytropic models with different polytropic indices n , and some models based on realistic equations of state. For weakly relativistic uniform density models and polytropes with small values of n , we can recover the growth times predicted from Newtonian theory when standard multipole formulae for the gravitational radiation are used. However, for more compact models, we find that relativistic linear perturbation theory predicts a weakening of the instability compared to the Newtonian results. When turning to polytropic equations of state, we find that for certain ranges of the polytropic index n , the r mode disappears, and instead of a growth, the time-evolutions show a rapid decay of the amplitude. This is clearly at variance with the Newtonian predictions. It is, however, fully consistent with our previous results obtained in the low-frequency approximation.     

A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment.
The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to  ∼3000 km s⁻¹  ) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.     

Stability of magnetic flux tubes in close binary stars

Photometric and Doppler imaging observations of active binaries indicate the existence of starspots at preferred longitudes (position angles with respect to the companion star). We investigate the stability of magnetic flux tubes in the convection zone of close, fast‐rotating binary stars and explore whether the observed preferred longitudes could be caused by tidal forces and the deformation of the active star. We assume a synchronized binary system with spin axes perpendicular to the orbital plane and a rotation period of a few days. The tidal force and the deviation from spherical structure are considered in lowest‐order perturbation theory. The magnetic field is in the form of toroidal magnetic flux rings, which are stored in mechanical equilibrium within the stably stratified overshoot region beneath the convection zone until the field has grown sufficiently strong for the undulatory instability to initiate the formation of rising loops. Frequencies and geometry of stable as well as growth rates of unstable eigenmodes are determined by linear stability analysis. Particular consideration is given to the question whether the effects of tidal forces and perturbations of the stellar structure can force a rising flux loop to enter the convection zone at specific longitudes.     

On the quasi-periodic oscillations in magnetars

We study torsional Alfvén oscillations of magnetars, that is neutron stars with a strong magnetic field. We consider the poloidal and toroidal components of the magnetic field and a wide range of equilibrium stellar models. We use a new coordinate system  ( X , Y )  , where     and     and a ₁ is the radial component of the magnetic field. In this coordinate system, the one+two-dimensional evolution equation describing the quasi-periodic oscillations (QPOs), see Sotani et al., is reduced to a one+one-dimensional equation where the perturbations propagate only along the y -axis. We solve the one+one-dimensional equation for different boundary conditions and the open magnetic field lines, that is magnetic field lines that reach the surface and there match up with the exterior dipole magnetic field as well as closed magnetic lines, i.e. magnetic lines that never reach the stellar surface. For the open field lines, we find two families of QPO frequencies: a family of 'lower' QPO frequencies which is located near the x -axis and a family of 'upper' frequencies located near the y -axis. According to Levin, the fundamental frequencies of these two families can be interpreted as the turning point of the continuous spectrum. We find that the upper frequencies are multiples of the lower ones by a constant equalling  2 n + 1  . For the closed lines, the corresponding factor is   n + 1  . By using these relations, we can explain both the lower and the higher observed frequencies in SGR 1806−20 and SGR 1900+14.     

Torsional oscillations of magnetized relativistic stars

Strong magnetic fields in relativistic stars can be a cause of crust fracturing, resulting in the excitation of global torsional oscillations. Such oscillations could become observable in gravitational waves or in high-energy radiation, thus becoming a tool for probing the equation of state of relativistic stars. As the eigenfrequency of torsional oscillation modes is affected by the presence of a strong magnetic field, we study torsional modes in magnetized relativistic stars. We derive the linearized perturbation equations that govern torsional oscillations coupled to the oscillations of a magnetic field, when variations in the metric are neglected (Cowling approximation). The oscillations are described by a single two-dimensional wave equation, which can be solved as a boundary-value problem to obtain eigenfrequencies. We find that, in the non-magnetized case, typical oscillation periods of the fundamental     torsional modes can be nearly a factor of 2 larger for relativistic stars than previously computed in the Newtonian limit. For magnetized stars, we show that the influence of the magnetic field is highly dependent on the assumed magnetic field configuration, and simple estimates obtained previously in the literature cannot be used for identifying normal modes observationally.     

Toroidal versus poloidal magnetic fields in Sun-like stars: a rotation threshold

From a set of stellar spectropolarimetric observations, we report the detection of surface magnetic fields in a sample of four solar-type stars, namely HD 73350, HD 76151, HD 146233 (18 Sco) and HD 190771. Assuming that the observed variability of polarimetric signal is controlled by stellar rotation, we establish the rotation periods of our targets, with values ranging from 8.8 d (for HD 190771) to 22.7 d (for HD 146233). Apart from rotation, fundamental parameters of the selected objects are very close to the Sun's, making this sample a practical basis to investigate the specific impact of rotation on magnetic properties of Sun-like stars.
We reconstruct the large-scale magnetic geometry of the targets as a low-order  (ℓ < 10)  spherical harmonic expansion of the surface magnetic field. From the set of magnetic maps, we draw two main conclusions. (i) The magnetic energy of the large-scale field increases with rotation rate. The increase in chromospheric emission with the mean magnetic field is flatter than observed in the Sun. Since the chromospheric flux is also sensitive to magnetic elements smaller than those contributing to the polarimetric signal, this observation suggests that a larger fraction of the surface magnetic energy is stored in large scales as rotation increases. (ii) Whereas the magnetic field is mostly poloidal for low rotation rates, more rapid rotators host a large-scale toroidal component in their surface field. From our observations, we infer that a rotation period lower than ≈12 d is necessary for the toroidal magnetic energy to dominate over the poloidal component.     

Physical characteristics of critically distorted or almost spherically shaped rotating magnetic polytropic models

We use the so-called complex plane iterative technique (CIT) to the computation of polytropic stars distorted by rotation (either rigid or differential) and magnetic fields (both toroidal and poloidal). We give emphasis on computing(i) critically rotating configurations, and (ii) configurations that – dueto the counterbalancing of the effects of rotation and poloidal magnetic field with the effects of toroidal magnetic field – obtain an almost spherical shape. This revised version was published online in July 2006 with corrections to the Cover Date.