The photometric evolution of FU Orionis objects: disc instability and wind–envelope interaction

We present the results of a photometric monitoring campaign of three well-studied FU Orionis systems (FU Orionis, V1057 Cygni and V1515 Cygni) undertaken at Maidanak Observatory between 1981 and 2003. When combined with photometric data in the literature, this data base provides a valuable resource for searching for short time-scale variability – both periodic and aperiodic – as well as for studying the secular evolution of these systems. In the case of V1057 Cyg (which is the system exhibiting the largest changes in brightness since it went into outburst) we compare the photometric data with time-dependent models. We show that prior to the end of the 'plateau' stage in 1996, the evolution of V1057 Cyg in the V –( B − V ) colour–magnitude diagram is well represented by disc instability models in which the outburst is triggered by some agent – such as an orbiting planet – in the inner disc. Following the end of the plateau phase in 1996, the dimming and irregular variations are consistent with occultation of the source by a variable dust screen, which has previously been interpreted in terms of dust condensation events in the observed disc wind. Here we instead suggest that this effect results from the interaction between the wind and an infalling dusty envelope, the existence of this envelope having been previously invoked in order to explain the mid-infrared emission of FU Orionis systems. We discuss how this model may explain some of the photometric and spectroscopic characteristics of FU Orionis systems in general.     

Inclusion of binaries in evolutionary population synthesis

Using evolutionary population synthesis we present integrated colours, integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as the Lick/IDS) system, for an extensive set of instantaneous-burst binary stellar populations with and without binary interactions. The ages of the populations are in the range 1–15 Gyr and the metallicities are in the range 0.0001–0.03. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the integrated U – B , B – V , V – R and R – I colours and all Lick/IDS spectral absorption indices (except for H_β) substantially smaller. In other words, binary evolution makes a population appear bluer. This effect raises the derived age and metallicity of the population.
We calculate several sets of additional solar-metallicity binary stellar populations to explore the influence of the binary evolution algorithm input parameters (the common-envelope ejection efficiency and the stellar wind mass-loss rate) on the resulting integrated colours. We also look at the dependence on the choice of distribution functions used to generate the initial binary population. The results show that variations in the choice of input model parameters and distributions can significantly affect the results. However, comparing the discrepancies that exist between the colours of various models, we find that the differences are less than those produced between the models with and those without binary interactions. Therefore it is very necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.     

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.     

Resolving the pulsations of the subdwarf B star KPD 2109 + 4401

We present the results of extensive time-series photometry of the pulsating subdwarf B star KPD 2109 + 4401. Our data set consists of 29 data runs with a total length of 182.6 h over 31 d, collected at five observatories in 2004. These data comprised high signal-to-noise ratio observations acquired with larger telescopes and wider time-coverage observations obtained with smaller telescopes. They are sufficient to resolve the pulsation structure to 0.4 μHz and are the most extensive data set for this star to date. With these data, we identify eight pulsation frequencies extending from 4701 to 5481 μHz, corresponding to periods of 182–213 s. The pulsation frequencies and their amplitudes are examined over several time-scales with some frequencies showing amplitude variability.     

The identification of the optical/IR counterpart to the ROSAT SMC X-ray transient RX J0117.6−7330

We report observations which identify the optical/IR counterpart to the ROSAT X-ray transient RX J0117.6−7330. The counterpart is suggested to be a B1–B2 star (luminosity class III–V) showing an IR excess and strong Balmer emission lines. The distance derived from reddening and systemic velocity measurements is consistent with the distance derived from X-ray measurements and puts the source in the Small Magellanic Cloud (SMC).     

WIYN Open Cluster Study – XVI. Optical/infrared photometry and comparisons with theoretical isochrones

We present combined optical/near-infrared photometry ( BVIK ) for six open clusters – M35, M37, NGC 1817, NGC 2477, NGC 2420 and M67. The open clusters span an age range from 150 Myr to 4 Gyr and have metal abundances from  [Fe/H]=−0.27  to +0.09 dex. We have utilized these data to test the robustness of theoretical main sequences constructed by several groups as denoted by the following designations – Padova, Baraffe, Y², Geneva and Siess. The comparisons of the models with the observations have been performed in the  [ M_V , ( B − V )₀], [ M_V , ( V − I )₀]  and  [ M_V , ( V − K )₀]  colour–magnitude diagrams as well as the distance-independent  [( V − K )₀, ( B − V )₀]  and  [( V − K )₀, ( V − I )₀]  two-colour diagrams. We conclude that none of the theoretical models reproduces the observational data in a consistent manner over the magnitude and colour range of the unevolved main sequence. In particular, there are significant zero-point and shape differences between the models and the observations. We speculate that the crux of the problem lies in the precise mismatch between theoretical and observational colour–temperature relations. These results underscore the importance of pursuing the study of stellar structure and stellar modelling with even greater intensity.     

The influence of binary interactions on infrared passbands of populations

In our evolutionary population synthesis models, samples of binaries are reproduced by a ' patched ' Monte Carlo simulation and the stellar masses, integrated   J ,  H ,  K ,  L ,  L 2  and M magnitudes, mass-to-light ratios and broad colours involving infrared bands are presented for an extensive set of instantaneous-burst binary stellar populations. In addition, the fluctuations in the integrated colours, which have been given by Zhang et al., are reduced.
By comparing the results for binary stellar populations with (Model A) and without (Model B) binary interactions, we show that the inclusion of binary interactions makes the stellar mass of a binary stellar population smaller (  ∼3.6–4.5  per cent during the past 15 Gyr), magnitudes greater (except   U , ∼ 0.18 mag  at the most), colours bluer (∼0.15 mag for   V − K   at the most) and mass-to-light ratios greater (∼0.06 for K band) except those in the U and B passbands at higher metallicities. Binary interactions make the V magnitude less sensitive to age, and R and I magnitudes more sensitive to metallicity.
Given an age, the absolute values of the differences in the stellar mass, magnitudes and mass-to-light ratios (except those in the U and B bands) between Models A and B reach a maximum at   Z = 0.0001  , i.e. the effects of binary interactions on these parameters reach a maximum, while the differences in some colours reach a maximum at   Z ∼ 0.01–0.0004  . In contrast, the absolute value of the difference in the stellar mass is minimal at   Z = 0.03  ; those in the   U ,  B ,  V   magnitudes and the mass-to-light ratios in the U and B bands reach a minimum at   Z ∼ 0.01–0.004  .     

Resumption of mass accretion in RS Oph

The latest outburst of the recurrent nova RS Oph occurred in 2006 February. Photometric data presented here show evidence of the resumption of optical flickering, indicating re-establishment of accretion by day 241 of the outburst. Magnitude variations of up to 0.32 mag in V band and 0.14 mag in B band on time-scales of 600–7000 s are detected. Over the two-week observational period, we also detect a 0.5 mag decline in the mean brightness, from   V ≈ 11.4  to 11.9, and record   B ≈ 12.9 mag  . Limits on the mass accretion rate of     are calculated, which span the range of accretion rates modelled for direct wind accretion and Roche lobe overflow mechanisms. The current accretion rates make it difficult for thermonuclear runaway models to explain the observed recurrence interval, and this implies average accretion rates are typically higher than seen immediately post-outburst.     

Do the physical properties of Ap binaries depend on their orbital elements?

We reveal sufficient evidence that the physical characteristics of Ap stars are related to binarity. The Ap star peculiarity [represented by the  Δ( V 1- G )  value and magnetic field strength] diminishes with eccentricity, and it may also increase with orbital period ( P _orb). This pattern, however, does not hold for large orbital periods. A striking gap that occurs in the orbital period distribution of Ap binaries at 160–600 d might well mark a discontinuity in the above-mentioned behaviour. There is also an interesting indication that the Ap star eccentricities are relatively lower than those of corresponding B9–A2 normal binaries for   P _orb>10 d  . All this gives serious support to the pioneering idea of Abt & Snowden concerning a possible interplay between the magnetism of Ap stars and their binarity. Nevertheless, we argue instead in favour of another mechanism, namely that it is binarity that affects magnetism and not the opposite, and suggest the presence of a new magnetohydrodynamical mechanism induced by the stellar companion and stretching to surprisingly large P _orb.     

The outer-halo globular cluster NGC 6229 – IV. Variable stars

We present time-series B , V photometry of NGC 6229, obtained with the purpose of providing the first extensive CCD variability study of this cluster. As a result, we were able to obtain periods for all NGC 6229 variables, with the exception of five stars from the candidate list of Borissova et al. located very close to the cluster centre. Two stars suspected to be variables by Carney, Fullton and Trammell are first-overtone RR Lyrae (RRc) stars, whereas seven of the 12 candidates of Borissova et al. are confirmed variables – three of them being fundamental RR Lyrae (RRab) pulsators, two first-overtone pulsators, one an eclipsing binary (most likely an Algol system) and one a bright star whose variability status could not be satisfactorily determined. A new image subtraction method (ISM) suggested by Alard has been employed which, together with the Welch–Stetson technique, has allowed us to discover twelve new RR Lyrae variables in the cluster, for which ephemerides are provided. Ten of these are RRabs, whereas the other two are RRcs. As originally suggested by Mayer, NGC 6229 is clearly an Oosterhoff type I globular cluster. We also confirm that V8 is a Population II Cepheid of the W Virginis type, and suspect that this is the case for V22 as well. The physical properties of the NGC 6229 RR Lyrae population are contrasted with those of M3 (NGC 5272) using several different methods, including a standard period-shift analysis. Possible differences between these two clusters are discussed.     

Multiwavelength study of the transient X-ray binary IGR J01583+6713

We have investigated multiband optical photometric variability and stability of the Hα line profile of the transient X-ray binary IGR J01583+6713. We set an upper limit of 0.05 mag on photometric variations in the V band over a time-scale of three months. The Hα line is found to consist of non-Gaussian profile and quite stable for a duration of two months. We have identified the spectral type of the companion star to be B2 IVe while the distance to the source is estimated to be ∼4.0 kpc. Along with the optical observations, we have also carried out analysis of X-ray data from three short observations of the source, two with the Swift –XRT and one with the RXTE –PCA. We have detected a variation in the absorption column density, from a value of  22.0 × 10²² cm⁻²  immediately after the outburst down to  2.6 × 10²² cm⁻²  four months afterwards. In the quiescent state, the X-ray absorption is consistent with the optical reddening measurement of   E ( B − V ) = 1.46  mag. From one of the Swift observations, during which the X-ray intensity was higher, we have a possible pulse detection with a period of 469.2 s. For a Be X-ray binary, this indicates an orbital period in the range of 216–561 d for this binary system.     

Variable stars in the field of open cluster NGC 6819

We report on the discovery of 25 variable stars plus 13 suspected variables found in the field of the open cluster NGC 6819. The stars were identified from time-series photometric data obtained on the Isaac Newton Telescope, La Palma, during two observing runs covering the 19 nights between 1999 June  22–30  and 1999 July  22–31  . The variables found include 12 eclipsing binaries with an additional three suspected, nine BY Draconis systems, plus four variables of other types, including one star believed to be a Cepheid. Three of the 15 eclipsing binaries are believed to be cluster members. Details of a further 10 suspected variable stars are also included.     

Two new pulsating hot subdwarf stars from the Edinburgh-Cape survey

We report the discovery of very rapid pulsations in two hot subdwarf stars from the Edinburgh-Cape blue object survey. The short periods, small amplitudes and multiperiodicity establish these stars as members of the class of rapidly-pulsating sdB stars. The spectrograms of both stars, however, show relatively strong He  ii 4686 and they are therefore more properly classified as sdOB. The light curve of EC 01541−1409 is dominated by two strong (∼1 per cent) variations with frequencies near 7114 and 7870 μHz (periods near 140.6 and 127.1 s), though at least five frequencies are present with amplitudes above about 0.002 mag. The light curve of EC 22221−3152 appears to be generated by at least 10 frequencies in the range 5670–11850 μHz (about 175–85 s) with amplitudes between about 0.01 and 0.001 mag, including the first overtone of the strongest variation. Somewhat surprisingly, this number of frequencies is detectable in observing runs as short as 3 h, probably due to the fact that the detected frequencies are well-separated.     

WR 7a: a V Sagittae or a qWR star?

The star WR 7a, also known as SPH 2, has a spectrum that resembles that of V Sagittae stars although no O  vi emission has been reported. The Temporal Variance Spectrum – TVS – analysis of our data shows weak but strongly variable emission of O  vi lines which is below the noise level in the intensity spectrum.
Contrary to what is seen in V Sagittae stars, optical photometric monitoring shows very little, if any, flickering. We found evidence of periodic variability. The most likely photometric period is   P _phot= 0.227(±14) d  , while radial velocities suggest a period of   P _spec= 0.204(±13) d  . One-day aliases of these periods can not be ruled out. We call attention to similarities with HD 45166 and DI Cru (= WR 46), where multiple periods are present. They may be associated to the binary motion or to non-radial oscillations.
In contrast to a previous conclusion by Pereira et al., we show that WR 7a contains hydrogen. The spectrum of the primary star seems to be detectable as the N  v 4604 Å  absorption line is visible. If so, it means that the wind is optically thin in the continuum and that it is likely to be a helium main sequence star.
Given the similarity to HD 45166, we suggests that WR 7a may be a qWR – quasi Wolf–Rayet – star. Its classification is WN4h/CE in the Smith, Shara & Moffat three-dimensional classification system.     

V473 Lyrae: a Blazhko-effect Cepheid?

The amplitude spectra of the light curves of RR Lyrae Blazhko stars can often be typified by symmetrical triplets. It is shown that the same applies to V473 Lyr. Furthermore, the star is probably a radial pulsator – as is thought to be the case for the RR Lyrae stars – thus ruling out rotational splitting of a purely non-radial mode.     

Atmospheric analysis of the M/L and M/T dwarf binary systems LHS 102 and Gliese 229

We present  0.9–2.5 μm  spectroscopy with   R ∼800  and  1.12–1.22 μm  spectroscopy with   R ∼5800  for the M dwarfs Gl 229A and LHS 102A, and for the L dwarf LHS 102B. We also report IZJHKL ' photometry for both components of the LHS 102 system, and L ' photometry for Gl 229A. The data are combined with previously published spectroscopy and photometry to produce flux distributions for each component of the kinematically old disc M/L dwarf binary system LHS 102 and the kinematically young disc M/T dwarf binary system Gliese 229. The data are analysed using synthetic spectra generated by the latest 'AMES-dusty' and 'AMES-cond' models by Allard & Hauschildt. Although the models are not able to reproduce the overall slope of the infrared flux distribution of the L dwarf, most likely because of the treatment of dust in the photosphere, the data for the M dwarfs and the T dwarf are well matched. We find that the Gl 229 system is metal-poor despite having kinematics of the young disc, and that the LHS 102 system has solar metallicity. The observed luminosities and derived temperatures and gravities are consistent with evolutionary model predictions if the Gl 229 system is very young  (age∼30 Myr)  with masses (A,B) of (0.38,≳0.007) M_⊙, and the LHS 102 system is older, aged  1–10 Gyr  with masses (A,B) of (0.19,0.07) M_⊙.     

uvbyβ photometry of the short-period binary VV Ursae Majoris

We present well-sampled uvby light curves, supplemented by a few β filter measurements, of the Algol binary VV UMa. The light curves are analysed using two different codes to derive the orbital and absolute stellar parameters of this binary. We find reasonably good fits to the light curves and determine the stellar effective temperatures T _eff,1≃9000–9600 K , and T _eff,2≃5300–5600 K with a mass ratio q ≃0.35 . From the light-curve fits we discard the possibility of an anomalous gravity-darkening exponent for the secondary star of this system, as previously suggested.
We find evidence of short-term, small-amplitude variations in the brightness of the system. Two periodicities of about 1.10 and 0.51 h seem to be present in the data for at least two different nights, even within the secondary eclipse. This suggests that VV UMa may be a new Algol binary with a low-amplitude variable primary star, but new data collected during longer observing runs are necessary to confirm the pulsating nature of the brightness variations.     

The Na i D resonance lines in main-sequence late-type stars

We study the sodium D lines (D1: 5895.92 Å; D2: 5889.95 Å) in late-type dwarf stars. The stars have spectral types between F6 and M5.5 ( B − V between 0.457 and 1.807) and metallicity between  [Fe/H]=−0.82  and 0.6. We obtained medium-resolution echelle spectra using the 2.15-m telescope at the Argentinian observatory Complejo Astronómico El Leoncito (CASLEO). The observations have been performed periodically since 1999. The spectra were calibrated in wavelength and in flux. A definition of the pseudo-continuum level is found for all our observations. We also define a continuum level for calibration purposes. The equivalent width of the D lines is computed in detail for all our spectra and related to the colour index ( B − V ) of the stars. When possible, we perform a careful comparison with previous studies. Finally, we construct a spectral index  ( R '_D)  as the ratio between the flux in the D lines and the bolometric flux. We find that, once corrected for the photospheric contribution, this index can be used as a chromospheric activity indicator in stars with a high level of activity. Additionally, we find that combining some of our results, we obtain a method to calibrate in flux stars of unknown colour.     

The 'off' state of GX 339–4

We report BeppoSAX and optical observations of the black hole candidate GX 339–4 during its X-ray 'off' state in 1999. The broad-band (0.8–50 keV) X-ray emission can be fitted by a single power law with spectral index, α ∼1.6. The observed luminosity is 6.6×10³³ erg s⁻¹ in the 0.5–10 keV band, which is at the higher end of the flux distribution of black hole soft X-ray transients in quiescence, comparable to that seen in GS 2023+338 and 4U 1630–47. An optical observation just before the BeppoSAX observation shows the source to be very faint at these wavelengths as well ( B =20.1, V =19.2). By comparing with previously reported 'off' and low states (LS), we conclude that the 'off' state is actually an extension of the LS, i.e. an LS at lower intensities. We propose that accretion models such as the advection-dominated accretion flows are able to explain the observed properties in such a state.