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.     

Numerical simulations of protostellar encounters — I. Star–disc encounters

It appears that most stars are born in clusters, and that at birth most stars have circumstellar discs which are comparable in size to the separations between the stars. Interactions between neighbouring stars and discs are therefore likely to play a key role in determining disc lifetimes, stellar masses, and the separations and eccentricities of binary orbits. Such interactions may also cause fragmentation of the discs, thereby triggering the formation of additional stars.   We have carried out a series of simulations of star–disc interactions using an SPH code which treats self-gravity, hydrodynamic and viscous forces. We find that interactions between discs and stars provide a mechanism for removing energy from, or adding energy to, the orbits of the stars, and for truncating the discs. However, capture during such encounters is unlikely to be an important binary formation mechanism.   A more significant consequence of such encounters is that they can trigger fragmentation of the disc, via tidally and compressionally induced gravitational instabilities, leading to the formation of additional stars and substellar objects. When the disc spins and stellar orbits are randomly oriented, encounters lead to the formation of new companions to the original star in 20 per cent of encounters. If most encounters are prograde and coplanar, as suggested by simulations of dynamically triggered star formation, then new companions are formed in approximately 50 per cent of encounters.     

Pulsar spin–velocity alignment from single and binary neutron star progenitors

The role of binary progenitors of neutron stars (NSs) in the apparent distribution of space velocities and spin–velocity alignment observed in young pulsars is studied. We performed a Monte Carlo synthesis of pulsar populations originated from single and binary stars with different assumptions about the NS natal kick (kick–spin alignment, kick amplitude and kick reduction in electron-capture supernovae in binary progenitors with initial main-sequence masses from the range  8–11 M_⊙  which experienced mass exchange due to Roche lobe overflow). The calculated spin–velocity alignment in pulsars is compared with data inferred from radio polarization measurements. The observed space velocity of pulsars is found to be mostly affected by the natal kick velocity form and its amplitude; the fraction of binaries is not important here for reasonably large kicks. The natal kick–spin alignment is found to strongly affect the spin–velocity correlation of pulsars. Comparison with the observed pulsar spin–velocity angles favours a sizeable fraction of binary progenitors and kick–spin angles  ∼5°–20°  .     

The dynamical stability of W Ursae Majoris-type systems

Theoretical study indicates that a contact binary system would merge into a rapidly rotating single star due to tidal instability when the spin angular momentum of the system is more than a third of its orbital angular momentum. Assuming that W Ursae Majoris (W UMa) contact binary systems rigorously comply with the Roche geometry and the dynamical stability limit is at a contact degree of about 70 per cent, we obtain that W UMa systems might suffer Darwin's instability when their mass ratios are in a region of about 0.076–0.078 and merge into the fast-rotating stars. This suggests that the W UMa systems with mass ratio   q ≤ 0.076  cannot be observed. Meanwhile, we find that the observed W UMa systems with a mass ratio of about 0.077, corresponding to a contact degree of about 86 per cent would suffer tidal instability and merge into the single fast-rotating stars. This suggests that the dynamical stability limit for the observed W UMa systems is higher than the theoretical value, implying that the observed systems have probably suffered the loss of angular momentum due to gravitational wave radiation (GR) or magnetic stellar wind (MSW).     

Binary stars and the fundamental initial mass function

We have carried out Monte Carlo simulations in which we generate a random pairing of objects drawn from a pre-assumed single-star power-law initial mass function (IMF), which we call the fundamental IMF. We show how the mass functions of primary stars and secondary stars and the mass function of the total mass of systems (if we could resolve them) differ from the underlying fundamental IMF for different slopes of this IMF. We also compare our results with the observed IMF, the binary frequency and the binary mass-ratio distributions for field stars and conclude that the fundamental IMF of subsolar mass stars could be steeper than is currently believed. In other words, the low-mass turn-over of the observed ('apparent') IMF could be spurious, if the main-sequence binary fraction of field stars is close to 100 per cent (perhaps owing to invisible companions).     

On the binarity of Herbig Ae/Be stars

We present high-resolution spectro-astrometry of a sample of 28 Herbig Ae/Be and three F-type pre-main-sequence stars. The spectro-astrometry, which is essentially the study of unresolved features in long-slit spectra, is shown from both empirical and simulated data to be capable of detecting binary companions that are fainter by up to 6 mag at separations larger than ∼0.1 arcsec. The nine targets that were previously known to be binary are all detected. In addition, we report the discovery of six new binaries and present five further possible binaries. The resulting binary fraction is 68 ± 11 per cent. This overall binary fraction is the largest reported for any observed sample of Herbig Ae/Be stars, presumably because of the exquisite sensitivity of spectro-astrometry for detecting binary systems. The data hint that the binary frequency of the Herbig Be stars is larger than that of the Herbig Ae stars. The Appendix presents model simulations to assess the capabilities of spectro-astrometry and reinforces the empirical findings. Most spectro-astrometric signatures in this sample of Herbig Ae/Be stars can be explained by the presence of a binary system. Two objects, HD 87643 and Z CMa, display evidence for asymmetric outflows. Finally, the position angles of the binary systems have been compared with available orientations of the circumprimary disc and these appear to be coplanar. The alignment between the circumprimary discs and the binary systems strongly suggests that the formation of binaries with intermediate-mass primaries is due to fragmentation as the alternative, stellar capture, does not naturally predict aligned discs. The alignment extends to the most massive B-type stars in our sample. This leads us to conclude that formation mechanisms that do result in massive stars, but predict random angles between the binaries and the circumprimary discs, such as stellar collisions, are also ruled out for the same reason.     

Be phenomenon in open clusters: results from a survey of emission-line stars in young open clusters

Emission-line stars in young open clusters are identified to study their properties, as a function of age, spectral type and evolutionary state. 207 open star clusters were observed using the slitless spectroscopy method and 157 emission stars were identified in 42 clusters. We have found 54 new emission-line stars in 24 open clusters, out of which 19 clusters are found to house emission stars for the first time. About 20 per cent clusters harbour emission stars. The fraction of clusters housing emission stars is maximum in both the 0–10 and 20–30 Myr age bin (∼40 per cent each). Most of the emission stars in our survey belong to Classical Be class (∼92 per cent) while a few are Herbig Be stars (∼6 per cent) and Herbig Ae stars (∼2 per cent). The youngest clusters to have Classical Be stars are IC 1590, NGC 637 and 1624 (all 4 Myr old) while NGC 6756 (125–150 Myr) is the oldest cluster to have Classical Be stars. The Classical Be stars are located all along the main sequence (MS) in the optical colour–magnitude diagrams (CMDs) of clusters of all ages, which indicates that the Be phenomenon is unlikely due to core contraction near the turn-off. The distribution of Classical Be stars as a function of spectral type shows peaks at B1–B2 and B6–B7 spectral types. The Be star fraction [N(Be)/N(B+Be)] is found to be less than 10 per cent for most of the clusters and NGC 2345 is found to have the largest fraction (∼26 per cent). Our results indicate there could be two mechanisms responsible for the Classical Be phenomenon. Some are born Classical Be stars (fast rotators), as indicated by their presence in clusters younger than 10 Myr. Some stars evolve to Classical Be stars, within the MS lifetime, as indicated by the enhancement in the fraction of clusters with Classical Be stars in the 20–30 Myr age bin.     

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.     

Stellar-mass black hole binaries as ultraluminous X-ray sources

Ultraluminous X-ray sources (ULXs) with   L _x > 10³⁹ erg s⁻¹  have been discovered in great numbers in external galaxies with ROSAT , Chandra and XMM-Newton . The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g. systems containing intermediate-mass BHs  (100–1000 M_⊙)  . We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass BHs. To this end, we have applied a unique set of binary evolution models for BH X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses  ≳10 M_⊙  the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, although the latter phases last only ∼5 per cent of the main-sequence phase. We show that with only a modest violation of the Eddington limit, e.g. a factor of ∼10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass BH binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is  ∼3 × 10⁻⁷ yr⁻¹  normalized to a core-collapse supernova rate of 0.01 yr⁻¹.     

The strongest cosmic magnets: soft gamma-ray repeaters and anomalous X-ray pulsars

Two classes of X-ray pulsars, the anomalous X-ray pulsars and the soft gamma-ray repeaters, have been recognized in the last decade as the most promising candidates for being magnetars: isolated neutron stars powered by magnetic energy. I review the observational properties of these objects, focussing on the most recent results, and their interpretation in the magnetar model. Alternative explanations, in particular those based on accretion from residual disks, are also considered. The possible relations between these sources and other classes of neutron stars and astrophysical objects are also discussed.     

Magnetic and spin evolution of neutron stars in close binaries

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

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.     

On the infant weight loss of low- to intermediate-mass star clusters

Star clusters are born in a highly compact configuration, typically with radii of less than about 1 pc roughly independently of mass. Since the star formation efficiency is less than 50 per cent by observation and because the residual gas is removed from the embedded cluster, the cluster must expand. In the process of doing so it only retains a fraction f _st of its stars. To date there are no observational constraints for f _st, although N -body calculations by Kroupa, Aarseth & Hurley suggest it to be about 20–30 per cent for Orion-type clusters. Here we use the data compiled by Testi et al., Testi, Palla & Natta and Testi, Palla & Natta for clusters around young Ae/Be stars and by de Wit et al. and de Wit et al. around young O stars and the study of de Zeeuw et al. of OB associations and combine these measurements with the expected number of stars in clusters with primary Ae/Be and O stars, respectively, using the empirical correlation between maximal stellar mass and star cluster mass of Weidner & Kroupa. We find that   f _st < 50  per cent with a decrease to higher cluster masses/more massive primaries. The interpretation would be that cluster formation is very disruptive. It appears that clusters with a birth stellar mass in the range  10–10³ M_⊙  keep at most 50 per cent of their stars.     

Automated classification of variable stars for All-Sky Automated Survey 1–2 data

With the advent of surveys generating multi-epoch photometry and the discovery of large numbers of variable stars, the classification of these stars has to be automatic. We have developed such a classification procedure for about 1700 stars from the variable star catalogue of the All-Sky Automated Survey 1–2 (ASAS 1–2) by selecting the periodic stars and by applying an unsupervised Bayesian classifier using parameters obtained through a Fourier decomposition of the light curve. For irregular light curves we used the period and moments of the magnitude distribution for the classification. In the case of ASAS 1–2, 83 per cent of variable objects are red giants. A general relation between the period and amplitude is found for a large fraction of those stars. The selection led to 302 periodic and 1429 semiperiodic stars, which are classified in six major groups: eclipsing binaries, 'sinusoidal curves', Cepheids, small amplitude red variables, SR and Mira stars. The type classification error level is estimated to be about 7 per cent.     

Aquila X-1: a low-inclination soft X-ray transient

We have obtained I -band photometry of the neutron star X-ray transient Aql X-1 during quiescence. We find a periodicity at 2.487 cycles d⁻¹, which we interpret as twice the orbital frequency (19.30±0.05 h). Folding the data on the orbital period, we model the light-curve variations as the ellipsoidal modulation of the secondary star. We determine the binary inclination to be 20°–30° (90 per cent confidence) and also determine the 95 per cent upper limits to the radial velocity semi-amplitude and rotational broadening of the secondary star to be 117 and 50 km s⁻¹, respectively.     

Optical properties of Small Magellanic Cloud X-ray binaries

This work represent the first major study of the optical and infrared characteristics of the mass donor companions to the X-ray pulsars in the Small Magellanic Cloud (SMC). In this work several new counterparts have been identified, and possible ones confirmed, as companions to X-ray pulsars in the SMC giving a total of 34 such objects now identified. In addition this work presents three new binary periods and confirms two X-ray periods using optical data for objects in this group. This homogeneous sample has been studied as a group to determine important general characteristics that may offer an insight into the evolution of such systems. In particular, the spectral class distribution shows a much greater agreement with those of isolated Be stars, and appears to be in some disagreement with the galactic population of Be stars in Be/X-ray binaries. Studies of the long-term optical modulation of the Be star companions reveal an extremely variable group of objects, a fact which will almost certainly make a major contribution to the pronounced X-ray variability. The spatial distribution of these systems within the SMC is investigated and strongly suggests a link between massive star formation and the H  i density distribution. Finally, studies of the circumstellar disc characteristics reveal a strong link with optical variability offering important clues into the long-term stability of such discs.     

The first magnetic maps of a pre-main-sequence binary star system – HD 155555

We present the first maps of the surface magnetic fields of a pre-main-sequence binary system. Spectropolarimetric observations of the young, 18 Myr, HD 155555 (V824 Ara, G5IV+K0IV) system were obtained at the Anglo-Australian Telescope in 2004 and 2007. Both data sets are analysed using a new binary Zeeman–Doppler imaging (ZDI) code. This allows us to simultaneously model the contribution of each component to the observed circularly polarized spectra. Stellar brightness maps are also produced for HD 155555 and compared to previous Doppler images.
Our radial magnetic maps reveal a complex surface magnetic topology with mixed polarities at all latitudes. We find rings of azimuthal field on both stars, most of which are found to be non-axisymmetric with the stellar rotational axis. We also examine the field strength and the relative fraction of magnetic energy stored in the radial and azimuthal field components at both epochs. A marked weakening of the field strength of the secondary star is observed between the 2004 and 2007 epochs. This is accompanied by an apparent shift in the location of magnetic energy from the azimuthal to radial field. We suggest that this could be indicative of a magnetic activity cycle. We use the radial magnetic maps to extrapolate the coronal field (by assuming a potential field) for each star individually – at present ignoring any possible interaction. The secondary star is found to exhibit an extreme tilt (≈75°) of its large-scale magnetic field to that of its rotation axis for both epochs. The field complexity that is apparent in the surface maps persists out to a significant fraction of the binary separation. Any interaction between the fields of the two stars is therefore likely to be complex also. Modelling this would require a full binary field extrapolation.     

Gamma-ray bursts from stellar remnants: probing the Universe at high redshift

A gamma-ray burst (GRB) releases an amount of energy similar to that of a supernova explosion, which combined with its rapid variability suggests an origin related to neutron stars or black holes. Since these compact stellar remnants form from the most massive stars not long after their birth, GRBs should trace the star formation rate in the Universe; we show that the GRB flux distribution is consistent with this. Because of the strong evolution of the star formation rate with redshift, it follows that the dimmest known bursts have z  ∼ 6, much above the value usually quoted and beyond the most distant quasars. This explains the absence of bright galaxies in well-studied GRB error boxes. The increased distances imply a peak luminosity of 8.3 × 10⁵¹ erg s⁻¹ and a rate density of 0.025 per million years per galaxy. These values are 20 times higher and 150 times lower, respectively, than are implied by fits with non-evolving GRB rates. This means either that GRBs are caused by a much rarer phenomenon than mergers of binary neutron stars, or that their gamma-ray emission is often invisible to us due to beaming. Precise burst locations from optical transients will discriminate between the various models for GRBs from stellar deaths, because the distance between progenitor birth place and burst varies greatly among them. The dimmest GRBs are then the most distant known objects, and may probe the Universe at an age when the first stars were forming.     

Magnetic fields and rotation in white dwarfs and neutron stars

Recent spectropolarimetric observations of Ap and Bp stars with improved sensitivity have suggested that most Ap and Bp stars are magnetic with dipolar fields of at least a few hundred gauss. These new estimates suggest that the range of magnetic fluxes found for the majority of magnetic white dwarfs is similar to that of main-sequence Ap–Bp stars, thus strengthening the empirical evidence for an evolutionary link between magnetism on the main sequence and magnetism in white dwarfs. We draw parallels between the magnetic white dwarfs and the magnetic neutron stars and argue that the observed range of magnetic fields in isolated neutron stars  ( B_p ∼ 10¹¹–10¹⁵ G)  could also be explained if their mainly O-type progenitors have effective dipolar fields in the range of a few gauss to a few kilogauss, assuming approximate magnetic flux conservation with the upper limit being consistent with the recent measurement of a field of   B_p ∼ 1100 G  for θ Orion C.
In the magnetic field–rotation diagram, the magnetic white dwarfs can be divided into three groups of different origin: a significant group of strongly magnetized slow rotators  ( P _rot∼ 50 –100 yr)  that have originated from single-star evolution, a group of strongly magnetized fast rotators  ( P _rot∼ 700 s)  , typified by EUVE J0317–853, that have originated from a merger, and a group of modest rotators ( P _rot∼ hours–days) of mixed origin (single-star and CV-type binary evolution). We propose that the neutron stars may similarly divide into distinct classes at birth , and suggest that the magnetars may be the counterparts of the slowly rotating high-field magnetic white dwarfs.     

Repeating microlensing events in the OGLE data

Microlensing events are usually selected among single-peaked non-repeating light curves in order to avoid confusion with variable stars. However, a microlensing event may exhibit a second microlensing brightening episode when the source or/and the lens is a binary system. A careful analysis of these repeating events provides an independent way to study the statistics of wide binary stars and to detect extrasolar planets. Previous theoretical studies predicted that 0.5–2 per cent of events should repeat due to wide binary lenses. We present a systematic search for such events in about 4000 light curves of microlensing candidates detected by the Optical Gravitational Lensing Experiment (OGLE) towards the Galactic bulge from 1992 to 2007. The search reveals a total of 19 repeating candidates, with six clearly due to a wide binary lens. As a by-product, we find that 64 events (∼2 per cent of the total OGLE-III sample) have been misclassified as microlensing; these misclassified events are mostly nova or other types of eruptive stars. The number and importance of repeating events will increase considerably when the next-generation wide-field microlensing experiments become fully operational in the future.