Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://star.arm.ac.uk/preprints/2007/514.pdf Äàòà èçìåíåíèÿ: Mon Jan 7 14:02:12 2008 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 02:48:49 2012 Êîäèðîâêà: ISO8859-5 Ïîèñêîâûå ñëîâà: arp 220

Mon. Not. R. Astron. Soc. 000, 1-14 (2002)

Printed 7 December 2007

A (MN L TEX style file v2.2)

C h a r a c t e r is a t io n o f t h e m a g n e t ic fi e ld o f t h e H e r b ig B e s t a r H D 2 0 0 7 7 5
E. Alecian1 2 , C. Catala2 , G.A. Wade1 , J.-F. Donati3 , P. Petit3 , J.D. Landstreet4, ? T. Bohm3, J.-C. Bouret5, S. Bagnulo6, C. Folsom6 , J. Grunhut 1 and J. Silvester 1
2 3 4 5 6

1

Dept. of Physics, Royal Military College of Canada, PO Box 17000, Stn Forces, Kingston K7K 7B4, Canada Observatoire de Paris, LESIA, 5, place Jules Janssen, F-92195 Meudon Principal CEDEX, France Laboratoire d'Astrophysique, Observatoire Midi-Pyrenees, 14 avenue Edouard Belin, F-31400 Toulouse, France ÄÄ Dept. of Physics & Astronomy, University of Western Ontario, London N6A 3K7, Canada Laboratoire d'Astrophysique de Marseille, Traverse du Siphon, BP8-13376 Marseille Cedex 12, France Armagh Observatory, College Hill, Armagh BT61 9DG, Northern Ireland, UK

Accepted . Received ; in original form

ABSTRACT

The origin of the magnetic fields observed in some intermediate mass and high mass main sequence stars is still a matter of vigorous debate. The favoured hypothesis is a fossil field origin, in which the observed fields are the condensed remnants of magnetic fields present in the original molecular cloud from which the stars formed. According to this theory a few percent of the PMS Herbig Ae/Be star should be magnetic with a magnetic topology similar to that of main sequence intermediate-mass stars. After our recent discovery of four magnetic Herbig stars, we have decided to study in detail one of them, HD 200775, to determine if its magnetic topology is similar to that of the main sequence magnetic stars. With this aim, we monitored this star in Stokes I and V over more than two years, using the new spectropolarimeters ESPaDOnS at CFHT, and Narval at TBL. By analysing the intensity spectrum we find that HD 200775 is a double-lined spectroscopic binary system, whose secondary seems similar, in temperature, to the primary. We have carefully compared the observed spectrum to a synthetic one, and we found no evidence of abundance anomalies in its spectrum. We infer the luminosity ratio of the components from the Stokes I profiles. Then, using the temperature and luminosity of HD 200775 found in the literature, we estimate the age, the mass and the radius of both components from their HR diagram positions. From our measurements of the radial velocities of both stars we determine the ephemeris and the orbital parameters of the system. A Stokes V Zeeman signature is clearly visible in most of the Least Square Deconvolution profiles and varies on a timescale on the order of one day. We have fitted the 30 profiles simultaneously, using a 2 minimisation method, with a centered and a decentered-dipole model. The best-fit model is obtained with a reduced 2 = 1.0 and provides a rotation period of 4.3281Á0.0010 d, an inclination angle of 60Á11, and a magnetic obliquity angle = 125Á8. The polar strength of the magnetic dipole field is 1000 Á 150 G, which is decentered by 0.05 Á 0.04 R from the center of the star. The derived magnetic field model is qualitatively identical to those commonly observed in the Ap/Bp stars. Our determination of the inclination of the rotation axis leads to a radius of the primary which is smaller than that derived from the HR diagram position. This can be explained by a larger intrinsic luminosity of the secondary relative to the primary, due to a larger circumstellar extinction of the secondary relative to the primary. Key words: Stars : magnetic field - Stars: pre-main-sequence - Stars: binaries: spectroscopic - Star: individual: HD 200775 - Instrumentation : spectropolarimetry.



Based on observations obtained at the Canada-France-Hawaii Telescope

(CFHT) which is operated by the National Research Council of Canada,


2

E. Alecian et al.

Table 1. Log of the observations. Columns 1 and 2 give the UT date and the Heliocentric Julian Date of the observations. Column 3 gives the total exposure time. Column 4 gives the peak signal to noise ratio (at 730 nm per 1.8 km/s spectral pixel) in the spectra and column 5 gives the signal to noise ratio in the LSD Stokes V profiles. Column 6 gives the longitudinal magnetic field, column 7 gives the rotation phase derived in Section 6.2, and columns 8 and 9 give the radial velocities of both components of the system. Column 9 gives the instrument used. Date UT Time 24/09/04 22/05/05 23/05/05 24/05/05 24/05/05 25/05/05 25/05/05 20/06/05 22/06/05 25/06/05 19/07/05 20/07/05 26/08/05 09/06/06 10/06/06 11/06/06 12/06/06 13/06/06 13/06/06 14/06/06 15/06/06 16/06/06 10/12/06 10/12/06 11/12/06 11/12/06 04/03/07 24/04/07 26/04/07 26/06/07 10:33 14:35 12:22 12:09 13:11 12:15 14:02 13:23 13:43 13:42 11:43 10:23 10:00 10:39 14:31 14:03 10:30 10:08 10:52 15:18 15:13 15:14 4:27 6:16 5:09 5:32 16:00 4:09 4:09 12:58 HJ D (2 450 000+) 3272. 3513. 3514. 3515. 3515. 3516. 3516. 3542. 3544. 3547. 3570. 3571. 3608. 3895. 3897. 3898. 3898. 3899. 3899. 3901. 3902. 3903. 4079. 4079. 4080. 4080. 4164. 4214. 4216. 4278. 9398 1072 0150 0065 0490 0101 0511 0596 0732 0730 9915 9358 9206 9447 1059 0867 9386 9231 9538 1387 1352 1365 6850 7610 7141 7302 1633 6708 6713 0424 texp (s) 1200 3600 3600 3600 2400 3600 2400 800 800 800 1200 1200 1600 1200 2160 2160 2160 2400 2400 800 1840 1720 1600 1600 1200 1200 1320 2700 2500 2400 S/N 480 540 420 560 470 520 300 260 240 280 370 530 610 550 780 710 760 790 810 540 760 740 330 220 470 470 600 410 180 620 S/N (LSD) 1810 1940 1470 2090 1690 1880 1660 840 540 920 1330 1830 2280 1970 2840 2590 2720 2780 2930 2000 2820 2700 1150 770 1860 1800 2130 1460 660 2280 B (G) -299Á89 91Á84 73Á104 -262Á78 -376Á98 -411Á78 -306Á88 -583Á180 229Á199 -378Á188 37Á118 -293Á90 156Á73 -404Á82 -383Á57 -24Á69 161Á63 -193Á60 -97Á59 -269Á79 -146Á58 180Á64 -93Á132 -198Á188 320Á83 34Á85 -198Á58 17Á95 -341Á180 -233Á52 phase 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 90 39 60 83 84 06 07 08 54 24 76 98 53 84 11 33 53 76 77 04 27 50 29 31 53 53 81 48 94 12 vradA (km.s-1 ) 8.2Á0. -19.3Á1. -19.8Á1. -19.9Á1. -20.0Á1. -19.9Á1. -20.3Á1. -21.3Á0. -21.5Á0. -21.2Á1. -22.6Á0. -21.6Á1. -23.3Á1. -20.2Á1. -19.8Á1. -20.4Á1. -20.2Á0. -20.1Á0. -20.2Á0. -20.4Á1. -19.9Á1. -21.5Á1. -14.2Á0. -14.0Á1. -13.7Á0. -13.9Á0. -9.9Á0. -7.3Á2. -7.5Á2. -1.8Á1. 5 1 0 0 0 0 0 9 9 0 9 9 0 0 0 0 9 9 8 1 1 3 9 0 9 9 8 1 7 0 vradB (km.s-1 ) -21. 1. 2. 0. 0. 0. -0. 2. 1. 3. 4. 2. 5. 1. 0. 4. 3. 1. 0. 2. 8. 9. -0. -2. -2. -3. -6. -10. -9. -9. 1Á1.0 2Á3.6 9Á3.5 0Á3.4 8Á3.4 9Á3.3 9Á3.4 3Á3.1 5Á3.0 0Á3.4 7Á3.2 9Á6.8 0Á3.4 9Á3.2 0Á3.2 4Á3.2 2Á3.2 1Á2.9 7Á2.8 6Á3.8 6Á4.1 3Á4.7 5Á3.3 7Á3.5 9Á3.4 4Á3.2 1Á3.0 2Á8.2 3Á10.0 1Á3.4 Instrument ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD ESPaD Narval Narval ESPaD O O O O O O O O O O O O O O O O O O O O O O O O O O O nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS nS

OnS

1 INTRODUCTION Some main sequence A, B and O stars host strong ( kG) organised magnetic fields. The origin of the magnetic fields of these intermediate and high mass stars is still a matter of debate. The favoured theory is the fossil field hypothesis. This theory assumes that the magnetic fields observed in these main sequence stars are relics of the magnetic fields which existed in the molecular clouds from which the stars formed. This theory implies that the remnant interstellar magnetic field should subsist throughout all the processes of formation encountered by the star, from the gravitational collapse to the main sequence phase, without being regenerated. According to this theory, some pre-main sequence (PMS) stars of intermediate mass should host magnetic fields. For a long time, this theory conflicted with the general belief that all stars pass through a completely convective phase during the PMS phase. The turbulent diffusion produced by the convection would dissipate the magnetic field during the Hayashi phase. However Palla & Stahler (1993) calculated the birthline (the locus in the HR diagram where stars become optically visible and start the PMS phase), showing that the Hayashi phase is considerably reduced for

stars between 1.5 M and 2 M and disappears completely for stars above 2 M . A fundamental conclusion of this work was that the magnetic field of the intermediate PMS stars can potentially survive the PMS phase. According to the fossil field hypothesis, some Herbig Ae/Be stars (PMS stars of intermediate mass) should be magnetic. Many authors have tried to detect a magnetic field in these stars. Catala et al. (1993) observed AB Aur using a Zeeman polarimeter placed before the Coude spectrograph of the CFHT. They found no circuÄ lar polarisation signal in the FeII 5018 Õ line and they obtained an upper limit around 1 kG. Catala et al. (1999) tried again to detect a magnetic field in AB Aur using the MUSICOS instrument temporarily installed at the CFHT. They observed no Zeeman signature and obtained a lower detection limit around 300 G. Donati et al. (1997) observed a large sample of cool and hot stars, including two Herbig stars, using the UCLES spectrograph of the AAT and the SemelPol polarimeter. Applying the Least Squares Deconvolution method, they detected a magnetic field of around 50 G in the Herbig Ae star HD 104237, but they obtained no detection in another Herbig Ae/Be star, HD 100546. Hubrig et al. (2004, 2007) claim detections in HD 139614 and HD 144432. Finally recent search for magnetic field in the field Herbig Ae/Be stars by Wade et al. (2007) has been carried out using the FORS1 spectropolarimeter at the ESO VLT. They identified two possible magnetic stars (HD 101412

the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii E-mail: evelyne.alecian@rmc.ca


Characterisation of the magnetic field of the Herbig Be star HD 200775
and BF Ori), of which one has been confirmed (HD 101412), with higher resolution data (Wade et al. 2007 in preparation). Recently, the new generation instrument ESPaDOnS was installed at the Cassegrain focus of the CFHT (Donati et al., in preparation). Thanks to the high efficiency and qualities of this high resolution spectropolarimeter, magnetic fields in several Herbig stars have been discovered. During the technical run of ESPaDOnS a magnetic field was discovered in the Herbig Ae/Be star HD 200775 (Donati et al., in preparation). Then, during scientific nights of the first semester of ESPaDOnS, three other stars were detected as magnetic : HD 72106, V380 Ori and HD 190073 (Wade et al. 2005; Catala et al. 2006). These recent discoveries bring new arguments to support the fossil field theory. However further exploration is necessary. According to the fossil field hypothesis, the structure and the intensity of the magnetic field of these PMS stars should be similar to those of the main sequence magnetic stars. In order to verify this, we have observed HD 200775, with ESPaDOnS, during many nights. Using the temporal variations of the Stokes V profile we were able to model the structure and the geometry of the magnetic field. The next section presents the observations and data reduction procedures. In Sect. 3, 4 and 5 we present the study of the intensity spectrum and in Section 6 the study of the polarization spectrum. Conclusions are given in Sect. 7.

3

tion between a profile (which is the same for all lines) and a mask including all choosen lines of the spectrum. We therefore apply a deconvolution to the observed spectrum using the mask, in order to obtain the average photospheric profiles of Stokes I and V . In this procedure, each line is weighted by its signal to noise ratio, its depth in the unbroadened model and its Lande factor. The mask Ä was first computed using Kurucz ATLAS 9 models (Kurucz 1993) with T eff = 19000 K and log g = 3.5 suitable for the star (Table 3). We excluded from this mask hydrogen Balmer lines, strong resonance lines and lines whose Lande factor is unknown. Then we Ä cleaned the mask by selecting the lines with a depth less than 0.4, in order to eliminate the lines contaminated by strong emission. We have also modified the line depths in order to take into account the relative depth of the lines of the observed spectrum. The final mask contains only 37 lines, but is sufficent to compute the LSD Stokes I and V profiles with a good S/N ratio (Table 1). To explore the sensitivity of the LSD profiles to the detailed line list, line masks were constructed using weak lines including He lines, and using weak lines excluding He lines. The various LSD profiles were analysed in the manner described in Section 3. No significant differences were found. The null N profile has been computed in the same way, from the null spectrum, and will be used, in Sect. 6, to check that the signature in the Stokes V profile is real. The LSD average line profiles were computed on a velocity grid with 3.6 km.s-1 sampling. The resulting relative noise in the LSD Stokes V profiles is given in the 5th column of Table 1.

2 OBSERVATIONS AND DATA REDUCTION Our data were obtained using the high resolution spectropolarimeter ESPaDOnS installed on the 3.6 m Canada-France-Hawaii Telescope (Donati et al., in preparation) during many scientific runs. Table 1 presents the log of the observations. One of our spectra has also been obtained using the instrument Narval, which is a copy of ESPaDOnS, installed on the 2 m Bernard Lyot Telescope (TBL) at the Pic du Midi observatory in France. We used the instruments ESPaDOnS and Narval in polarimetric mode, and we obtained spectra with resolving power of 65000. Each exposure was divided in 4 sub-exposures of equal time in order to compute the optimal extraction of the polarisation spectra (Donati et al. 1997, Donati et al., in preparation). We recorded only circular polarisation, as the Zeeman signature expected in linear polarisation is about one order of magnitude lower than circular polarisation. The data were reduced using the "Libre ESpRIT" package especially developed for ESPaDOnS and Narval, and installed at the CFHT and at the TBL (Donati et al. 1997, Donati et al., in preparation). After reduction, we obtained the intensity Stokes I and the circular polarisation Stokes V spectra of the star observed, both normalised to the continuum intensity of HD 200775. A null spectrum (N ) is also computed in order to diagnose spurious polarisation signatures, and to help to verify that the signatures in the Stokes V spectrum are real (Donati et al. 1997). The data of May 2005 were affected by a 1.3 mag loss compared to the data obtained more recently with ESPaDOnS. This problem, which was due to damage to the external jacket of the optical fibres, was fixed prior to the July run. In order to compensate for this damage and to obtain a satisfactory S/N ratio, we exposed longer (see Table 1). We then applied the Least Squares Deconvolution procedure to all spectra (Donati et al. 1997). This method assumes that all selected lines of the intensity spectrum have a profile of similar shape. Hence, this supposes that all lines are broadened in the same way. We can therefore consider that the observed spectrum is a convolu3 BINARITY In the litterature, many authors refer to the binarity of HD 200775. Corporon & Lagrange (1999) detected radial velocity variations of photospheric lines of the spectrum of HD 200775, but they were not able to determine the orbital period. Miroshnichenko et al. (1998) analysed spectroscopic data of this star obtained by many authors over 20 years. They noted a cyclic variation of the equivalent width of H with a 1345 day period, and they suggested that a companion star could be the trigger of the H activity. Pogodin et al. (2004) discussed again this theory using new data and observational evidence of radial velocity variations of some photospheric lines. They plotted two radial velocity curves, obtained on the one hand from the H wings, and on the other hand from photospheric lines. They fitted them with a synthetic radial velocity curve in order to obtain the ephemeris and the orbital parameters of the double system. Finally Monnier et al. (2006) were able to separate both stars of HD 200775 using interferometry, and determine orbital parameters for the system. Figure 1 shows the LSD average Stokes I profiles (thin full line) of all our observations, obtained using the cleaned mask as explained in Section 2. We can see the presence of a second component in each profile : on the blue for the Sept. 2004 profile and on the red for the others. Furthermore, Fig. 2 shows some individual absorption lines of the intensity spectrum, observed in Sept. 2004 and May 2005, where the secondary component is visible. Between Sept. 2004 and May 2005 the primary component of the profile and the photospheric lines shifted towards the blue, while the secondary component moved in the other direction. We therefore assign this second component of the average profile to a companion star. HD 200775 is therefore an SB2 system whose orbital period is likely greater than one year. In the following we call HD 200775A (with the shaper, deeper lines) the primary component and HD 200775B (with the broader, shallower lines) the secondary component.


4

E. Alecian et al.

Figure 1. LSD Stokes I profiles of all spectra (thin full line) superimposed with their fit (thick dashed-line) described in Section 3

Figure 2. Photospheric lines of SiII 4128Õ, SII 4153Õ, SII 5032Õ, and SII 5454Õ, oberved on September 23rd 2004 (redshifted one) and on August 25th 2005 (blueshifted one). The fit (thick dashed line) with a convolution of a gaussian and a rotation profile is superimposed to the observed profiles (thin full line).

In order to study the variations of the radial velocities of our spectra, we performed a simultaneous Least-Squares fit to the 30 LSD I profiles1 . Each profile is fitted with the sum of two func1

tions, each function modeling the line profile of one components. Each one of these two functions is the convolution of a rotation function (for which the projected rotational velocity v sin i is a free parameter in the fit) and a gaussian whose width is fixed and computed from the spectral resolution and the inferred macroturbulent velocity (Gray 1992). We adopted an isotropic macroturbulent velocity of 15 km.s-1 in order to fit the wings of the profiles of individual lines in the spectrum (both strong and weak lines), and of the LSD I profiles. The free parameters of the fitting procedure are the centroids, depths, and projected rotational velocities (v sin i) of both components. The centroids of both functions can vary from one profile to another., whereas the depths and v sin i of both components cannot. This fitting procedure therefore assumes that the depths and v sin i of both components do not vary with time, which we confirm (within the error bars) by fitting each profile separately. Figure 3 shows an example of such a fit for the Sept. 2004 and 22 May 2005 profiles, as well as the profiles of the isolated components calculated by subtracting the synthetic profile of the other component from the observed LSD Stokes I profile. This automatic fitting procedure enables us to measure the v sin i and radial velocities of both components. We obtain projected rotational velocities of 26Á 2 km.s-1 and 59 Á 5 km.s-1 for the primary and the secondary components, respectively. Assuming an orbital period greater than one year (consistent with our data) and considering a projected rotational velocity of the primary of 26 km.s-1 , we conclude that the system is not synchronised. The inferred macroturbulent velocity is not determined very precisely: acceptable values range from about 10-30 km.s-1 . However, the fitting procedure is largely insensitive to these uncertainties - changing the macroturbulence within the uncertainties results

Although a potentially more effective approach would be to perform spectral disentangling, our data set was not suitable for such an analysis due to

the sparse coverage of the orbital period and the small velocity separation of the lines of the components (O. Kochukhov, private communication).


Characterisation of the magnetic field of the Herbig Be star HD 200775

5

Figure 3. Upper panels: LSD Stokes I profiles from spectra obtained on September 23, 2004 and September 21, 2005 (full line). The fitted line profiles are the superimposed dashed lines. Lower panels: Individual intensity profiles of the primary (narrow red profile) and secondary (large green profile) components of the system.

Figure 4. Upper panel : Equivalent width of H in function of time. H appears in two different orders in the spectrum. Therefore the values measured on both orders are plotted. Lower panel : Radial velocities of the primary (red filled circle) and the secondary (green filled diamond) components of HD 200775. The full lines are the fitting curves. The arrows indicate the periastron passage.

in only very small changes to the inferred v sin is, and to no significant change in the radial velocities. We observe a strong variation of the radial velocity of the photospheric lines of the primary star between Sept. 2004 and May 2005 from +7 km.s-1 to -20 km.s-1 . Then it is slightly varying, and since August 2005 it has been increasing slowly. In Fig. 4 (lower panel), we plot the radial velocity of both components as a function of time, fitted with the radial velocity curves of an eccentric binary system, using a 2 minimization method. We obtain the following parameters for the system: orbital period P = 1412 Á 54 d, periastron epoch T 0 = 2448991 Á 152 d, systemic radial velocity = -7.9 Á 0.9 km.s-1 , eccentricity e = 0.32 Á 0.06, periastron longitude = 216 Á 12 , and radial velocity amplitude of both components KP = 20.9 Á 2.5 km.s-1 and KS = 17.0 Á 2.5 km.s-1 , MP S leading to the mass ratio of the system q = MS = KP = 0.81 Á 0.22. K The orbital parameters that we found are similar to those obtained by Pogodin et al. (2004) from the bisector velocities of H. In order to compare our work to theirs, we also measured the bisector velocities of the H line, as well as the other emission lines of the intensity spectrum. We observe that all of them closely follow the radial velocity of the secondary star (see Fig. 5). We therefore conclude that the region from which these emission lines originate is linked to the secondary component of the system. These results are in contradiction with those of Pogodin et al. who found that the bisector velocity follows the radial velocity of the primary star. Whether this contradiction is related to a real change of behavior of the binary system and its environment, or to problems of measurement of bisector and/or photospheric radial velocities in the data of Pogodin et al., remains to be determined. Assuming that the variations of the bisector velocities of H originate from the orbital movement of the secondary, we included them into the fit of the radial velocities. We found no significant differences of the inferred orbital parameters compared to those obtained from the fit of solely the radial velocities.

Figure 5. Radial velocities of the primary (red filled circle) and the secondary components (green filled diamond) of HD 200775 superimposed with the bisector velocities of the H line (black open square), measured at a level of 1.5 - 2.0Fc . Note that the bisector velocities trace the radial velocities of the secondary, and clearly do not follow those of the primary.

Using our determination of the orbital parameters we can derive the semi-major axis of the orbit, and the masses of the primary and the secondary, as functions of the inclination of the orbit: a sin iorb , MP sin3 iorb , and MS sin3 iorb (see Table 2). Using our determination of MP and MS (Table 3) from stellar evolutionary models, we can therefore estimate the orbital inclination of the system. We find sin iorb (P) = 0.70 Á 0.16 from MP and MP sin3 iorb , and sin iorb (S) = 0.78 Á 0.17 from MS and MS sin3 iorb . We determine the final value of sin iorb by taking a weighted avarage of sin iorb (P) and sin iorb (S). We obtain sin iorb = 0.74 Á 0.17, leading to iorb = 48+17 . - 13 Then we use this value to derive the semi-major axis of the orbit : a = 6.7 Á 1.9 AU. Using the Hipparcos parallax of HD 200775


6

E. Alecian et al.

Table 2. Orbital parameters of the system found in this work, compared to previous work (Pogodin et al. 2004; Monnier et al. 2006). Parameter Radial Velocity This work Results of RV fit 1412 Á 54 2448991 Á 152 -7.9 Á 0.9 0.32 Á 0.06 216 Á 12 20.9 Á 2.5 17.0 Á 2.5 Derived parameters 4.9 Á 0.8 2.73 Á 0.46 2.21 Á 0.38 3.6 Á 1.6 4.4 Á 1.9 48+17 - 13 6.7 Á 1.9 16 Á 9 Bisector velocity (Pogodin et al. 2004) 1341 Á 23 2449149 Á 87 5.0 Á 0.6 0.29 Á 0.07 203 Á 22 11.2 Á 0.7 Interferometry (Monnier et al. 2006) 1377 Á 25 2449152 Á 90 0.30 Á 0.06 224 Á 16

Period (days) T 0 (HDJ) (km/s) e ( ) KP (km/s) KS (km/s)

a sin iorb (AU) aP sin iorb (AU) aS sin iorb (AU) MP sin3 iorb ( M ) MS sin3 iorb ( M ) iorb ( ) a (AU) a (mas)

1.38 Á 0.11

65 Á 8 15.14 Á 0.70

( = 2.33 Á 0.62 mas), this corresponds to a projected separation of the components a = 16 Á 9 mas. The values of the orbital parameters that we derived are summarised in Table 2 and compared to the orbital elements of Pogodin et al. (2004), as well as to the orbital parameters determined from interferometric data by Monnier et al. (2006). Most of our measurements of the orbital parameters agree well with both works. Our determination of KS is larger than that of Pogodin et al. which is likely due to our lack of data during the minimum separation of both stars. We need to continue to observe the system in order to confirm that value.

mainly due to its higher rotational velocity (v sin i = 59 km.s-1 , Section 3) and slightly due to a fainter luminosity. Furthermore, all lines observed in the spectrum of the secondary are also observed in the spectrum of the primary star. For all these reasons, we suspect that the temperature of the secondary star is similar to that of the primary component. Hernandez et al. (2004) determined the spectral type of Ä HD 200775 mainly by comparing the strength of atomic absorption lines to those of standard stars and found log(T eff ) = 4.27. We compared our spectra with a synthetic spectrum of a double star whose components have the same effective temperature and surface gravity : T eff = 19000 K and log g = 3.5. First, we calculated separate spectra of both components using TLUSTY non-LTE atmosphere models and the SYNSPEC code (Hubeny 1988; Hubeny & Lanz 1992, 1995). Then we computed the spectrum of the double star, with the BINMAG1 program (O. Kochukhov, private communication), using a macroturbulent velocity of 15 km/s, our measured v sin i, and radial velocities, and our estimation of the luminosity ratio LS /LP = 0.67 (see below). We conclude that the temperatures of the synthetic spectrum are in good agreement with the observed spectrum of HD 200775. Figure 6 shows a portion of the observed spectrum superimposed on the synthetic one. In addition, we also compared our spectra with a synthetic spectrum of a single star of T eff = 19000 K and log g = 3.5. We found that the synthetic spectrum of a double star fits better the depth and the width of the lines than the spectrum of a single star. That adds additional support to our hypothesis that both components of the system have similar temperature. In order to quantify the temperature range of the stars, we have calculated synthetic spectra by varying the effective temperature of both stars around 19000K. We found that these spectra are able to fit our observations within Á2000 K around 19000K. The luminosity of Hernandez et al. (2004) is the sum of the Ä luminosity of the primary (LP ) and the secondary (LS ) components of the system. We can estimate the luminosity ratio rL = LS /LP from the Stokes I profiles (normalised to the continuum of the binary). Section 3 describes the fitting procedure of these profiles. The result of the fitting procedure is the function:

4 FUNDAMENTAL PARAMETERS OF THE PRIMARY AND THE SECONDARY We used the fundamental parameters of HD 200775 of Hernandez Ä et al. (2004). These authors studied the spectra of 75 stars, most of which are classified as Herbig Ae/Be stars. They determined the spectral type of each star as well as the luminosity with particular attention to the total-to-selective extinction (RV ) used. In the case of classical main sequence stars the interstellar extinction law gives RV = 3.1. However, the circumstellar matter of Herbig stars may be dominated by dust and can lead to a higher reddening of the star. Using U BV R photometric data of Herbst & Shevchenko (1999), Hernandez et al. (2004) found that RV = 5 fits better most stars Ä of their sample than the standard value RV = 3.1. The knowledge of RV is fundamental to determine the luminosity: using RV = 3.1, they found log(L/L ) = 3.73 Á 0.27, while with RV = 5, they found log(L/L ) = 4.17 Á 0.27. For the purposes of this analysis, we adopt the most likely value of RV = 5, as HD 200775, illuminating the reflection nebulae NGC 7023, is still largely surrounded by dust and gas (Fuente et al. 1998), which certainly have an impact on the extinction law for this star. The observed luminosity, derived by Hernandez et al. (2004), Ä corresponds the total of the luminosities of both stars. A priori, we do not know the luminosity ratio of the components. However, we can see in the spectrum that the lines of the secondary component are less deep than those of the primary. We suppose that this is


Characterisation of the magnetic field of the Herbig Be star HD 200775

7

Figure 6. ESPaDOnS spectrum (full line) of HD 200775 between 438 and 450 nm superimposed with the synthetic spectrum of a double star (dashed-line) calculated using TLUSTY non-LTE atmosphere models and the SYNSPEC code (Hubeny 1988; Hubeny & Lanz 1992, 1995).

Table 3. Fundamental parameters of HD 200775A and HD 200775B. Column 7 gives the age of the system. Star HD 200775A HD 200775B T eff (K) 18600 Á 2000 18600 Á 2000 log(L/L ) 3.95 Á 0.30 3.77 Á 0.30 M/M 10.7 Á 2.5 9.3 Á 2.1 R/R


log g (cm.s-2 ) 3.4 Á 1.2 3.6 Á 1.2

age (Myr) 0.1 Á 0.05

v sin i (km.s-1 ) 26 Á 2 59 Á 5

10.4 Á 4.9 8.3 Á 3.9

I = 1 - fP - fS = 1 -

1 rL fP - fS , 1 + rL 1 + rL

(1)

where fP and fS are the separate Stokes I profile shapes of the pri mary and the secondary stars. The ratio fS dv/ fP dv is therefore equal to rL WS /WP , where WP and WS are the equivalent widths of the separate Stokes I profiles of both stars. These profiles have been computed using the same mask, and therefore the same lines. Assuming a common origin for both stars of the system and that they are not peculiar, they should have approximately the same chemical composition. The intrinsic equivalent widths should therefore be identical and the ratio fS dv/ fP dv be equal to rL . We computed this ratio for all our observations. The mean and the standard deviation of these 30 values gives a luminosity ratio LS /LP equal to 0.67 Á 0.05. We therefore derive the luminosity of the stars using : log LP L LS log L = = L - log(1 + rL ) = 3.95 Á 0.30 L L 1 log - log(1 + ) = 3.77 Á 0.30 L rL log (2) (3)

where, L = LP + LS is the observed luminosity of HD 200775 (Hernandez et al. 2004). Ä We used the derived luminosities and the temperature of HD 200775 to place both stars in the HR diagram and to compare their

position with evolutionary tracks of different masses calculated with the CESAM stellar evolutionary code (Morel 1997). In this way we obtained the mass and the radius of both stars, as well as the age of the system, assuming a common origin for both stars (the age of the system is the intersection of both age ranges determined for each star separatly) (Fig. 7). The value of the fundamental parameters are summarized in Table 3. The mass ratio is therefore q = 1.1 Á 0.5, which is consistent with the value determined from the orbit analysis (q = 0.81 Á 0.22, Sec. 3). Note that in both cases the error bars on the mass ratio are very large. This is partly due to our poor coverage of the orbital phase : our data cover only the half of the orbital period. In particular we do not have data where the maximum radial velocities are predicted. Therefore we need to get more data during the next years to improve it. The large error bars of the orbital parameters are also due to the non-unique solution of the fit of the Stokes I profiles, leading to very large error bars on the radial velocities and the luminosity ratio. To better constrain the mass of the system, photometric observations of both separate components would be necessary. It would provide us the luminosity of each star and therefore we would no longer be dependent on our estimation of the luminosity ratio. Furthermore, we note that using the orbit analysis we find that the secondary is more massive than the primary, whereas in this


8

E. Alecian et al.

Figure 7. Evolutionnary tracks (full lines) of different mass : from 7 M to 16 M , plotted in a HR diagram. The dashed lines are the 0.01, 0.025, 0.04 Myr, 0.065 Myr, 0.09 Myr, 0.115 Myr, 0.14 Myr, 0.17 and 0.2 Myr isochrones. Crosses represent the error bars in temperature and luminosity of HD 200775A (red square) and HD 200775B (green triangle).

Figure 8. Stokes I (bottom) and V (top) mean profiles of HD 200775 observed on Sept. 23rd 2004 (full line) and May 23rd 2005 (dashed line). The null profile N is plotted in the middle, and the 3 error bars of V and N are plotted next to the respective profiles.

analysis we find that it is less massive than the primary. This inconsistency could result from an inaccurate estimation of the luminosity of both components, determined by assuming a similar reddening for both stars. This uncertainty, which we are unable to take into account in our error bars, could lead to an underestimation of the luminosity of the secondary relative to the primary and therefore to an overestimation of the mass of the secondary relative to the primary. This point is discussed in detail in Sec. 7.

5 PROPERTIES OF THE INTENSITY SPECTRUM 5.1 Abundances anomalies Magnetic intermediate mass stars on the main sequence show strong photospheric abundance peculiarities. As some Herbig Ae/Be stars must be the evolutionary progenitors of these chemically peculiar Ap/Bp stars, we searched for abundances anomalies in our spectra, as well as time profile variations due to abundance spots. We chose in the same spectral region lines of different chemical species whose depths predicted by the synthetic spectra are identical. Then in the observed spectra, we compared the depths of these lines. We used this method, insensitive to the veiling phenomenon, to try to find over- or under-abundant species. We found no systematic differences between observed and calculated equivalent widths of numerous chemical species. We conclude that there are no strong abundance anomalies, in the limit of Á0.40 dex in abundance, in the spectrum of HD 200775. In particular, we observe no peculiarity in the helium lines nor any variability of the equivalent width of various species from one night to another and over many nights from Sept 2004 to April 2007. These results are different from the behavior of most of the hot magnetic B stars on the main sequence, which are frequently He-rich or He-weak, and which show abundances starspots varying with the rotational phase of the star. This suggests that at the young age of HD 200775, abundance anomalies have not had sufficient time to develop, or that they are limited by ongoing mixing due to accretion and mass loss. 5.2 The emission lines We observe two kinds of emission lines in the spectrum of HD 200775 : broad lines and narrow lines. Amongst the broad

emission lines we find mainly the Balmer lines and the OI 7772 Õ lines, and some Fe II and Si II with FWHM around 130 km.s-1 . The quality of our data as well as the complex structure of the emission lines do not allow us to fit them with a gaussian function, and therefore to measure the radial velocities of these lines. However the bisector velocities of these lines follow the radial velocity of the secondary (Sec. 3). We note a decrease of the equivalent widths of these emission lines from Sept 2004 to April 2007 (see Fig. 4 upper panel, for H). We observe in Fig. 4 that the maximum equivalent width of H occurs close to the periastron passage of the orbit (indicated by the arrows) that we derive in Section 3. This is consistent with the hypothesis of Pogodin et al. (2004) who suggest that binarity could be the origin of these variations. In addition to these broad lines, we note a large number of narrow emission lines, about a hundred over the whole spectrum, with a FWHM around 19 km.s-1 . Some of these are clearly visible in Fig. 6. We do not observe any shift in radial velocity of these lines during our 2-year observations. Using the numerous observations of Carinae (Thackeray 1953, 1962, 1967; Hamann et al. 1994) we were able to identify all of them and to conclude that they are nebular lines (lines coming from the nebulosity illuminated by HD 200775). Appendix A gives the identification of all lines.

6 MAGNETIC FIELD ANALYSIS Figure 8 shows the Stokes I and V LSD profiles obtained on Sept. 2004 and May 23rd 2005. First we note that the null profile N (Sec. 2) is totally flat indicating that the signature in the V profile is real. Secondly the V profile has shifted from the blue to the red, as the primary component of the I profile, between Sept. 2004 and May 2005. This behaviour is observed throughout the entire data set. For this reason we attribute the observed magnetic field solely to the primary component. To investigate the possible presence of a magnetic field in the secondary, we calculated synthetic Stokes V profiles for a star similar to the secondary one, with a v sin i of 59 km.s-1 , and considering the same magnetic proprieties of the primary, described in Section 6.2. The maximum signal in circular polarisation predicted is around 10-3 , which is the same order of magnitude as the error bars of our observed V profiles. Therefore, if the secondary component of the system hosts a magnetic field of similar topology and


Characterisation of the magnetic field of the Herbig Be star HD 200775
similar intensity to that of the primary, it would be difficult to detect it with our data. On the other hand, if our polarisation spectra contain a signal of the secondary star corresponding to a similar field intensity, it is negligible with respect to the signal of the primary star.

9

6.1 Phased longitudinal field variation Figure 10 shows the variations of the 30 Stokes V profiles, normalised to the continuum intensity of the binary, observed from 2004 to 2007. The longitudinal magnetic field B , which is the projected magnetic field onto the line of sight, can be obtained from both LSD Stokes I and V profiles, as explained by Donati et al. (1997) and Wade et al. (2000). With this aim we compute the Stokes I and V profiles of the primary (IP and VP ) as follows. Using Eq. (1) and the fit of the Stokes I profiles of the binary (Sec. 3), we subtract from the observed Stokes I profile the synthetic profile of the secondary component. Then we renormalised to the continuum of the primary:
IP = 1 - (1 + rL )(1 - I - fS ).

Figure 9. Longitudinal magnetic field plotted in function of the rotation phase. The dashed line is the best sinusoidal fit.

(4)

Table 4. Ranges and minimum bins of parameters explored in the fit of the Stokes V profiles Parameters P (days) T 0 (HJD) i ( ) ( ) Bd (G) ddip (R ) Min 4.0 2453514 0 0 0 -0.4 Max 5.0 2453519 90 180 2000 0.4 bin 0.0001 0.02 1 1 10 0.01

As we assume that the observed magnetic field is only from the primary, we only need to renormalise the observed Stokes V profile by the continuum of the primary, in order to get VP : VP = (1 + rL )V (5)

Using these new profiles and Eq. (1) of Wade et al. (2000), we measured the longitudinal magnetic field B of our data. We fit the variations of B with a sinusoidal function of 4 parameters : Ç Ç Ç Ç P, t0 , B, B0 the rotation period of the star, the reference Julian Date of the maximum of B , the semi-amplitude of the curve, , the shift of the sinusoidal curve with respect to B = 0,

V (v, , ) = -C g0 cb (, )

dI , dv

(6)

consistent with a dipole centered inside the star. The best fit superimposed on the data in Fig. 9, gives the following parameters : P = 4.328 Á 0.003 d, t0 = 2453515.8 Á 0.3 d, B = -309 Á 115 G and B0 = -139 Á 95 G, with a reduced 2 = 1.1. We used that value as a first estimation of the rotation period in the fitting procedure of the Stokes V profiles. This analysis reveals the basic properties of the magnetic field of HD 200775A: that it is organised on large scales; that it has an important global dipole component; and that the dipole geometry is stable on timescales of several years.

Ä where C = 4.67 ç 10-13 Õ-1 G-1 , g and 0 (Õ) are the mean Lande factor and wavelength of the lines used in the mask (Section 2), c is the speed of light, I is the local intensity profile, and v is the velocity. Then we integrated over the visible stellar surface using the limb darkening law with a parameter equal to 0.4 (Claret 2000). We obtain the synthetic Stokes V (v) profile, that we normalised to the intensity continuum, to compare to the observed V profiles. This model depends on five parameters : Ç P the rotation period Ç t0 , the reference Julian Date of the maximum of the surface magnetic intensity, used with P to compute the rotation phase Ç i, the inclination of the stellar rotation axis to the observers line-of-sight, Ç the magnetic obliquity angle, Ç Bd the dipole magnetic intensity, Ç ddip the displacement of the dipole from the center of the star, along the magnetic axis, in stellar radii (R ). We calculated a grid of V profiles for each date of observation (see Table 4 for details on the grid), varying the six parameters (and assuming for the initial value of P the solution obtained from modelling the longitudinal field variation). Then we applied a 2 minimisation to find the best model which matches simultaneously our 30 observed profiles. The best model that we found, with 2 = 1.0, corresponds to P = 4.3281 Á 0.0010 d, i = 60 Á 11 , = 125 Á 8 , Bd = 1000 Á 150 G, and ddi p = 0.05 Á 0.04 R , where the error bars correspond to 3 confidence. Figure 10 shows the synthetic Stokes V profiles superimposed on the observed ones. We see that this model acceptably reproduces most of the observed V profiles.

6.2 Fitting of Stokes V profiles To fit the Stokes V profiles directly, we use the decentred dipole oblique rotator model described by Landstreet (1970). Using the relations of Landstreet (1970), giving the intensity of the magnetic field inside a star, we calculate the longitudinal magnetic field b (, ), at each point (,) of the surface of the star of velocity v, in classical spherical coordinates in the observer's frame. We assume a gaussian local intensity profile of width and depth d. The width is calculated using the resolving power of the instrument and the macroturbulent velocity determined in Section 3. The depth is determined by fitting the Stokes I profiles of the primary component, determined in Section 6.1. We then calculated the local Stokes V profile at each point on the surface of the star using the weak magnetic field relation of Landi degl'Innocenti & Landi degl'Innocenti (1973) :


10

E. Alecian et al.

Figure 10. LSD V profiles (noisy black line) of the 30 spectra superimposed to the best oblique rotator model (smooth red line). The numbers close to the profiles are the rotation phase, and the little bars on the right of the profiles are the mean error bars in V . The profiles are sorted by increasing date, as in Fig. 1.

The small value of the dipole decentring parameter is nearly consistent with zero. The available data are therefore consistent with a dipolar magnetic field at the centre of the star In order to compare this result to the fit of the longitudinal field variations, we estimate the intensity Bd and the obliquity angle of the magnetic dipole component using the fit of the longitudinal field variations. Using the well-known equations relating i, and Bd to the longitudinal field extrema (Borra & Landstreet 1980), and adopting i = 60 Á 11 , we find = 128 Á 13 and Bd = 1500+1100 G. - 700 These values agree within the error bars with those determined with the fit of the Stokes V parameters, which confirms the topology of the magnetic field found using this fit. Table 5 summarises the parameters of the adopted dipole model. Because the Stokes I and V profiles are influenced in similar ways by uncertainties in the procedure that we have used to recover the individual LSD profiles of the primary, the magnetic modeling is only weakly sensitive to the details of the procedures described in Sect. 4 and Sect. 6. From a variety of experiments we find that neither the assumed luminosity ratio, nor the details of the fitting procedure yielding the v sin is, profile depths and radial velocities, have any significant direct influence on the derived magnetic models. Finally, we point out that the magnetic modelling provides a strong constraint on the inclination of the rotation axis i. In particular, for smaller values of i, we are unable to reproduce the Stokes V profile shapes.

Table 5. Magnetic dipole model of HD 200775A P (days) T 0 (HJD) i ( ) ( ) Bd (G) ddip (R ) 4.3281 Á 0.0010 2453515.8 Á 0.8 60 Á 11 125 Á 8 1000 Á 150 0.05 Á 0.04

7 DISCUSSION 7.1 Consequences of the magnetic model on the fundamental parameters of the star Using our determinations of the v sin i (from the LSD I profiles) and the period and inclination (from the magnetic model), we can estimate the radius of HD 200775A. We find R = 2.6 Á 0.5R ,

which is significantly smaller than the radius found from stellar evolutionary models (Table 3). In particular, this value suggest a star close to the main sequence, contrary to the stellar models. The model value is determined from the stellar effective temperature (which we consider to be well-determined) and the inferred luminosity of the primary star. The primary's luminosity is itself determined based on the ratio of equivalent widths of the LSD profiles, and from the dereddened photometry of the combined system supported by Hernandez et al. (2004), which assumed equal reddening Ä for each of the two components. The inferred luminosities could be affected by two possibly incorrect assumptions inherent in this analysis. First, based on our derivation of nearly identical effective temperatures for the two components, we have assumed that the intrinsic equivalent widths of their LSD profiles are equal. We then used the observed equivalent widths of the two components, determined from the LSD profiles, to infer their relative luminosities. However, if the temperatures of the primary and secondary differ by the maximum allowed by the error bars (about 4000 K), the combined effects of differ4 ences in excitation and ioniosation, and the T eff dependence of the flux contribution, could combine to generate important relative differences (either larger or smaller depending on the effective temperatures) in the contributions of the spectra of the two components to the observed LSD profiles. The second assumption that could influence the inferred luminosities is that the two components are affected by similar redden-


Characterisation of the magnetic field of the Herbig Be star HD 200775
ing. In their interferometric study of the HD 200775 system, Monnier et al. (2006) resolved the circumstellar disc of the secondary and found that the secondary was brighter than the primary in H band by a factor of 6.5. If this IR flux difference results from the presence of CS material around the secondary that is not present around the primary (which would be consistent with our spectroscopic observations, indicating that the H emission is associated with the secondary), it would imply that the secondary is affected by additional reddening that has not been taken into account in our analysis. The increased extinction of the secondary would result in an underestimated luminosity of this star relative to the primary. Consequently, the primary's luminosity would be overestimated, leading to a systematic overestimation of its inferred radius. It would also lead to an overestimation of the mass of the primary relative to the secondary and therefore to a mass ratio q = MP /MS greater than one (as observed), which is inconsistent with the orbit determination of q = 0.81 Á 0.22. Neither of these potentially important sources of systematic error can be investigated in detail with the available data, and ultimately we are led to conclude that the primary's radius could be anywhere from the ZAMS radius to the radius reported in Table 3. Likewise the masses and ages of the stars derived in Sec. 4 from the luminosity and temperature must be considered with caution. Individual photometric observations of each of the stars is required to get independent measurements of the luminosities of both stars, and therefore accurate values of their masses, radii and ages. 7.2 Consequences of the magnetic model on the origin of the system Assuming i = 60 implies that the primary would be close to the main sequence, and therefore would have a lower luminosity than the secondary. According to the luminosity ratio in H band of Monnier et al. (2006), the secondary is redder that the primary and therefore seems to be surrounded by denser circumstellar matter than the secondary, which is consistent with our spectroscopic data. This would indicate that the secondary is in a younger evolutionary state than the primary. Furthermore the eccentricity of the orbit found in Sec. 3 is very large (Table 2). These properties may well suggest that the two components are at two very different evolutionary stages, which could imply that the system formed by capture instead of forming initially as a simple close binary.

11

8 CONCLUSIONS In the framework of the understanding of the origin of the magnetic field of the intermediate and high mass stars, we have begun to acquire spectropolarimetric observations of their potential progenitors, the Herbig Ae/Be stars, including HD 200775. A magnetic field has been discovered in the latter and in order to determine its topology, we monitored it during many successive nights and over more than 2 years, using the high resolution spectropolarimeters ESPaDOnS and Narval. First we inspected the intensity spectrum and found that this star is a double-lined spectroscopic binary, whose secondary seems similar, in temperature and luminosity, to the primary. We measured accurate values of the projected rotational velocity of both stars: 26 Á 2 km.s-1 and 59 Á 5 km.s-1 . We fitted the radial velocity curves of both components and found an ephemeris similar to that of Pogodin et al. (2004). The amplitude of the radial velocity curves leads to a poorly constrained mass ratio. Observations of the

separate components, in order to better determine the luminosity ratio of the system, are required to determine accurately the luminosities and hence the masses of the components. According to our determination of the orbital period and of the masses of both stars, we estimate the separation between both components of the system around 16 mas (the distance of the system is estimated around 430 pc by van den Ancker et al. (1997)). Only interferometric observations of the system would be able to separate the two components and to estimate the luminosity ratio. We have also shown that the broad emission observed in the intensity spectrum is linked to the secondary component of the system. We found no evidence of abundance anomalies in the spectrum. Finally, many narrow emission lines have been observed over the whole spectrum, which we identify as nebular lines. The position and velocity variations of the Stokes V profile shows that the detected Zeeman signature corresponds solely to the primary component. However, if the secondary hosts a magnetic field of similar topology and intensity to the primary, the signal in polarisation would be negligible with respect to the signal of the primary. On the other hand, this implies that our observations are relatively insensitive to the presence of a field in the secondary. We modeled the temporal variations of the Stokes V profiles in two different ways. First we measured the longitudinal magnetic field of the star and we fitted a sinusoidal curve, as predicted by the dipole oblique rotator model. Then we used the period found as a first estimation of the rotation period of the star in the Stokes V profiles fitting procedure. We considered an oblique rotator model and we applied 2 minimisation to match simultaneously the 30 Stokes V profiles that we observed over more than one year. We find that this simple decentered dipole model is sufficient to fit most of the profiles. In this way we provide the magnetic field topology: a dipole of intensity 1000 G, displaced by 0.05 Á 0.04 R from the center of the star towards the positive magnetic pole, whose the rotation axis is inclined by 60 Á 11 with respect to the observer line-of-sight, and magnetic axis is inclined by 125 Á 8 with respect to the rotation axis. HD 200775A rotates with a period of 4.3281 Á 0.0010 d. We therefore show that the magnetic field of this star is approximately dipolar, strongly inclined with respect to the rotation axis, with a polar intensity of 1000 G, and stable over more than two years. These characteristics are similar to the topology of the magnetic fields of the Ap/Bp stars (Borra & Landstreet 1980; Bohlender et al. 1987). We have concluded that significant uncertainties exist related to the luminosity, and therefore the mass, radius and age of the magnetic primary star. To clarify these uncertainties and to move forward, additional photometric observations of the individual components are required. In addition, a more complete spectroscopic coverage of the orbital cycle is needed, with the hope of ultimately performing spectral disentangling to yield the individual spectra of the primary and secondary.

ACKNOWLEDGMENTS We warmly thank Georges Alecian for fruitful discussions. We are very grateful to O. Kochukhov, who provided his BINMAG1 code, and who tried to apply his spectral disentangling method to our data. We thank G. Mathys, the referee, for his judicious comments, which led to major improvements in the paper. EA is supported by the Marie Curie FP6 program. GAW acknowledges support from


12

E. Alecian et al.
Wade G. A., Bagnulo S., Drouin D., Landstreet J. D., Monin D., 2007, MNRAS, 376, 1145 Wade G. A., Donati J.-F., Landstreet J. D., Shorlin S. L. S., 2000, MNRAS, 313, 851 Wade G. A., Drouin D., Bagnulo S., Landstreet J. D., Mason E., Silvester J., Alecian E., Bohm T., Bouret J.-C., Catala C., Donati ? J.-F., 2005, A&A, 442, L31 Whitcomb S. E., Gatley I., Hildebrand R. H., Keene J., Sellgren K., Werner M. W., 1981, ApJ, 246, 416

the Natural Science and Engineering Research Council of Canada (NSERC) and the DND Academic Research Programme (ARP).

RE F E RE NCE S Bohlender D. A., Landstreet J. D., Brown D. N., Thompson I. B., 1987, ApJ, 323, 325 Borra E. F., Landstreet J. D., 1980, ApJS, 42, 421 Catala C., Alecian E., Donati J.-F., Wade G., Landstreet J., Bom ? T., Bouret J.-C., Bagnulo S., Folsom C., Silvester J., 2006, A&A, p. submitted Catala C., Bohm T., Donati J.-F., Semel M., 1993, A&A, 278, 187 Catala C., Donati J. F., Bohm T., Landstreet J., Henrichs H. F., ? Unruh Y., Hao J., Collier Cameron A., et al. 1999, A&A, 345, 884 Claret A., 2000, A&A, 363, 1081 Corporon P., Lagrange A.-M., 1999, A&AS, 136, 429 Donati J.-F., Semel M., Carter B. D., Rees D. E., Collier Cameron A., 1997, MNRAS, 291, 658 Fuente A., Martin-Pintado J., Rodriguez-Franco A., MoriartySchieven G. D., 1998, A&A, 339, 575 Gray D. F., 1992, The observation and analysis of stellar photospheres. Cambridge Astrophysics Series, Cambridge: Cambridge University Press, 1992, 2nd ed., ISBN 0521403200. Hamann F., Depoy D. L., Johansson S., Elias J., 1994, ApJ, 422, 626 Herbst W., Shevchenko V. S., 1999, AJ, 118, 1043 Hernandez J., Calvet N., Briceno C., Hartmann L., Berlind P., Ä ~ 2004, AJ, 127, 1682 Hubeny I., 1988, Comput. Phys. Comm., 52, 103 Hubeny I., Lanz T., 1992, A&A, 262, 501 Hubeny I., Lanz T., 1995, ApJ, 439, 875 Hubrig S., Pogodin M. A., Yudin R. V., Scholler M., Schnerr R. S., ? 2007, A&A, 463, 1039 Hubrig S., Scholler M., Yudin R. V., 2004, A&A, 428, L1 ? Kurucz R., 1993, Opacities for Stellar Atmospheres: [-3.5],[4.0],[-4.5]. Kurucz CD-ROM No. 7. Cambridge, Mass.: Smithsonian Astrophysical Observatory, 1993., 7 Landi degl'Innocenti E., Landi degl'Innocenti M., 1973, Sol. Phys., 31, 299 Landstreet J. D., 1970, ApJ, 159, 1001 Miroshnichenko A. S., Mulliss C. L., Bjorkman K. S., Morrison N. D., Glagolevskij Y. V., Chountonov G. A., 1998, PASP, 110, 883 Monnier J. D., Berger J.-P., Millan-Gabet R., Traub W. A., Schloerb F. P., Pedretti E., Benisty M., Carleton N. P., Haguenauer P., Kern P., Labeye P., Lacasse M. G., Malbet F., Perraut K., Pearlman M., Zhao M., 2006, ApJ, 647, 444 Morel P., 1997, A&AS, 124, 597 Palla F., Stahler S. W., 1993, ApJ, 418, 414 Pogodin M. A., Miroshnichenko A. S., Tarasov A. E., Mitskevich M. P., Chountonov G. A., Klochkova V. G., Yushkin M. V., Manset N., Bjorkman K. S., Morrison N. D., Wisniewski J. P., 2004, A&A, 417, 715 Stift M. J., 1975, MNRAS, 172, 133 Thackeray A. D., 1953, MNRAS, 113, 211 Thackeray A. D., 1962, MNRAS, 124, 251 Thackeray A. D., 1967, MNRAS, 135, 51 van den Ancker M. E., The P. S., Tjin A Djie H. R. E., Catala C., de Winter D., Blondel P. F. C., Waters L. B. F. M., 1997, A&A, 324, L33

APPENDIX A: NEBULAR LINES


Characterisation of the magnetic field of the Herbig Be star HD 200775
Table A1: Identification of the nebular lines in the spectrum of HD 200775. Column 1 gives the wavelength observed in our spectra, column 2 gives the equivalent width of the line, column 3 gives the ion, column 4 gives the multiplet number, column 5 gives the rest wavelength of the line, and column 6 gives the reference. I: Thackeray (1953), II: Thackeray (1962), III: Thackeray (1967), IV: Hamann et al. (1994). ( Õ) 4114. 4177. 4210. 4243. 4244. 4276. 4287. 4305. 4319. 4326. 4346. 4352. 4358. 4359. 4372. 4382. 4413. 4416. 4451. 4457. 4474. 4488. 4492. 4414. 4528. 4639. 4664. 4727. 4773. 4814. 4874. 4889. 4898. 4905. 4947. 4950. 4973. 5005. 5020. 5040. 5043. 5107. 5111. 5157. 5158. 5163. 5168. 5181. 5199. 5219. W (mÕ) 8.5 10.9 P CI 2 42.6 11.7 40.7 85.5 9.3 19.7 5.4 10.0 19.5 22.9 70.6 9.4 5.2 56.4 51.1 30.2 29.7 20.0 12.1 7.6 7.4 6.6 12.9 4.5 24.9 17.4 54.2 16.3 33.6 9.7 29.5 7.5 12.1 15.1 13.7 13.6 26.7 6.0 7.6 17.7 17.9 53.5 20.8 9.8 10.0 8.0 17.7 Ion Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Ni I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I Fe I F F F F F F F F F F F F F F e e e e e e e e e e e e e e I I I I I I I I I I I I I I Identification 23F 21F 23F 21F 21F 21F 7F 21F 21F 3F 21F 21F 21F 7F 21F 6F 7F 6F 7F 6F 7F 6F 6F 6F 6F 4F 4F 4F 4F 20F 20F 4F 20F 20F 20F 20F 20F 20F 20F 18F 19F 18F 19F 35F 42 18F 19F 0 ( Õ) 4114.48 4177.21 4211.10 4243.98 4244.81 4276.83 4287.40 4305.90 4319.62 4326.28 4246.85 4352.78 4358.37 4359.34 4372.43 4382.75 4413.78 4416.27 4452.95 4457.95 4474.91 4488.75 4492.64 4414.90 4528.39 4639.68 4664.45 4728.07 4774.74 4814.49 4874.49 4889.63 4905. 4947. 4950. 4973. 5005. 5020. 5043. 5107. 5111. 5158. 5158. 5163. 5169. 5181. 35 38 74 39 52 24 53 95 63 00 81 94 03 97 R II II II II II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I f. I I I I I

13

34 07 99 83 68 69 25 75 47 13 72 64 23 18 28 57 62 11 95 80 75 60 49 74 23 51 29 91 56 38 33 46 44 19 23 58 22 37 08 81 81 75 46 83 61 79 85 78 00 89

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I

3

I I I I I I I

I I I I I I I

F e II

5220.06

III

2 3

Inverse Pcygni profile, the equivalent width have therefore not been determined Blended with 5168.73 unidentified


14

E. Alecian et al.
Table A1: continued. ( Õ) 5261. 5268. 5273. 5282. 5296. 5333. 5347. 5376. 5412. 5414. 5432. 5477. 5527. 5556. 5580. 5673. 5746. 5835. 6364. 6371. 6666. 6729. 6808. 6813. 7154. 7171. 7377. 7387. 7411. 7452. 7495. 7512. 7731. 7999. 8125. 8184. 8216. 8308. 8616. 8891. 8229. 8287. 8694. 8715. 8768. 8819. 9122. 9186. 9203. 9217. 9243. W (mÕ) 64.4 15.6 4.3 5.9 12.6 39.7 4.7 32.8 9.1 4.6 16.6 9.4 19.3 2.1 4.8 8.1 22.0 6.6 9.0 24.1 13.4 4.0 6.0 6.7 51.8 18.2 73.5 15.9 31.4 24.5 11.3 25.6 8.3 129.4 104.0 2.7 97.8 60.1 65.1 26.0 64.6 56.3 31.8 11.1 23.4 17.2 21.9 16.6 42.2 86.7 177.7 Ion F e II F e II F e II F F F F F F F F F F F F F e e e e e e e e e e e e e I I I I I I I I I I I I I I I I I I I I I I I I I I Identification 19F 18F 18F 19F 19F 18F 19F 17F 18F 34F 17F 18F 39F F 34F F 2 2F 31F 31F 8F 14F 14F 2F 14F 2F 14F 0 ( Õ) 5261.61 5268.88 5273.38 5296. 5333. 5347. 5376. 5412. 5433. 5477. 5527. 5556. 5580. 5673. 5746. 5835. 84 65 67 47 64 15 25 33 31 82 22 96 44 R II II II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I f. I I I I I I I I I I I I I I I I

46 70 18 94 67 48 50 28 49 80 99 05 17 18 65 01 78 25 94 07 60 64 99 41 93 74 63 91 41 30 36 91 41 80 05 62 04 22 65 62 43 53 94 46 89 24 64 94 80 98 94

Si Ni Fe Fe Ni Fe Fe Ni Fe Ni Fe Fe Fe Fe Cr Cr

II II II II II II II II II II II II II II II II

1F 1F 2 1F 13F 13F 1F 6 42F

6371.36 6668.8 6729.85 6809.21 6813.73 7155.16 7172.00 7377.90 7388.18 7411.60 7452.54 7495.62 7513.16 7731.68 7999.85 8125.22 8216. 8308. 8616. 8891. 8229. 28 39 96 88 81

II II I I I V V V V V V V V V V V V V V V I

NI Cr I Fe I Fe I Cr I

I I I I

SI F e II

8694.70 8715.84

II II

F e II F e II F e II

9122.91 9187.15 9204.05

IV IV IV

A This paper has been typeset from a TEX/ LTEX file prepared by the author.