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.     

Oscillations of rapidly rotating superfluid stars

Using time evolutions of the relevant linearized equations, we study non-axisymmetric oscillations of rapidly rotating and superfluid neutron stars. We consider perturbations of Newtonian axisymmetric background configurations and account for the presence of superfluid components via the standard two-fluid model. Within the Cowling approximation, we are able to carry out evolutions for uniformly rotating stars up to the mass-shedding limit. This leads to the first detailed analysis of superfluid neutron star oscillations in the fast rotation regime, where the star is significantly deformed by the centrifugal force. For simplicity, we focus on background models where the two fluids (superfluid neutrons and protons) corotate, are in β-equilibrium and co-exist throughout the volume of the star. We construct sequences of rotating stars for two analytical model equations of state. These models represent relatively simple generalizations of single fluid, polytropic stars. We study the effects of entrainment, rotation and symmetry energy on non-radial oscillations of these models. Our results show that entrainment and symmetry energy can have a significant effect on the rotational splitting of non-axisymmetric modes. In particular, the symmetry energy modifies the inertial mode frequencies considerably in the regime of fast rotation.     

Magnetic field evolution in neutron stars

Neutron stars contain persistent, ordered magnetic fields that are the strongest known in the Universe. However, their magnetic fluxes are similar to those in magnetic A and B stars and white dwarfs, suggesting that flux conservation during gravitational collapse may play an important role in establishing the field, although it might also be modified substantially by early convection, differential rotation, and magnetic instabilities. The equilibrium field configuration, established within hours (at most) of the formation of the star, is likely to be roughly axisymmetric, involving both poloidal and toroidal components. The stable stratification of the neutron star matter (due to its radial composition gradient) probably plays a crucial role in holding this magnetic structure inside the star. The field can evolve on long time scales by processes that overcome the stable stratification, such as weak interactions changing the relative abundances and ambipolar diffusion of charged particles with respect to neutrons. These processes become more effective for stronger magnetic fields, thus naturally explaining the magnetic energy dissipation expected in magnetars, at the same time as the longer-lived, weaker fields in classical and millisecond pulsars. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-linear hydrodynamical evolution of rotating relativistic stars: numerical methods and code tests

We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, non-linear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes (based on approximate Riemann solvers) and compare their accuracy in preserving uniformly rotating stationary initial configurations in long-term evolutions. Among these four schemes, we find that the third-order piecewise parabolic method scheme is superior in maintaining the initial rotation law in long-term evolutions, especially near the surface of the star. It is further shown that HRSC schemes are suitable for the evolution of perturbed neutron stars and for the accurate identification (via Fourier transforms) of normal modes of oscillation. This is demonstrated for radial and quadrupolar pulsations in the non-rotating limit, where we find good agreement with frequencies obtained with a linear perturbation code. The code can be used for studying small-amplitude or non-linear pulsations of differentially rotating neutron stars, while our present results serve as testbed computations for three-dimensional general-relativistic evolution codes.     

Burial of the polar magnetic field of an accreting neutron star – II. Hydromagnetic stability of axisymmetric equilibria

The theory of polar magnetic burial in accreting neutron stars predicts that a mountain of accreted material accumulates at the magnetic poles of the star, and that, as the mountain spreads equatorward, it is confined by, and compresses, the equatorial magnetic field. Here, we extend previous, axisymmetric, Grad–Shafranov calculations of the hydromagnetic structure of a magnetic mountain up to accreted masses as high as   M _a= 6 × 10⁻⁴ M_⊙  , by importing the output from previous calculations (which were limited by numerical problems and the formation of closed bubbles to   M _a < 10⁻⁴ M_⊙  ) into the time-dependent, ideal-magnetohydrodynamic code zeus-3d and loading additional mass on to the star dynamically. The rise of buoyant magnetic bubbles through the accreted layer is observed in these experiments. We also investigate the stability of the resulting hydromagnetic equilibria by perturbing them in zeus-3d . Surprisingly, it is observed that the equilibria are marginally stable for all   M _a≤ 6 × 10⁻⁴ M_⊙  ; the mountain oscillates persistently when perturbed, in a combination of Alfvén and acoustic modes, without appreciable damping or growth, and is therefore not disrupted (apart from a transient Parker instability initially, which expels <1 per cent of the mass and magnetic flux).     

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)     

Axisymmetric modes of rotating relativistic stars in the Cowling approximation

Axisymmetric pulsations of rotating neutron stars can be excited in several scenarios, such as core collapse, crust- and core-quakes or binary mergers, and could become detectable in either gravitational waves or high-energy radiation. Here, we present a comprehensive study of all low-order axisymmetric modes of uniformly and rapidly rotating relativistic stars. Initial stationary configurations are appropriately perturbed and are numerically evolved using an axisymmetric, non-linear relativistic hydrodynamics code, assuming time-independence of the gravitational field (Cowling approximation). The simulations are performed using a high-resolution shock-capturing finite-difference scheme accurate enough to maintain the initial rotation law for a large number of rotational periods, even for stars at the mass-shedding limit. Through Fourier transforms of the time evolution of selected fluid variables, we compute the frequencies of quasi-radial and non-radial modes with spherical harmonic indices l =0 , 1, 2 and 3, for a sequence of rotating stars from the non-rotating limit to the mass-shedding limit. The frequencies of the axisymmetric modes are affected significantly by rotation only when the rotation rate exceeds about 50 per cent of the maximum allowed. As expected, at large rotation rates, apparent mode crossings between different modes appear. In addition to the above modes, several axisymmetric inertial modes are also excited in our numerical evolutions.     

Extrasolar planets and the rotation and axisymmetric mass-loss of evolved stars

I examine the implications of the recently found extrasolar planets on the planet-induced axisymmetric mass-loss model for the formation of elliptical planetary nebulae (PNe). This model attributes the low departure from spherical mass-loss of upper asymptotic giant branch (AGB) stars to envelope rotation which results from deposition of orbital angular momentum of the planets. Since about half of all PNe are elliptical, i.e., have low equatorial to polar density contrast, it was predicted that about 50 per cent of all Sun-like stars have Jupiter-like planets around them, i.e., a mass about equal to that of Jupiter, M _J, or more massive. In the light of the new findings that only 5 per cent of Sun-like stars have such planets, and a newly proposed mechanism for axisymmetric mass-loss, the cool magnetic spots model, I revise this prediction. I predict that indeed ∼50 per cent of PN progenitors do have close planets around them, but the planets can have much lower masses, as low as ∼0.01 M _J, in order to spin-up the envelopes of AGB stars efficiently. To support this claim, I follow the angular momentum evolution of single stars with main-sequence mass in the range of 1.3–2.4 M_⊙ , as they evolve to the post-AGB phase. I find that single stars rotate much too slowly to possess any significant non-spherical mass-loss as they reach the upper AGB. It seems, therefore, that planets, in some cases even Earth-like planets, are sufficient to spin-up the envelope of these AGB stars for them to form elliptical PNe. The prediction that on average several such planets orbit each star, as in the Solar system, still holds.     

Models of rapidly rotating neutron stars: remnants of accretion-induced collapse

Equilibrium models of differentially rotating nascent neutron stars are constructed, which represent the result of the accretion-induced collapse of rapidly rotating white dwarfs. The models are built in a two-step procedure: (1) a rapidly rotating pre-collapse white dwarf model is constructed; (2) a stationary axisymmetric neutron star having the same total mass and angular momentum distribution as the white dwarf is constructed. The resulting collapsed objects consist of a high-density central core of size roughly 20 km, surrounded by a massive accretion torus extending over 1000 km from the rotation axis. The ratio of the rotational kinetic energy to the gravitational potential energy of these neutron stars ranges from 0.13 to 0.26, suggesting that some of these objects may have a non-axisymmetric dynamical instability that could emit a significant amount of gravitational radiation.     

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.     

A new numerical scheme for structures of rotating magnetic stars

We have developed a new numerical scheme for obtaining structures of rapidly rotating stars with strong magnetic fields. In our scheme, both poloidal and toroidal magnetic fields can be treated for stars with compressibility and infinite conductivity. By introducing the vector potential and its integral representation, we can treat the boundary condition for the magnetic fields across the surface properly. We show structures and distributions of magnetic fields as well as the distributions of the currents of rotating magnetic polytropic stars with polytropic index   N = 1.5  . The shapes of magnetic stars are oblate as long as the magnetic vector potential decreases as 1/ r when   r →∞  . For extremely strong magnetic fields, equilibrium configurations can be of toroidal shapes.     

On the excitation mechanism in roAp stars

We investigate a model for the excitation of high-order oscillations in roAp stars. In this model we assume that the strong concentration of magnetic field about the magnetic poles is enough to suppress convection. Thus the model considered is composed of two polar regions, in which convection is presumed to be suppressed totally, and an equatorial region, where the convection is unaffected. This model is generated by building pairs of locally spherically symmetrical equilibria to represent the polar and equatorial regions of the star, which are patched together below the base of the convection zone. Gravitational settling of heavy elements is taken into account by choosing appropriate chemical composition profiles for both the polar and equatorial regions. Our results indicate that the composite model is unstable against axisymmetric non-radial high-order modes of pulsation that are aligned with the magnetic poles. The oscillations are excited by the κ mechanism acting principally in the hydrogen ionization zones of the polar regions. The effect of the lateral inhomogeneity on the second frequency differences is also investigated; we find that the perturbation to them by the inhomogeneity is of the same order as the second differences themselves, thereby hindering potential attempts to use such differences to identify the degrees of the modes in a straightforward way.     

Towards multiple-star population synthesis

The multiplicities of stars, and some other properties, were collected recently by Eggleton & Tokovinin, for the set of 4559 stars with Hipparcos magnitude brighter than 6.0 (4558 excluding the Sun). In this paper I give a numerical recipe for constructing, by a Monte Carlo technique, a theoretical ensemble of multiple stars that resembles the observed sample. Only multiplicities up to eight are allowed; the observed set contains only multiplicities up to seven. In addition, recipes are suggested for dealing with the selection effects and observational uncertainties that attend the determination of multiplicity. These recipes imply, for example, that to achieve the observed average multiplicity of 1.53, it would be necessary to suppose that the real population has an average multiplicity slightly over 2.0.
This numerical model may be useful for (i) comparison with the results of star and star cluster formation theory, (ii) population synthesis that does not ignore multiplicity above 2 and (iii) initial conditions for dynamical cluster simulations.     

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.     

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.     

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.     

Relativistic models of magnetars: the twisted torus magnetic field configuration

We find general relativistic solutions of equilibrium magnetic field configurations in magnetars, extending previous results of Colaiuda et al. Our method is based on the solution of the relativistic Grad–Shafranov equation, to which Maxwell's equations can be reduced. We obtain equilibrium solutions with the toroidal magnetic field component confined into a finite region inside the star, and the poloidal component extending to the exterior. These so-called twisted torus configurations have been found to be the final outcome of dynamical simulations in the framework of Newtonian gravity, and appear to be more stable than other configurations. The solutions include higher-order multipoles, which are coupled to the dominant dipolar field. We use arguments of minimal energy to constrain the ratio of the toroidal to the poloidal field.     

High-velocity stars from decay of small stellar systems

In this paper we present numerical results on the decay of small stellar systems under different initial conditions (multiplicity 3 ≤  N  ≤ 10, and various mass spectra, initial velocities and initial configurations). The numerical treatment uses the CHAIN1 code (Mikkola &38; Aarseth). Particular attention is paid to the distribution of high-velocity escapers: we define these as stars with velocity above 30 km s⁻¹. These numerical experiments confirm that small N -body systems are dynamically unstable and produce cascades of escapers in the process of their decay. It is shown that the fraction of stars that escape from small dense stellar systems with an escape velocity greater than 30 km s⁻¹ is ∼1 per cent for all systems treated here. This relatively small fraction must be considered in relation to the rate of star formation in the Galaxy in small groups: this could explain some moderately high-velocity stars observed in the Galactic disc and possibly some young stars with relatively high metallicity in the thick disc.     

The circumstellar envelopes of Be stars: viscous disc dynamics

We have constructed models of axisymmetric, circumstellar envelopes for Be star discs by successfully combining two numerical codes: a non-local thermodynamic equilibrium (non-LTE) radiative transfer code which calculates the level populations and disc temperature distribution self-consistently, and a hydrodynamical code. The output of one code is used as input to the other code, and hence evolving density and thermal structures may be examined. The temperatures, disc density and velocity distribution are used to investigate the outflowing viscous disc model for Be stars. We find that these simulations place constraints on the power-law density decrease in the disc with increasing distance from the star. We find that the power-law index for the line-forming region of the disc lies between 3 and 3.5 with a small dispersion.     

The first magnetic maps of a pre-main-sequence binary star system – HD 155555

We present the first maps of the surface magnetic fields of a pre-main-sequence binary system. Spectropolarimetric observations of the young, 18 Myr, HD 155555 (V824 Ara, G5IV+K0IV) system were obtained at the Anglo-Australian Telescope in 2004 and 2007. Both data sets are analysed using a new binary Zeeman–Doppler imaging (ZDI) code. This allows us to simultaneously model the contribution of each component to the observed circularly polarized spectra. Stellar brightness maps are also produced for HD 155555 and compared to previous Doppler images.
Our radial magnetic maps reveal a complex surface magnetic topology with mixed polarities at all latitudes. We find rings of azimuthal field on both stars, most of which are found to be non-axisymmetric with the stellar rotational axis. We also examine the field strength and the relative fraction of magnetic energy stored in the radial and azimuthal field components at both epochs. A marked weakening of the field strength of the secondary star is observed between the 2004 and 2007 epochs. This is accompanied by an apparent shift in the location of magnetic energy from the azimuthal to radial field. We suggest that this could be indicative of a magnetic activity cycle. We use the radial magnetic maps to extrapolate the coronal field (by assuming a potential field) for each star individually – at present ignoring any possible interaction. The secondary star is found to exhibit an extreme tilt (≈75°) of its large-scale magnetic field to that of its rotation axis for both epochs. The field complexity that is apparent in the surface maps persists out to a significant fraction of the binary separation. Any interaction between the fields of the two stars is therefore likely to be complex also. Modelling this would require a full binary field extrapolation.