arXiv:1111.2186v1 [astro-ph.SR] 9 Nov 2011 aino ik id crto,adjt.T ae these has date, of behaviour number To controlled representative Magnetically a jets. . for Ae/Be and studied Herbig been accretion, combi- interact not wind, a have fields fields disk, displaying these of environment, how Herbig nation circumstellar and im- in the is generated fields It with magnetic are stars. the stars higher how Ae/Be for understand clear to stars), less Tauri portant is T picture classical the low- (the but stars of sequence properties pre-main suc- observational mass outflows many and reproduce accretion driven cessfully process. magnetically formation of in- important Models intermediate-mass are the fields of magnetic gredients Ae/Be that indicated Herbig have several stars of observations Spectropolarimetric Introduction 1. evtosa h .6mNri pia eecp nL Palm La on Telescope Optical Nordic ob- 2.56-m SOFIN the and at servations 085.C-0137(A)) programme (ESO Observatory ⋆ srnm Astrophysics & Astronomy oebr1,2011 10, November ae nosrain banda h uoenSouthern European the at obtained observations on Based hrceiigtemgei ed fteHri eB stars Ae/Be Herbig the of fields magnetic the Characterising ope neato ewe h tla antcfil and field words. magnetic Key stellar the between interaction complex a 3 4 5 o D168.Temgei edgoer a ieyb descr be likely can geometry field magnetic strength The 176386. HD for fasrn Gfil,wihwsudtce norpeiu low- previous our compon in circumstellar undetected the in was detected which was field field, magnetic kG weak strong a of 6 Results. Methods. 8 7 2 1 Aims. poor. very remains iecrual oaie pcr r nepee nterms in interpreted are spectra techniqu moment polarised the geometry. circularly using field measured side magnetic were 480 the MWC for on Herb obtained fund constraints stellar the place Using we of periods. curves, observations rotational low-resolution their determine 2 FORS and SOFI 41 the with discuss 480 MWC We for m obtained were spectropolarimetric spectra in polarimetric obtained spectra low-resolution 9 Context. 2011 10, November version: online Preprint .Hubrig S. eateto srnm,Uiest fMcia,AnArbor Ann Michigan, of University Astronomy, of Department eateto hoeia hsc n srpyis Masar Astrophysics, and Physics Theoretical of Department bevtr n lntru fJ aia V Palisa, J. of Planetarium and Observatory NFOsraoi srnmc iPlro izadlParl del Piazza Palermo, di Astronomico INAF-Osservatorio sa etnIsiueo hl,SitPtrbr Branch, Saint-Petersburg Chile, of Institute Newton Isaac Russia 196140, Saint-Petersburg, Observatory, Pulkovo uoenSuhr bevtr,Karl-Schwarzschild-Str. Observatory, Southern European ebi-ntttfu srpyi osa AP,A e St der An e-mail: (AIP), Potsdam f¨ur Astrophysik Leibniz-Institut eatmnod ´sc srnmı,Fcla eCiencia de Astronom´ıa, Facultad F´ısica y de Departamento epeettemgei rpriso orHri eB star Ae/Be Herbig four of properties magnetic the present We D908 D109,H 736 n W 480 MWC and 176386, HD 150193, HD 97048, HD u erhfrproiiisrsle in resulted periodicities for search Our u nweg ftepeec n h oeo antcfilsi fields magnetic of role the and presence the of knowledge Our B [email protected] u eut o h he ebgA/esasH 74,H 1501 HD 97048, HD stars Ae/Be Herbig three the for results Our d fsvrlhnrdGus h ogtdnlmgei-edm magnetic-field longitudinal The Gauss. hundred several of tr:pemi eune—sas topee tr:magn stars: — atmospheres stars: — sequence pre-main stars: 1 .Sch¨oller M. , aucitn.hubrig17407 no. manuscript 2 .Ilyin I. , 1 .R Cowley R. C. , .V Yudin V. R. B-TcnclUiest,Otaa zc Republic Czech Ostrava, University, Technical - SB ˇ P 0 = ABSTRACT . 9 o D97048, HD for d 693 a. 7 mna aaeesadtelniuia antcfil phas field magnetic longitudinal the and parameters amental , falniuia antcfield magnetic longitudinal a of d ihFR one nteVT e high-resolution New VLT. the on mounted 2 FORS with ode 8 3 n .Cur´e M. and , .Mikul´aˇsek Z. , ,878Grhn e M¨unchen, Germany bei Garching 85748 2, antcfil esrmnsfrteHri eB stars Ae/Be Herbig the for measurements magnetic-field look sequence to pre-main detail of in type targets this stars. encompassing of from patterns number for come sufficient disk- only a the however, studying understanding can, mag- in interaction stellar Progress a magnetosphere required. and is some plasma field photosphere, to accreting netic stellar midplane the disk the between the ac- near coupling from the latitudes lifted Since stellar be higher stars. to has Ae/Be material Herbig creting understanding of to magnetospheres importance configura- utmost field the of Magnetic obviously d. 42.1 are of tions period modulus rotation field uncov- the magnetic line UVES over mean analysed velocities, and of radial asymmetries, been line amount widths, have widths, equivalent small star in variations The this ering for 2011). 101412 acquired HD 2010, star spectra al. Ae Herbig et magnetic vari- (Hubrig strongly photometric and the spectroscopic of the ability both in found been nso h Ca the of ents he ihrslto iclryplrsdSFNspectra SOFIN polarised circularly high-resolution Three Russia h icmtla environment. circumstellar the rwre1,142Ptdm Germany Potsdam, 14482 16, ernwarte ,Uiesddd Valpara´ıso, Chile de Universidad s, pcrgahisalda h odcOtclTelescope. Optical Nordic the at installed spectrograph N kUiest,Bn,CehRepublic Czech Brno, University, yk mno1 03 aem,Italy Palermo, 90134 1, amento gA/esasH 74,H 513 n D176386 HD and 150193, HD 97048, HD stars Ae/Be ig nti ok epeetsre fma longitudinal mean of series present we work, this In hr aeeghsit ewe ih-adleft-hand and right- between shifts wavelength where e bdb ete ioewt oa antcfield magnetic polar a with dipole centred a by ibed I41914,USA 48109-1042, MI , eouinplrmti bevtoswt OS1 A 1. FORS with observations polarimetric resolution hthv o enpeiul tde ndetail. in studied previously been not have that s aueet fMC40rva h presence the reveal 480 MWC of easurements nemdaems r-ansqec stars pre-main-sequence intermediate-mass n 3 n D168 r ae nmulti- on based are 176386 HD and 93, P ii tcfil tr:vrals general variables: stars: — field etic 9 1 = 4 & ie n h Na the and lines H&K , 5 .Stelzer B. , . 1 o D109,and 150193, HD for d 317 h B z i 6 . .A Pogodin A. M. , i ie,indicating lines, D P ⋆ 0 =

c S 2011 ESO . 9 d 899 7 , 8 , e 1 S. Hubrig et al.: Magnetic fields of Herbig Ae/Be stars

Table 1. Target stars with multi-epoch spectropolarimet- Table 2. Magnetic field measurements of Herbig Ae/Be ric observations. stars with FORS 1/2.

Object Other V Spectral Name Identifier Type MJD Phase hBziall [G] hBzihyd [G] MWC 480 HD 31648 7.7 A3Ve HD 97048 CU Cha 8.5 A0pshe HD 97048 HD 150193 V2307 Oph 8.9 A1Ve 54609.137 0.083 164±42 188±47 HD 176386 CD−37 13023 7.3 B9V 55279.131 0.349 −92±41 −135±58 55311.109 0.467 −163±53 −195±59 Notes. Spectral types and visual magnitudes are taken from 55312.035 0.802 −6±74 12±91 SIMBAD. 55320.094 0.426 −44±54 −62±62 55322.119 0.346 −123±42 −62±48 HD97048, HD150193, and HD176386 obtained at low reso- 55324.105 0.210 4±59 9±63 lution with the multi-mode instrument FORS 2 at the VLT, 55326.987 0.367 −124±34 −156±44 55334.039 0.537 −154±53 −193±61 and describe their magnetic field geometries. Magnetic 55335.119 0.094 127±120 135±49 fields of the order of 120–250G were detected in these stars 55345.987 0.768 54±59 26±63 for the first time during our visitor run with FORS 1 in May 55347.180 0.488 −186±61 −231±64 2008 (Hubrig et al. 2009). The follow-up program developed 55348.070 0.772 −36±38 −12±64 to search for magnetic field variations over the rotation cy- 55349.102 0.260 −7±74 65±78 cle was allocated for observations with FORS 2 in 2010. HD 150193 The previous spectropolarimetric observations of the star 54609.092 0.118 −144±32 −252±48 MWC480 (=HD31648) with FORS1 revealed the presence 55321.364 0.523 −124±60 −184±62 of a weak longitudinal magnetic field hBzi = 87 ± 22 G 55312.319 0.648 −97±84 −239±102 and distinct circular polarisation signatures in spectral lines 55320.179 0.622 −148±69 −147±80 originating in the circumstellar (CS) environment (Hubrig 55322.295 0.230 −109±48 −111±64 et al. 2006, 2007). We report here on three new high- 55323.311 0.003 −39±34 −43±46 resolution polarimetric spectra obtained with the SOFIN 55324.285 0.744 −107±43 −124±64 spectrograph installed at the Nordic Optical Telescope, 55325.307 0.520 −189±45 −228±59 which confirm the presence of a magnetic field in this star 55326.341 0.306 −142±53 −176±70 in both photospheric and CS lines. 55327.340 0.066 −19±36 −15±48 55328.356 0.838 −62±74 −46±99 55334.314 0.366 −143±32 −132±74 2. Magnetic field measurements and period 55341.279 0.661 −217±58 −235±74 55348.125 0.864 −120±96 −306±130 determination 55354.296 0.555 −201±55 −229±82 The targets for which we present new spectropolarimet- HD 176386 ric measurements are listed in Table 1. Multi-epoch po- 54610.272 0.174 −119±33 −121±35 larimetric spectra of the Herbig Ae/Be stars HD 97048, 55312.353 0.145 −10±69 −42±88 HD150193, and HD176386 were obtained with FORS21 on 55322.313 0.224 −41±65 −44±89 Antu (UT1) from 2010 March 23 to 2010 June 7 in service 55323.247 0.263 −62±62 −126±64 mode. Using a slit width of 0′′. 4, the spectral resolving power 55324.308 0.444 −190±62 −179±65 of FORS2 achieved with the GRISM 600B is about 2000. 55325.331 0.582 −201±63 −225±71 A detailed description of the assessment of the longitudinal 55326.365 0.732 −115±54 −119±68 magnetic-field measurements using FORS 2 was presented 55327.244 0.709 −6±64 −49±79 55328.375 0.967 99±59 111±79 in our previous papers (e.g., Hubrig et al. 2004a, 2004b, and 55334.339 0.602 −65±49 −50±63 references therein). The mean longitudinal magnetic field, 55335.264 0.631 −102±43 −144±54 hBzi, was derived using 55341.317 0.364 −101±65 −104±71 2 55349.128 0.052 115±50 121±68 V geff eλ 1 dI 55354.226 0.724 −130±41 −162±62 = − hBzi , (1) I 4πmec2 I dλ 55355.134 0.733 −25±38 −30±61 where V is the Stokes parameter that measures the circular polarisation, I is the intensity in the unpolarised spectrum, g is the effective Land´efactor, e is the electron charge, eff Cassegrain focus of the Nordic Optical Telescope (NOT). λ is the wavelength, me the electron mass, c the speed of We used a 2K Loral CCD detector to register 40 echelle light, dI/dλ is the derivative of Stokes I, and hBzi is the ˚ mean longitudinal magnetic field. orders partially covering the range from 3500 to 10000 A with a length of the spectral orders of about 140 A˚ at Three spectropolarimetric observations of MWC480 ˚ with S/N ≥ 200 were obtained on 2009 December 30 and 5500 A. The polarimeter is located in front of the entrance 2010 December 14, and 2010 December 25 with the low- slit of the spectrograph and consists of a fixed calcite resolution camera (R = λ/∆λ ≈ 30000) of the echelle spec- beam splitter aligned along the slit and a rotating super- trograph SOFIN (Tuominen et al. 1999) mounted at the achromatic quarter-wave plate. Two spectra polarised in opposite senses are recorded simultaneously for each echelle 1 The spectropolarimetric capabilities of FORS 1 were moved order providing sufficient separation by the cross-dispersion to FORS 2 in 2009. prism. Usually two sub-exposures with the quarter-wave

2 S. Hubrig et al.: Magnetic fields of Herbig Ae/Be stars

HD 97048 The frequency analysis was performed on the longi- 300 tudinal magnetic field measurements using a non-linear 200 least squares fit to the multiple harmonics utilizing the 100 Levenberg-Marquardt method (Press et al. 1992) with an 0 optional possibility of pre-whitening the trial harmonics. -100 To detect the most probable period, we calculated the fre- Mean longitudinal field [G] -200 quency spectrum and for each trial frequency we performed -300 400 a statistical F-test of the null hypothesis for the absence of -0.25 0 0.25 0.5 0.75 1 1.25 300 200 Phase periodicity (Seber 1977). The resulting F-statistics can be 100 thought of as the total sum including covariances of the 0

Residuals [G] -100 ratio of harmonic amplitudes to their standard deviations, -200 i.e. a signal-to-noise ratio. The derived ephemeris for the -300 -0.25 0 0.25 0.5 0.75 1 1.25 detected periods are Phase pos ext HD 150193 HD 97048 : hBzi = 200 MJD55318.41923 ± 0.01372+ 0.69334 ± 0.00039E 100 HD 150193 : hB ipos ext = 0 z

-100 MJD55318.05538 ± 0.03072+ 1.31697 ± 0.00013E pos ext -200 HD 176386 : hBzi = Mean longitudinal field [G]

-300 MJD55325.70578 ± 0.02185+ 0.89920 ± 0.00008E

200 -0.25 0 0.25 0.5 0.75 1 1.25 Phase The logbook of the new FORS2 and the old FORS1 100 spectropolarimetric observations is presented in Table 2. 0

Residuals [G] In the first column, we indicate the MJD value at mid -100

-200 exposure. The phases of the measurements of the mag-

-0.25 0 0.25 0.5 0.75 1 1.25 netic field are listed in Column 2. In Columns 3 and 4, we Phase present the longitudinal magnetic field hBziall measured us- ing the whole spectrum and the longitudinal magnetic field HD 176386 hBzi using only the hydrogen lines. All quoted errors 300 hyd

200 are 1σ uncertainties. The FORS 1 measurements published

100 by Hubrig et al. (2009) are presented in the first line, re-

0 spectively. Phase diagrams of the data folded with the de-

-100 termined periods are presented in Fig. 1. The quality of

Mean longitudinal field [G] 2 -200 our fits is described by a reduced χ -value that is 0.71 for -300 HD98048, 0.32 for HD150193, and 0.92 for HD176386. 300 -0.25 0 0.25 0.5 0.75 1 1.25 200 All three Herbig Ae/Be stars exhibit a single-wave vari- Phase 100 ation in the longitudinal magnetic field during the stel- 0 lar rotation cycle. These observations are usually consid- Residuals [G] -100 -200 ered as evidence of a dominant dipolar contribution to the -300 magnetic field topology. Assuming that the studied Herbig -0.25 0 0.25 0.5 0.75 1 1.25 Phase Ae/Be stars are oblique dipole rotators, the magnetic dipole max min axis tilt β is constrained by hBzi / hBzi = cos(i + Fig. 1. Phase diagrams with the best sinusoidal fit for β)/ cos(i − β), where the inclination angle i is derived from the longitudinal magnetic field measurements of HD 97048 resolved observations of their disks. In Table 3, we show for (top), HD 150193 (middle), and HD 176386 (bottom). The each star in rows 1 and 2 the mean value hBzi and the am- residuals (observed – calculated) are shown in the lower plitude of the field variation AhBzi. For the star HD97048, panel. The deviations are mostly of the same order as the Lagage et al. (2006) used VISIR observations of the mid- error bars, and no systematic trends are obvious, which infrared emission of polycyclic hydrocarbons at the surface justifies a single sinusoid as a fit function. of the disk to determine i = 42.8 ± 2.5◦. Using this value, we obtain a magnetic obliquity of β = 84.8 ± 5.0◦ and the dipole strength Bd = 720 ± 120 G. The disk inclination of plate angles separated by 90◦ are necessary to derive circu- HD 150193, i = 38 ± 9◦, was determined by Fukagawa et larly polarised spectra. The spectra are reduced with the 4A al. (2003) using Subaru near-infrared imaging. With this software package (Ilyin 2000). Bias subtraction, master flat- value, we find a magnetic obliquity β = 42 ± 10◦ and the field correction, scattered light subtraction, and weighted dipole strength Bd = 590 ± 40 G. extraction of spectral orders comprise the standard steps For HD 176386, the inclination angle i is not known be- of the image processing. A ThAr spectral lamp is used for cause both IR (Siebenmorgen et al. 2000) and far-UV spec- wavelength calibration, taken before and after each target troscopy (Martin-Za¨ıdi et al. 2008) suggest that the circum- exposure to minimise temporal variations in the spectro- stellar matter has the form of an envelope rather than the graph. Wavelength shifts between right- and left-hand side form of a disk. To constrain the magnetic field geometry of circularly polarised spectra were interpreted in terms of a this star, we determined that v sin i = 165 ± 8kms−1 from longitudinal magnetic field hBzi, using the moment tech- our NTT/EMMI spectrum obtained for this star in 1995 nique described by Mathys (1994). July 8. In Fig. 2, we present the line profile of Mg ii 4481

3 S. Hubrig et al.: Magnetic fields of Herbig Ae/Be stars

Table 3. Magnetic field models constrained by the FORS 1/2 spectropolarimetric observations.

Object HD 97048 HD 150193 HD 176386

hBzi [G] −14.8±14.3 −106.9±7.4 −36.2±15.9 AhBzi [G] 151.6±19.6 74.8±9.8 127.5±24.7 v sin i [km s−1] 165±8 R [R⊙] 3.0±0.5 −1 veq [km s ] 169±28 i [◦] 42.8±2.5 38±9 78±46 β [◦] 84.8±5.0 41.8±10.1 37±106 Bd [G] 720±120 590±40 700±650

Table 4. Magnetic field measurements of MWC 480 with SOFIN.

MJD hBzi [G] 53296.350 87±22 55195.086 −952±177 55544.090 458±149 55555.982 189±204

the radius R =3 ± 0.5 R⊙. With this radius and the period −1 P =0.899 d, we obtain veq = 169 ± 28kms , which leads to the poorly constrained i = 78 ± 46◦, magnetic obliquity ◦ β = 37 ± 106 , and the dipole strength Bd = 700 ± 650 G. The parameters of the magnetic field dipole models for all three stars are listed in the last two rows of Table 3. Fig. 2. Determination of the v sin i value for HD 176386. ii The logbook of our new SOFIN spectropolarimetric ob- The line profile of Mg 4481 in the NTT/EMMI spectrum servations of MWC480 is presented in Table 4. In the first is fitted by the corresponding synthetic spectrum assuming column, we list the MJD values for the middle of the spec- T = 9900K and log g =3.7. eff tropolarimetric observations and in Col. 2 we present the mean longitudinal magnetic field hBzi. In the first line, we list the earlier FORS 1 measurement published by Hubrig et al. (2006). All quoted errors are 1σ uncertainties. Since only four magnetic field measurements are available, no search for a periodicity can be undertaken. On the other hand, the rotation period can be estimated using the known in- clination angle i = 36 ± 1◦ (Pi´etu et al. 2006), the ra- dius R = 1.67 R⊙ (Blondel & Tijn A Djie 2006), and v sin i = 90 ± 2kms−1 estimated on our SOFIN spectra. In Fig. 3, we present the line profile of Mg ii 4481 over- plotted by the synthetic spectrum assuming Teff = 8500K and log g = 3.5. With this value for v sin i, we obtain −1 veq = 153 ± 8kms . By substituting these values into the relation P = 50.6 R/veq, it follows that P =0.531±0.029d. This period is probably very close to the true rotation pe- riod, since Grady et al. (2010) report that their FUSE ob- Fig. 3. Determination of the v sin i value for MWC 480. servations integrated over the individual FUSE orbits indi- ii The line profile of Mg 4481 in the SOFIN spectrum is cate a marginal detection of variability on the timescale of fitted by the corresponding synthetic spectrum assuming 0.5 d. Teff = 8500K and log g =3.5. 3. Magnetic field and CS environment of MWC 480 overplotted by the corresponding synthetic spectrum as- In our earlier studies based on the low-resolution po- suming Teff = 9900K and log g = 3.7, and solar abun- larimetric FORS1 spectra, we reported the detection of dances. We note that the literature values for Teff widely distinct Zeeman features in several Herbig Ae/Be stars differ from each other, from Teff = 10715K as given by van at the positions of the Ca ii H&K lines (Hubrig et al. den Ancker et al. (1997) to Teff = 9527K derived by Meyer 2004c, 2006, 2007, 2009). The Ca ii H&K lines frequently & Wilking (2006). The equatorial rotation velocity can be display variable multi-component complex structures in calculated using the relation veq = 50.6 R/P , where R is both the Stokes V and Stokes I spectra and it has been sug- the stellar radius in solar units and P is the rotation pe- gested that they originate at the base of the stellar wind, riod in days. Taking into account the uncertainties in atmo- as well as in the gaseous accretion flow. The star MWC480 spheric parameter determinations for HD 176386, we adopt exhibits a significant spectral variability, which implies that

4 S. Hubrig et al.: Magnetic fields of Herbig Ae/Be stars

1.5 15

1.4

2010-12-25 10 1.3 N 1.2 2010-12-14 5 1.1

Normalized Flux 1.0 0 7.72 7.74 7.76 7.78 7.8 7.82 7.84 7.86 7.88 Hp [mag] 0.9 2009-12-30 Fig. 6. Available Hipparcos measurements for MWC 480 0.8 with outliers observed during HJD 8674.75-8674.94. -400 -200 0 200 400 Velocity [km/s] magnetic field geometry, then we expect to measure a mag- netic field strength of −76G detected in the CS compo- Fig. 4. The behaviour of the Fe ii line at λ5018.4 at the nent of the Ca ii K line at a distance of ∼2.3 stellar radii. three observational epochs 2009 December 30 (bottom), This consideration is fully in line with the generally as- 2010 December 14 (middle), and 2010 December 25 (top). sumed scenario of magnetospheric accretion in Herbig Ae and stars. The stellar magnetic field truncates the accretion disk at a few stellar radii and gas accretes along there is a strong interaction between the star and its CS en- magnetic channels from the protoplanetary disk to the star. vironment. The temporal behaviour of the complex variable The disk is truncated at a radius RT, which depends upon structure of Hα, He i, and the Na i D line was monitored the strength of the surface magnetic field, the mass accre- in several previous studies (e.g., Beskrovnaya & Pogodin tion radius, the , and the stellar radius (e.g. 2004; Kozlova et al. 2007), which showed that the stellar Bouvier et al. 2007). For a dipole configuration of the mag- netic field geometry, a magnetic field strength of −120G wind is non-homogeneous and consists of several compo- i nents, differing from each other in terms of their velocities. detected in the CS component of the Na D is expected to Furthermore, structural variations in the stellar wind have be measured at a distance of ∼2.0 stellar radii. been found to correlate with changes in the accretion pro- cess in the envelope. 4. Discussion In our high-resolution SOFIN spectra, the contamina- tion of the spectra by the CS material is clearly visible. In Using FORS 2 low-resolution multi-epoch polarimetric Fig. 4, we show the presence of a sharp blue component in spectra, we have determined rotation/magnetic periods for the Mult. 42 Fe ii line at λ5018.4 during the first and the the three Herbig Ae/Be stars HD97048, HD150193, and second observational epochs in 2009 and 2010. During the HD176386. In the framework of the rigid rotator model last observation, the wings of the Fe ii line appear in emis- usually assumed for magnetic stars on the upper main- sion. In addition, two other lines, λλ5169.0 and 4923.9 of sequence, the period of the magnetic field variation cor- the same multiplet, exhibit a similar behaviour. The shift responds to the period. Our study indicates of this sharp blue component, reaches ∼153kms−1 during that dipole models provide a satisfactory fit to the acquired the first SOFIN observation in 2009, and is of the same or- magnetic data. der as the shifts measured for the blue components of the On the other hand, our search for rotation-modulated Na i D, the Ca ii H&K lines, and the cores of the hydrogen photometric variability using the available photometric lines. databases has been unsuccessful: no clear variation was In agreement with previous studies, the Na i D lines are found for the three stars in either the ASAS database highly variable, but the Ca ii H&K lines also show a com- (Pojma´nski 2002) or Hipparcos photometry. The star plex structure consisting of several components and pro- MWC 480 was not observed by ASAS, and only Hipparcos nounced variability. The behaviour of the Na i and Ca ii photometry is available for a period search. It displayed doublets at three different observation epochs is presented moderate light changes of up to 0.098 mag on timescales in Fig. 5. Owing to strong blending, it is impossible to from tens of hours to several days. The behaviour of measure the shift between the line profiles in right- and MWC480 resembles the characteristics typical of other left-hand side circularly polarised spectra for all line com- photometrically monitored Herbig Ae/Be stars, for which ponents. The compilation of measurements of the compo- light variations are likely of stochastic nature and caused nent positions in the Ca ii and Na i doublet lines and our by the pre-main-sequence disk accretion phenomena (e.g. attempts to determine the longitudinal magnetic field are Rucinski et al. 2010). Using archive Hipparcos data, we presented in Table A.1. have found an indication of a variability with a period of Weak magnetic fields of the order of several tens of 1.6 d. However, it seems that six outliers appear during the Gauss were detected in the CS components of the Na i and period HJD 8674.75–8674.94. This group of outliers is lo- Ca ii doublet lines and can probably be considered as CS cated on the right side of the histogram presented in Fig. 6. magnetic fields. If we assume a dipole configuration for the Excluding these outliers from the analysis, we found no

5 S. Hubrig et al.: Magnetic fields of Herbig Ae/Be stars

2.0 2010-12-25 2.5 2010-12-25 2010-12-25

1.5 2.0 2010-12-14 2010-12-14 2010-12-14 1.5

1.0 2009-12-30 1.5 1.0 Normalized Flux Normalized Flux Normalized Flux

2009-12-30 2009-12-30 1.0 0.5 0.5

5885 5890 5895 5900 3930 3932 3934 3936 3964 3966 3968 3970 3972 3974 Wavelength [A]Ê Wavelength [A]Ê Wavelength [A]Ê

Fig. 5. The variability and complex structure of Na i D (left) and Ca ii K&H (middle and right) in the spectra of MWC480 at the same observational epochs as in Fig. 4.

7.72 0.06 7.74

7.76 0.05

7.78 0.04

7.8 0.03 Hp [mag] 7.82 amplitude [mag] 0.02 7.84

7.86 0.01 MWC 480

7.88 0 7900 8000 8100 8200 8300 8400 8500 8600 8700 0 0.5 1 1.5 2 HJD−2440000 frequency [d−1]

0.04

−0.15 V380 Ori P = 4.31726 d 0.035 rot

0.03 −0.1

0.025 −0.05

0.02 0

amplitude [mag] 0.015 0.05

0.01 detrended magnitude [mag] 0.1 0.005 0.15 0 0.5 1 1.5 2 2.5 3 −1 frequency f [d ] 0 0.2 0.4 0.6 0.8 1 rotational phase Fig. 7. Hipparcos observations of MWC480 and the corre- sponding frequency periodogram. The six outliers are indi- Fig. 8. Frequency periodogram from Hipparcos observa- cated by an ellipse. tions of V380Ori. The red arrow in this figure denotes the rotation period determined by Alecian et al. (2009). The lower panel presents the corresponding phase diagram. No indication for periodicity. In Fig. 7, we present Hipparcos periodicity is detected. observations of MWC480 and the corresponding frequency periodogram. Kozlova et al. (2007) studied the spectral variability of of the red-shifted component in the Ca ii H&K lines during MWC 480 in the region of Hα. They found an equivalent the second and third epochs of the SOFIN observations. width variability on a timescale of about 1200 d and sug- Our experience shows that the photometric periodic gest that the most plausible reason for the detected cyclic variability of Herbig Ae stars is generally difficult to de- variability is the reconstruction of the inner structure of tect. As an example, we show in Fig. 8 Hipparcos photo- the CS gas envelope, caused by the presence of a low mass metric data for V380 Ori, for which the rotation period companion or planet around the star. It is not yet clear how P = 4.31276 ± 0.00042d was determined from magnetic and whether this periodicity is related to the appearance field measurements (Alecian et al. 2009). As we show in

6 S. Hubrig et al.: Magnetic fields of Herbig Ae/Be stars this figure, no periodicity can be detected in the photomet- Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. ric data. 1992, Numerical Recipes, 2nd ed. (Cambridge University Press: Until recently, rotation periods had been known for Cambridge) Rucinski, S. M., Zwintz, K., Hareter, M., et al. 2010, A&A, 522, A113 only three Herbig Ae stars, V380 Ori, HD101412, and Seber, G. A. F. 1977, Linear Regression Analysis (Wiley: New York) HD135544B. They had been determined from direct mea- Siebenmorgen, R., Prusti, T., Natta, A., & M¨uller, T. G. 2000, A&A, surements of longitudinal magnetic fields (for V380 Ori and 361, 258 HD101412), photometry (for HD101412), or radial veloc- Szeifert, T., Hubrig, S., Sch¨oller, M., et al. 2010, A&A, 509, L7 Tuominen, I., Ilyin, I., & Petrov, P. 1999, in “Astrophysics with the ities (for HD 135544B) (Alecian et al. 2009; Hubrig et al. NOT”, eds. H. Karttunen & V. Piirola, University of Turku, 2011; M¨uller et al. 2011). Tuorla Observatory, 47 The longitudinal magnetic fields of HD 97048, van den Ancker, M. E., The, P. S., Tjin A Djie, H. R. E., et al. 1997, HD150193, and HD176386 are not very strong com- A&A, 324, L33 pared to the kG fields detected in V380Ori, HD101412, and Z CMa (Alecian et al. 2009; Hubrig et al. 2010, 2011; Appendix A: Magnetic field measurements of Szeifert et al. 2010). Grady et al. (2010) placed limits on ii i the field strength for MWC 480 assuming that this star is components identified in the Ca and Na +5 accreting magnetospherically. Their result, Baver =3−2 kG, doublet line profiles. appears comparable with our measurements presented in During the first epoch of our SOFIN observations in this work. December 2009, the Ca doublet lines were found to have The detection of a weak longitudinal magnetic field in three components. Two of them are shifted to the blue by the blue-shifted CS components of the Na i and the Ca ii 146.4 km s−1 and 83.1 km s−1, respectively, and the third doublet lines in the polarimetric spectra of MWC 480 im- one is shifted to the red by 13.7 km s−1. This velocity cor- plies that for Herbig Ae/Be stars the Doppler-shifted spec- responds to the of the photospheric spec- tropolarimetric contributions from photospheric and cir- tral lines. During the epochs 2 and 3, an additional strong cumstellar environmental diagnostic lines should be investi- broad component appears on the red side. The shift mea- gated using high-resolution polarimetric observations over sured during the second epoch is almost 100kms−1. Such the rotation period. These studies will allow us to apply the a strong broad red-shifted component is also discovered in technique of Doppler Zeeman tomography to determine the the higher number Balmer lines, providing potential evi- correspondence between the magnetic field structure and dence of a lower mass companion. In any case, future studies the radial density and temperature profiles. should focus on multi-epoch spectroscopic observations to Acknowledgements. This research has made use of the SIMBAD monitor the radial velocities and the line profile variations database, operated at CDS, Strasbourg, France. MAP and RVY ac- in this star and prove the presence of such a companion. knowledge the support obtained by RFBR grant No 07-02-00535a and Wavelength shifts for MWC 480 and the corresponding ve- Sci.Schole No 6110.2008.2. locities and magnetic field measurements of its components identified in the Ca ii and Na i doublet line profiles are pre- sented in Table A.1. The blue CS components in the Na i D References lines appear to be clearly separated at the first epoch al- lowing us to measure a weak longitudinal magnetic field of Alecian, E., Wade, G. A., Catala, C., et al. 2009, MNRAS, 400, 354 i Beskrovnaya, N. G., & Pogodin, M. A. 2004, A&A, 414, 955 the order of ∼−100G. The Na D lines themselves are in Blondel, P. F. C., & Tjin A Djie, H. R. E. 2006, A&A, 456, 1045 emission. They have overlapping sharp absorption compo- Bouvier, J., Alencar, S. H. P., Harries, T. J., et al. 2007, Protostars nents. During the second and third epochs of observations, and Planets V, B. Reipurth, D. Jewitt, and K. Keil (eds.), p. 479 the blue CS components indicate a double structure and (University of Arizona Press: Tucson) are located closer to the emission Na i D lines. Fedele, D., van den Ancker, M. E., Acke, B., et al. 2008, A&A, 491, 809 As for the large shifts in the blue components of the Fukagawa, M., Tamura, M., Itoh, Y., et al. 2003, ApJ, 590, L49 Na i D lines and the Ca ii H&K lines detected in the SOFIN Grady, C. A., Hamaguchi, K., Schneider, G., et al. 2010, ApJ, 719, spectrum in 2009, the cores of the hydrogen lines from Hβ 1565 to H8 in our FORS1 spectra obtained in 2004 are also Hubrig, S., Kurtz, D. W., Bagnulo, S., et al. 2004a, A&A, 415, 661 −1 Hubrig, S., Szeifert, T., Sch¨oller, M., et al. 2004b, A&A, 415, 685 shifted to the blue by ∼150kms . In Fig. A.1, we show Hubrig, S., Sch¨oller, M., & Yudin, R. V. 2004c, A&A, 428, L1 the behaviour of the Hγ and Hδ line profiles together with Hubrig, S., Yudin, R. V., Sch¨oller, M., & Pogodin, M. A. 2006, A&A, the calculated synthetic spectra for atmosphere model pa- 446, 1089 rameters Teff = 8500K and log g = 3.5. Interestingly, this Hubrig, S., Pogodin, M. A., Yudin, R. V., et al. 2007, A&A, 463, 1039 Hubrig, S., Stelzer, B., Sch¨oller, M., et al. 2009, A&A, 502, 283 pattern is the opposite of that seen for the Herbig Ae star Hubrig, S., Sch¨oller, M., Savanov, I., et al. 2010, Astron. Nach., 331, HD 101412, where the higher order Balmer lines display red 361 shifts that decrease towards the lower order Balmer lines. Hubrig, S., Mikulasek, Z., Gonz´alez, J. F., et al. 2011, A&A, 525, L4 It is possible that the opposite trend for hydrogen line- Ilyin, I. 2000, Ph.D. Thesis, University of Oulu, Finland core velocities in both stars can be explained by different Kozlova, O. V., Alekseev, I. Yu., & Shakhovskoi, D. N. 2007, Astrophysics, 50, 467 viewing angles, i.e. we view MWC 480 mostly in the wind ◦ Lagage, P-O., Doucet, C., Pantin, E., et al. 2006, Science, 314, 621 direction at an angle i = 36±1 (Pi´etu et al. 2006), whereas Martin-Za¨ıdi, C., Deleuil, M., Le Bourlot, J., et al. 2008, A&A, 484, HD 101412 is viewed in the direction of the accretion flow 225 with i = 80 ± 7◦ (Fedele et al. 2008). Mathys, G. 1994, A&AS, 108, 547 M¨uller, A., van den Ancker, M. E., Launhardt, R., et al. 2011, A&A, 530, A85 Meyer, M. R., & Wilking, B. A. 2006, PASP, 121, 350 Pi´etu, V., Dutrey, A., Guilloteau, S., et al. 2006, A&A, 460, L43 Pojma´nski, G. 2002, Acta Astron., 52, 397

7 S. Hubrig et al.: Magnetic fields of Herbig Ae/Be stars

Table A.1. Wavelength shifts for MWC 480 with corresponding velocities and magnetic field measurements of compo- nents identified in the Ca ii and Na i doublet line profiles. Numbers with colons indicate the most uncertain values due to blends.

Epoch 1 Epoch 2 Epoch 3 Component ∆λ velocity shift hBzi ∆λ velocity shift hBzi ∆λ velocity shift hBzi [A]˚ [km s−1] [G] [A]˚ [km s−1] [G] [A]˚ [km s−1] [G] Ca ii 3933.7 1st −1.92 −146.4 −0.86 −65.6 17 2nd −1.09 −83.1 −24 −0.84 −64.1 3rd 0.18 13.7 −30 0.18 13.7 0.18 13.7 4th 1.31: 99.9: Ca ii 3968.5 1st −1.92 −146.4 −76 −0.85 −64.3 56 2nd −1.09 −83.1 −20 3rd 0.17 13.0 0.17 13.0 4th 1.06: 80.1: 0.67: 50.7: Hǫ 0.19 14.4 −706 Na i 5890.0 1st −3.01 −152.8 −116 −1.43 −72.8 2nd −0.90 −45.8 −1.31 −66.7 3rd 0.28 14.3 −16 0.28 14.3 −21 0.28 14.3 Na i 5895.9 1st −2.99 −152.1 −90 −1.44 −72.8 2nd −0.90 −45.8 −1.3 −66.2 3rd 0.30 15.2 −120 0.30 15.2 0.30 15.2

Fig. A.1. The Hγ (upper plot) and Hδ (lower plot) line profiles in the spectrum of MWC480 observed with FORS1 in 2004 with an overplotted synthetic spectrum assuming Teff = 8500K and log g =3.5.

8