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.     

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

Boundary layer on the surface of a neutron star

In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take   M _star= 1.4 × 10⁻² M_⊙  and for the gas temperature   T = 5 × 10⁶ K  . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 10⁸ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.     

Excitation of oscillation modes by tides in close binaries: constraints on stellar and orbital parameters

The parameter space favourable for the resonant excitation of free oscillation modes by dynamic tides in close binary components is explored using qualitative considerations to estimate the order of magnitude of the tidal force and the frequency range covered by the tidally induced oscillations. The investigation is valid for slowly rotating stars with masses in the interval between 2 and  20 M_⊙  , and an evolutionary stage ranging from the beginning to the end of the main sequence. Oscillation modes with eigenfrequencies of the order of five times the inverse of the dynamical time-scale  τ_dyn  of the star, i.e. the lowest-order p -modes, the f -mode and the lowest-order g ⁺-modes, are found to be outside the favourable parameter space since their resonant excitation requires orbital eccentricities that are too high for the binary to stay detached when the components pass through the periastron of their relative orbit. Resonances between dynamic tides and g ⁺-modes with frequencies of the order of half of the inverse of the dynamical time-scale of the star on the other hand are found to be favourable for orbital periods up to  ∼200τ_dyn  , provided that the binary mass ratio q is larger than 1/3, and the orbital eccentricity e is larger than ∼0.25. This favourable range comes down to orbital periods of up to 5–12 d in the case of  2–20 M_⊙  zero-age main-sequence binary components, and orbital periods of up to 21–70 d in the case of terminal main-sequence binary components.     

The bottom magnetic field and magnetosphere evolution of neutron star in low-mass X-ray binary

The accretion-induced neutron star (NS) magnetic field evolution is studied through considering the accretion flow to drag the field lines aside and dilute the polar-field strength, and as a result the equatorial field strength increases, which is buried inside the crust on account of the accretion-induced global compression of star crust. The main conclusions of model are as follows: (i) the polar field decays with increase in the accreted mass; (ii) the bottom magnetic field strength of about 10⁸ G can occur when the NS magnetosphere radius approaches the star radius, and it depends on the accretion rate as     ; and (iii) the NS magnetosphere radius decreases with accretion until it reaches the star radius, and its evolution is little influenced by the initial field and the accretion rate after accreting  ∼0.01 M_⊙  , which implies that the magnetosphere radii of NSs in low-mass X-ray binaries would be homogeneous if they accreted the comparable masses. As an extension, the physical effects of the possible strong magnetic zone in the X-ray NSs and recycled pulsars are discussed. Moreover, the strong magnetic fields in the binary pulsars PSR 1831−00 and PSR 1718−19 after accreting about  0.5 M_⊙  in the binary-accretion phase,  8.7 × 10¹⁰  and  1.28 × 10¹² G  , respectively, can be explained through considering the incomplete frozen flow in the polar zone. As an expectation of the model, the existence of the low magnetic field  (∼3 × 10⁷ G)  NSs or millisecond pulsars is suggested.     

Optical spectroscopy and photometry of SAX J1808.4−3658 in outburst

We present phase resolved optical spectroscopy and photometry of V4580 Sagittarii, the optical counterpart to the accretion powered millisecond pulsar SAX J1808.4−3658, obtained during the 2008 September/October outburst. Doppler tomography of the N  iii λ4640.64 Bowen blend emission line reveals a focused spot of emission at a location consistent with the secondary star. The velocity of this emission occurs at  324 ± 15 km s⁻¹  ; applying a ' K -correction', we find the velocity of the secondary star projected on to the line of sight to be  370 ± 40 km s⁻¹  . Based on existing pulse timing measurements, this constrains the mass ratio of the system to be  0.044^+0.005_−0.004  , and the mass function for the pulsar to be  0.44^+0.16_−0.13 M_⊙  . Combining this mass function with various inclination estimates from other authors, we find no evidence to suggest that the neutron star in SAX J1808.4−3658 is more massive than the canonical value of  1.4 M_⊙  . Our optical light curves exhibit a possible superhump modulation, expected for a system with such a low mass ratio. The equivalent width of the Ca  ii H and K interstellar absorption lines suggest that the distance to the source is ∼2.5 kpc. This is consistent with previous distance estimates based on type-I X-ray bursts which assume cosmic abundances of hydrogen, but lower than more recent estimates which assume helium-rich bursts.     

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.     

The correlations between the twin kHz quasi-periodic oscillation frequencies of low-mass X-ray binaries

We analyzed the recently published kHz quasi-period oscillaiton (QPO) data in the neutron star low-mass X-ray binaries (LMXBs), in order to investigate the different correlations of the twin-peak kHz QPOs in bright Z sources and in the less luminous Atoll sources. We find a power-law relation  ν₁∼ν ^b ₂  between the upper and the lower kHz QPOs with different indices: b ≃ 1.5 for the Atoll source 4U 1728-34 and b ≃ 1.9 for the Z source Sco X-1. The implications of our results for the theoretical models for kHz QPOs are discussed.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

The γ-ray burst phenomenon treated as the collapse of a QED magnetized vacuum bubble: analogy with sonoluminescence

We treat the phenomenon of a γ -ray burst as the non-linear collapse of a magnetic cavity surrounding a neutron star with extremely large magnetic field B ∼10¹⁵–10¹⁶ G due to the process of bubble shape instability in the resonant MHD field of an accreting plasma or on a neutron star surface. The QED effect of vacuum polarizability by a strong magnetic field is taken into a consideration. We develop an analogy with the phenomenon of sonoluminescence in which the gas bubble is located in a surrounding liquid with a driven sound intensity.     

The influence of Hall drift on the magnetic fields of neutron stars

We consider the evolution of magnetic fields under the influence of Hall drift and Ohmic decay. The governing equation is solved numerically, in a spherical shell with   r _i / r _o = 0.75  . Starting with simple free-decay modes as initial conditions, we then consider the subsequent evolution. The Hall effect induces so-called helicoidal oscillations, in which energy is redistributed among the different modes. We find that the amplitude of these oscillations can be quite substantial, with some of the higher harmonics becoming comparable with the original field. Nevertheless, this transfer of energy to the higher harmonics is not sufficient to accelerate significantly the decay of the original field, at least not at the   R _B = O (100)  parameter values accessible to us, where this Hall parameter   R _B   measures the ratio of the Ohmic time-scale to the Hall time-scale. We do find clear evidence though of increasingly fine structures developing for increasingly large   R _B   , suggesting that perhaps this Hall-induced cascade to ever-shorter length-scales is eventually sufficiently vigorous to enhance the decay of the original field. Finally, the implications for the evolution of neutron star magnetic fields are discussed.     

An iterative filter to reconstruct planetary transit signals in the presence of stellar variability

We report the identification, from a photometric, astrometric and spectroscopic study, of a massive white dwarf member of the nearby, approximately solar metallicity, Coma Berenices open star cluster (Melotte 111). We find the optical to near-infrared energy distribution of WD 1216+260 to be entirely consistent with that of an isolated DA and determine the effective temperature and surface gravity of this object to be   T _eff= 15 739⁺¹⁹⁷₋₁₉₆ K  and  log  g = 8.46^+0.03_−0.02  . We set tight limits on the mass of a putative cool companion,   M ≳ 0.036 M_⊙  (spatially unresolved) and   M ≳ 0.034 M_⊙  (spatially resolved and   a ≲ 2500 au  ). Based on the predictions of CO core, thick H layer evolutionary models we determine the mass and cooling time of WD 1216+260 to be   M _WD= 0.90 ± 0.04 M_⊙  and  τ_cool= 363⁺⁴⁶₋₄₁ Myr  , respectively. For an adopted cluster age of  τ= 500 ± 100 Myr  we infer the mass of its progenitor star to be   M _init= 4.77^+5.37_−0.97 M_⊙  . We briefly discuss this result in the context of the form of the stellar initial mass–final mass relation.     

Kelvin–Helmholtz driven propeller in AE Aquarii: A unified model for thermal and non-thermal flares

In this paper, an attempt is made to integrate the propeller ejection of material by the fast rotating white dwarf in AE Aquarii with the highly transient thermal and non-thermal emission in a single unifying model. It has been shown that the violent interaction between the fast rotating magnetosphere and a clumpy fragmented stream, in AE Aquarii specifically, may result in the growth of unstable modes of the Kelvin–Helmholtz instability and associated turbulence over length scales comparable to the stream radius on time-scales  τ_K-H ∼  t _dyn (∼ 600 s  ). For all conversion efficiencies of magnetohydrodynamic (MHD) power to mechanical energy ε≥ 0.1, these instabilities result in the effective azimuthal acceleration of the gas parcels to the escape velocity over time-scales   t _acc≤ 1000 s (∼ t _dyn)  . Further, it has been shown that the turbulence in the flow will cascade down to the dissipative level over time-scales  τ_cas∼ 3 h  . If released through dissipative shocks, this reservoir can drive a luminosity of   L ∼ 10³³ erg s⁻¹  , which can significantly contribute to the total emission when blobs collide in the exit stream, resulting in shock heating and associated flares. During the propeller process, particles can also be accelerated to high energies, which may be the driving mechanism behind the non-thermal radio to mid-infrared emission. The confluence of these ejected magnetized clouds may result in radio remnant surrounding AE Aquarii, which is optically thin between frequencies ν≥ 100 MHz–1 GHz.     

Optical spectroscopy of the quiescent counterpart to EXO 0748−676

We present phase resolved optical spectroscopy and X-ray timing of the neutron star X-ray binary EXO 0748−676 after the source returned to quiescence in the autumn of 2008. The X-ray light curve displays eclipses consistent in orbital period, orbital phase and duration with the predictions and measurements before the return to quiescence. Hα and He  i emission lines are present in the optical spectra and show the signature of the orbit of the binary companion, placing a lower limit on the radial velocity semi-amplitude of   K ₂ > 405 km s⁻¹  . Both the flux in the continuum and the emission lines show orbital modulations, indicating that we observe the hemisphere of the binary companion that is being irradiated by the neutron star. Effects due to this irradiation preclude a direct measurement of the radial velocity semi-amplitude of the binary companion; in fact, no stellar absorption lines are seen in the spectrum. Nevertheless, our observations place a stringent lower limit on the neutron star mass of   M ₁ > 1.27 M_⊙  . For the canonical neutron star mass of   M ₁= 1.4 M_⊙  , the mass ratio is constrained to  0.075 < q < 0.105  .     

Periodicities in the high-mass X-ray binary system RX J0146.9+6121/LS I+61°235

The high-mass X-ray binary RX J0146.9+6121, with optical counterpart LS I+61°235 (V831 Cas), is an intriguing system on the outskirts of the open cluster NGC 663. It contains the slowest Be type X-ray pulsar known with a pulse period of around 1400 s and, primarily from the study of variation in the emission line profile of Hα, it is known to have a Be decretion disc with a one-armed density wave period of approximately 1240 d. Here we present the results of an extensive photometric campaign, supplemented with optical spectroscopy, aimed at measuring short time-scale periodicities. We find three significant periodicities in the photometric data at, in order of statistical significance, 0.34, 0.67 and 0.10 d. We give arguments to support the interpretation that the 0.34 and 0.10 d periods could be due to stellar oscillations of the B-type primary star and that the 0.67 d period is the spin period of the Be star with a spin axis inclination of  23⁺¹⁰₋₈  degrees. We measured a systemic velocity of  −37.0 ± 4.3 km s⁻¹  confirming that LS I+61°235 has a high probability of membership in the young cluster NGC 663 from which the system's age can be estimated as 20–25 Myr. From archival RXTE All Sky Monitor (ASM) data we further find 'super' X-ray outbursts roughly every 450 d. If these super outbursts are caused by the alignment of the compact star with the one-armed decretion disc enhancement, then the orbital period is approximately 330 d.     

Explosive reconnection in magnetars

The X-ray activity of anomalous X-ray pulsars and soft γ-ray repeaters may result from the heating of their magnetic corona by direct currents dissipated by magnetic reconnection. We investigate the possibility that X-ray flares and bursts observed from anomalous X-ray pulsars and soft γ-ray repeaters result from magnetospheric reconnection events initiated by development of the tearing mode in magnetically dominated relativistic plasma. We formulate equations of resistive force-free electrodynamics, discuss the relation of the latter to ideal electrodynamics, and give examples of both ideal and resistive equilibria. Resistive force-free current layers are unstable towards the development of small-scale current sheets where resistive effects become important. Thin current sheets are found to be unstable due to the development of the resistive force-free tearing mode. The growth rate of the tearing mode is intermediate between the short Alfvén time-scale  τ_A  and a long resistive time-scale  τ_R: Γ∼ 1/(τ_Rτ_A)^1/2  , similar to the case of non-relativistic non-force-free plasma. We propose that growth of the tearing mode is related to the typical rise time of flares, ∼10 ms. Finally, we discuss how reconnection may explain other magnetar phenomena and ways to test the model.     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

What is the nature of RX J0720.4–3125?

RX J0720.4–3125 has recently been identified as a pulsating soft X-ray source in the ROSAT all-sky survey with a period of 8.391 s. Its spectrum is well characterized by a blackbody with a temperature of 8 × 10⁵ K. We propose that the radiation from this object is thermal emission from a cooling neutron star. For this blackbody temperature we can obtain a robust estimate of the object's age of ∼ 3 × 10⁵ yr, yielding a polar field ∼ 10¹⁴ G for magnetic dipole spin-down and a value of P compatible with current observations.     

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.     

The boundary layer of VW Hydri in quiescence

In this short paper, we suggest that the missing boundary layer luminosity of dwarf novae in quiescence is released mainly in the ultraviolet (UV) as the second component commonly identified in the far-UV as the 'accretion belt'. We present the well-studied SU UMa-type system VW Hyi in detail as a prototype for such a scenario. We consider detailed multiwavelength observations and in particular the recent Far Ultraviolet Spectroscopic Explorer ( FUSE ) observations of VW Hyi in quiescence, which confirm the presence of a second component (the 'accretion belt') in the far-UV spectrum. The temperature  (≈ 50 000 K)  and rotational velocity  (≈ 3000 km s ⁻¹)  of this second component are entirely consistent with the optically thick region  (τ≈ 1)  located just at the outer edge of optically thin boundary layer in the simulations of Popham. This second component contributes about 20 per cent of the far-UV flux. Using geometrical assumptions and taking into account the X-ray luminosity, we find that the total boundary layer luminosity sums up to   L _BL= (0.53 ± 0.25) L _disc  , while the theory (Kluźniak) predicts, for the rotation rate of VW Hyi's white dwarf,   L _BL≈ (0.76 ± 0.03) L _disc  . About one-fifth of the boundary layer energy is emitted in the X-ray and the remaining is emitted in the UV. This scenario is consistent with the recent simultaneous X-ray and UV observations of VW Hyi by Pandel, Córdova & Howell, from which we deduce here that the viscosity in the boundary layer region must be of the order of  ν≈ 10¹³–10¹⁴ cm² s ⁻¹  , depending on the white dwarf mass and the size of the boundary layer.