An optical spectroscopic HR diagram for low-mass stars and brown dwarfs in Orion

The masses and temperatures of young low-mass stars and brown dwarfs in star-forming regions are not yet well established because of uncertainties in the age of individual objects and the spectral type–temperature scale appropriate for objects with ages of only a few Myr. Using multi-object optical spectroscopy, 45 low-mass stars and brown dwarfs in the Trapezium Cluster in Orion have been classified and 44 of these confirmed as bona fide cluster members. The spectral types obtained have been converted to effective temperatures using a temperature scale intermediate between those of dwarfs and giants, which is suitable for young pre-main-sequence objects. The objects have been placed on a Hertzsprung–Russell (HR) diagram overlaid with theoretical isochrones. The low-mass stars and the higher mass substellar objects are found to be clustered around the 1 Myr isochrone, while many of the lower mass substellar objects are located well above this isochrone. An average age of 1 Myr is found for the majority of the objects. Assuming coevality of the sources and an average age of 1 Myr, the masses of the objects have been estimated and range from  0.018 to 0.44 M_⊙  . The spectra also allow an investigation of the surface gravity of the objects by measurement of the sodium doublet equivalent width. With one possible exception, all objects have low gravities, in line with young ages, and the Na indices for the Trapezium objects lie systematically below those of young stars and brown dwarfs in Chamaeleon, suggesting that the 820 nm Na index may provide a sensitive means of estimating ages in young clusters.     

High-resolution spectra of very low-mass stars

We present the results of high-resolution (1–0.4 Å) optical spectroscopy of a sample of very low-mass stars. These data are used to examine the kinematics of the stars at the bottom of the hydrogen-burning main sequence. No evidence is found for a significant difference between the kinematics of the stars in our sample with I  −  K  > 3.5 ( M _bol ≳ 12.8) and those of more massive M dwarfs ( M _bol ≈ 7–10). A spectral atlas at high (0.4-Å) resolution for M8–M9+ stars is provided, and the equivalent widths of Cs  I , Rb  I and Hα lines present in our spectra are examined. We analyse our data to search for the presence of rapid rotation, and find that the brown dwarf LP 944-20 is a member of the class of 'inactive, rapid rotators'. Such objects seem to be common at and below the hydrogen-burning main sequence. It seems that in low-mass/low-temperature dwarf objects either the mechanism that heats the chromosphere, or the mechanism that generates magnetic fields, is greatly suppressed.     

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.     

Radial and rotational velocities of young brown dwarfs and very low-mass stars in the Upper Scorpius OB association and the ρ Ophiuchi cloud core

We present the results of a radial velocity (RV) survey of 14 brown dwarfs (BDs) and very low-mass (VLM) stars in the Upper Scorpius OB association (UScoOB) and three BD candidates in the ρ Ophiuchi dark cloud core. We obtained high-resolution echelle spectra at the Very Large Telescope using the Ultraviolet and Visual Echelle Spectrograph (UVES) at two different epochs for each object, and measured the shifts in their RVs to identify candidates for binary/multiple systems in the sample. The average time separation of the RV measurements is 21.6 d, and our survey is sensitive to the binaries with separation <0.1 au. We found that four out of 17 objects (or 24⁺¹⁶₋₁₃ per cent by fraction) show a significant RV change in 4–33 d time-scale, and are considered as binary/multiple 'candidates'. We found no double-lined spectroscopic binaries in our sample, based on the shape of cross-correlation curves. The RV dispersion of the objects in UScoOB is found to be very similar to that of the BD and VLM stars in Chamaeleon I (Cha I). We also found the distribution of the mean rotational velocities ( v   sin   i ) of the UScoOB objects is similar to that of the Cha I, but the dispersion of v   sin   i is much larger than that of the Cha I objects.     

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.     

Optical spectroscopic classification and membership of young M dwarfs in star-forming regions

The spectral type is a key parameter in calibrating the temperature which is required to estimate the mass of young stars and brown dwarfs. We describe an approach developed to classify low-mass stars and brown dwarfs in the Trapezium Cluster using red optical spectra, which can be applied to other star-forming regions. The classification uses two methods for greater accuracy: the use of narrow-band spectral indices which rely on the variation of the strength of molecular lines with spectral type and a comparison with other previously classified young, low-mass objects in the Chamaeleon I star-forming region. We have investigated and compared many different molecular indices and have identified a small number of indices which work well for classifying M-type objects in nebular regions. The indices are calibrated for young, pre-main-sequence objects whose spectra are affected by their lower surface gravities compared with those on the main sequence. Spectral types obtained are essentially independent of both reddening and nebular emission lines.
Confirmation of candidate young stars and brown dwarfs as bona fide cluster members may be accomplished with moderate resolution spectra in the optical region by an analysis of the strength of the gravity-sensitive Na doublet. It has been established that this feature is much weaker in these very young objects than in field dwarfs. A sodium spectral index is used to estimate the surface gravity and to demonstrate quantitatively the difference between young (1–2 Myr) objects, and dwarf and giant field stars.     

Spectral analysis of water vapour in cool stars

M-star spectra, at wavelengths beyond 1.35 μm, are dominated by water vapour, yet terrestrial water vapour makes it notoriously difficult to obtain accurate measurement from ground-based observations. We have used the short-wavelength spectrometer on the Infrared Space Observatory at four wavelength settings to cover the  2.5–3.0 μm  region for a range of M stars. The observations show a good match with previous ground-based observations and with synthetic spectra based on the Partridge & Schwenke line list, although not with the SCAN line list. We have used a least-squared minimization technique to systematically find best-fitting parameters for the sample of stars. The temperatures that we find indicate a relatively hot temperature scale for M dwarfs. We consider that this could be a consequence of problems with the Partridge & Schwenke line list which leads to synthetic spectra predicting water bands that are too strong for a given temperature. Such problems need to be solved in the next generation of water vapour line lists, which will extend the calculation of water vapour to higher energy levels with the good convergence necessary for reliable modelling of hot water vapour. Then water bands can assume their natural role as the primary tool for the spectroscopic analysis of M stars.     

The distance to the Orion Nebula cluster

The distance to the Orion Nebula cluster (ONC) is estimated using the rotational properties of its low-mass pre-main-sequence (PMS) stars. Rotation periods, projected equatorial velocities and distance-dependent radius estimates are used to form an observational sin  i distribution (where i is the axial inclination), which is modelled to obtain the distance estimate. A distance of  440 ± 34 pc  is found from a sample of 74 PMS stars with spectral types between G6 and M2, but this falls to  392 ± 32 pc  when PMS stars with accretion discs are excluded on the basis of their near-infrared excess. Since the radii of accreting stars are more uncertain and probably systematically underestimated, then this closer distance is preferred. The quoted uncertainties include statistical errors and uncertainties due to a number of systematic effects including binarity and inclination bias. This method is geometric and independent of stellar evolution models, though does rely on the assumption of random axial orientations and the Cohen & Kuhi effective temperature scale for PMS stars. The new distance is consistent with, although lower and more precise, than most previous ONC distance estimates. A closer ONC distance implies smaller luminosities and an increased age based on the positions of PMS stars in the Hertzsprung–Russell diagram.     

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.     

The discovery of a low-mass, pre-main-sequence stellar association around γ Velorum

We report the serendipitous discovery of a population of low-mass, pre-main-sequence (PMS) stars in the direction of the Wolf–Rayet/O-star binary system γ ²  Vel and the Vela OB2 association. We argue that γ ²  Vel and the low-mass stars are truly associated and approximately coeval, and that both are at distances between 360 and 490 pc, disagreeing at the 2 σ level with the recent Hipparcos parallax of γ ²  Vel, but consistent with older distance estimates. Our results clearly have implications for the physical parameters of the γ ²  Vel system, but also offer an exciting opportunity to investigate the influence of high-mass stars on the mass function and circumstellar disc lifetimes of their lower mass PMS siblings.     

Metallicity and photospheric abundances in field K and M dwarfs

Elemental abundances in late-type stars are of interest in several ways: they determine the location of the stars in the HR diagram and therefore their ages, as well as the atmospheric structure in their middle and upper photospheres. Especially in the case of chromospherically active late-type stars the question arises to what degree the upper photosphere is influenced by the nearby chromosphere. Analysing S/N ∼ 200 and Δλ/λ ∼ 20 000 data, we found a mean metallicity index [M/H] = −0.2 for programme K and M field stars based on an analysis of spectra in the region 5500–9000 Å. We also found that the Ca  I 6162-Å transition is a potential surface gravity indicator for K-type stars. For the chromospheric activity interval 4.4 < log  F _{Mg II}  < 6.6 we did not find any chromospheric activity impact on photospheric and upper photospheric transitions. With the derived metallicity, we confirmed the Li abundance from our previous paper and thus its dependence on the Mg  II chromospheric activity index. The nature of the spectrum for the active M-type star Gl 896A is explained by pure rotation of 14 km s⁻¹. As far as the lithium–rotation relation is concerned, the spectrum of Gl 517 is rotationally broadened as well, by 12 km s⁻¹, and the Li abundance is the second highest in our sample of stars. However, there is no link between very high Li abundance, 2.2 dex, in the K dwarf star Gl 5 and stellar rotation.     

Observations and modelling of spectral energy distributions of carbon stars with optically thin envelopes

We present broad-band photometry in the optical, near-infrared and submillimetre, and mid-infrared spectrophotometry of a selection of carbon stars with optically thin envelopes. Most of the observations were carried out simultaneously.   Beside the emission feature at 11.3 μ m due to silicon carbide grains in the circumstellar environment, many of our mid-infrared spectra show an emission feature at 8.6 μ m. All the observed spectral energy distributions exhibit a very large far-infrared flux excess. Both these features are indeed common to many carbon stars surrounded by optically thin envelopes.   We have modelled the observed spectral energy distributions by means of a full radiative transfer treatment, paying particular attention to the features quoted above. The peak at 8.6 μ m is usually ascribed to the presence of hydrogenated amorphous carbon grains. We find also that the feature at 8.6 μ m might be reproduced by assuming that the stars have a circumstellar environment formed of both carbon- and oxygen-rich dust grains, although this is in contrast with what one should expect in a carbon-rich environment. The far-infrared flux excess is usually explained by the presence of a cool detached dust shell. Following this hypothesis, our models suggest a time-scale for the modulation of the mass-loss rate of the order of some 10³ yr.     

Rotation, planets, and the 'second parameter' of the horizontal branch

We analyse the angular momentum evolution from the red giant branch (RGB) to the horizontal branch (HB) and along the HB. Using rotation velocities for stars in the globular cluster M13, we find that the required angular momentum for the fast rotators is up to 1–3 orders of magnitude (depending on some assumptions) larger than that of the Sun. Planets of masses up to 5 times Jupiter's mass and up to an initial orbital separation of ~2 au are sufficient to spin-up the RGB progenitors of most of these fast rotators. Other stars have been spun-up by brown dwarfs or low-mass main-sequence stars. Our results show that the fast rotating HB stars have been probably spun-up by planets, brown dwarfs or low-mass main-sequence stars while they evolved on the RGB. We argue that the angular momentum considerations presented in this paper further support the 'planet second parameter' model. In this model, the 'second parameter' process, which determines the distribution of stars on the HB, is interaction with low-mass companions, in most cases with gas-giant planets, and in a minority of cases with brown dwarfs or low-mass main-sequence stars. The masses and initial orbital separations of the planets (or brown dwarfs or low-mass main-sequence stars) form a rich spectrum of different physical parameters, which manifests itself in the rich varieties of HB morphologies observed in the different globular clusters.     

Brown dwarfs and the cataclysmic variable period minimum

Using improved, up-to-date stellar input physics tested against observations of low-mass stars and brown dwarfs, we calculate the secular evolution of low-mass donor cataclysmic variables (CVs), including those that form with a brown-dwarf donor. Our models confirm the mismatch between the calculated minimum period ( P _min70 min) and the observed short-period cut-off (80 min) in the CV period histogram. We find that tidal and rotational corrections applied to the one-dimensional stellar structure equations have no significant effect on the period minimum. Theoretical period distributions synthesized from our model sequences always show an accumulation of systems at the minimum period, a feature absent from the observed distribution. We suggest that non-magnetic CVs become unobservable as they are effectively trapped in permanent quiescence before they reach P _min, and that small-number statistics may hide the period spike for magnetic CVs.     

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.     

100-yr mass-loss modulations on the asymptotic giant branch

We analyse the differences in infrared circumstellar dust emission between oxygen-rich Mira and non-Mira stars, and find that they are statistically significant. In particular, we find that these stars segregate in the K–[12] versus [12]–[25] colour–colour diagram, and have distinct properties of the IRAS LRS spectra, including the peak position of the silicate emission feature. We show that the infrared emission from the majority of non-Mira stars cannot be explained within the context of standard steady-state outflow models.
The models can be altered to fit the data for non-Mira stars by postulating non-standard optical properties for silicate grains, or by assuming that the dust temperature at the inner envelope radius is significantly lower (300–400 K) than typical silicate grain condensation temperatures (800–1000 K) . We argue that the latter is more probable and provide detailed model fits to the IRAS LRS spectra for 342 stars. These fits imply that two-thirds of non-Mira stars and one-third of Mira stars do not have hot dust (>500 K) in their envelopes.
The absence of hot dust can be interpreted as a recent (∼100 yr) decrease in the mass-loss rate. The distribution of best-fitting model parameters agrees with this interpretation and strongly suggests that the mass loss resumes on similar time-scales. Such a possibility appears to be supported by a number of spatially resolved observations (e.g. recent Hubble Space Telescope images of the multiple shells in the Egg Nebula) and is consistent with new dynamical models for mass loss on the asymptotic giant branch.     

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 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.     

High-resolution spectroscopy of ROSAT-discovered weak-line T Tauri stars near Lupus

We present high-resolution optical echelle spectroscopy for a large fraction of the Li-rich late-type stars recently discovered in the vicinity of the Lupus dark clouds. Our results confirm the high Li  i   λ 6708 equivalent widths previously estimated from medium-resolution spectra, thus adding strength to the conclusion that the large majority of these stars are still in the pre-main-sequence phase of their evolution, contrary to claims from other authors that many of them might be zero-age main-sequence stars. We present a statistical approach to derive a mean distance for the sample, and find that it is consistent with, or slightly lower than, the Hipparcos distance of the Lupus star-forming region. The radial velocities measured for part of these stars are consistent with those observed for the Lupus star-forming region, while stars outside the dark clouds show a mean difference of the order of 3 km s⁻¹. The projected rotational velocities show a lack of slow rotators, which is interpreted as a consequence of the X-ray selection of the sample. The Li-rich stars in Lupus studied in this work yield a fairly 'clean' sample of very young stars, while in other star-forming regions a larger fraction of older zero-age main-sequence stars has been found among ROSAT -discovered Li-rich stars. We argue that this fact reflects the relation of these stars with the Gould Belt.