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.     

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.     

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.     

The continuing saga of Sakurai's object (V4334 Sgr): dust production and helium line emission

We report further UKIRT spectroscopic observations of Sakurai's object (V4334 Sgr) made in 1999 April/May in the 1–4.75 μm range, and find that the emission is dominated by amorphous carbon at T _d~600 K. The estimated maximum grain size is 0.6 μm, and the mass lower limit is 1.7±0.2×10⁻⁸ M_⊙ to 8.9±0.6×10⁻⁷ M_⊙ for distances of 1.1–8 kpc. For 3.8 kpc the mass is 2.0±0.1×10⁻⁷ M_⊙.
We also report strong He  i emission at 1.083 μm, in contrast to the strong absorption in this line in 1998. We conclude that the excitation is collisional, and is probably caused by a wind, consistent with the P Cygni profile observed by Eyres et al. in 1998.     

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.     

An extremely slow nova?

The binary companion to the peculiar F supergiant HD 172481 is shown to be a Mira variable with a pulsation period of 312 d. Its characteristics are within the normal range found for solitary Miras of that period, although its pulsation amplitude and mass-loss rate ̇ ∼3×10⁻⁶ M_⊙ yr⁻¹ are higher than average. Reasons are given for suspecting that the F supergiant, which has L ∼10⁴ L_⊙, is a white dwarf burning hydrogen accreted from its companion.     

Cygnus X-3 with ISO: investigating the wind

We have observed the energetic binary Cygnus X-3 in both quiescent and flaring states between 4 and 16 μm using the ISO satellite. We find that the quiescent source shows the thermal free–free spectrum typical of a hot, fast stellar wind, such as from a massive helium star. The quiescent mass-loss rate arising from a spherically symmetric, non-accelerating wind is found to be in the range (0.4–2.9)×10⁻⁴ M_⊙ yr⁻¹, consistent with other infrared and radio observations, but considerably larger than the 10⁻⁵ M_⊙ yr⁻¹ deduced from both the orbital change and the X-ray column density. There is rapid, large-amplitude flaring at 4.5 and 11.5 μm at the same time as enhanced radio and X-ray activity, with the infrared spectrum apparently becoming flatter in the flaring state. We believe that non-thermal processes are operating, perhaps along with enhanced thermal emission.     

Spectroscopic and photometric observations of the eclipsing star UV Leo

High-resolution spectroscopic observations around the Hα line and BVRI photometry of the eclipsing short-period RS CVn star UV Leo are presented. The simultaneous light-curve solution and radial velocity-curve solution led to the following values of the global parameters of the binary: temperatures   T ₁= 6000 ± 100 K  and   T ₂= 5970 ± 20 K  ; masses   M ₁= 0.976 ± 0.067 M_⊙  and   M ₂= 0.931 ± 0.052 M_⊙  ; separation   a = 3.716 ± 0.048 R_⊙  ; orbital inclination     ; radii   R ₁= 1.115 ± 0.052 R_⊙  and   R ₂= 1.078 ± 0.051 R_⊙  ; equatorial velocities   V ₁= 98.8 ± 2.3 km s⁻¹  and   V ₂= 89.6 ± 2.7 km s⁻¹  . These results lead to the conclusion that the two components of UV Leo are slightly oversized for their masses and lie within the main-sequence band on the mass–radius diagram, close to the isochrone 9 × 10¹⁰ yr.     

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  .     

The component masses of the cataclysmic variable V347 Puppis

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variable V347 Pup (=LB 1800). There is evidence of irradiation on the inner hemisphere of the secondary star, which we correct for using a model to give a secondary-star radial velocity of   K _R= 198 ± 5 km s⁻¹  . The rotational velocity of the secondary star in V347 Pup is found to be   v sin  i = 131 ± 5 km s⁻¹  and the system inclination is   i = 840 ± 23  . From these parameters we obtain masses of   M ₁= 0.63 ± 0.04 M_⊙  for the white dwarf primary and   M ₂= 0.52 ± 0.06 M_⊙  for the M0.5V secondary star, giving a mass ratio of   q = 0.83 ± 0.05  . On the basis of the component masses, and the spectral type and radius of the secondary star in V347 Pup, we find tentative evidence for an evolved companion. V347 Pup shows many of the characteristics of the SW Sex stars, exhibiting single-peaked emission lines, high-velocity S-wave components and phase-offsets in the radial velocity curve. We find spiral arms in the accretion disc of V347 Pup and measure the disc radius to be close to the maximum allowed in a pressureless disc.     

KPD 0422 + 5421: a new short-period subdwarf B/white dwarf binary

The subdwarf B (sdB) star KPD 0422 + 5421 was discovered to be a single-lined spectroscopic binary with a period of P  = 0.090 1795 ± (3 × 10⁻⁷) d (2 h 10 min). The U B light curves display an ellipsoidal modulation with amplitudes of ≈ 0.02 mag. The sdB star contributes nearly all of the observed flux. This and the absence of any reflection effect suggest that the unseen companion star is small (i.e. R _comp ≈ 0.01 R_⊙) and therefore degenerate. We modelled the U B light curves and derived i  = 78.05° ± 0.50° and a mass ratio of q  =  M _comp/ M _sdB = 0.87 ± 0.15. The sdB star fills 69 per cent of its Roche lobe. These quantities may be combined with the mass function of the companion [ f ( M ) = 0.126 ± 0.028 M_⊙] to derive M _sdB = 0.72 ± 0.26 M_⊙ and M _comp = 0.62 ±  0.18 M_⊙. We used model spectra to derive the effective temperature, surface gravity and helium abundance of the sdB star. We found T _eff = 25 000 ± 1500 K, log g  = 5.4 ± 0.1 and [He/H] = −1.0. With a period of 2 h 10 min, KPD 0422 + 5421 has one of the shortest known orbital periods of a detached binary. This system is also one of only a few known binaries that contain a subdwarf B star and a white dwarf. Thus KPD 0422 + 5421 represents a relatively unobserved, and short-lived, stage of binary star evolution.     

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.     

The evolution of low-metallicity asymptotic giant branch stars and the formation of carbon-enhanced metal-poor stars

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities   Z = 10⁻⁴  and 10⁻⁸. We determine which stars undergo an episode of flash-driven mixing, where protons are ingested into the intershell convection zone, as they enter the thermally pulsing AGB phase and which undergo third dredge-up. We find that flash-driven mixing does not occur above a metallicity of   Z = 10⁻⁵  for any mass of star and that stars above  2 M_⊙  do not experience this phenomenon at any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of hydrogen-burning region, occurs in the mass range  2 M_⊙  to around  4 M_⊙  . We suggest that CI may be a weak version of the flash-driven mechanism. We also investigate the effects of convective overshooting on the behaviour of these objects. Our models struggle to explain the frequency of Carbon-Enhanced Metal-Poor (CEMP) stars that have both significant carbon and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI or just third dredge-up. Nitrogen can be enhanced through hot bottom burning and the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio may be a good indicator of the mass of the primary AGB stars.     

The local star formation rate and radio luminosity density

We present a new determination of the local volume-averaged star formation rate from the 1.4-GHz luminosity function of star forming galaxies. Our sample, taken from the   B ≤12  Revised Shapley–Ames catalogue (231 normal spiral galaxies over an effective area of 7.1 sr) has ≃100 per cent complete radio detections and is insensitive to dust obscuration and cirrus contamination. After removal of known active galaxies, the best-fitting Schechter function has a faint-end slope of  −1.27±0.07  in agreement with the local H α luminosity function, characteristic luminosity   L ∗=(2.6±0.7)×10²² W Hz⁻¹  and density   φ ∗=(4.8±1.1)×10⁻⁴ Mpc⁻³.  The inferred local radio luminosity density of  (1.73±0.37±0.03)×10¹⁹ W Hz⁻¹ Mpc⁻³  (Poisson noise, large-scale structure fluctuations) implies a volume-averaged star formation rate ∼2 times larger than the Gallego et al. H α estimate, i.e.   ρ _1.4 GHz=(2.10±0.45±0.04)×10⁻² M_⊙ yr⁻¹ Mpc⁻³  for a Salpeter initial mass function from  0.1–125 M_⊙  and Hubble constant of 50 km s⁻¹ Mpc⁻¹. We demonstrate that the Balmer decrement is a highly unreliable extinction estimator, and argue that optical–ultraviolet (UV) star formation rates (SFRs) are easily underestimated, particularly at high redshift.     

Type IIP supernova SN 2004et: a multiwavelength study in X-ray, optical and radio

We present X-ray, broad-band optical and low-frequency radio observations of the bright type IIP supernova SN 2004et. The Chandra X-ray Observatory observed the supernova at three epochs, and the optical coverage spans a period of ∼470 d since explosion. The X-ray emission softens with time, and we characterize the X-ray luminosity evolution as   L _X∝ t ^−0.4  . We use the observed X-ray luminosity to estimate a mass-loss rate for the progenitor star of  ∼2 × 10⁻⁶ M_⊙ yr⁻¹  . The optical light curve shows a pronounced plateau lasting for about 110 d. Temporal evolution of photospheric radius and colour temperature during the plateau phase is determined by making blackbody fits. We estimate the ejected mass of ⁵⁶Ni to be  0.06 ± 0.03 M_⊙  . Using the expressions of Litvinova & Nadëzhin we estimate an explosion energy of  (0.98 ± 0.25) × 10⁵¹ erg  . We also present a single epoch radio observation of SN 2004et. We compare this with the predictions of the model proposed by Chevalier, Fransson & Nymark. These multiwavelength studies suggest a main-sequence progenitor mass of  ∼20 M_⊙  for SN 2004et.     

Masses for the Local Group and the Milky Way

We use the very large Millennium Simulation of the concordance Λ cold dark matter cosmogony to calibrate the bias and error distribution of Timing Argument estimators of the masses of the Local Group and of the Milky Way. From a large number of isolated spiral–spiral pairs similar to the Milky Way/Andromeda system, we find the interquartile range of the ratio of timing mass to true mass to be a factor of 1.8, while the 5 and 95 per cent points of the distribution of this ratio are separated by a factor of 5.7. Here, we define true mass as the sum of the 'virial' masses, M ₂₀₀, of the two dominant galaxies. For present best values of the distance and approach velocity of Andromeda, this leads to a median likelihood estimate of the true mass of the Local Group of  5.27 × 10¹² M_⊙  or  log  M _LG/M_⊙= 12.72  , with an interquartile range of [12.58, 12.83] and a 5–95 per cent range of [12.26, 13.01]. Thus, a 95 per cent lower confidence limit on the true mass of the Local Group is  1.81 × 10¹² M_⊙  . A timing estimate of the Milky Way's mass based on the large recession velocity observed for the distant satellite Leo I works equally well, although with larger systematic uncertainties. It gives an estimated virial mass for the Milky Way of  2.43 × 10¹² M_⊙  with a 95 per cent lower confidence limit of  0.80 × 10¹² M_⊙  .     

Characterizing the pulsations of the ZZ Ceti star KUV 02464+3239

We present the results on period search and modelling of the cool DAV star KUV 02464+3239. Our observations resolved the multiperiodic pulsational behaviour of the star. In agreement with its position near the red edge of the DAV instability strip, it shows large amplitude, long-period pulsation modes, and has a strongly non-sinusoidal light curve. We determined six frequencies as normal modes and revealed remarkable short-term amplitude variations. A rigorous test was performed for the possible source of amplitude variation: beating of modes, effect of noise, unresolved frequencies or rotational triplets. Among the best-fitting models resulting from a grid search, we selected three that gave   l = 1  solutions for the largest amplitude modes. These models had masses of 0.645, 0.650 and  0.680  M _⊙  . The three 'favoured' models have   M _H  between  2.5 × 10⁻⁵ and 6.3 × 10⁻⁶  M _*  and give 14.2–14.8 mas seismological parallax. The  0.645  M _⊙  (11 400 K) model also matches the spectroscopic  log  g   and   T _eff  within 1σ. We investigated the possibility of mode trapping and concluded that while it can explain high amplitude modes, it is not required.     

Quantifying the coexistence of massive black holes and dense nuclear star clusters

In large spheroidal stellar systems, such as elliptical galaxies, one invariably finds a  10⁶–10⁹ M_⊙  supermassive black hole at their centre. In contrast, within dwarf elliptical galaxies one predominantly observes a  10⁵–10⁷ M_⊙  nuclear star cluster. To date, few galaxies have been found with both types of nuclei coexisting and even less have had the masses determined for both central components. Here, we identify one dozen galaxies housing nuclear star clusters and supermassive black holes whose masses have been measured. This doubles the known number of such hermaphrodite nuclei – which are expected to be fruitful sources of gravitational radiation. Over the host spheroid (stellar) mass range  10⁸–10¹¹ M_⊙  , we find that a galaxy's nucleus-to-spheroid (baryon) mass ratio is not a constant value but decreases from a few per cent to ∼0.3 per cent such that  log[( M _BH+ M _NC)/ M _sph]=−(0.39 ± 0.07) log[ M _sph/10¹⁰ M_⊙]− (2.18 ± 0.07)  . Once dry merging commences and the nuclear star clusters disappear, this ratio is expected to become a constant value.
As a byproduct of our investigation, we have found that the projected flux from resolved nuclear star clusters is well approximated with Sérsic functions having a range of indices from ∼0.5 to ∼3, the latter index describing the Milky Way's nuclear star cluster.     

Mass transfer in eccentric binary stars

The concept of Roche lobe overflow is fundamental to the theory of interacting binaries. Based on potential theory, it is dependent on all the relevant material corotating in a single frame of reference. Therefore if the mass losing star is asynchronous with the orbital motion or the orbit is eccentric, the simple theory no longer applies and no exact analytical treatment has been found. We use an analytic approximation whose predictions are largely justified by smoothed particle hydrodynamic simulations (SPH). We present SPH simulations of binary systems with the same semi-major axis   a = 5.55 R_⊙  , masses   M ₁= 1 M_⊙, M ₂= 2 M_⊙  and radius   R ₁= 0.89 R_⊙  for the primary star but with different eccentricities   e = 0.4, 0.5, 0.6  and 0.7. In each case the secondary star is treated as a point mass. When   e = 0.4  no mass is lost from the primary while at   e = 0.7  catastrophic mass transfer, partly through the L₂ point, takes place near periastron. This would probably lead to common-envelope evolution if star 1 were a giant or to coalescence for a main-sequence star. In between, at   e ≥ 0.5  , some mass is lost through the L₁ point from the primary close to periastron. However, rather than being all accreted by the secondary, some of the stream appears to leave the system. Our results indicate that the radius of the Roche lobe is similar to circular binaries when calculated for the separation and angular velocity at periastron. Part of the mass loss occurs through the L₂ point.     

The radial velocity of the companion star in the low-mass X-ray binary 2S 0921–630: limits on the mass of the compact object

In this paper we report on optical spectroscopic observations of the low-mass X-ray binary 2S 0921–630 obtained with the Very Large Telescope. We found sinusoidal radial velocity variations of the companion star with a semi-amplitude of  99.1 ± 3.1 km s⁻¹  modulated on a period of 9.006 ± 0.007 d, consistent with the orbital period found previously for this source, and a systemic velocity of  44.4 ± 2.4 km s⁻¹  . Owing to X-ray irradiation, the centre of light measured by the absorption lines from the companion star is probably shifted with respect to the centre of mass. We try to correct for this using the so-called K -correction. Conservatively applying the maximum correction possible and using the previously measured rotational velocity of the companion star, we find a lower limit to the mass of the compact object in 2S 0921–630 of   M_X sin³ i > 1.90 ± 0.25 M_⊙  (1σ errors). The inclination in this system is well constrained since partial eclipses have been observed in X-ray and optical bands. For inclinations in the range  60° < i < 90°  we find  1.90 ± 0.25 < M_X < 2.9 ± 0.4 M_⊙  . However, using this maximum K -correction we find that the ratio between the mass of the companion star and that of the compact object, q , is 1.32 ± 0.37, implying super-Eddington mass-transfer rates; however, evidence for that has not been found in 2S 0921–630. We conclude that the compact object in 2S 0921–630 is either a (massive) neutron star or a low-mass black hole.