New Praesepe white dwarfs and the initial mass–final mass relation

We report the spectroscopic confirmation of four further white dwarf members of Praesepe. This brings the total number of confirmed white dwarf members to 11, making this the second largest collection of these objects in an open cluster identified to date. This number is consistent with the high-mass end of the initial mass function of Praesepe being Salpeter in form. Furthermore, it suggests that the bulk of Praesepe white dwarfs did not gain a substantial recoil kick velocity from possible asymmetries in their loss of mass during the asymptotic giant branch phase of evolution. By comparing our estimates of the effective temperatures and the surface gravities of WD0833+194, WD0840+190, WD0840+205 and WD0843+184 to modern theoretical evolutionary tracks, we have derived their masses to be in the range  0.72–0.76 M_⊙  and their cooling ages ∼300 Myr. For an assumed cluster age of 625 ± 50 Myr, the inferred progenitor masses are between 3.3 and  3.5 M_⊙  . Examining these new data in the context of the initial mass–final mass relation, we find that it can be adequately represented by a linear function  ( a ₀= 0.289 ± 0.051,  a ₁= 0.133 ± 0.015)  over the initial mass range 2.7–6  M_⊙  . Assuming an extrapolation of this relation to larger initial masses is valid and adopting a maximum white dwarf mass of  1.3 M_⊙  , our results support a minimum mass for core-collapse supernovae progenitors in the range  ∼6.8–8.6 M_⊙  .     

Praesepe – two merging clusters?

A membership catalogue for Praesepe was compiled and split into four mass bins. A contour plot indicates the presence of a subcluster some 3 pc from the centre of the cluster, of approximately 30 M_⊙. A tidally truncated King profile was fitted to the remainder of the cluster and the tidal radius is found to be 12.1 pc; the mass of the cluster (excluding the subcluster) is 630 M_⊙. From the calculated velocity dispersions we find that the cluster appears to have too much kinetic energy and should be rapidly disintegrating. X-ray data suggest that there may be an age spread between the main core stars and the subcluster stars. This leads us to the conclusion that Praesepe is two merging clusters.     

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.     

The lithium depletion boundary in NGC 2547 as a test of pre-main-sequence evolutionary models

Intermediate resolution spectroscopy from the European Southern Observatory Very Large Telescope is analysed for 63 photometrically selected low-mass  (0.08–0.30 M_⊙)  candidates of the open cluster NGC 2547. We have confirmed membership for most of these stars using radial velocities, and found that lithium remains undepleted for cluster stars with   I > 17.54 ± 0.14  and   K_s > 14.86 ± 0.12  . From these results, several pre–main-sequence evolutionary models give almost model independent ages of 34–36 Myr, with a precision of 10 per cent. These ages are only slightly larger than the ages of 25–35(±5) Myr obtained using the same models to fit isochrones to higher mass stars descending towards the zero-age main-sequence, both in empirically calibrated and theoretical colour–magnitude diagrams. This agreement between age determinations in different mass ranges is an excellent test of the current generation of low-mass pre–main-sequence stellar models and lends confidence to ages determined with either method between 30 and 120 Myr.     

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.     

Predictions of the extent of self-enrichment in oxygen of giant metal-poor H ii regions

In general, H  ii regions do not show clear signs of self-enrichment in products from massive stars  ( M ≥ 8 M_⊙)  . In order to explore why I modelled the contamination with Wolf–Rayet star ejecta of metal-poor  ( Z = 0.001)  H  ii regions, ionized either by a  10⁶ M_⊙  cluster of coeval stars (cluster 1) or by a cluster resulting from continuous star formation at a rate of  1 M_⊙ yr⁻¹  (cluster 2). The clusters have   Z = 0.001  and a Salpeter initial mass function from 0.1 to  120 M_⊙  . Independent one-dimensional constant density simulations of the emission-line spectra of unenriched H  ii regions were computed at the discrete ages 1, 2, 3, 4 and 5 Myr, with the photoionization code cloudy , using as input, radiative and mechanical stellar feedbacks predicted by the evolutionary synthesis code starburst99 . Each H  ii region was placed at the outer radius of the adiabatically expanding superbubble of Mac Low & McCray. For models with thermal and ionization balance time-scales of less than 1 Myr, and with oxygen emission-line ratios in agreement with observations, the volume of the superbubble and the H  ii region was uniformly and instantaneously polluted with stellar ejecta predicted by starburst99 . I obtained a maximum oxygen abundance enhancement of 0.025 dex, with cluster 1, at 4 Myr. It would be unobservable.     

Do accretion discs regulate the rotation of young stars?

We present a photometric study of I -band variability in the young cluster IC 348. The main purpose of the study was to identify periodic stars. In all, we find 50 periodic stars, of which 32 were previously unknown. For the first time in IC 348, we discover periods in significant numbers of lower-mass stars  ( M < 0.25 M_⊙)  and classical T Tauri stars. This increased sensitivity to periodicities is a result of the enhanced depth and temporal density of our observations, compared with previous studies. The period distribution is at first glance similar to that seen in the Orion nebula cluster (ONC), with the higher-mass stars  ( M > 0.25 M_⊙)  showing a bi-modal period distribution concentrated around periods of 2 and 8 d, and the lower-mass stars showing a uni-modal distribution, heavily biased towards fast rotators. Closer inspection of the period distribution shows that the higher-mass stars show a significant dearth of fast rotators, compared to the ONC, whilst the low-mass stars are rotating significantly faster than those in Orion. We find no correlation between rotation period and K – L colour or Hα equivalent width.
We also present a discussion of our own IC 348 data in the context of previously published period distributions for the ONC, the Orion flanking fields and NGC 2264. We find that the previously claimed correlation between infrared excess and rotation period in the ONC might, in fact, result from a correlation between infrared excess and mass. We also find a marked difference in period distributions between NGC 2264 and IC 348, which presents a serious challenge to the disc-locking paradigm, given the similarity in ages and disc fractions between the two clusters.     

The low-mass initial mass function in the young cluster NGC 6611

NGC 6611 is the massive young cluster (2–3 Myr) that ionizes the Eagle Nebula. We present very deep photometric observations of the central region of NGC 6611 obtained with the Hubble Space Telescope and the following filters: ACS/WFC F775W and F850LP and NIC2 F110W and F160W, loosely equivalent to ground-based IZJH filters. This survey reaches down to   I ∼ 26 mag  . We construct the initial mass function (IMF) from  ∼1.5 M_⊙  well into the brown dwarf regime (down to  ∼0.02 M_⊙  ). We have detected 30–35 brown dwarf candidates in this sample. The low-mass IMF is combined with a higher-mass IMF constructed from the ground-based catalogue from Oliveira et al. We compare the final IMF with those of well-studied star-forming regions: we find that the IMF of NGC 6611 more closely resembles that of the low-mass star-forming region in Taurus than that of the more massive Orion Nebula Cluster. We conclude that there seems to be no severe environmental effect in the IMF due to the proximity of the massive stars in NGC 6611.     

Rotation of the cluster of galaxies A2107

We present indications of rotation in the galaxy cluster A2107 by a method that searches for the maximum gradient in the velocity field in a flat disc-like model of a cluster. Galaxies from cumulative subsamples containing more and more distant members from the cluster centre are projected on to an axis passing through the centre and we apply a linear regression model on the projected distances x and the line-of-sight velocities V . The axis with the maximum linear correlation coefficient   r _max= max[ r ( V , x )]  defines the direction of the maximum velocity gradient, and consequently it presents the major axis of the apparently elliptical cluster. Because the effects of rotation are subtle, we put strong emphasis on the estimation of the uncertainties of the results by implementing different bootstrap techniques. We have found that rotational effects are more strongly expressed at distances of 0.26–0.54 Mpc from the cluster centre. The total virial mass of the cluster is  (3.2 ± 0.6) × 10¹⁴ M_⊙  , while the virial mass, corrected for rotation, is  (2.8 ± 0.5) × 10¹⁴ M_⊙  .     

Early-A stars from IPHAS, and their distribution in and around the Cyg OB2 association

Stellar photometry derived from the INT/WFC Photometric Hα Survey (IPHAS) of the Northern Galactic plane can be used to identify large, reliable samples of A0–A5 stars. For every A-type star, so identified, it is also possible to derive individual reddening and distance estimates, under the assumption that most selected objects are on or near the main sequence, at a mean absolute r ' magnitude of 1.5–1.6. This study presents the method for obtaining such samples and shows that the known reddenings and distances to the open clusters NGC 7510 and NGC 7790 are successfully recovered. A sample of over 1000 A-type stars is then obtained from IPHAS data in the magnitude range  13.5 < r ' < 20  from the region of sky including the massive northern OB association Cyg OB2. An analysis of these data reveals a concentration of ∼200 A stars over an area about a degree across, offset mainly to the south of the known 1–3 Myr old OB stars in Cyg OB2: their dereddened r ' magnitudes fall in the range 11.8–12.5. These are consistent with a ∼7 Myr old stellar population at distance modulus DM = 10.8, or with an age of ∼5 Myr at DM = 11.2. The number of A-type stars found in this clustering alone is consistent with a lower limit to the cluster mass of  ∼10⁴ M_⊙  .     

Using the minimum spanning tree to trace mass segregation

We present a new method to detect and quantify mass segregation in star clusters. It compares the minimum spanning tree (MST) of massive stars with that of random stars. If mass segregation is present, the MST length of the most massive stars will be shorter than that of random stars. This difference can be quantified (with an associated significance) to measure the degree of mass segregation. We test the method on simulated clusters in both 2D and 3D and show that the method works as expected.
We apply the method to the Orion Nebula Cluster (ONC) and show that the method is able to detect the mass segregation in the Trapezium with a 'mass segregation ratio (MSR)'  Λ_MSR= 8.0 ± 3.5  (where  Λ_MSR= 1  is no mass segregation) down to  16 M_⊙  , and also that the ONC is mass segregated at a lower level  (∼2.0 ± 0.5)  down to  5 M_⊙  . Below  5 M_⊙  we find no evidence for any further mass segregation in the ONC.     

On the first generation of stars

We argue that the first stars may have spanned the conventional mass range rather than be identified with the very massive objects  (∼100–10³ M_⊙)  favoured by numerical simulations. Specifically, we find that magnetic field generation processes acting in the first protostellar systems suffice to produce fields that exceed the threshold for magneto-rotational instability (MRI) to operate, and thereby allow the MRI dynamo to generate equipartition-amplitude magnetic fields on protostellar mass scales below  ∼50 M_⊙  . Such fields allow primordial star formation to occur at essentially any metallicity by regulating angular momentum transfer, fragmentation, accretion and feedback in much the same way as occurs in conventional molecular clouds.     

On the origin of high-velocity runaway stars

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary  50–100 M_⊙  star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of  ≳200–400 km s⁻¹  (typical of pulsars), while  3–4 M_⊙  stars can attain velocities of  ≳300–400 km s⁻¹  (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.     

Orbital periods of the binary sdB stars PG 0940+068 and PG 1247+554

We have used the radial velocity variations of two sdB stars previously reported to be binaries to establish their orbital periods. They are PG 0940+068 ( P =8.33 d) and PG 1247+554 ( P =0.599 d). The minimum masses of the unseen companions, assuming a mass of 0.5 M_⊙ for the sdB stars, are 0.090±0.003 M. for PG 1247+554 and 0.63±0.02 M_⊙ for PG 0940+068. The nature of the companions is not constrained further by our data.     

Long Secondary Periods in variable red giants

We present a study of a sample of Large Magellanic Cloud red giants exhibiting Long Secondary Periods (LSPs). We use radial velocities obtained from VLT spectral observations and MACHO and OGLE light curves to examine properties of the stars and to evaluate models for the cause of LSPs. This sample is much larger than the combined previous studies of Hinkle et al. and Wood, Olivier & Kawaler.
Binary and pulsation models have enjoyed much support in recent years. Assuming stellar pulsation, we calculate from the velocity curves that the typical fractional radius change over an LSP cycle is greater than 30 per cent. This should lead to large changes in T _eff that are not observed. Also, the small light amplitude of these stars seems inconsistent with the radius amplitude. We conclude that pulsation is not a likely explanation for the LSPs. The main alternative, physical movement of the star – binary motion – also has severe problems. If the velocity variations are due to binary motion, the distribution of the angle of periastron in our large sample of stars has a probability of  1.4 × 10⁻³  that it comes from randomly aligned binary orbits. In addition, we calculate a typical companion mass of  0.09 M_⊙  . Less than 1 per cent of low-mass main-sequence stars have companions near this mass  (0.06–0.12 M_⊙)  whereas ∼25–50 per cent of low-mass red giants end up with LSPs. We are unable to find a suitable model for the LSPs and conclude by listing their known properties.     

The abundance of brown dwarfs

The amount of mass contained in low-mass objects is investigated anew. Instead of using a mass–luminosity relation to convert a luminosity function to a mass function, I predict the mass–luminosity relation from assumed mass functions and the luminosity functions of Jahreiss & Wielen and Gould, Bahcall & Flynn. Comparison of the resulting mass–luminosity relations with data for binary stars constrains the permissible mass functions. If the mass function is assumed to be a power law, the best-fitting slope lies either side of the critical slope, α =−2, below which the mass in low-mass objects is divergent, depending on the luminosity function adopted. If these power-law mass functions are truncated at 0.001 M_⊙, the contribution to the local density from stars lies between 0.013 and 0.10 M_⊙ pc⁻³ depending on the mass at which the mass function is normalized and the adopted value of α . Recent dynamical estimates of the local mass density rule out stellar mass densities above ∼0.05 M_⊙ pc⁻³. Hence, power laws steeper than α =−2 that extend down to 0.001 M_⊙ are allowed only if one adopts an implausible normalization of the mass function. If the mass function is generalized from a power law to a low-order polynomial in log( M ), the mass in stars with M <0.1 M_⊙ is either negligible or strongly divergent, depending on the order of the polynomial adopted.     

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.     

Modelling the binary progenitor of Supernova 1993J

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code, we investigate potential progenitors for the Type IIb Supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al.) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient, the observations can be reproduced by a system consisting of a  15 M_⊙  primary and a  14 M_⊙  secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a  17 M_⊙  primary and a  16 M_⊙  secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing.     

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.     

Cosmological origin of the lowest metallicity halo stars

We explore the predictions of the standard hierarchical clustering scenario of galaxy formation, regarding the numbers and metallicities of PopIII stars that are likely to be found within our Galaxy today. By PopIII we refer to stars formed at large redshift ( z >4), with low metallicities ([ Z /Z_⊙]<−2.5) and in small systems (total mass ≲ 2×10⁸ M_⊙) that are extremely sensitive to stellar feedback, and which through a prescribed merging history end up becoming part of the Milky Way today. An analytic, extended Press–Schechter formalism is used to obtain the mass functions of haloes which will host PopIII stars at a given redshift, and which will end up in Milky Way sized systems today. Each of these is modelled as a mini-galaxy, with a detailed treatment of the dark halo structure, angular momentum distribution, final gas temperature and disc instabilities, all of which determine the fraction of the baryons that are subject to star formation. The use of new primordial metallicity stellar evolutionary models allows us to trace the history of the stars formed, and give accurate estimates of their expected numbers today and their location in L /L_⊙ versus T /K Hertzsprung–Russell (HR) diagrams. A first comparison with observational data suggests that the initial mass function (IMF) of the first stars was increasingly high-mass weighted towards high redshifts, levelling off at z ≳9 at a characteristic stellar mass scale m _s=10–15 M_⊙.