The mass of the neutron star in the binary millisecond pulsar PSR J 1012 + 5307

We have measured the radial velocity variation of the white dwarf secondary in the binary system containing the millisecond pulsar PSR J 1012 + 5307. Combined with the orbital parameters of the radio pulsar, we infer a mass ratio q (≡ M ₁/ M ₂) = 10.5 ± 0.5. Our optical spectroscopy has also allowed us to determine the mass of the white dwarf companion by fitting the spectrum to a grid of DA model atmospheres: we estimate M ₂ = 0.16 ± 0.02 M⊙, and hence the mass of the neutron star is 1.64 ± 0.22 M⊙, where the error is dominated by that of M ₂. The orbital inclination is 52 ± 4°. For an initial neutron star mass of ∼ 1.4 M⊙, only a few tenths of a solar mass at most has been successfully accreted over the lifetime of the progenitor low-mass X-ray binary. If the initial mass of the secondary was ∼ 1 M⊙, our result suggests that the mass transfer may have been non-conservative.     

Full evolution of low-mass white dwarfs with helium and oxygen cores

We study the full evolution of low-mass white dwarfs with helium and oxygen cores. We revisit the age dichotomy observed in many white dwarf companions to millisecond pulsar on the basis of white dwarf configurations derived from binary evolution computations. We evolve 11 dwarf sequences for helium cores with final masses of 0.1604, 0.1869, 0.2026, 0.2495, 0.3056, 0.3333, 0.3515, 0.3844, 0.3986, 0.4160 and  0.4481 M_⊙  . In addition, we compute the evolution of five sequences for oxygen cores with final masses of 0.3515, 0.3844, 0.3986, 0.4160 and  0.4481 M_⊙  . A metallicity of   Z = 0.02  is assumed. Gravitational settling, chemical and thermal diffusion are accounted for during the white dwarf regime. Our study reinforces the result that diffusion processes are a key ingredient in explaining the observed age and envelope dichotomy in low-mass helium-core white dwarfs, a conclusion we arrived at earlier on the basis of a simplified treatment for the binary evolution of progenitor stars. We determine the mass threshold where the age dichotomy occurs. For the oxygen white dwarf sequences, we report the occurrence of diffusion-induced, hydrogen-shell flashes, which, as in the case of their helium counterparts, strongly influence the late stages of white dwarf cooling. Finally, we present our results as a set of white dwarf mass–radius relations for helium and oxygen cores.     

PSR J1012+5307: younger than it looks?

Lorimer et al. have recently reported that the spin-down age (∼7 × 10⁹ yr) of the low-mass binary pulsar PSR J1012+5307 is much higher than the cooling age (3 × 10⁸ yr) of its white dwarf companion. The proposed solutions for this discrepancy are outlined and discussed. In particular, the revised cooling age estimate proposed by Alberts et al. agrees with data from other low-mass binary pulsar systems if a transition to the 'classical' cooling regime occurs between ∼0.14 and ∼0.28 M_⊙. If this transition is excluded, PSR J1012+5307 seems to have finished its accretion phase far from the spin-up line.     

The impact of element diffusion on the formation and evolution of helium white dwarf stars

The aim of this work is to investigate the effect of element diffusion on the evolution of helium white dwarfs. To this end, we couple the multicomponent flow equations that describe gravitational settling, chemical and thermal diffusion to an evolutionary code. We compute the evolution of a set of helium white dwarf models with masses ranging from 0.169 to 0.406 M_⊙. In particular, several low-mass white dwarfs have been found in binary systems as companion to millisecond pulsars. In these systems, pulsar emission is activated by mass transfer episodes so that, if we place the zero-age point at the end of such mass transfer, then the pulsar and the white dwarf ages should be equal. Interestingly enough, available models of helium white dwarfs neglect element diffusion. Using such models, good agreement has been found between the ages of the components of the PSR J1012+5307 system. However, recent observations of the PSR B1855+09 system cast doubts on the correctness of such models, which predict a white dwarf age twice as long as the spin-down age of the pulsar. In this work, we find that element diffusion induces thermonuclear hydrogen shell flashes for models in the mass interval 0.18≲ M /M_⊙ ≲ 0.41 . We show, in particular, that the occurrence of these diffusion-induced flashes eventually leads to white dwarf models with hydrogen envelope masses too small to support any further nuclear burning, thus implying much shorter cooling ages than in the case when diffusion is neglected. In particular, excellent agreement is found between the ages of PSR B1855+09 system components, solving the age discrepancy from first principles.     

The case and fate of HD 75767 – neutron star or supernova?

We report the discovery of the nearby  ( d = 24 pc)  HD 75767 as an eight billion year old quadruple system consisting of a distant M dwarf pair, HD 75767 C–D, in orbit around the known short-period   P = 10.25 d  single-lined binary HD 75767 A–B, the primary of which is a solar-like G star. On the reasonable assumption of synchronous orbital rotation as well as rotational and orbital coplanarity for the inner pair, we get   M _B= 0.96 M_⊙  for the unseen HD 75767 B, that is, the case of a massive white dwarf. Upon future evolution, mass transfer towards HD 75767 B will render the   M _A= 0.96 M_⊙  G-type primary, now a turnoff star, to become a helium white dwarf of   M _A∼ 0.33 M_⊙  . Depending on the mass accretion rate, accretion efficiency and composition of the massive white dwarf, this in turn may result in a collapse of HD 75767 B with the formation of a millisecond pulsar, i.e. the creation of a low-mass binary pulsar (LMBP), or, instead, a Type Ia supernova explosion and the complete disruption of HD 75767 B. Irrespective of which scenario applies, we point to the importance of the distant M dwarfs as the likely agents for the formation of the inner, short-period HD 75767 A–B pair, and hence a path that particularly avoids preceding phases of common envelope evolution.     

Merged binary white dwarf evolution: rapidly accreting carbon–oxygen white dwarfs and the progeny of extreme helium stars

We have examined the evolution of merged low-mass double white dwarfs that become luminous helium stars. We have approximated the merging process by the rapid accretion of matter, consisting mostly of helium, on to a carbon–oxygen (CO) white dwarf. After a certain mass is accumulated, a helium shell flash occurs, the radius and luminosity increase and the star becomes a yellow giant. Mass accretion is stopped artificially when the total mass reaches a pre-determined value. When the mass above the helium-burning shell becomes small enough, the star evolves blueward almost horizontally in the Hertzsprung–Russell diagram. The theoretical models for the merger of a 0.6-M_⊙ CO white dwarf with a 0.3-M_⊙ He white dwarf agree very well with the observed locations of extreme helium stars in the  log  T _eff–log  g   diagram, with their observed rates of blueward evolution, and with luminosities and masses obtained from their pulsations. Together with predicted merger rates for  CO+He  white dwarf pairs, the evolutionary time-scales are roughly consistent with the observed numbers of extreme helium stars. Predicted surface carbon and oxygen abundances can be consistent with the observed values if carbon and oxygen produced in the helium shell during a previous asymptotic giant branch phase are assumed to exist in the helium zone of the initial CO white dwarfs. These results establish the  CO+He  white dwarf merger as the best, if not only, viable model for the creation of extreme helium stars and, by association, the majority of R Coronae Borealis stars.     

The possible companions of young radio pulsars

We discuss the formation of pulsars with massive companions in eccentric orbits. We demonstrate that the probability for a non-recycled radio pulsar to have a white dwarf as a companion is comparable to that of having an old neutron star as a companion. Special emphasis is given to PSR B1820−11 and PSR B2303+46. Based on population synthesis calculations we argue that PSR B1820−11 and PSR B2303+46 could very well be accompanied by white dwarfs with mass ≳1.1 M_⊙. For PSR B1820−11, however, we cannot exclude the possibility that its companion is a main-sequence star with a mass between ∼0.7 M_⊙ and ∼5 M_⊙.     

Stellar forensics — I. Cooling curves

The presence of low-mass, degenerate secondaries in millisecond pulsar binaries offers the opportunity to determine an age for the binary system independently of the rotational properties of the pulsar. To this end, we present here a detailed calculation of the evolution of a grid of low-mass (< 0.05 M⊙) helium core white dwarfs. We investigate the effects of different hydrogen layer masses and provide results for well-known optical bandpasses. We supplement the OPAL opacity calculations with our own calculations for low effective temperatures ( T _eff < 6000 K) and also provide fitting formulae for the gravity as a function of mass and effective temperature. In Paper II we apply these results to individual cases.     

Diffusion in helium white dwarf stars

This paper is aimed at exploring the effects of diffusion on the structure and evolution of low-mass helium white dwarfs. To this end, we solve the multicomponent flow equations describing gravitational settling and chemical and thermal diffusion. The diffusion calculations are coupled to an evolutionary code in order to follow the cooling of low-mass, helium core white dwarf models having envelopes made up of a mixture of hydrogen and helium, as recently suggested by detailed evolutionary calculations for white dwarf progenitors in binary systems. We find that diffusion causes hydrogen to float and the other elements to sink over time-scales shorter than evolutionary time-scales. This produces a noticeable change in the structure of the outer layers, making the star inflate. Thus, in order to compute accurately the mass–radius relation for low-mass helium white dwarfs we need to account for the diffusion processes during (at least) the white dwarf stages of the evolution of these objects. This should be particularly important when studying the general characteristics of binary systems containing a helium white dwarf and a pulsar.
In addition, we present an analytic, approximate model for the outer layers of the white dwarf aimed at interpreting the physical reasons for the change in the surface gravity for low-mass white dwarfs induced by diffusion.     

Cooling curves and initial models for low-mass white dwarfs (<0.25 M⊙) with helium cores

We present a detailed calculation of the evolution of low-mass (<0.25 M_⊙) helium white dwarfs. These white dwarfs (the optical companions to binary millisecond pulsars) are formed via long-term, low-mass binary evolution. After detachment from the Roche lobe, the hot helium cores have a rather thick hydrogen layer with mass between 0.01 and 0.06 M_⊙. As a result of mixing between the core and outer envelope, the surface hydrogen content ( X _surf) is 0.5–0.35 , depending on the initial value of the heavy element Z and the initial secondary mass. We found that the majority of our computed models experience one or two hydrogen shell flashes. We found that the mass of the helium dwarf in which the hydrogen shell flash occurs depends on the chemical composition. The minimum helium white dwarf mass in which a hydrogen flash takes place is 0.213 M_⊙ ( Z =0.003), 0.198 M_⊙ ( Z =0.01), 0.192 M_⊙ ( Z =0.02) or 0.183 M_⊙ ( Z =0.03). The duration of the flashes (independent of chemical composition) is between a few ×10⁶ and a few ×10⁷ yr. In several flashes the white dwarf radius will increase so much that it forces the model to fill its Roche lobe again. Our calculations show that the cooling history of the helium white dwarf depends dramatically on the thickness of the hydrogen layer. We show that the transition from a cooling white dwarf with a temporarily stable hydrogen-burning shell to a cooling white dwarf in which almost all residual hydrogen is lost in a few thermal flashes (via Roche lobe overflow) occurs between 0.183 and 0.213 M_⊙ (depending on the heavy element value).     

Scenarios for the formation of binary and millisecond pulsars – A critical assessment

The evolution of high-and low-mass X-ray binaries (HMXB and LMXB) into different types of binary radio pulsars, the ‘high-mass binary pulsars’(HMBP) and ‘low-mass binary pulsars’ (LMBP) is discussed. The HMXB evolve either into Thorne-Zytkow objects or into short-period binaries consisting of a helium star plus a neutron star (or a black hole), resembling Cygnus X-3. The latter systems evolve (with or without a second common-envelope phase) into close binary pulsars, in which the companion of the pulsar may be a massive white dwarf, a neutron star or a black hole ( some final systems may also consist of two black holes). A considerable fraction of the systems may also be disrupted in the second supernova explosion. We discuss the possible reasons why the observed numbers of double neutron stars and of systems like Cyg X-3 are several orders of magnitude lower than theoretically predicted. It is argued that the observed systems form the tip of an iceberg of much larger populations of unobserved systems, some of which may become observable in the future. As to the LMBP, we consider in some detail the origins of systems with orbital periods in the range 1–20 days. We show that to explain their existence, losses of orbital angular momentum (e.g., by magnetic braking) and in a number of cases: also of mass, have to be taken into account. The masses of the low-mass white dwarf companions in these systems can be predicted accurately. We notice a clear correlation between spin period and orbital period for these systems, as well as a clear correlation between pulsar magnetic field strength and orbital period. These relations strongly suggest that increased amounts of mass accreted by the neutron stars lead to increased decay of their magnetic fields: we suggest a simple way to understand the observed value of the ‘bottom’ field strengths of a few times 10⁸ G. Furthermore, we find that the LMBP-systems in which the pulsar has a strong magnetic field (> 10¹¹ G) have an about two orders of magnitude larger birth rate (i.e., about 4 × 10^-4 yr^-1 in the Galaxy) than the systems with millisecond pulsars (which have B < 10⁹ G). Using the observational fact that neutron stars receive a velocity kick of ∼450 km/s at birth, we find that some 90% of the potential progenitor systems of the strong-field LMBP must have been disrupted in the Supernovae in which their neutron stars were formed. Hence, the formation rate of the progenitors of the strong-field LMBP is of the same order as the galactic supernova rate (4 × 10^-3 yr^-1). This implies that a large fraction of all Supernovae take place in binaries with a close low-mass (< 2.3 M⊙) companion.     

The DODO survey – II. A Gemini direct imaging search for substellar and planetary mass companions around nearby equatorial and Northern hemisphere white dwarfs

The aim of the Degenerate Objects around Degenerate Objects (DODO) survey is to search for very low-mass brown dwarfs and extrasolar planets in wide orbits around white dwarfs via direct imaging. The direct detection of such companions would allow the spectroscopic investigation of objects with temperatures much lower  (<500 K)  than the coolest brown dwarfs currently observed. These ultra-low-mass substellar objects would have spectral types >T8.5, and so could belong to the proposed Y dwarf spectral sequence. The detection of a planet around a white dwarf would prove that such objects can survive the final stages of stellar evolution and place constraints on the frequency of planetary systems around their progenitors (with masses between 1.5 and 8   M_⊙  , i.e. early B to mid-F). This paper presents the results of a multi epoch J band common proper motion survey of 23 nearby equatorial and Northern hemisphere white dwarfs. We rule out the presence of any common proper motion companions, with limiting masses determined from the completeness limit of each observation, to 18 white dwarfs. For the remaining five targets, the motion of the white dwarf is not sufficiently separated from the non-moving background objects in each field. These targets require additional observations to conclusively rule out the presence of any common proper motion companions. From our completeness limits, we tentatively suggest that  ≲5 per cent  of white dwarfs have substellar companions with   T _eff≳ 500 K  between projected physical separations of 60–200 au.     

Discovery of 10 pulsars in an Arecibo drift-scan survey

We present the results of a 430-MHz survey for pulsars conducted during the upgrade to the 305-m Arecibo radio telescope. Our survey covered a total of 1147 deg² of sky using a drift-scan technique. We detected 33 pulsars, 10 of which were not known prior to the survey observations. The highlight of the new discoveries is PSR J0407+1607, which has a spin period of 25.7 ms, a characteristic age of 1.5 Gyr and is in a 1.8-yr orbit about a low-mass  (>0.2 M_⊙)  companion. The long orbital period and small eccentricity  ( e = 0.0009)  make the binary system an important new addition to the ensemble of binary pulsars suitable to test for violations of the strong equivalence principle. We also report on our initially unsuccessful attempts to detect optically the companion to J0407+1607, which imply that its absolute visual magnitude is >12.1. If, as expected on evolutionary grounds, the companion is an He white dwarf, our non-detection implies a cooling age of least 1 Gyr.     

The low-mass white dwarf companion to β Crateris

We present FUSE H Lyman series spectroscopy of the hot white dwarf companion to the 4th magnitude A1 III star β  Crt, which shows that it has an unusually low mass,     , and has almost certainly evolved through binary interaction. This system could be a long-sought remnant of Algol-type evolution, although radial velocity measurements appear to show that the pair are not close. Instead, micro-variations in the proper motion of β  Crt as measured by Hipparcos suggest that the period could be as high as ∼10 yr. However, a low-mass white dwarf in a system with a period ≳3 yr is difficult to explain by conventional models for binary evolution. We speculate on alternative models for the evolution of this system which involve an eccentric binary or multiple components.     

Non-adiabatic linear pulsation models for low-mass helium stars

Non-adiabatic linear pulsation models have been calculated for low-mass stars with effective temperatures between 16 000 and 35 000 K, and with surface gravities in the range 3⊙, X =0.00, Z =0.02. It is shown that the Z -bump instability persists to low masses ( M ∼0.4 M_⊙) but is suppressed either by a reduction in metallicity Z or by a selective enhancement of the carbon abundance. An unexpected result is the discovery that Z -bump instability persists at hydrogen abundances X >0.3, although the position of the red edge is sensitive to X . We have found that non-radial pulsations are also excited in the same instability region as radial pulsations.
The implications of these results for individual low-mass helium stars are discussed. It is concluded that Z -bump driven pulsations (radial and/or non-radial) may be excited in some helium-rich subdwarf B stars, representing a possible major extension to the class of variable stars represented by the prototype V652 Her.     

The enigma of the oldest 'nova': the central star and nebula of CK Vul

CK Vul is classified as, amongst others, the slowest known nova, a hibernating nova or a very late thermal pulse object. Following its eruption in ad 1670, the star remained visible for 2 yr. A 15-arcsec nebula was discovered in the 1980s, but the star itself has not been detected since the eruption. We here present radio images which reveal a 0.1-arcsec radio source with a flux of 1.5 mJy at 5 GHz. Deep Hα images show a bipolar nebula with a longest extension of 70 arcsec, with the previously known compact nebula at its waist. The emission-line ratios show that the gas is shock-ionized, at velocities  >100 km s⁻¹  . Dust emission yields an envelope mass of  ∼5 × 10⁻² M_⊙  . Echelle spectra indicate outflow velocities up to 360 km s⁻¹. From a comparison of images obtained in 1991 and 2004 we find evidence for expansion of the nebula, consistent with an origin in the 1670 explosion; the measured expansion is centred on the radio source. No optical or infrared counterpart is found at the position of the radio source. The radio emission is interpreted as thermal free–free emission from gas with   T _e∼ 10⁴ K  . The radio source may be due to a remnant circumbinary disc, similar to those seen in some binary post-AGB stars. We discuss possible classifications of this unique outburst, including that of a sub-Chandrasekhar mass supernova, a nova eruption on a cool, low-mass white dwarf or a thermal pulse induced by accretion from a circumbinary disc.     

The ages and colours of cool helium-core white dwarf stars

The purpose of this work is to explore the evolution of helium-core white dwarf stars in a self-consistent way with the predictions of detailed non-grey model atmospheres and element diffusion. To this end, we consider helium-core white dwarf models with stellar masses of 0.406, 0.360, 0.327, 0.292, 0.242, 0.196 and 0.169 M_⊙ and follow their evolution from the end of mass-loss episodes, during their pre-white dwarf evolution, down to very low surface luminosities.
We find that when the effective temperature decreases below 4000 K, the emergent spectrum of these stars becomes bluer within time-scales of astrophysical interest. In particular, we analyse the evolution of our models in the colour–colour and in the colour–magnitude diagrams and find that helium-core white dwarfs with masses ranging from ∼0.18 to 0.3 M_⊙ can reach the turn-off in their colours and become blue again within cooling times much less than 15 Gyr and then remain brighter than M _V ≈16.5 . In view of these results, many low-mass helium white dwarfs could have had enough time to evolve to the domain of collision-induced absorption from molecular hydrogen, showing blue colours.     

Models for the formation of binary and millisecond radio pulsars

The peculiar combination of a relatively short pulse period and a relatively weak surface dipole magnetic field strength of binary radio pulsars finds a consistent explanation in terms of (i) decay of the surface dipole component of neutron-star magnetic fields on a timescale of (2–5) × 10⁶ yr, in combination with (ii) spin-up of the rotation of the neutron star during a subsequent mass-transfer phase. The four known binary radio pulsars appear to fall into two different categories. Two of them, PSR 0655 + 64 and PSR 1913 + 16, have short orbital periods (<25 h) and high mass functions, indicating companion masses 0.7M⊙ (∼1 (± 0.3) M⊙ and 1.4 M⊙, respectively). The other two, PSR 0820 + 02 and PSR 1953 + 29, have long orbital periods (117d), nearly circular orbits, and low, almost identical mass functions of about 3×10^-3 M⊙, suggesting companion masses of about 0.3M⊙. It is pointed out that these two classes of systems are expected to be formed by the later evolution of binaries consisting of a neutron star and a normal companion star, in which the companion was (considerably) more massive than the neutron star, or less massive than the neutron star, respectively. In the first case the companion of the neutron star in the final system will be a massive white dwarf, in a circular orbit, or a neutron star in an eccentric orbit. In the second case the final companion to the neutron star will be a low-mass (∼ 0.3 M⊙) helium white dwarf in a wide and nearly circular orbit. In systems of the second type the neutron star was most probably formed by the accretion-induced collapse of a white dwarf. This explains in a natural way why PSR 1953 + 29 has a millisecond rotation period and PSR 0820 + 02 has not. Among the binary models proposed for the formation of the 1.5-millisecond pulsar, the only ones that appear to be viable are those in which the companion disappeared by coalescence with the neutron star. In such models the companion may have been a red dwarf of mass 0.03M⊙, a neutron star, or a massive (>0.7M⊙) white dwarf. Only in the last-mentioned case is a position of the pulsar close to the galactic plane a natural consequence. In the first-mentioned case the progenitor system most probably was a cataclysmic-variable binary in which the white dwarf collapsed by accretion.     

Metallicity measurements using atomic lines in M and K dwarf stars

We report the first survey of chemical abundances in M and K dwarf stars using atomic absorption lines in high-resolution spectra. We have measured Fe and Ti abundances in 35 M and K dwarf stars using equivalent widths measured from  λ/Δλ≈ 33 000  spectra. Our analysis takes advantage of recent improvements in model atmospheres of low-temperature dwarf stars. The stars have temperatures between 3300 and 4700 K, with most cooler than 4100 K. They cover an iron abundance range of  −2.44 < [Fe/H] < +0.16  . Our measurements show [Ti/Fe] decreasing with increasing [Fe/H], a trend similar to that measured for warmer stars, where abundance analysis techniques have been tested more thoroughly. This study is a step towards the observational calibration of procedures to estimate the metallicity of low-mass dwarf stars using photometric and low-resolution spectral indices.     

Modelling the formation of galaxy clusters in MOND

We present time-resolved spectroscopy and photometry of the cataclysmic variable (CV) SDSS J133941.11+484727.5 (SDSS 1339) which has been discovered in the Sloan Digital Sky Survey (SDSS) Data Release 4. The orbital period determined from radial velocity studies is 82.524(24) min, close to the observed period minimum. The optical spectrum of SDSS 1339 is dominated to 90 per cent by emission from the white dwarf (WD). The spectrum can be successfully reproduced by a three-component model (white dwarf, disc, secondary) with   T _WD=12 500 K  for a fixed  log   g = 8.0, d = 170 pc  , and a spectral type of the secondary later than M8. The mass-transfer rate corresponding to the optical luminosity of the accretion disc is very low,  ≃ 1.7 × 10⁻¹³ M_⊙ yr⁻¹  . Optical photometry reveals a coherent variability at 641 s with an amplitude of 0.025 mag, which we interpret as non-radial pulsations of the white dwarf. In addition, a long-period photometric variation with a period of either 320 or 344 min and an amplitude of 0.025 mag is detected, which bears no apparent relation with the orbital period of the system. Similar long-period photometric signals have been found in the CVs SDSS J123813.73−033933.0, SDSS J204817.85−061044.8, GW Lib and FS Aur, but so far no working model for this behaviour is available.