Chemical abundances of magnetic and non-magnetic Herbig Ae/Be stars
The photospheres of about 10-20% of main sequence A- and B-type stars exhibit a wide range of chemical peculiarities, often associated with the presence of a magnetic field. It is not exactly known at which stage of stellar evolution these chemical peculiarities develop. To investigate this issue, in this paper we study the photospheric compositions of a sample of Herbig Ae and Be stars, which are considered to be the pre-main sequence progenitors of A and B stars. We have performed a detailed abundance analysis of 20 Herbig stars (three of which have confirmed magnetic fields), and one dusty young star. We have found that half the stars in our sample show Boo chemical peculiarities to varying degrees, only one star shows weak Ap/Bp peculiarities, and all the remaining stars are chemically normal. The incidence of Boo chemical peculiarities we find in Herbig stars is much higher than what is seen on the main sequence. We argue that a selective accretion model for Boo star formation is a natural explanation for the remarkably large number of Boo stars in our sample. We also find that the magnetic Herbig stars do not exhibit a range of chemical compositions remarkably different from the normal stars: one magnetic star displays Boo chemical peculiarities (HD 101412), one displays weak Ap/Bp peculiarities (V380 Ori A), and one (HD 190073) is chemically normal. This is completely different from what is seen on the main sequence, where all magnetic A and cool B stars show Ap/Bp chemical peculiarities, and this is consistent with the idea that the magnetic field precedes the formation of the chemical peculiarities typical of Ap and Bp stars.
keywords:stars: magnetic fields, stars: abundances, stars: chemically peculiar, stars: pre-main-sequence
In main sequence A and B stars a wide range of peculiar chemical abundances are observed. While much study has been devoted to these peculiar objects, a lot remains to be learnt about the physical processes behind these chemical peculiarities. An avenue of study that has not yet been properly investigated is to examine the chemistry in the pre-main sequence progenitors of main sequence peculiar stars.
Herbig Ae and Be (HAeBe) stars are pre-main sequence A and B stars, and thus evolve into the wide range of different main sequence A and B stars. Observationally, Herbig stars are generally identified by optical emission lines, infrared excess, and are usually associated with nebulous regions (Vieira et al., 2003). Chemically, these stars are thought to possess approximately solar abundances (Acke & Waelkens, 2004), similar to most young, nearby main sequence stars.
A small number of HAeBe stars have recently been found to have strong, globally ordered magnetic fields (Donati et al., 1997; Wade et al., 2005, 2007; Catala et al., 2007; Alecian et al., 2008a, b). The strength and morphology of these magnetic fields are very similar to those seen in the magnetic chemically peculiar Ap and Bp stars (Alecian et al., 2008a; Folsom et al., 2008; Alecian et al., 2009), making these magnetic HAeBe stars strong candidates for the progenitors of Ap and Bp stars (Wade et al., 2005). This discovery in particular raises the question of whether any sign of chemical peculiarity can be found on the pre-main sequence.
In main sequence A and cooler B stars, magnetic fields are always seen together with the characteristic chemical peculiarities of Ap and Bp stars (Aurière et al., 2007). Hotter magnetic B stars also usually show chemical peculiarities, such as He-weak or He-strong stars. Consequently, if magnetic HAeBe stars evolve into Ap and Bp stars, then at some point they must develop chemical peculiarities. Detecting such peculiarities on the pre-main sequence would provide new constraints on the timescales on which, and the conditions under which, these peculiarities evolve.
A few detections of chemical peculiarities in very young stars near or on the zero-age main sequence provide tantalising hints that chemical peculiarities may be common in magnetic Herbig Ae/Be stars. Some notable cases of this are V380 Ori (Alecian et al., 2009), HD 72106 A (Folsom et al., 2008), NGC 6611 W601 (Alecian et al., 2008b), and perhaps NGC 2244-334 (Bagnulo et al., 2004) which is a very young main sequence star. Recent modelling by Vick et al. (2011) of chemical diffusion in the presence of modest mass loss suggest it is possible for chemical peculiarities to form during the pre-main sequence phase. Interestingly, Cowley et al. (2010) found a different form of chemical peculiarities, Boötis peculiarities, in the Herbig Ae star HD 101412.
Boötis stars are mostly main sequence A stars, with strong underabundances of many metals, particularly iron peak elements. Lighter elements, specifically C, N, O, S, and in some cases Na, have normal abundances, while iron peak elements are usually depleted by 1 dex (e.g. Venn & Lambert, 1990; Andrievsky et al., 2002; Heiter, 2002). This is in contrast to almost all other chemically peculiar A and B stars, which are characterised by strong overabundances of iron peak elements. Boo stars have not been found to have magnetic fields, unlike Ap and Bp stars (Bohlender & Landstreet, 1990). Boo stars have a distribution of values that is the same as for normal A and B stars (e.g. Abt & Morrell, 1995), unlike Ap or Am stars.
The cause of the peculiar abundances seen in Boo stars remains unknown, though a number of hypotheses have been suggested. The most popular hypothesis is that the peculiarities are a result of a selective accretion process (Venn & Lambert, 1990; Waters et al., 1992). In this scenario, gas depleted in heavier elements is accreted preferentially, building up a layer of relative underabundance in heavier elements at the surface of the star. A proposed mechanism for this suggests that heavier elements are bound into dust grains, which are then preferentially driven away from the star by its radiation, while gas, which is depleted in heavier elements as a consequence of dust formation, is more readily accreted (e.g. Andrievsky & Paunzen, 2000). Thus, because of the lack of significant convective mixing in the atmospheres of A-type stars, a layer that is underabundant in heavier elements is quickly built up at the surface of the star.
A variation on this hypothesis is that the Boo peculiarities result from selective accretion of gas in the interstellar medium, rather than pre-existing circumstellar material (Kamp & Paunzen, 2002). In this scenario the star passes through a diffuse interstellar cloud, in which heavier elements are already preferentially bound into dust grains. The star then accretes the metal poor gas, while driving away the dust by radiation pressure.
Other hypotheses have been put forward to explore the origins of Boo stars. For example it has been suggested that Boo stars are actually spectroscopic binaries (Faraggiana & Bonifacio, 2005), thus the apparent under abundances would be due to continuum emission from a secondary, but this hypothesis is largely discounted (Stütz & Paunzen, 2006). Radiatively driven atomic diffusion, which is the most widely accepted mechanism for forming peculiarities in most chemically peculiar A and B stars, has been considered for Boo stars with the addition of mass loss (Michaud & Charland, 1986). However, any turbulent or rotational mixing impedes the process sufficiently to prevent it from forming Boo peculiarities (Charbonneau, 1993). Since Boo stars rotate with similar to normal A stars, they potentially can have significant rotational mixing. Thus atomic diffusion is generally considered to be insufficient to cause Boo peculiarities.
Numerical models of Boo star atmospheres including the accretion of material depleted in heavier elements, together with atomic diffusion, have been made by Turcotte & Charbonneau (1993) and Turcotte (2002). They find that those models can produce Boo peculiarities quickly (0.1 Myr) for gas accretion rates of yr, but the peculiarities in these models do not last long after accretion has stopped (1 Myr). The peculiarities in these models appear to persist despite rotational mixing (Turcotte & Charbonneau, 1993), but larger surface convection zones require larger accretion rates to produce the same peculiarities (Turcotte, 2002).
In order to help address the question of Boo star formation, and to investigate when Ap/Bp stars develop chemical peculiarities, we have analysed high resolution spectra of 20 Herbig Ae/Be stars, and one dusty young star. We have determined atmospheric parameters and performed a detailed abundance analysis for these stars, with the goal of detecting chemical peculiarities. Three of these stars have confirmed magnetic field detections, and thus may be the progenitors of Ap/Bp stars. The other 18 stars have been carefully searched for magnetic fields, but no confirmed magnetic fields have been found.
Observations for this study were obtained with the Echelle Spectropolarimetric Device for the Observation of Stars (ESPaDOnS) instrument at the Canada France Hawaii Telescope (CFHT). ESPaDOnS is a high-resolution échelle spectropolarimeter, providing nearly continuous wavelength coverage from 3700 to 10500 Å at a resolution of . All observations were obtained in spectropolarimetric mode, which provides circularly polarised Stokes spectra, as well as total intensity Stokes spectra. Data reduction was performed with the Libre-ESpRIT (Donati et al., 1997) package, which is optimised for the ESPaDOnS instrument, and performs calibrations and optimal spectrum extraction.
The observations were obtained over several years, between 2004 and 2006 as part of an extended campaign investigating the presence of magnetic fields in Herbig Ae and Be stars. This study is discussed as a whole by Alecian et al. (2012a), and individual results are reported by Wade et al. (2005), Catala et al. (2007), Alecian et al. (2008a), Alecian et al. (2008b), Folsom et al. (2008), and Alecian et al. (2009). The observations that were analysed for this paper are reported in Table 1. Additional observations from Alecian et al. (2012a) were used to check for variability in spectral lines. The high resolution and high S/N necessary for the detection of magnetic fields in metallic lines also makes these observations very well suited to a detailed abundance analysis, thus we have an excellent dataset for our study.
The stars analysed in this study are only a subset of the full set of targets observed by Alecian et al. (2012a). The stars selected for our study were chosen to cover a range of and values. The stars were also selected to generally have only modest amounts of emission in their optical spectra, which allows for a more complete and accurate abundance analysis. Preference was given to stars that were not clearly spectroscopic binaries, and to observations with a peak S/N over 200. Stars with a confirmed magnetic field and a below 15000 K were all included (except for HD 72106 A, which was studied in detail by Folsom et al., 2008).
Observations of one star (HD 101412) were obtained from the Anglo-Australian Telescope (AAT) with the University College London Echelle Spectrograph (UCLES) together with SEMPOL polarimeter. This instrument consists of a bench mounted cross-dispersed échelle spectrograph (UCLES) fibre fed from a Cassegrain mounted polarimeter unit (SEMPOL), and is fundamentally a similar instrument to ESPaDOnS. Details of the instrument can be found in Semel et al. (1993) and Donati et al. (1997); Donati et al. (2003), and the observing run that the observations of HD 101412 were obtained in is described by Marsden et al. (2011). Data reduction and optimal spectral extraction were performed with a generic version of ESpRIT (Donati et al., 1997).
Supplemental archival observations from the FORS1 instrument at the Very Large Telescope (VLT) were used for examining Balmer lines in many stars, and are listed in Table 2. FORS1 is a low resolution spectropolarimeter, which produces an observed spectrum entirely in one order. This is useful as it reduces the possibility of normalisation errors across Balmer lines. The observations were obtained from the study of Wade et al. (2007), and a detailed description of the observations and data reduction techniques can be found therein. While Wade et al. (2007) focused on detecting magnetic fields, we only concern ourselves with Balmer line profile shapes.
3 Fundamental parameters
Temperature and surface gravity were first measured by fitting Balmer lines, and then confirmed and, if possible, refined by enforcing local thermodynamic equilibrium (LTE) ionisation and excitation balances in the metallic line spectra of the stars during the modelling procedure.
Balmer line fitting was performed on FORS1 spectra when available, since it is easier to normalise broad spectral features when they are recorded entirely in one spectral order. When FORS1 spectra were not available, ESPaDOnS spectra were used and careful attention was paid to the normalisation of Balmer lines. Continuum points well outside the wings on both sides of a Balmer line were carefully chosen, a low order polynomial was then fit through these points, and the polynomial was checked to ensure it varied only a small amount over the Balmer line. In the ESPaDOnS spectra, careful data reduction and the absence of instrumental artifacts produced an excellent correspondence between spectral orders in regions where the orders overlap. These overlap regions were used to evaluate and minimise errors in normalisation across spectral orders. For stars with both FORS1 and ESPaDOnS spectra available the normalised Balmer lines were compared to ensure the accuracy of the normalisation procedure for ESPaDOnS spectra. A good agreement between the normalised ESPaDOnS and FORS1 Balmer lines was consistently found. In all cases potential normalisation errors are included in the quoted uncertainties of physical parameters.
Careful attention was paid to avoid contamination of the Balmer lines by emission. In our observations the cores of most Balmer lines were partially infilled by emission. H was often heavily contaminated by emission, and in many stars was avoided entirely. Despite this, the wings of H and H generally remain free from emission. The wavelength range of possible contamination in the higher Balmer lines can usually be assessed by the more prominent emission in H. As a consequence of these complications, the fitting procedure was performed by hand and rather conservative uncertainties were adopted. Synthetic Balmer lines were computed from ATLAS9 (Kurucz, 1993) model atmospheres using solar chemical abundances.
Problematically, for many of our stars there is a substantial degeneracy in effective temperature and surface gravity when fitting only the wings of Balmer lines. As a result there is a large covariance between and , and our uncertainties often represent long ellipses in the two dimensional parameter space. Typically these degenerate regions have a slope of +250 K in for +0.1 dex in , though the slope generally becomes steeper towards hotter temperatures.
When good quality fits to the observations were obtainable, fitting of and for the metallic spectra was included in our determination. The fitting procedure is described in Sect. 4. This method relies on the ionisation and excitation balance in the star assuming LTE, but a number of potential problems can occur. Lines with any emission infilling must be avoided. Poor quality atomic data, undiagnosed line blends and non-LTE effects may also be problematic. This also requires lines with a wide range of excitation energies and multiple ionisation states. Consequently, this technique could not successfully be applied to all stars in our sample. In a number of cases the excitation potential could not provide a useful constraint, but the ionisation potential was still useful. In these cases the constraint provided by Balmer line fits, combined with the ionisation balance could still provide fairly precise and values.
In order to place the stars of our sample on the H-R diagram, we needed to determine accurate luminosities. Photometric Johnson magnitudes and colours were derived from the Hipparcos magnitudes using the conversion from ESA (1997). For stars without Hipparcos observations, Johnson magnitudes and colours were taken from Vieira et al. (2003), Herbst & Shevchenko (1999), and de Winter et al. (2001). These magnitudes and colours are consistent with those used by Alecian et al. (2012a).
Distance were derived using Hipparcos parallaxes, when available, from the new reduction by van Leeuwen (2007). When no parallax was available, the star was associated with an OB association or star forming region, and a literature distance to that was taken. Associations were taken from Alecian et al. (2012a), which agree with Vieira et al. (2003). Uncertainties for these distances were taken to be the spatial extent of the associations, which were estimated from their angular size on the sky, assuming the associations are approximately spherical. For HD 169142 no association could reliably be made, and thus a literature photometric distance from Sylvester et al. (1996) was used, and generous uncertainties were assumed. Literature sources for individual distances, along with the distances themselves, can be found in Table, 3.
Using these distances, bolometric luminosities were calculated. was calculated using the intrinsic colours of Kenyon & Hartmann (1995), and total extinction was calculated using an of 5, as suggested for HAeBe stars by Hernández et al. (2004). The bolometric correction of Balona (1994) was used.
Using our luminosities and temperatures, we placed the stars on the H-R diagram, as shown in Fig. 1. Pre-main sequence evolutionary tracks and isochrones were calculated with the cesam (Morel, 1997) evolutionary code, assuming solar metallicities. The birth line, the locus of points on the H-R diagram where a star becomes optically visible, was taken from Palla & Stahler (1993), for an accretion rate of yr. By comparison to these evolutionary tracks and isochrones, we determined stellar masses, and ages with respect to the birthline. These values, along with fractional pre-main sequence ages () and radii, are presented in Table 3. The uncertainties on mass and age for a star were based on the range of evolutionary tracks and isochrones that intersect the ellipse on the H-R diagram described by the uncertainties in the star’s luminosity and . Note that the choice of a birthline may introduce a further systematic uncertainty into our ages. For example, using the birthline of Behrend & Maeder (2001), with a mass dependent accretion rate, would generally increase our ages by 0.5 Myr (but with the actual increase varying with mass). While we do not present uncertainties on , the values should be considered approximate.
|ID||(K)||distance (pc)||age (Myr)||Magnetic|
|V380 Ori A||M|
|V380 Ori B||N|
4 Abundance analysis
The abundance analysis and determination of and microturbulence was performed by directly fitting synthetic spectra to the observations. Model spectra were calculated using the Zeeman spectrum synthesis code (Landstreet, 1988; Wade et al., 2001). This code performs polarised radiative transfer under the assumption of LTE. Optimisations to the code for stars with negligible magnetic fields have been included. A Levenberg-Marquardt minimisation procedure (e.g. Press et al., 1992) was used to determine best fit parameters.
The optimisations to the Zeeman spectrum synthesis routine exploit symmetries in a non-magnetic chemically homogeneous model of a star to dramatically reduce the amount of computation required. As in the original version of Zeeman, the visible disk of the star is divided into a number of surface elements. In this version it is assumed that the local emergent spectra only vary with projected radial distance from the centre of the disk of the star. Thus radiative transfer is only performed for a set of surface elements with different radial distances, and the results are reused for different angular positions, which minimises the number of times radiative transfer must be performed. The optimisations also assume the local line absorption (Voigt) and anomalous dispersion (Faraday-Voigt) profiles only vary vertically through the atmosphere, and not with position on the stellar disk, substantially reducing the number of times these profiles must be calculated. In this version of the code, limb darkening is still calculated directly, by performing radiative transfer at different radial positions on the disk of the star. Rotational broadening is also calculated directly by summing local emergent spectra that have been Doppler shifted across the disk of the star.
Detailed comparisons of the optimised code to the original version have been performed. For the same input parameters, with no magnetic field, identical spectra are produced. The original code has been checked against other spectrum synthesis codes in detail by Wade et al. (2001), with a very good agreement found. Similarly, the optimised version of Zeeman has been checked against Synth3 (Kochukhov, 2007) and again a good agreement was found.
The Levenberg-Marquardt minimisation procedure (Press et al., 1992) provides a robust and efficient routine for iteratively fitting synthetic spectra to observations. Best fit models from this routine have been compared to fits obtained by hand for a wide variety of stars, and good agreement between the two methods has been consistently found.
Input atomic data were taken from the Vienna Atomic Line Database (VALD) (Kupka et al., 1999), using an ‘extract stellar’ request and solar chemical abundances. Updated line lists with peculiar abundances were extracted as necessary. Model atmospheres were computed using the ATLAS9 code (Kurucz, 1993) with solar abundances. Model atmospheres were calculated on a grid in steps of 250 K in and 0.5 in , then interpolated linearly when more precise values were desired. Interpolation on this grid of models introduces less then 1% relative error to the line depths in the resulting spectrum.
Spectra were initially fit simultaneously for chemical abundances, , microturbulence, and radial velocity, using the best fit and determined from the Balmer lines. Then, if a good constraint could be obtained, the fits were repeated including and as free parameters. If this produced a good fit with well constrained parameters, then these and values were used instead of the initial Balmer line fit values. If only one of and was well constrained from the metallic line fit, then the fit was repeated using the Balmer line value for the unconstrained parameter. We considered a parameter to be well constrained if the fitting routine reliably converged to the same value for different initial conditions, and if that value was consistent with the Balmer line profiles. The fits were generally performed on 5 independent regions 500 Å long. The approximate wavelength ranges usually were 4400-4800, 4900-5500, 5500-6000, 6000-6500, and 6600-7600 Å, and occasionally 4150-4280 Å, with significant gaps due to non-photospheric features. These windows typically contained several hundred spectral lines (the number varying greatly with wavelength and ), however many of these lines are very weak and hence do not provide significant constraints on the best fit abundances or atmospheric parameters. The exact wavelength ranges varied between stars, due to varying regions of emission and different Balmer line widths.
Microturbulence and were determined by minimisation, simultaneously with the other stellar parameters, using the entire spectral 500 Å window. This method relies on resolved, rotationally broadened spectral lines for constraining , and a range of both weak and strong lines for constraining microturbulence.
We verified this method of determining and by analysing well studied stars, and comparing to literature results. In this study we have analysed Cet, and compared our results with the very precise study by Fossati et al. (2009). Fossati et al. (2009) used photometry and spectral energy distributions, as well as Balmer lines, and excitation and ionisation balances when determining their parameters. They also used different software tools to perform their modelling, thus their results are truly independent. Our results are fully consistent with theirs, as discussed in Sect. A.1, which demonstrates the accuracy of our methodology.
For this method to be successful, lines contaminated by emission must be avoided. While a line entirely in emission is easily identified, lines with small amounts of emission infilling are harder to identify. When multiple observations were available, variability in the emission could be used to identify lines with emission infilling. If only one observation was available, special attention was paid to lines with low excitation potentials and large oscillator strengths, such lines being more likely to contain emission. If inconsistent fits were obtained between normal and low excitation potential lines, and no other clear explanation for the inconsistency could be found (such as an error in or the atomic data), the low excitation potential lines were considered to likely contain emission and were excluded from the final fit. Additionally, the shapes of line profiles were examined and, if they could not be reproduced by the synthetic spectra, they were considered to be contaminated by a non-photospheric component. The presence of veiling in our spectra is assessed in detail in Section 6.4, and we find it is not present at a significant level.
The use of multiple spectral windows is valuable, as it allows us to take the average and standard deviation of abundances determined across the windows. The standard deviation in particular provides a reasonable estimate of the uncertainty of derived parameters. Since line-to-line abundances scatter is introduced by errors in the atomic data as well as errors in the model atmosphere, these effects are included in the standard deviation. The standard deviation may underestimate the impact of such errors, but provides a robust and easily understandable uncertainty estimate. The use of multiple spectral windows also allows us to verify the atmospheric parameters derived from metal line fits. If the parameters differ substantially from window to window, they are likely poorly constrained in some if not all windows, in which case we can rely more heavily on the parameters from Balmer line fitting.
For chemical elements with fewer than four lines providing good constraints on the abundance in any spectral window (including all cases of elements with lines in only one or two spectral windows) the uncertainties were estimated by eye rather than using a standard deviation. These uncertainties were chosen to include line-to-line scatter, uncertainties in the atmospheric parameters, and potential normalisation errors.
Examples of typical best fit synthetic spectra compared to the observed spectrum of HD 139614 in the 5000-5050 Å and 6100-6200 Å windows are presented in Fig. 2. Final averaged best fit abundances and atmospheric parameters are presented in Table 4, together with uncertainties. Elements with abundances based on less than four lines are indicated with an asterisk in the table.
We find 11 stars with varying degrees of Boo peculiarity, one star (V380 Ori A) that appears to be a weak Ap/Bp star, and 9 stars that are chemically normal (although results for V380 Ori B are very uncertain). The properties of the individual stars are discussed in detail in the Appendix. The final best fit abundances are plotted relative to solar abundances from Grevesse et al. (2005) in Figures 6, 7, 8, and 9. Illustrations of our fits for all of the stars in our study can be found in Figures 10 to 29.
|HD 17081||HD 31293||HD 31648||HD 36112||HD 68695||HD 139614||HD 141569||Solar|
|( Cet)||(AB Aur)|
|HD 142666||HD 144432||HD 163296||HD 169142||HD 176386||HD 179218||HD 244604||Solar|
In order to investigate trends in chemical abundance relative to other stellar parameters, we constructed a peculiarity index to describe whether a star is chemically normal, more like a Boo star, or more like an Ap/Bp star. This index is designed to be positive for stars with overabundances of iron peak elements, zero for chemically normal stars, and negative for stars underabundant in iron peak elements. The peculiarity index is the difference between the average of the Cr, Fe, and Ni abundances relative to solar and the average of the C, N, and O abundances relative to solar. Specifically, the index is calculated as:
where the abundances are relative to solar: . These particular elements were chosen because they tend to have peculiar abundances in chemically peculiar stars (for Cr, Fe, and Ni) or normal abundances in chemically normal stars (for C, N and O). These elements were also chosen because reliable abundances for them are available for most stars in our sample. In the few cases in which we could not determine an abundance for an element used in the peculiarity index, that element was left out of the average.
We plotted the index against effective temperature, mass, age, and fractional pre-main sequence age, as shown in Fig. 3. No clear trends were found, with a large scatter in chemical abundance at most temperatures, masses, and ages. This peculiarity index was also compared to stellar radius, , and microturbulence, but again no clear trends were found. Thus we conclude that the Boo peculiarities that we see are not restricted to a narrow mass range, and that they are present throughout the pre-main sequence.
6.1 Boo peculiarities
The most striking result of this study is the discovery of a large number of HAeBe stars with moderate Boo type chemical peculiarities. This trend seems to be present in both stars with and without detected magnetic fields, as shown by HD 101412. We find 7 stars showing fairly strong Boo peculiarities, and another 4 stars with weaker Boo peculiarities, while 8 stars are chemically normal and 1 is rather uncertain but appears to be normal. With the derived abundances, we see a continuous distribution of Boo peculiarities, from fairly strong to nearly undetectable. This incidence of roughly 50% is much higher than the incidence of Boo stars on the main sequence of roughly 2% (Paunzen, 2001). While a detection rate of 50% may be an overestimate (since our sample is biased towards stars with limited emission in their optical spectra), the true incidence of Boo peculiarities in HAeBe stars certainly is much higher than observed on the main sequence, suggesting an important connection with the pre-main sequence evolutionary phase.
The sample of stars in our study is biased towards stars with weaker circumstellar emission or absorption in their spectra, since these stars allow for a more complete analysis of their photospheric spectra. However, the stars were drawn from the larger sample of Alecian et al. (2012a) which did not have this bias. We found 65% of the stars in their sample had modest amounts of circumstellar contamination, and would have been suitable for our photospheric analysis. Thus, as a lower limit, 33% of HAeBe stars should show Boo peculiarities, assuming that all the unsuitable stars in the sample of Alecian et al. (2012a) are chemically normal, and using our derived incidence of peculiarities for the suitable stars. The true incidence is likely higher than this lower limit.
There are almost certainly more undetected Boo stars in the sample of Alecian et al. (2012a) that we did not analyse, simply because analysing their complete sample was beyond the scope of this study. Indeed Alecian et al. (2012a) mention that HD 34282 appears to have weak lines of iron peak elements, but a normal O line. To roughly estimate the number of potential Boo stars in the remaining sample, we compared synthetic spectra computed using literature values for their atmospheric parameters and assuming solar abundances with the ESPaDOnS observations. In doing so we identified 5 additional stars that have suspiciously weak metal lines for their literature temperatures: HD 37806, HD 98922, HD 144668, HD 76534, and VX Cas. If all these stars do indeed show Boo peculiarities, 35% of the stars observed with ESPaDOnS by Alecian et al. (2012a) would be Boo stars, which is consistent with our results.
This large incidence of Boo type peculiarities in HAeBe stars is qualitatively consistent with a selective accretion model of Boo stars. HAeBe stars have recently undergone accretion, and may still be accreting, so it is reasonable that a process which depends on accretion is seen more frequently in these stars. Thus we consider our results supportive of the selective accretion model of Boo stars.
The large number of pre-main sequence Boo stars we find suggests that the majority of Boo stars may develop their peculiarities on the pre-main sequence by accreting their circumstellar material. However, we cannot rule out the existence more of evolved Boo stars that only develop peculiarities on the main sequence, by accreting from diffuse interstellar clouds. Due to the relatively short lifetime of Boo peculiarities after accretion has halted, the presence of such peculiarities later in the main sequence (e.g. Paunzen et al., 2002) may require a different source of accreted material from that of pre-main sequence Boo stars.
Turcotte & Charbonneau (1993) modelled the diffusion of chemical elements in the atmosphere of a star (with K and ) that was accreting gas depleted in iron peak elements. They found that Boo peculiarities can be generated quickly (0.1 Myr), but also dissipate quickly after accretion has stopped (1 Myr). This implies that, in the context of selective accretion, there is a good chance that HAeBe stars with Boo peculiarities are still accreting, or were very recently accreting.
The range in strengths of the Boo peculiarities seen in the stars in our sample might be interpreted as the consequence of different accretion rates in the stars. However, Turcotte & Charbonneau (1993) do not find a strong impact on surface abundances from different accretion rates, as long as they are sufficient to overwhelm diffusion and rotational mixing.
Alternately, the differences in the strengths of peculiarities could reflect the abundances of the accreted gas. These differences being due, for example, to the efficiency with which iron peak elements are bound into dust grains, or to the degree to which dust is blown away from the star while gas is accreted.
In this context, the chemically normal stars are most likely not accreting significant amounts of material, and thus any Boo peculiarities have dissipated. However, it is possible that accretion is ongoing but the selection process in the circumstellar material has broken down. For example, dust grain formation may not be proceeding efficiently, and thus the accreted gas would not be depleted in iron peak elements
Unfortunately, deriving accretion rates for HAeBe stars is not a simple process. For example, Donehew & Brittain (2011) measured accretion rates based on Balmer jump excess fluxes. They obtained values between a few and yr. These values rely on models from Muzerolle et al. (2004), which are constructed assuming magnetospheric accretion. Given the absence of detectable magnetic fields in most of the stars in our study, magnetospheric accretion may not be an accurate model (see Wade et al., 2007), and thus these accretion rates may not be accurate either. Generally, accretion rates derived for HAeBe stars are based on rather simple models, and thus it is not clear how reliable those accretion rates are. Eleven of the stars in our sample have measured accretion rates by Donehew & Brittain (2011). Supplementing these results with accretion rates from Mendigutía et al. (2011), who used a similar methodology, and from Garcia Lopez et al. (2006), who based their accretion rates on Br emission, provides literature values for 15 stars in our sample. Worryingly, the literature accretion rates differ for five of these stars by between 0.5 and 1 dex. Whether this apparent disagreement is due to intrinsic stellar variability, differences in methodology, or simply large uncertainties is not clear (Mendigutía et al., 2011). We compare these accretion rates with our peculiarity index (defined in Sect. 5), stellar mass, and fractional pre-main sequence age in Fig. 4. In cases where multiple accretion rates are available the average value was taken. V380 Ori was left out of the comparison due to its binarity. We see no correlation between accretion rate and chemical peculiarity or mass. There is possibly a weak correlation with fractional pre-main sequence age, but due to the large uncertainties on both the accretion rates and the fractional ages it is not entirely clear.
Our results provide further evidence against the hypothesis that Boo stars are formed by radiative diffusion with mass loss, as proposed by Michaud & Charland (1986). In their calculations, Michaud & Charland (1986) find that it takes or years to build up Boo peculiarities by diffusion in the presence of mass loss. This is much greater than the ages of the Boo-like stars in our sample, which are only a few years old. Thus the time scale for this mechanism is dramatically incompatible with our observations. This is consistent with the generally accepted view that diffusion with mass loss is insufficient for producing Boo peculiarities (e.g. Charbonneau, 1993) .
6.2 Magnetic stars
The one possible Bp star in our sample is V380 Ori A. This star shows modest overabundances of Fe, Mn, and Ni of 0.6 dex, and Si is also enhanced relative to solar. C, Ne, and S are consistent with solar, as are O and Mg at 2. This pattern is characteristic of Ap/Bp stars. Arguably, results are marginally consistent with the star simply having a high metallicity (0.4 dex above solar). However, since Fe, Mn and Ni are more enhanced than 0.4 dex, and C, O, Ne, Mg and S are less enhanced than this, we consider it more likely that the star is weakly chemically peculiar. If the star is a Bp star, the star has begun developing Ap/Bp chemical peculiarities, but they have not yet reached the degree seen in most main sequence Ap and Bp stars.
Atomic diffusion producing chemical peculiarities was modelled in V380 Ori A by Vick et al. (2011), using a model that includes stellar evolution, as well as mass loss and atomic diffusion. They used the stellar parameters determined by Alecian et al. (2009), who find a somewhat lower mass than we do (mostly due to assuming rather than 5, thus using less extinction). Nevertheless, the abundances they derive should be roughly representative. Vick et al. (2011) find Mn, Fe and Ni to be overabundant by 0.5 dex, which agrees nicely with our results. They find roughly solar abundances for C, O and S, which agrees with our results. The overabundance of P that they find is also roughly consistent with our results. However, the underabundant Mg, and solar Si abundance that they find are inconsistent with our results. Thus, while this model is imperfect, we find it does provide a reasonable approximation of the majority of the observed peculiarities in V380 Ori A.
Interestingly, HD 190073 and HD 101412 are both magnetic but do not show Ap/Bp chemical peculiarities. HD 190073 is chemically normal, while HD 101412 shows Boo type chemical peculiarities. Thus, unlike on the main sequence, not all pre-main sequence magnetic A and B stars are chemically peculiar. Atomic diffusion theory suggests that enough time must elapse after accretion and deep surface convection have halted for significant peculiarities to develop. The dissipation of a large convective envelope can be modelled easily enough, and in most models dissipates early in pre-main sequence life of A and B stars. However, determining when accretion has dropped sufficiently for atomic diffusion to dominate is much harder. Indeed it may depend on the circumstellar environment of the particular star, and thus not be reliably predictable.
The Boo peculiarities of HD 101412 presumably simply reflect the high incidence of Boo peculiarities in HAeBe stars in general. Thus, we suspect that these peculiarities are not a consequence of the magnetic field of HD 101412.
Modelling by Vick et al. (2011) suggests that HD 190073 has not yet had time to develop Ap/Bp chemical peculiarities. In their models, a star of this mass requires 2 Myr to develop peculiarities by atomic diffusion, while we find HD 190073 is only Myr old. HD 101412 has almost the same mass, and we find an age of only Myr, thus by the models of Vick et al. (2011) it is also too young to have developed Ap/Bp peculiarities. The younger age of HD 101412 is consistent with the picture of its Boo type chemical peculiarities being due to selective accretion, since accretion is more likely to be ongoing earlier in a stars life.
It is interesting that the two stars in this sample with the largest amount of emission in their optical spectra are V380 Ori A and HD 190073, both of which are magnetic. In particular, the emission mostly appears in lines that would be strongly in absorption for a purely photospheric spectrum of the star. HD 101412, the third magnetic star, displays a fairly modest amount of emission. Additionally, the measured accretion rates for V380 Ori and HD 190073 are both substantially higher than for other stars in our sample ( yr for V380 Ori, yr for HD 190073; Donehew & Brittain, 2011; Mendigutía et al., 2011). The impact of the secondary in V380 Ori does not appear to be considered in the measurement of the star’s accretion rate, and thus it may be somewhat inaccurate. In classical T Tauri stars, the strong emission in their spectra is thought to result from magnetically channelled accretion flows falling onto the star. We speculate that the strong emission in V380 Ori A and HD 190073 might be produced by a similar process. However, it is not clear whether this is the case, both because of the small sample size, and because the non-magnetic stars in the sample are biased towards lower emission stars (since large amounts of emission limit the abundance analysis). Magnetic fields are not detected in individual emission lines for any of our stars, thus we cannot use those lines to assess the presence of magnetospheric accretion directly.
The three stars in our sample with confirmed magnetic fields also have the lowest values (for more analysis see Alecian et al., 2012b). In this respect these three magnetic HAeBe star are similar to Ap/Bp stars, which rotate slower than normal non-magnetic A and B stars. The low of the magnetic HAeBe stars HD 200775A ( km s, Alecian et al., 2008a) and of NGC 2244 OI 201 ( km s, Alecian et al., 2008b), and the moderate of HD 72106 A ( km s, Folsom et al., 2008) support this observation. However, the very high of NGC 6611 W610 ( km s, Alecian et al., 2008b) indicates that not all magnetic HAeBe stars rotate slowly. Nevertheless, the large number of low magnetic HAeBe stars suggests that magnetic braking may have occurred in many of these stars (see Alecian et al., 2012b).
Several other very young magnetic A and B stars are known to have chemical peculiarities. Folsom et al. (2008) examined the magnetic star HD 72106 A, and found Ap/Bp peculiarities. They find that the star has a mass of 2.4 . While they were unable to firmly conclude that HD 72106 A is on the pre-main sequence, they do find that the system is between 6 and 13 Myr from the birthline, based on the H-R diagram position of the secondary in the binary system. The pattern of abundances seen in HD 72106 A is similar to that seen in V380 Ori A, but more extreme. Overabundances of iron peak elements and a weak overabundance of Si is seen in both stars. However, the strong He underabundance seen in HD 72106 A is not seen in V380 Ori A. The differences in abundances between the stars could be due to V380 Ori A being hotter and more massive HD 72106 A, or due to HD 72106 A being more evolved.
Bagnulo et al. (2004) investigated the main sequence magnetic star NGC 2244-334, which is a member of the open cluster NGC 2244, and found Ap/Bp peculiarities. They found that the star has a mass of 4 and commented that, as a cluster member, the star is likely 2 Myr old as measured from the ZAMS. This is a similar mass to V380 Ori A, though NGC 2244-334 is somewhat older. The peculiarities they found are very similar to those of HD 72106 A, with overabundances in iron peak elements of 1-2 dex, a somewhat weaker overabundance of Si, and a strongly underabundant He. This is again similar to what we see in V380 Ori A, though much stronger, and with a He underabundance.
Alecian et al. (2008b) studied the magnetic Herbig Be star NGC 6611 W601, and noted peculiarly strong He and Si lines in its spectrum. This suggest the star likely represents a pre-main sequence He-strong star, but a detailed abundance analysis has yet to be performed. While this pattern of peculiarities would be different from that of V380 Ori A, the star has a much higher temperature of K and a much higher mass of (Alecian et al., 2008b). For a main sequence magnetic star in this temperature and mass range, one would expect to see He-strong peculiarities, rather than Ap/Bp peculiarities. Alecian et al. (2008b) estimate the age of NGC 6611 W601 to be less than 0.06 Myr from the birthline of Palla & Stahler (1993, for an accretion rate of yr), based on its H-R diagram position. While this is very young, there are no models from Vick et al. (2011) for a comparable mass. We do note, however, that Vick et al. (2011) find peculiarities developing faster for more massive models.
Thus, while the peculiarities seen in V380 Ori A are weak, the pattern of peculiarities is consistent with those seen in other young intermediate mass magnetic stars. This is consistent with the hypothesis that the weak peculiarities we see in this star will develop into the stronger peculiarities seen in most main sequence Ap/Bp stars.
6.3 Chemically normal stars
We find no hints of any chemical peculiarities in the stars of our sample other than Boo peculiarities and the weak Bp peculiarities in V380 Ori A. Thus we find no evidence for Am or HgMn stars on the pre-main sequence. Our sample is not large enough to definitively exclude the existence of Am and HgMn stars at the same incidence as on the main sequence. However, if the roughly 10-20% incidence (e.g. Smith, 1996) of these stars on the main sequence is extended to the pre-main sequence, we would expect to have seen two to four such stars. Thus it is possible that these peculiarities do not form on the pre-main sequence, unlike for Ap/Bp stars. The atomic diffusion thought to give rise to these peculiarities may be disrupted by weak accretion more readily than the diffusion in the presence of a magnetic field that produces Ap/Bp stars. The presence of a magnetic field may stabilise a stellar atmosphere and increase diffusion velocities, potentially making diffusion in the presence of a magnetic field more robust (Michaud, 1970; Michaud et al., 1981). Alternately, diffusion in non-magnetic stars may simply require more time than the pre-main sequence lifetime of these stars to produce significant peculiarities.
The large majority of young A and B stars analysed, both on the main sequence and on the pre-main sequence (e.g. Acke & Waelkens, 2004), have been found to be chemically normal. A very precise study of A and B star abundances by Fossati et al. (2009) investigated what ‘normal’ abundances are for A and B stars. They found that solar abundances are a good approximation, though there were small deviations for specific elements, which appeared to depend on the temperature of the star. Studies of nearby open clusters support this conclusion (e.g. Fossati et al., 2008; Villanova et al., 2009; Fossati et al., 2011).
Acke & Waelkens (2004) performed an abundance analysis of 24 pre-main sequence and very young main sequence stars, based on equivalent widths. A comparison of individual results for the 7 stars we have in common with their study is discussed in Sect. A. Our results are generally consistent, although we are able to determine abundances for a wider range of elements, and in many cases provide more precise abundances. The one notable difference in our results is for HD 139614. For this star our more accurate allows us to classify the star as Boo, while Acke & Waelkens (2004) could only comment that it appeared to be metal weak. In total Acke & Waelkens (2004) found only one Herbig star that was also a clear Boo star (HD 100546), and they noted that AB Aur shows hints of Boo peculiarity. In part this difference in detection rates may due to our more accurate and values, as the case of HD 139614 demonstrates. Our higher precision, due to an improved methodology and higher S/N observations, is also a significant part of this difference. Finally, both our studies deal with a fairly small number of stars, and thus the true incidence of Boo peculiarities in HAeBe stars remains somewhat uncertain.
6.4 Possible Veiling
One concern when analysing spectra of T Tauri stars is veiling. Veiling is additional continuum emission, usually from an accretion shock (or other circumstellar material), which makes spectral lines appear weaker than they would in a purely photospheric spectrum. In HAeBe stars, veiling can be seen in the UV (e.g. Donehew & Brittain, 2011), but usually is considered to be unimportant in the optical (e.g. Böhm & Catala, 1993; Hernández et al., 2004; Cowley et al., 2010; Donehew & Brittain, 2011). HAeBe stars are intrinsically much more luminous than T Tauri stars, thus much more emission from an accretion shock would be needed to have the same veiling effect. Measured accretion rates for HAeBe stars are usually not much larger than for T Tauri stars, thus it is unlikely that HAeBe stars would have much larger amounts of emission from accretion shocks. Consequently, we would not expect veiling to have much impact on HAeBe star spectra.
If veiling is present in our observations it could make a star appear to be metal poor when it is not, thus it is important to check for veiling. If it is present in our spectra, veiling should have some wavelength dependence. We have compared abundances derived from different wavelength regions for individual stars. No clear correlation of abundance with wavelength was found, suggesting that veiling is unimportant in these stars. Veiling should affect all lines in a wavelength region similarly. Thus, if it is present veiling should affect most elements similarly, since the lines for most elements are widely distributed in our spectra. In particular, it could not produce underabundances in iron peak elements while keeping solar abundances of C and O that we see in many of the stars in our sample. Therefore veiling should not affect the relative abundances of C, N, and O to Fe, even if it were present at a significant level in our spectra. We conclude it is unlikely that veiling has a significant impact on our derived abundances, and even if veiling is present it is not obvious how it could produce peculiar abundance patterns.
Arguably, the pattern of Boo underabundances we observe could be due to weak emission infilling in all the lines of some elements, however this is very unlikely. We have searched carefully for emission in lines, both by looking for unexpected variability and by looking for lines with shapes inconsistent with rotational broadening. Thus we have excluded most if not all lines with emission. More importantly, lines of an element with lower excitation potentials are more likely to have stronger emission, thus emission would not be uniform across all lines. If emission infilling was occurring and otherwise undetected, we would observe it as an inability to simultaneously fit a wide range of lines of one element with one abundance. We do not find this, and consequently can safely conclude that the Boo peculiarities we detect are real, and not simply due to emission infilling.
We have performed a detailed abundance analysis of 21 stars, 20 of which are Herbig Ae and Be stars, 1 of which ( Cet) is not a Herbig star, but may still be on the pre-main sequence. Three of these stars have confirmed magnetic fields, while magnetic fields do not appear to be present at detectable levels in the remaining 18.
We find 9 stars that are chemically normal (one of which is rather uncertain), 11 that show Boötis chemical peculiarities, and one that is weakly Ap/Bp. This is a remarkably large fraction of stars that display Boo type peculiarities. On the main sequence, only 2% of stars display Boo peculiarities (Paunzen, 2001). We interpret this as evidence in favour of a selective accretion hypothesis for the development of Boo peculiarities. Since Herbig stars are either accreting or have recently finished accreting, a process which depends on accretion should be seen more frequently in these stars.
Of the three magnetic stars, HD 101412 displays Boo peculiarities, HD 190073 is chemically normal, and V380 Ori is weakly Bp. Thus it appears that magnetic Herbig stars can display Ap/Bp type peculiarities, but not all magnetic Herbig stars do. This is in contrast to main sequence magnetic A and late B stars, which always display Ap/Bp type peculiarities (Aurière et al., 2007). The presence of Boo peculiarities in a magnetic star appears to simply reflect the high incidence of these peculiarities in all Herbig stars.
GAW and JDL are supported by Natural Science and Engineering Research Council (NSERC Canada) Discovery Grants, and GAW is supported by a Department of National Defence (Canada) ARP grant.
- Abt & Morrell (1995) Abt H. A., Morrell N. I., 1995, ApJS, 99, 135
- Acke & Waelkens (2004) Acke B., Waelkens C., 2004, A&A, 427, 1009
- Alecian et al. (2008a) Alecian E., Catala C., Wade G. A., Donati J.-F., Petit P., Landstreet J. D., Böhm T., Bouret J.-C., Bagnulo S., Folsom C., Grunhut J., Silvester J., 2008a, MNRAS, 385, 391
- Alecian et al. (2008b) Alecian E., Wade G. A., Catala C., Bagnulo S., Boehm T., Bohlender D., Bouret J.-C., Donati J.-F., Folsom C. P., Grunhut J., Landstreet J. D., 2008b, A&A, 481, L99
- Alecian et al. (2009) Alecian E., Wade G. A., Catala C., Bagnulo S., Böhm T., Bouret J., Donati J., Folsom C. P., Grunhut J., Landstreet J. D., 2009, MNRAS, 400, 354
- Alecian et al. (2012a) Alecian E., Wade G. A., Catala C., Grunhut J. H., Boehm T., Bouret J.-C., Flood J., Folsom C. P., Landstreet J. D., Marsden S. C., Waite I. A., 2012a, MNRAS, submitted
- Alecian et al. (2012b) Alecian E., Wade G. A., Catala C., Grunhut J. H., Landstreet J. D., Boehm T., Bouret J.-C., Flood J., Folsom C. P., 2012b, MNRAS, submitted
- Andrievsky et al. (2002) Andrievsky S. M., Chernyshova I. V., Paunzen E., Weiss W. W., Korotin S. A., Beletsky Y. V., Handler G., Heiter U., Korotina L., Stütz C., Weber M., 2002, A&A, 396, 641
- Andrievsky & Paunzen (2000) Andrievsky S. M., Paunzen E., 2000, MNRAS, 313, 547
- Aurière et al. (2007) Aurière M., Wade G. A., Silvester J., Lignières F., Bagnulo S., Bale K., Dintrans B., Donati J. F., et al. 2007, A&A, 475, 1053
- Bagnulo et al. (2004) Bagnulo S., Hensberge H., Landstreet J. D., Szeifert T., Wade G. A., 2004, A&A, 416, 1149
- Bagnulo et al. (2011) Bagnulo S., Landstreet J. D., Fossati L., Kochukhov O., 2011, A&A, submitted
- Balona (1994) Balona L. A., 1994, MNRAS, 268, 119
- Behrend & Maeder (2001) Behrend R., Maeder A., 2001, A&A, 373, 190
- Bohlender & Landstreet (1990) Bohlender D. A., Landstreet J. D., 1990, MNRAS, 247, 606
- Böhm & Catala (1993) Böhm T., Catala C., 1993, A&AS, 101, 629
- Brown et al. (1994) Brown A. G. A., de Geus E. J., de Zeeuw P. T., 1994, A&A, 289, 101
- Catala et al. (2007) Catala C., Alecian E., Donati J.-F., Wade G. A., Landstreet J. D., Böhm T., Bouret J.-C., Bagnulo S., Folsom C., Silvester J., 2007, A&A, 462, 293
- Charbonneau (1993) Charbonneau P., 1993, ApJ, 405, 720
- Cowley et al. (2010) Cowley C. R., Hubrig S., González J. F., Savanov I., 2010, A&A, 523, A65
- de Winter et al. (2001) de Winter D., van den Ancker M. E., Maira A., Thé P. S., Djie H. R. E. T. A., Redondo I., Eiroa C., Molster F. J., 2001, A&A, 380, 609
- de Zeeuw et al. (1999) de Zeeuw P. T., Hoogerwerf R., de Bruijne J. H. J., Brown A. G. A., Blaauw A., 1999, AJ, 117, 354
- Dolan & Mathieu (2001) Dolan C. J., Mathieu R. D., 2001, AJ, 121, 2124
- Donati et al. (2003) Donati J.-F., Collier Cameron A., Semel M., Hussain G. A. J., Petit P., Carter B. D., Marsden S. C., Mengel M., López Ariste A., Jeffers S. V., Rees D. E., 2003, MNRAS, 345, 1145
- Donati et al. (1997) Donati J.-F., Semel M., Carter B. D., Rees D. E., Collier Cameron A., 1997, MNRAS, 291, 658
- Donehew & Brittain (2011) Donehew B., Brittain S., 2011, AJ, 141, 46
- ESA (1997) ESA, 1997, The Hipparcos and Tycho Catalogues, ESA SP-1200
- Faraggiana & Bonifacio (2005) Faraggiana R., Bonifacio P., 2005, A&A, 436, 697
- Folsom et al. (2008) Folsom C. P., Wade G. A., Kochukhov O., Alecian E., Catala C., Bagnulo S., Böhm T., Bouret J., Donati J., Grunhut J., Hanes D. A., Landstreet J. D., 2008, MNRAS, 391, 901
- Fossati et al. (2008) Fossati L., Bagnulo S., Landstreet J., Wade G., Kochukhov O., Monier R., Weiss W., Gebran M., 2008, A&A, 483, 891
- Fossati et al. (2011) Fossati L., Folsom C. P., Bagnulo S., Grunhut J. H., Kochukhov O., Landstreet J. D., Paladini C., Wade G. A., 2011, MNRAS, 413, 1132
- Fossati et al. (2009) Fossati L., Ryabchikova T., Bagnulo S., Alecian E., Grunhut J., Kochukhov O., Wade G., 2009, A&A, 503, 945
- Garcia Lopez et al. (2006) Garcia Lopez R., Natta A., Testi L., Habart E., 2006, A&A, 459, 837
- Grevesse et al. (2005) Grevesse N., Asplund M., Sauval A. J., 2005, in Alecian G., Richard O., Vauclair S., eds, Element Stratification in Stars: 40 Years of Atomic Diffusion Vol. 17 of EAS Pub. Ser., The New Solar Chemical Composition. p. 21
- Guimarães et al. (2006) Guimarães M. M., Alencar S. H. P., Corradi W. J. B., Vieira S. L. A., 2006, A&A, 457, 581
- Heiter (2002) Heiter U., 2002, A&A, 381, 959
- Herbst & Shevchenko (1999) Herbst W., Shevchenko V. S., 1999, AJ, 118, 1043
- Hernández et al. (2004) Hernández J., Calvet N., Briceño C., Hartmann L., Berlind P., 2004, AJ, 127, 1682
- Hubrig et al. (2011) Hubrig S., Mikulášek Z., González J. F., Schöller M., Ilyin I., Curé M., Zejda M., Cowley C. R., Elkin V. G., Pogodin M. A., Yudin R. V., 2011, A&A, 525, L4
- Hubrig et al. (2007) Hubrig S., Pogodin M. A., Yudin R. V., Schöller M., Schnerr R. S., 2007, A&A, 463, 1039
- Hubrig et al. (2004) Hubrig S., Schöller M., Yudin R. V., 2004, A&A, 428, L1
- Hubrig et al. (2009) Hubrig S., Stelzer B., Schöller M., Grady C., Schütz O., Pogodin M. A., Curé M., Hamaguchi K., Yudin R. V., 2009, A&A, 502, 283
- Hubrig et al. (2006) Hubrig S., Yudin R. V., Schöller M., Pogodin M. A., 2006, A&A, 446, 1089
- Kamp & Paunzen (2002) Kamp I., Paunzen E., 2002, MNRAS, 335, L45
- Kenyon & Hartmann (1995) Kenyon S. J., Hartmann L., 1995, ApJS, 101, 117
- Kochukhov (2007) Kochukhov O. P., 2007, in I. I. Romanyuk, D. O. Kudryavtsev, O. M. Neizvestnaya, & V. M. Shapoval ed., Physics of Magnetic Stars Spectrum synthesis for magnetic, chemically stratified stellar atmospheres. pp 109–118
- Kupka et al. (1999) Kupka F., Piskunov N., Ryabchikova T. A., Stempels H. C., Weiss W. W., 1999, A&AS, 138, 119
- Kurucz (1993) Kurucz R., 1993, CDROM Model Distribution, Smithsonian Astrophys. Obs.
- Landstreet (1988) Landstreet J. D., 1988, ApJ, 326, 967
- Leinert et al. (1997) Leinert C., Richichi A., Haas M., 1997, A&A, 318, 472
- Malfait et al. (1998) Malfait K., Bogaert E., Waelkens C., 1998, A&A, 331, 211
- Marsden et al. (2011) Marsden S. C., Jardine M. M., Ramírez Vélez J. C., Alecian E., Brown C. J., Carter B. D., Donati J.-F., Dunstone N., Hart R., Semel M., Waite I. A., 2011, MNRASA, 413, 1922
- Mendigutía et al. (2011) Mendigutía I., Calvet N., Montesinos B., Mora A., Muzerolle J., Eiroa C., Oudmaijer R. D., Merín B., 2011, A&A, 535, A99
- Michaud (1970) Michaud G., 1970, ApJ, 160, 641
- Michaud & Charland (1986) Michaud G., Charland Y., 1986, ApJ, 311, 326
- Michaud et al. (1981) Michaud G., Charland Y., Megessier C., 1981, A&A, 103, 244
- Miroshnichenko et al. (2001) Miroshnichenko A. S., Bjorkman K. S., Chentsov E. L., Klochkova V. G., Gray R. O., García-Lario P., Perea Calderón J. V., 2001, A&A, 377, 854
- Morel (1997) Morel P., 1997, A&AS, 124, 597
- Muzerolle et al. (2004) Muzerolle J., D’Alessio P., Calvet N., Hartmann L., 2004, ApJ, 617, 406
- Palla & Stahler (1993) Palla F., Stahler S. W., 1993, ApJ, 418, 414
- Paunzen (2001) Paunzen E., 2001, A&A, 373, 633
- Paunzen et al. (2002) Paunzen E., Iliev I. K., Kamp I., Barzova I. S., 2002, MNRAS, 336, 1030
- Press et al. (1992) Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P., 1992, Numerical Recipes in FORTRAN, 2nd edn. Cambridge University Press, Cambridge
- Ripepi et al. (2006) Ripepi V., Bernabei S., Marconi M., Palla F., et al. 2006, A&A, 449, 335
- Semel et al. (1993) Semel M., Donati J.-F., Rees D. E., 1993, A&A, 278, 231
- Smith (1996) Smith K. C., 1996, Ap&SS, 237, 77
- Stütz & Paunzen (2006) Stütz C., Paunzen E., 2006, A&A, 458, L17
- Sylvester et al. (1996) Sylvester R. J., Skinner C. J., Barlow M. J., Mannings V., 1996, MNRAS, 279, 915
- Turcotte (2002) Turcotte S., 2002, ApJ, 573, L129
- Turcotte & Charbonneau (1993) Turcotte S., Charbonneau P., 1993, ApJ, 413, 376
- van Leeuwen (2007) van Leeuwen F., 2007, A&A, 474, 653
- Venn & Lambert (1990) Venn K. A., Lambert D. L., 1990, ApJ, 363, 234
- Vick et al. (2011) Vick M., Michaud G., Richer J., Richard O., 2011, A&A, 526, A37
- Vieira et al. (2003) Vieira S. L. A., Corradi W. J. B., Alencar S. H. P., Mendes L. T. S., Torres C. A. O., Quast G. R., Guimarães M. M., da Silva L., 2003, AJ, 126, 2971
- Villanova et al. (2009) Villanova S., Carraro G., Saviane I., 2009, A&A, 504, 845
- Wade et al. (2007) Wade G. A., Bagnulo S., Drouin D., Landstreet J. D., Monin D., 2007, MNRAS, 376, 1145
- Wade et al. (2001) Wade G. A., Bagnulo S., Kochukhov O., Landstreet J. D., Piskunov N., Stift M. J., 2001, A&A, 374, 265
- Wade et al. (2005) Wade G. A., Drouin D., Bagnulo S., Landstreet J. D., Mason E., Silvester J., Alecian E., Böhm T., Bouret J.-C., Catala C., Donati J.-F., 2005, A&A, 442, L31
- Waters et al. (1992) Waters L. B. F. M., Trams N. R., Waelkens C., 1992, A&A, 262, L37
Appendix A Results for individual stars
a.1 Cet (HD 17081)
Cet (HD 17081) is a cool B star that was studied in detail by Fossati et al. (2009). They performed a very precise, detailed abundance analysis, and considered the star to be a good example of a ‘normal’ cool B star. Cet is not a Herbig Be star by the strict definition, since it has no strong emission in its optical spectrum. However, the star does display an infrared excess (Malfait et al., 1998), as well as some unexplained line profile variations (Fossati et al., 2009), and an H-R diagram position to the right of the main sequence. Thus the star is certainly very young, and possibly still on the pre-main sequence, despite the lack of clear emission. Based on the H-R diagram position of the star, it appears to be less than 0.5 Myr from the birth line, and almost certainly is less than 1 Myr old. Malfait et al. (1998) suggested this star was a Vega-type star, a classification which is consistent with our observations.
In order to check the accuracy of our methodology, we analysed an observation of Cet and compared our results with the results of Fossati et al. (2009). Our results are completely consistent with those of Fossati et al. (2009), for temperature, gravity, , and abundances. Our atmospheric parameters and abundances are compared with those of Fossati et al. (2009) in Table 6, and the differences between our abundances and those of Fossati et al. (2009) are plotted in Fig. 5. The excellent agreement between our results and Fossati et al. (2009) demonstrates the accuracy of our methodology.
Cet was observed spectropolarimetricly by Wade et al. (2007), who did not detect a magnetic field. Alecian et al. (2012a) also found no magnetic field in Cet, and thus we consider the star to be non-magnetic.
|Fossati et al. (2009)||This work|
a.2 HD 31293 (AB Aur)
We derive K and for HD 31293 (AB Aur). This is consistent with Acke & Waelkens (2004), but our is inconsistent with their (remarkably high) value of . The and are both consistent with Böhm & Catala (1993). We find solar He and O abundances, C and N are 1 above solar, while the iron peak elements are consistently 1 below solar. The systematic underabundance of iron peak elements by 0.3 dex suggests the star may be a weak Boo star. This is strengthened somewhat by the observation that the iron peak elements are 0.5 dex below the average C, N, and O abundance. Based on these underabundances, we conclude that the star most likely is a weak Boo star, but it is a marginal case.
Acke & Waelkens (2004) derive abundances for O, N, and Si that are consistent with ours, but their Fe abundance of 1 dex below solar (based on 2 Fe i lines) is inconsistent with our value. This difference may be due to the higher Acke & Waelkens (2004) used. Alternately, it is possible their measurement was impacted by the large amount of emission in the spectrum of HD 31293. Both Fe abundances are consistent with a Boo classification of the star, but we prefer our larger abundance since it is based on larger number and a wider range of lines, and it is consistent with the neighbouring iron peak element abundances we derive.
a.3 Hd 31648
For HD 31648 we derive an of K and a of , which is consistent with Acke & Waelkens (2004). We find most abundances within 1 of solar, and all abundances except Mg and Ti within 2. The abundances are systematically slightly above solar and may reflect a slightly enhanced metallicity, or a slightly overestimated .
Our Ca and Fe abundances are consistent with Acke & Waelkens (2004), our C abundance is consistent at 2, as is our Y abundance if we assume the uncertainty on their value is comparable to ours. Our Si and Mg abundances are not consistent with Acke & Waelkens (2004), although their abundances are based on only 1 line and thus are likely uncertain. The large of this star causes significant blending in its spectrum. This could explain the discrepancies between our results, since the Acke & Waelkens (2004) abundances are based on equivalent width measurements, which can easily be influenced by unaccounted-for blending.
HD 31648 was not found to have a magnetic field by Alecian et al. (2012a). No magnetic field was found by Wade et al. (2007) either. Hubrig et al. (2006) claimed a weak detection of a magnetic field in their observation of HD 31648 but, in the re-reduction of the observation made by Hubrig et al. (2007) the field is below their detection threshold. Surprisingly, Hubrig et al. (2007) report a signal the H, H, H, and Ca H & K lines of their spectrum, which they interpret as a “circumstellar” magnetic field. This result has yet to be independently confirmed. A careful examination of the H, H, H, and Ca H & K lines in our spectra shows no polarisation signature. Bagnulo et al. (2011) re-reduced a wide range of spectropolarimetric observations from the FORS1 instrument at the VLT, including the observations of Hubrig et al. (2006). They find no evidence for a magnetic field in HD 31648. Considering the lack of evidence for a photospheric magnetic field, and the absence of a confirmation of a “circumstellar” magnetic field, we conclude that HD 31648 is non-magnetic.
a.4 Hd 36112
For HD 36112 we find K and , both of which are somewhat greater than the and of Acke & Waelkens (2004). Our uncertainty on may be overly optimistic, though it does reflect the scatter in derived from different spectral windows. The uncertainties on the calibration of become important at this level, and thus there may be some systematic uncertainty not included in our error bar. Most of our derived abundances are within 1 of solar, and all are within 2. There are no systematic enhancements or depletions, thus we conclude HD 36112 is chemically normal, with solar abundances.
Our abundances are consistent with those of Acke & Waelkens (2004) for C, Si, S, Ca, Sc, Cr, Mn, Fe, and Ni. Our abundance for Ti and Y differ by 2. These differences may be a consequence of our hotter temperature, or they may be due to the small number of lines Acke & Waelkens (2004) used for determining those abundances.
a.5 Hd 68695
For HD 68695 we derive K and . We find C, N, O, and S abundances that are 1 above solar, while iron peak abundances are 2 to 3 below solar (0.7 dex below solar), except for the rather uncertain Mn abundance. We conclude that this star has clear Boo peculiarities.
a.6 Hd 139614
HD 139614 was analysed by Acke & Waelkens (2004) using an of 8000 K. They found an overall deficiency of metals, but no clear evidence of the selective depletion characteristic of Boo stars. We find HD 139614 to have an of K, and a of , which is consistent with the results of K and reported by Guimarães et al. (2006), and marginally consistent with Acke & Waelkens (2004) (at 1.3 in and 2 in ). We consider our parameters (and Guimarães et al., 2006) to be more accurate, since Acke & Waelkens (2004) values were based on the spectral classification of Malfait et al. (1998), which could be influenced by chemical peculiarities in the star. In fact, for our parameters we find a clear Boo pattern of abundances, with solar C, N, and O, and underabundant iron peak elements.
No magnetic field was detected in HD 139614 by Alecian et al. (2012a). Hubrig et al. (2004) claimed a detection of a magnetic field in this star, which was not supported by Wade et al. (2005). Further observations by Hubrig et al. (2009) did not detect a magnetic field in the star, though they comment on a possible signal in the Ca H and K lines. In their re-reduction of FORS data, Bagnulo et al. (2011) find no evidence for a magnetic field in HD 139614, and conclude that the reported detection by Hubrig et al. (2004) is likely spurious. Consequently, we consider HD 139614 not to have a significant magnetic field.
a.7 Hd 141569
We find K and for HD 141569, which is consistent with the parameters found by Guimarães et al. (2006). We find Boo peculiarities in this star, with solar He, C and N, possibly overabundant O, and underabundant iron peak elements.
a.8 Hd 142666
For HD 142666 we find K and , which is consistent with Guimarães et al. (2006). We find weak Boo peculiarities in this star. C, O, Na and S have solar abundance, while the iron peak elements are consistently 0.5 dex below solar, except for the rather uncertain Co abundance.
a.9 Hd 144432
For HD 144432 we find K and , which is consistent with Guimarães et al. (2006). The abundances of all elements we derive are within 2 of solar, and most are within 1. There are no clear trends in abundances of different elements, thus we consider the star chemically normal, with nearly solar metallicity.
Alecian et al. (2012a) observed HD 144432 twice and found no magnetic field. Wade et al. (2007) also found no magnetic field in the star. Hubrig et al. (2007) claim a weak magnetic field detection in this star in one observation, but Hubrig et al. (2004) and Hubrig et al. (2009) did not detect a significant magnetic field in the star. Bagnulo et al. (2011) re-reduced the one observation of HD 144432 with a magnetic field detection from Hubrig et al. (2007). Bagnulo et al. (2011) do find a magnetic field in the observation at a level ( G). However, they caution that the impact of instrumental effects on the detection of weak magnetic fields is not well understood for the FORS1 instrument. They conclude that for weak magnetic fields, detections at less than should not be considered definite, and require supporting observations. Thus, while we cannot completely exclude the possibility of a magnetic field in HD 144432, based on the large majority of the evidence (five out of six observations) we conclude that HD 144432 likely has no significant magnetic field.
a.10 Hd 163296
We determine an of K and a of , which is consistent with Guimarães et al. (2006). We find the abundances of most elements to be within uncertainty of solar, though Mg and Cr are both overabundant at slightly more than 2. With largely solar abundance and no clear selective depletion or enhancement, we conclude that the star is chemically normal.
a.11 Hd 169142
For HD 169142 we find K, which is consistent with Guimarães et al. (2006), and , which is 2 greater than Guimarães et al. (2006) who find . If we were to adopt the of Guimarães et al. (2006), it would decrease our abundances by between 0.05 and 0.2 dex depending on the element, which is generally within our uncertainties. We find a clear Boo pattern of abundances in this star. C and O are solar, and S is only slightly underabundant, while iron peak elements are underabundant by between 0.5 and 1 dex. Interestingly Ba and Nd are roughly solar, though Eu may be underabundant.
a.12 Hd 176386
We derive an of K and a of for HD 176386. The abundances we derive for the star are consistent with solar, with the possible exception of Ca. Thus we conclude that the star is chemically normal.
No magnetic field was found for HD 176386 by Alecian et al. (2012a). Hubrig et al. (2009) claim a marginal detection at 3 in one observation. Bagnulo et al. (2011) re-reduced this observation and found no magnetic field. Thus, based on the majority of the evidence, we consider HD 176386 to be a non-magnetic star.
a.13 Hd 179218
For HD 179218 we derived K and , both of which are consistent with Guimarães et al. (2006). We find a solar He abundance, abundances of C, O, and Na that are marginally enhanced relative to solar, and N and S abundances that are uncertain, but apparently enhanced relative to solar. Iron peak elements are depleted relative to solar by 0.5 dex. We conclude that HD 179218 displays Boo peculiarities. The overabundances of C, N, O, and S may reflect an above solar intrinsic metallicity in the star, which the Boo peculiarities are superimposed upon. However, this is not certain, as He appears to have a nearly solar abundance.
a.14 Hd 244604
We found K and for HD 244604, which is consistent with the parameters from Acke & Waelkens (2004). We find most elements have solar abundances or are enhanced by around 1. There is no evidence for selective depletion or enhancement. We conclude that the star is chemically normal. Acke & Waelkens (