Magnetic, Chemical and Rotational Properties of the Herbig Ae/Be Binary System HD 72106

Magnetic, Chemical and Rotational Properties of the Herbig Ae/Be Binary System HD 72106


Recently, strong, globally-ordered magnetic fields have been detected in some Herbig Ae and Be (HAeBe) stars, suggesting a possible evolutionary connection to main sequence magnetic chemically peculiar Ap and Bp stars. We have undertaken a detailed study of the binary system HD 72106, which contains a B9 magnetic primary and a HAeBe secondary, using the ESPaDOnS spectropolarimeter mounted on the CFHT. A careful analysis of the very young primary reveals that it has an approximately dipolar magnetic field geometry, strong chemical peculiarities, and strong surface chemical abundance inhomogeneities. Thus the primary is very similar to an Ap/Bp star despite having completed less then 1.5% of its main sequence life, and possibly still being on the pre-main sequence. In contrast, a similar analysis of the secondary reveals solar chemical abundances and no magnetic field.

Stars: abundances – Stars: magnetic fields – Stars: chemically peculiar – Stars: pre-main sequence – Stars: individual: HD 72106

123 2007 \volume1 \volnumber1 \firstpage1


C.P. Folsom, G.A. Wade, et al.

1 Introduction

Strong, globally-ordered magnetic fields have been known to exist in Ap/Bp stars for quite some time. However, the origin and evolution of these magnetic fields, as well as their impact on the observed peculiar chemical abundances, remains unclear. Recently, magnetic fields have been discovered in some Herbig Ae and Be (HAeBe) stars (Donati et al., 1997; Hubrig et al., 2004; Wade et al., 2005; Wade et al., 2007). HAeBe stars are pre-main sequence intermediate-mass stars and are the direct evolutionary progenitors of A and B stars.

The discovery of magnetic fields in HAeBe stars is of particular importance because it suggests a possible evolutionary link between Ap/Bp stars and magnetic HAeBe stars. The observed magnetic fields in HAeBe stars appear to have similar strengths and geometries to those observed in Ap/Bp stars (Alecian et al., P27, these proceedings). Furthermore, magnetic fields are detected in HAeBe stars with approximately the same frequency as in main sequence A and B stars, 5-10% (Wade et al., 2007).

These recent discoveries prompted us to perform a detailed investigation of one particularly interesting case: the very young binary system HD 72106.

2 Hd 72106

HD 72106 is a double star system with a 0.8\arcsec  separation between the components. The brighter component (HD 72106A) was observed to possess a magnetic field by Wade et al. (2005), and the secondary was identified as being a HAeBe star by Vieira et al. (2003).

The HD 72106 system was observed by Hipparcos, and solved as a binary, producing a “good quality” solution. Very large proper motions in both RA: mas yr, and in DEC: mas yr were observed, with no evidence for relative motion between the components (ESA, 1997). The Hipparcos parallax for the system places the two stars at a distance of pc. The radial velocities we determine from our observations of both components are identical at km s. Thus the stars are at approximately the same point in three dimensional space, and are moving together in three dimensions. This strongly suggesting that the system is in fact a true binary.

Based on both Balmer line and metallic line fitting, we find   = K and   = for the primary, and   = K and   = for the secondary.

The two stars can be placed on the H-R diagram, as seen in Fig. 1. When compared with the evolutionary tracks of Palla & Stahler (1993), we find that both stars sit near the zero-age main sequence (ZAMS). Based on the position of the secondary (which we infer to be a pre-main sequence star based on emission at H) we can constrain the age of the system to be between 3 and 13 Myr. At the younger end of this range, the primary would still have of its pre-main sequence lifetime remaining. However, at the older limit, the primary would already be on the main sequence, with a fractional main sequence age () of 0.015. Thus the evolutionary status of the primary is unclear, although the star is certainly very young.

Figure 1: The pre-main sequence H-R diagram of the HD 72106 system. Evolutionary track (solid lines) and isochrones (dashed lines), as well as the ZAMS and the birthline (for an accretion rate of M yr) from Palla & Stahler (1993) are shown.

In this investigation we used a set of 20 high-resolution Stokes and spectra obtained with ESPaDOnS, the high-resolution echelle spectropolarimeter mounted on the CFHT.

3 Magnetic Fields

We applied the Least Squares Deconvolution (LSD) technique (Donati et al., 1997) to our observations of both HD 72106A & B. We find clear evidence for a magnetic field in nearly all our Stokes profiles of HD 72106A, with maximum longitudinal field of  G, but we see no evidence for a magnetic field in any of our observations of HD 72106B (maximum longitudinal field of  G).

In order to determine the magnetic field geometry of HD 72106, we need to know the rotation period of the star. We searched for periodicity in our longitudinal magnetic field observations using a modified Lomb-Scargle method. We also applied this technique to both our Stokes and LSD profiles themselves; searching each pixel in the line profile for periodic variability, and then averaging the results over all pixels. Additionally, the phasing of LSD profiles produced by candidate periods was investigated by eye, looking for variability physically consistent with rotation. The one period meeting all these criteria, days, was adopted as the rotation period of the star.

Using this rotation period, together with the inferred stellar radius and , we determined the inclination of the rotation axis to the line-of-sight of HD 72106A to be . We modeled the magnetic field of HD 72106A by fitting the observed Stokes LSD profiles with synthetic profiles computed by ZEEMAN2 (Landstreet, 1988; Wade et al., 2001). The mean Stokes LSD profile was fit, then the phased set of Stokes LSD profiles was fit using a dipole field model. From this procedure we derive the dipole polar strength G and obliquity , with error bars. An analysis of the longitudinal field measurements shows that they are also consistent with this model.

Figure 2: Longitudinal magnetic field measurements, phased according to the rotation period of days, together with the best fit sinusoid for the data.

4 Abundance Analysis

The spectra of both components of HD 72106 were fit, using synthetic spectra calculated with the ZEEMAN2 spectrum synthesis code. Spectra were fit with a Levenberg-Marquardt minimization method, and the results were extensively verified by eye. To approximate our magnetic observations, a 1 kG dipole magnetic field was included in the model of HD 72106A, and no field was used for HD 72106B. The final best fit values, reported in Fig. 3, are averages over seven spectral windows between 100 and 200 Å  in length, and the uncertainties are the standard deviations. When an abundance was based on only a few windows, a conservative estimate of the uncertainty was made by eye.

For HD 72106A we find remarkably strong over-abundances of Cr, Fe and Nd. Additionally, Si is overabundant by dex while He appears to be underabundant by dex. We find   = 41.0 0.3 km s  for this star. This pattern of strong peculiarities is distinctive of Bp stars.

In contrast, for HD 72106B we find abundances that are consistent with solar, at least at the 2 level. Small departures from solar abundance for C and Sc may exist, however they are consistent with the normal scatter for an A star. We find   = km s  and a microturbulence of km s.

Figure 3: Mean best fit abundances, relative to solar, for both HD 72106A & B. Points marked with only an arrow indicate the value is an upper limit only. Solar abundances are from Grevesse et al.(2005).

5 Doppler Imaging

We performed Doppler Imaging of HD 72106A using the INVERS12 code (Kochukhov et al., 2004). LSD profiles were used for the Doppler Imaging, rather then individual lines, due to the S/N of our spectra. LSD profiles for specific elements were constructed using line masks containing only lines of one element.

Surface abundance maps were calculated for Si, Ti, Cr, and Fe, using a inclination angle, and are shown in Fig. 4. Strongly inhomogeneous distributions can be seen for all four elements. Ti, Cr, and Si have similar distributions. The distribution of Ti appears to differ slightly from that of Cr and Fe in the southern hemisphere, but the inversion is poorly constrained in that region. The large spot seen in the northern hemisphere near phase 0 in all four images is close to the location of the positive magnetic pole.

Figure 4: Surface abundance maps of HD 72106A for Si, Ti, Cr and Fe.

6 Conclusions

HD 72106A is in many respects very much like an Ap/Bp star. It possess a predominantly dipolar magnetic field, strong chemical peculiarities, and strong surface abundance inhomogeneities. Despite this, the star is very young, having completed at most 1.5 % of its main sequence life and possibly still being on the pre-main sequence. In contrast, HD 72106B is chemically normal, and has no detectable magnetic field.

The HD 72106 system provides evidence that chemical peculiarities can form at the ZAMS, and possibly even on the pre-main sequence. Additionally, it provides support for an evolutionary link between magnetic HAeBe stars and Ap/Bp stars. Finally this system raises questions about when in a star’s evolution chemical peculiar first appear, in particular whether peculiarities can arise during the pre-main sequence or only early in the main sequence.

EA is funded by the Marie-Curie FP6 programme, while JDL and GAW acknowledge support from NSERC and DND-ARP (Canada).


  1. Donati, J.-F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A., 1997, MNRAS, 291, 658
  2. ESA 1997, The Hipparcos and Tycho Catalogues, ESA SP-1200
  3. Grevesse, N., Asplund, M., & Sauval, A. J., 2005, in the proceedings of ‘Element Stratification in Stars: 40 Years of Atomic Diffusion’, EAS 17, 21
  4. Hubrig, S., Schöller, M., & Yudin, R. V., 2004, A&A, 428, 1
  5. Kochukhov, O., Drake, N. A., Piskunov, N., & de la Reza, R., 2004, A&A, 424, 935
  6. Landstreet, J. D., 1988, ApJ, 326, 967
  7. Palla, F. & Stahler, S. W., 1993, ApJ 418, 414
  8. 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
  9. Wade, G. A., Bagnulo, S., Drouin, D., Landstreet, J. D., Monin, D., 2007, MNRAS, 376, 1145
  10. Wade, G. A,. Bagnulo, S., Kochukov, O., Landstreet, J. D., Piskunov, N., & Stift, M. J., 2001, A&A, 374, 265
  11. 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, 31
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description