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.     

Can a slowly rotating neutron star be a radio pulsar?

It is shown that the radius of curvature of magnetic field lines in the polar region of a rotating magnetized neutron star can be significantly less than the usual radius of curvature of the dipole magnetic field. The magnetic field in the polar cap is distorted by toroidal electric currents flowing in the neutron star crust. These currents close up the magnetospheric currents driven by the electron–positron plasma generation process in the pulsar magnetosphere. Owing to the decrease in the radius of curvature, electron–positron plasma generation becomes possible even for slowly rotating neutron stars, with   PB ^−2/3₁₂ < 10 s  , where P is the period of star rotation and   B ₁₂= B /10¹² G  is the magnitude of the magnetic field on the star surface.     

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)     

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.     

Magnetic and spin evolution of neutron stars in close binaries

The evolution of neutron stars in close binary systems with a low-mass companion is considered, assuming the magnetic field to be confined within the solid crust. We adopt the standard scenario for the evolution in a close binary system, in which the neutron star passes through four evolutionary phases ('isolated pulsar'–'propeller'– accretion from the wind of a companion – accretion resulting from Roche-lobe overflow). Calculations have been performed for a great variety of parameters characterizing the properties of both the neutron star and the low-mass companion. We find that neutron stars with more or less standard magnetic field and spin period that are processed in low-mass binaries can evolve to low-field rapidly rotating pulsars. Even if the main-sequence life of a companion is as long as 10¹⁰ yr, the neutron star can maintain a relatively strong magnetic field to the end of the accretion phase. The model that is considered can account well for the origin of millisecond pulsars.     

Low-frequency gravitational waves from cosmological compact binaries

We consider gravitational waves emitted by various populations of compact binaries at cosmological distances. We use population synthesis models to characterize the properties of double neutron stars, double black holes and double white dwarf binaries, and white dwarf–neutron star, white dwarf–black hole and black hole–neutron star systems.
We use the observationally determined cosmic star formation history to reconstruct the redshift distribution of these sources and their merging rate evolution.
The gravitational signals emitted by each source during its early spiralling in phase add randomly to produce a stochastic background in the low-frequency band with spectral strain amplitude between ~10⁻¹⁸ and ~5×10⁻¹⁷ Hz^−1/2 at frequencies in the interval ~5×10⁻⁶–5×10⁻⁵ Hz.
The overall signal, which at frequencies above 10⁻⁴ Hz is largely dominated by double white dwarf systems, might be detectable with LISA in the frequency range 1–10 mHz and acts like a confusion-limited noise component, which might limit the LISA sensitivity at frequencies above 1 mHz.     

Minimal models of cooling neutron stars with accreted envelopes

We study the 'minimal' cooling scenario of superfluid neutron stars with nucleon cores, where the direct Urca process is forbidden and enhanced cooling is produced by neutrino emission due to the Cooper pairing of neutrons. Extending our recent previous work, we include the effects of surface accreted envelopes of light elements. We employ the phenomenological density-dependent critical temperatures   T _cp(ρ)  and   T _cnt(ρ)  of singlet-state proton and triplet-state neutron pairing in a stellar core, as well as the critical temperature   T _cns(ρ)  of singlet-state neutron pairing in a stellar crust. We show that the presence of accreted envelopes simplifies the interpretation of observations of thermal radiation from isolated neutron stars in the scenario of our recent previous work and widens the class of models for nucleon superfluidity in neutron star interiors consistent with the observations.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

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.     

Resistive relaxation of a magnetically confined mountain on an accreting neutron star

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to  τ_I∼ 10⁵–10⁸ yr  (depending on the conductivity of the neutron star crust) for an accreted mass of   M _a= 1.2 × 10⁻⁴ M_⊙  . The magnetic dipole moment simultaneously re-emerges as the screening currents dissipate over  τ_I  . For non-axisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfvén time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease  τ_I  by one order of magnitude relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.     

Radiative transitions of the helium atom in highly magnetized neutron star atmospheres

Recent observations of thermally emitting isolated neutron stars revealed spectral features that could be interpreted as radiative transitions of He in a magnetized neutron star atmosphere. We present Hartree–Fock calculations of the polarization-dependent photoionization cross-sections of the He atom in strong magnetic fields ranging from 10¹² to 10¹⁴ G. Convenient fitting formulae for the cross-sections are given along with the related oscillator strengths for various bound–bound transitions. The effects of finite nucleus mass on the radiative absorption cross-sections are examined using perturbation theory.     

The EC 14026 stars – VIII. PG 1336−018: a pulsating sdB star in an HW Vir-type eclipsing binary

The sdB star PG 1336−018 is found to be a very short-period eclipsing binary system, remarkably similar to the previously unique system HW Vir. In addition, and unlike HW Vir, the sdB star in the PG 1336 system shows rapid oscillations of the type found in the recently discovered sdB pulsators, or EC 14026 stars. The orbital period, 0.101 0174 d, is one of the shortest known for a detached binary. Analysis of photoelectric and CCD photometry reveals pulsation periods near 184 and 141 s, with semi-amplitudes of ∼0.01 and ∼0.005 mag respectively. Both oscillations might have variable amplitude, and it is probable that other frequencies are present with amplitudes ∼0.003 mag or less. The 184- and 141-s pulsations are in the range of periods predicted by models for hot horizontal-branch stars. Analysis of medium-dispersion spectrograms yields T _eff=33 000±1000 K and log g =5.7±0.1 for the sdB primary star, a radial velocity semi-amplitude K ₁=78±3 km s⁻¹ and a system velocity γ=6±2 km s⁻¹. Spectrograms from the IUE Final Archive give T _eff=33 000±3000 K and E ( B − V )=0.05 for log g =6.0 models. The derived angular radius leads to a distance of 710±50 pc for the system, and an absolute magnitude for the sdB star of +4.1±0.2. A preliminary analysis of U , V and R light curves indicates the orbital inclination to be near 81° and the relative radii to be r ₁=0.19 and r ₂=0.205. Assuming the mass of the sdB primary to be 0.5 M⊙ leads to a mass ratio q =0.3 for the system, and indicates that the secondary is a late-type dwarf of type ∼M5. As with HW Vir, it is necessary to invoke small limb-darkening coefficients and high albedos for the secondary star to obtain reasonable fits to the observed light curves.     

Early evolution of newly born magnetars with a strong toroidal field

We present a state-of-the-art scenario for newly born magnetars as strong sources of gravitational waves (GWs) in the early days after formation. We address several aspects of the astrophysics of rapidly rotating, ultra-magnetized neutron stars (NSs), including early cooling before transition to superfluidity, the effects of the magnetic field on the equilibrium shape of NSs, the internal dynamical state of a fully degenerate, oblique rotator and the strength of the electromagnetic torque on the newly born NS. We show that our scenario is consistent with recent studies of supernova remnant surrounding Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) in the Galaxy that constrains the electromagnetic energy input from the central NS to be  ≤ 10⁵¹  erg. We further show that if this condition is met, then the GW signal from such sources is potentially detectable with the forthcoming generation of GW detectors up to Virgo cluster distances where an event rate ∼1 yr⁻¹ can be estimated. Finally, we point out that the decay of an internal magnetic field in the 10¹⁶ G range couples strongly with the NS cooling at very early stages, thus significantly slowing down both processes: the field can remain this strong for at least 10³ yr, during which the core temperature stays higher than several times 10⁸ K.     

On the time dependence of differential rotation in young late-type stars

A model for the angular momentum transfer within the convection zone of a rapidly rotating star is introduced and applied to the analysis of recent observations of temporal fluctuations of the differential rotation on the young late-type stars AB Doradus (AB Dor) and LQ Hydrae (LQ Hya). Under the hypothesis that the mean magnetic field produced by the stellar dynamo rules the angular momentum exchanges and that the angular velocity depends only on the distance s from the rotation axis and the time, the minimum azimuthal Maxwell stress  | B_sB _φ|  , averaged over the convection zone, is found to range from ∼0.04 to  ∼0.14 T²  . If the poloidal mean magnetic field   B _s   is of the order of 0.01 T, as indicated by the Zeeman–Doppler imaging maps of those stars, then the azimuthal mean field   B _φ  can reach an intensity of several teslas, which significantly exceeds equipartition with the turbulent kinetic energy. Such strong fields can account also for the orbital period modulation observed in cataclysmic variables and RS Canum Venaticorum systems with a main-sequence secondary component. Moreover, the model allows us to compute the kinetic energy dissipation rate during the maintenance of the differential rotation. Only in the case of the largest surface shear observed on LQ Hya may the dissipated power exceed the stellar luminosity, but the lack of a sufficient statistic on the occurrence of such episodes of large shear does not allow us to estimate their impact on the energy budget of the convection zone.     

Model pulsar magnetospheres: the perpendicular rotator

The simplest model illustrating the effect of the magnetospheric charge-current field on the structure of a pulsar magnetic field has the region within the light-cylinder filled with the GoldreichJulian charge density which corotates with the neutron star, but has no electric currents along the magnetic field lines. This model has previously been studied for the axisymmetric case, with the rotation and magnetic dipolar axes aligned. The analogous problem is now solved with the two axes mutually perpendicular, so that not only the material current arising from the rotating charges but also the displacement current contributes. Again, the constructed magnetic field B ₀ crosses the light-cylinder normally, and there is no energy flux to infinity. However, in a more realistic model there is a flow of current along B ₀, generating a field B ₁ which has a non-vanishing toroidal component at the light-cylinder, so yielding a finite integrated Poynting flux.     

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).     

Turbulent dynamo in asymptotic giant branch stars

Using recent results on the operation of turbulent dynamos, we show that a turbulent dynamo may amplify a large-scale magnetic field in the envelopes of asymptotic giant branch (AGB) stars. We propose that a slow rotation of the AGB envelope can fix the symmetry axis, leading to the formation of an axisymmetric magnetic field structure. Unlike solar-type αω dynamos, the rotation has only a small role in amplifying the toroidal component of the magnetic field; instead of an αω dynamo we propose an α ² ω . The magnetic field may reach a value of     , where B _e is the equipartition (between the turbulent and magnetic energy densities) magnetic field. The large-scale magnetic field is strong enough for the formation of magnetic cool spots on the AGB stellar surface. The spots may regulate dust formation, and hence the mass-loss rate, leading to axisymmetric mass loss and the formation of elliptical planetary nebulae (PNe). Despite its role in forming cool spots, the large-scale magnetic field is too weak to play a dynamic role and directly influence the wind from the AGB star, as required by some models. We discuss other possible problems in models where the magnetic field plays a dynamic role in shaping the AGB winds, and argue that they cannot explain the formation of non-spherical PNe.     

Cooling of quark stars in the colour superconductive phase: effect of photons from glueball decay

The cooling history of a quark star in the colour superconductive phase is investigated. Here we specifically focus on the two-flavour colour (2SC) phase where the novel process of photon generation via glueball (GLB) decay has already been investigated. The picture we present here can, in principle, be generalized to quark stars entering a superconductive phase where similar photon generation mechanisms are at play. As much as 10⁴⁵–10⁴⁷ erg of energy is provided by the GLB decay in the 2SC phase. The generated photons slowly diffuse out of the quark star, keeping it hot and radiating as a blackbody (with possibly a Wien spectrum in gamma-rays) for millions of years. We discuss hot radio-quiet isolated neutron stars in our picture (such as RX J185635–3754 and RX J0720.4–3125) and argue that their nearly blackbody spectra (with a few broad features) and their remarkably tiny hydrogen atmosphere are indications that these might be quark stars in the colour superconductive phase where some sort of photon generation mechanism (reminiscent of the GLB decay) has taken place. Fits to observed data of cooling compact stars favour models with superconductive gaps of  Δ_2SC∼ 15–35 MeV  and densities  ρ_2SC= (2.5–3.0) ×ρ_N  (ρ_N being the nuclear matter saturation density) for quark matter in the 2SC phase. If correct, our model combined with more observations of isolated compact stars could provide vital information to studies of quark matter and its exotic phases.     

A numerical study of the r-mode instability of rapidly rotating nascent neutron stars

The first results of numerical analysis of classical r-modes of rapidly rotating compressible stellar models are reported. The full set of linear perturbation equations of rotating stars in Newtonian gravity is solved numerically without the slow rotation approximation. A critical curve of gravitational wave emission induced instability, which restricts the rotational frequencies of hot young neutron stars, is obtained. Taking the standard cooling mechanisms of neutron stars into account, we also show the 'evolutionary curves' along which neutron stars are supposed to evolve as cooling and spinning down proceed. Rotational frequencies of 1.4-M_⊙ stars suffering from this instability decrease to around 100 Hz when the standard cooling mechanism of neutron stars is employed. This result confirms the results of other authors, who adopted the slow rotation approximation.     

The secondary stars in cataclysmic variables and low-mass X-ray binaries

We critically re-examine the available data on the spectral types, masses and radii of the secondary stars in cataclysmic variables (CVs) and low-mass X-ray binaries (LMXBs), using the new catalogue of Ritter &38; Kolb as a starting point. We find there are 55 reliable spectral type determinations and only 14 reliable mass determinations of CV secondary stars (10 and 5, respectively, in the case of LMXBs). We derive new spectral type–period, mass–radius, mass–period and radius–period relations, and compare them with theoretical predictions. We find that CV secondary stars with orbital periods shorter than 7–8 h are, as a group, indistinguishable from main-sequence stars in detached binaries. We find that it is not valid, however, to estimate the mass from the spectral type of the secondary star in CVs or LMXBs. We find that LMXB secondary stars show some evidence for evolution, with secondary stars which are slightly too large for their mass. We show how the masses and radii of the secondary stars in CVs can be used to test the validity of the disrupted magnetic braking model of CV evolution, but we find that the currently available data are not sufficiently accurate or numerous to allow such an analysis. As well as considering secondary star masses, we also discuss the masses of the white dwarfs in CVs, and find mean values of ⁻ M  = 0.69 ± 0.13 M_⊙ below the period gap, and ⁻ M  = 0.80 ± 0.22 M_⊙ above the period gap.