Arxiv:1109.3447V1 [Astro-Ph.SR] 15 Sep 2011 Ntefr Faoeamdsia Est Atr (E.G

arXiv:1109.3447v1 [astro-ph.SR] 15 Sep 2011 ntefr faoeamdsia est atr (e.g. pattern density spiral 1998; one-armed likely al. a disk, et of Vakili circumstellar form the the in in oscillations global with ylcvrain ntertobtentebu-adred- H and instance, the blue- For the of between wing quasi- ratio directly. shifted the associating radii in allowed variations stellar cyclic have few studies structure a disk interferometric of inner wave- ob- the scales unravel infrared be on to and also one visual can allowing at disks lengths, interferometry these into of with insights physics tained Po- unique and Further constrain structure disk. can structure. the 2011) inner density al. disk the et the feed Draper to to (e.g. way photosphere studies a larimetric the provide (e.g. might from which pulsations material 1998), equato- al. non-radial an et spectroscopic in revealed Rivinius Optical/infrared arranged have be observations to structure. believed disk-like circumstel- line is ionized rial which hydrogen of gas, presence with lar the associated indicating stars emission, B-type sequence) 3 4 8.-91A n ihteCAAarray. CHARA and the 084.C-0848(A) with IDs and programme 085.C-0911(A) under Observatory Paranal [email protected] umte oTeAtohsclJunlo 01 ue6 accep 6; June 2011, L on using Journal typeset Astrophysical Preprint The to Submitted lsia esasaemi-eune(rna main- near (or main-sequence are stars Be Classical etrfrHg nua eouinAtooyadDepartmen and Astronomy Resolution Angular High for Center eateto hsclSine,EbyRdl Aeronautic Embry-Riddle Sciences, Physical of Department A ITIUIN IEAIS N XIAINSRCUEIN STRUCTURE EXCITATION AND KINEMATICS, DISTRIBUTION, GAS 1 .Kraus S. ae nosrain aewt S eecpsa the at telescopes ESO with made observations on Based .Thureau N. o30m n ihhg pcrlrslto u to (up resolution spectral high with and m) 330 to ehiu nw rmrdoatooy o u rttre star target first our For few astronomy. a estimate of radio mass scales dynamical from on a known field interferometry, velocity optical technique hig in gas of time the combination first constrain The the to stars. us Be allows classical achieved two around kinematics and +-adcniumadBr and continuum H+K-band emtial hneutra ik etn ieetds oainla rotation disk different Testing ( rotation disk. Keplerian equatorial in thin geometrically a h yrgnPudsre P 42) oprn h pta orig (Br spatial Brackett the Comparing the that 14-22). (Pf find series we Pfund hydrogen the ujc headings: star. Subject Keplerian-rotat rotating the near-critical where the ou models, strong from disk a decretion predicting viscous models support compression wind scena with disk-formation inconsistent different are Discussing computation. transfer (38 rttpclB star Be prototypical tloeti ai ( radii stellocentric sn HR n LIna-nrrdsetoitreoer wit spectro-interferometry near-infrared VLTI and CHARA Using . 5 1 umte oTeAtohsclJunlo 01 ue6 acce 6; June 2011, on Journal Astrophysical The to Submitted 1 ± ..Monnier J.D. , eateto srnm,Uiest fMcia,98Denni 918 Michigan, of University Astronomy, of Department 1 ◦ 1. ,adtems ftecnrlsa (3 star central the of mass the and ), 2 5 tfle l 09 aco ta.2009; al. et Carciofi 2009; al. et Stefl ˇ ..tnBrummelaar ten T.A. , h HR ra,GogaSaeUiest,PO o 95 At 3965, Box P.O. University, State Georgia Array, CHARA The INTRODUCTION A T α 5 E ln msin( emission -line tl mltajv 11/10/09 v. emulateapj style X eateto hsc n srnm,Uiest fS.Andr St. of University Astronomy, and Physics of Department icmtla atr–sas msinln,B tr:idvda ( individual stars: – Be emission-line, stars: – matter circumstellar LSIA ESTARS BE CLASSICAL tr:fnaetlprmtr ehius interferometric techniques: – parameters fundamental stars: R cont ζ 1 .Che X. , ar n ptal eovdteBr the resolved spatially and Tauri v ( R < r ) ∝ Pf 1 γ γ r V/R .Schaefer G. , R < ,Pud(f1-7,adBle (H Balmer and 14-17), (Pf Pfund ), − ln emtywt erciia oaigselrpoopeean photosphere stellar rotating near-critical a with geometry -line 0 . 5 2 ± Br ..McAlister H.A. , variability) 0 tat,G 00-16 USA 30302-4106, GA Atlanta, γ . 1 n eieteds oiinage(140 angle position disk the derive and ) ∼ R ABSTRACT H β α lUiest,60S ld orsBv. atn ec L32 FL Beach Daytona Blvd., Morris Clyde S. 600 University, al 2 . .Touhami Y. , AI IOI AND MINORIS CANIS ,wihw a erdc iha T ieradiative line LTE an with reproduce can we which ), 5 ± fPyisadAtooy eri tt nvriy ..B P.O. University, State Georgia Astronomy, and Physics of t 0 . λ/ 2 hr a enacnrvryaotteapaac of appearance the observations spectro-interferometric about in controversy inversion a phase a recently been Furthermore, has 2003). there Rivinius de- & viscous (Porter and cretion compression, or wind pressure induced ram magnetically mecha- winds, essen- formation driven which radiatively disk is scenarios including the nism, explain different evidence to between such proposed decide been com- Obtaining have to velocity order in disk. expanding tial the an in versus ponent between rotational distinguish to purely evidence detailed a the the and field about velocity uncertainties gas cover- significant (PA) angle leaving of position age, baseline terms or in resolution limited stud- spectral were earlier spectro-interferometry most However, using ies 2010). (e.g. Carciofi to scenarios by allow review disk-formation see might several out which exhibit rule stars profile, effectively Be rotation growing classical near-Keplerian provided around have a disks series studies the Pfund These that 2011), and evidence al. 1994; 2010). 2011), et al. al. 2007, Delaa et et 2005; (Pott al. Quirrenbach al. et et (e.g. (Meilland Tycner Brackett Balmer 1998; spec- al. on the hydrogen et constraints the Vakili of for direct lines particular first in tral con- the kinematics, disk gas provided the spec- the have on in lines observations studies tral interferometric Besides geometry, 2010). tinuum al. et Schaefer 2 M ∆ e o ulcto n21,Spebr11. September 2011, on publication for ted ..Turner N.H. , λ ⊙ tdfrpbiaino 01 etme 11. September 2011, on publication for pted .A eodtre tr eosre the observed we star, target second a As ). γ o ulig n ro,M 80-00 USA 48109-1090, MI Arbor, Ann Building, son 200,w tde h a distribution gas the studied we 000), 12 = 3 msina ela ietastosfrom transitions nine as well as emission ..Gies D.R. , is ecnld htorconstraints our that conclude we rios, n iki elnse ihmaterial with replenished is disk ing sn h oiinvlct analysis position-velocity the using no h ieetln transitions, line different the of in , at,G 00-95 USA 30302-3965, GA lanta, etmti aeielnts(up lengths baseline hectometric h flwn eoiycmoet but component, velocity tflowing α β pta n pcrlresolution spectral and spatial h w,Soln,UK Scotland, ews, 2 ie rgnt rmdifferent from originate lines ) .Sturmann J. , tla ai n ooti,for obtain, to and radii stellar ai ioi,w oe the model we Minoris, Canis ζ 3 s efidta h ikis disk the that find we ws, ..Aufdenberg J.P. , TAURI H IK RUDTHE AROUND DISKS THE 1 ± 1 2 β . .Sturmann L. , 7 CMi, ◦ inclination ) 4 .Baron F. , ζ a)– Tau) 1,USA 114, d 2 x4106, ox 1 , 2 Bre ta.17;Bneue l 06 evla&Fouqué 2 & Kervella 2006; al. et Bonneau 1978; al. et (Barnes h Ddaeeso D749(0 79469 HD of diameters UD The (0 . Note of several classical Be stars, which triggered speculations 536 about secondary dynamical effects or the need for an ± oun5(P eoe h ubro onig bando t on obtained pointings of number the denotes (NP) 5 Column — . additional kinematical component beyond the canonical 0 . star+disk paradigm (Stefl et al. 2011). (0 43042 HD mas), 037 Here, we present near-infrared spectro-interferometric observations on the classical Be star β CMi with a high spectral resolution of R = 12, 000 in the hydrogen Brγ- line, enabling us to constrain the rotation profile directly. In addition, we observed the classical Be star ζ Tau for the first time in multiple hydrogen line transitions (Brγ and Pfund lines), providing direct information about the . agtDt ntuetSeta PDTTlsoeCalibrato Telescope DIT NP Spectral Instrument Date Target 62 β ζ

excitation structure within the disk. . 471 M 091-1ABRH- .7 T-T-T D71095 HD UT1-UT2-UT4 1 1 2.172 HR-K AMBER 2009-12-31 CMi a 000-1ABRM- . . T-T-T D71095 HD UT1-UT2-UT4 0.2 1 2.3 MR-K AMBER 2010-01-01 Tau In the following, we present our CHARA and VLTI ± 0 ± .

interferometric observations (Sect. 2). The observations (2 71095 HD mas), 05 0 . in spectral lines are then ﬁrst interpreted 73262 HD 73262 HD 73262 HD using a 43042 HD model- 97633 43042 S2-E2-W2 HD HD 79469, HD 43042 S1-E2-W1 HD S2-W1-W2 - - S1-E2-W1-W2 - S2-E2-W1-W1 4 71095 HD 976 HD 79469, S1-E1-W1-W2 - HD 1 25490, HD S1-E1-W1-W2 - 2 S1-E1-W1-W2 - 1 K - 2 K - 3 K 3 H35 K0-G1-A0 1 H35 CLIMB H35 CLIMB H35 2010-12-03 CLIMB 6 H35 MIRC 2010-12-02 MIRC 2010-11-30 1 MIRC 2010-12-14 MIRC 2010-11-03 2.172 HR-K MIRC 2010-11-02 2009-11-10 AMBER 2008-12-02 2010-04-23 2 a)adH 73 (0 97633 HD and mas) 027

independent photocenter analysis approach (Sect. 3). In conﬁguration [s] mode (UT) Sect. 4, we present our continuum and kinematical modeling on β CMi, followed by our discussion of the results on ζ Tau (Sect. 5). Finally, we summarize our ﬁndings in Sect. 6.

2. observations. VLTI and CHARA our of log Observation .

OBSERVATIONS 008). 009 ± 0 . 3 a) n D722(0 73262 HD and mas), 139 . 68 ± 0 esinesa.Frteclbaos easm h followin the assume we calibrators, the For star. science he . 6ms aebe siae yaeaigteetmtsfo t from estimates the averaging by estimated been have mas) 06 AL 1 TABLE . 47 ± 0 . 3ms euetedaee siae rmsacCl(Bonnea searchCal from estimates diameter the use we mas) 03

Fig. 1.— uv-coverage achieved with our CHARA/MIRC (H- band), CHARA/CLIMB (K-band), and VLTI/AMBER (K-band) interferometric observations on β CMi. r(s)

Our near-infrared H-band continuum observations 25490 HD For diameters: (UD) disk uniform g on β CMi were obtained using the CHARA array

(ten Brummelaar et al. 2005), which is operated by methods photometric independent hree

Georgia State University. The MIRC beam combiner 33 (Monnier et al. 2006b) allowed us to combine the light from four of the six CHARA 1 m-telescopes simultaneously, yielding baseline lengths of up to 330 m, with good baseline coverage (Fig. 1). The MIRC data cover the H-band with low spectral dispersion (R = 35) and was reduced using the University of Michigan MIRC data reduction pipeline (Monnier et al. 2007).

In order to investigate the K-band disk geometry of 2006). al. et u β CMi, we employed the CHARA/CLIMB 3-telescope beam combiner (Sturmann et al. 2010). Visibilities and closure phases were derived using the “redclimb” and “reduceir” software. Besides the statistical errors, we also add a calibration uncertainty of 0.05 for the derived Gas distribution, kinematics, and excitation structure in the disks around β CMi and ζ Tau 3

MR-mode covering the upper K-band (2.12 to 2.46 µm). The assumed systemic velocity for ζ Tau is +21.8 kms−1 (Duflot et al. 1995). In the ζ Tau data, we detected not only the Brγ transition (7-4), but also hydrogen Pfund transitions (Fig. 4, 1st row), including clear detection of Pf14 (19-5, 2.4477 µm) to Pf22 (27-5, 2.3591 µm). Higher Pfund transitions are also present in the spectrum, but cannot be clearly separated due to the wide, double-peaked profile of the individual lines. For ζ Tau, the derived CPs are too noisy to provide additional information and are therefore not included in our further analysis. Details about the observational setup for all interferometric observations are listed in Tab. 1. Each observation on a science star was accompanied by observations on in- Fig. 2.— Closure phases measured with the MIRC beam combiner on β CMi (blue data points), overplotted with the model terferometric calibrators, allowing us to monitor and cor- predictions from our best-fit H-band model (red data points; rect for the atmospheric and instrumental transfer func- Sect. 4.2). tion. 3. PHOTOCENTER ANALYSIS visibilities, which represents an empirical value for the typical scatter in the instrument transfer function. Differential phases provide unique information about Spectro-interferometric observations with medium small-scale (sub-mas) photocenter displacements be- (MR mode, R = 1500) and high spectral dispersion (HR tween the blue- and red-shifted line wings. These dis- mode, R = 12000) were obtained with the Very Large placements provide a very sensitive measure of the gas Telescope Interferometer (VLTI) of the European South- kinematics on scales of a few stellar radii. In a first anal- ern Observatory and the AMBER 3-telescope beam com- ysis step, we reconstruct the on-sky 2-D photocenter dis- biner instrument (Petrov et al. 2007). The AMBER ob- placement from the measured DPs by solving the system servations on β CMi were recorded using three 8.2m of linear equations unit telescopes (2009-12-31) and three 1.8 m auxiliary φ λ ~p = − i · , (1) telescopes (2010-04-23), respectively. The atmospheric 2π ~ piston was stabilized using the FINITO fringe tracker Bi (Le Bouquin et al. 2008), which allowed us to use long where φi is the differential phase measured on baseline detector integration times (DITs) of 1s and 6s and to i, B~i is the corresponding baseline vector, and λ is the record data with high spectral dispersion around the central wavelength (Le Bouquin et al. 2009). The de- vacuum hydrogen Brγ-line (λBrγ = 2.166078 µm). Due to rived photocenter plots for β CMi and ζ Tau are shown in the presence of residual phase jitter, the absolute vis- Figs. 3 and 4 (middle panel), respectively, and clearly re- ibility calibration of the AMBER data is not reliable, veal rotation-dominated velocity fields for both objects, while the important wavelength-differential observables as indicated by the linear alignment of the photocenter are not affected. Spectra and wavelength-differential vis- vectors for the different gas velocities and the opposite ibilities and phases (DPs, Fig. 3) were derived from the sign of the photocenter displacement for the blue- and AMBER data using the amdlib (V3.0) data reduction red-shifted emission. In order to associate the DP with software (Tatulli et al. 2007; Chelli et al. 2009). The the on-sky orientation, it is necessary to calibrate the wavelength calibration was done using atmospheric tel- sign of the DP measurements. For this purpose, we re- luric features close to the Brγ-line (yielding an accu- processed the ζ Tau VLTI/AMBER data set presented racy of approximately 1 spectral channel) and by ap- by Steflˇ et al. (2009) and calibrated our DP sign in order plying a heliocentric-barycentric system correction us- to match the published on-sky orientation. −1 ing heliocentric velocities of +5.97 kms (2009-12-31) For β CMi, we determine the position angles θ for the and −28.61 kms−1 (2010-04-23), respectively. For the disk rotation plane at the two epochs to be 138.3 ± 1.5◦ systemic velocity we assume +22.0 kms−1 (Duflot et al. (2009-12-31) and 141.7 ± 1.4◦ (2010-04-22). The photo- 1995). center vectors corresponding to our highest-SNR obser- From all CHARA and VLTI interferometric observa- vation (2009-12-31) show an interesting arc-like structure tions, we also derived closure phases (CPs). Both the in the red-shifted line wing (Fig. 3, bottom left), where H- and K-band continuum closure phases are consistent the photocenter vectors corresponding to low gas veloci- with zero on a 2σ-level, which leads us to conclude that ties are above the derived disk plane, while the high ve- the brightness distribution does not show significant indi- locities are displaced in the opposite direction. Although cations for deviations from centro-symmetry. The most the significance of this pattern is still only marginal in constraining CPs have been recorded with MIRC in the our data, we speculate that this pattern might result H-band, which are shown in Fig. 2. Our AMBER HR from opacity effects, with the more distant parts of the measurements from 2009-12-31 provide us also with a CP disk appearing fainter than the disk parts facing the ob- measurement in the Brγ-line of β CMi, while the derived server. Such an obscuration screen would displace the CPs from the 2010-04-23 dataset are rather noisy and photocenter perpendicular to the disk plane, where the are therefore not included for our model fits. amplitude of the displacement is stronger for lower gas ζ Tau was observed on 2010-01-01 using AMBER’s velocities, since the low-velocity emission is distributed 4

Fig. 3.— Upper panel: VLTI/AMBER spectra (1st row), visibilities (2nd row), and DPs (3rd row) measured on β CMi for the epochs 2009-12-31 (left) and 2010-04-22 (right). Middle panel: From the measured DPs, we derive for each spectral channel the 2-D photocenter displacement vector (including continuum and line contributions; East is plotted left and North is up). The ﬂux and DP data points as well as the derived photocenter vectors have been color-coded based on the Doppler velocity (see the upper panel to relate each color to a wavelength). The DPs corresponding to the astrometric solutions are shown in the middle panel as solid black lines. Bottom panel: Using the procedure outlined in Sect. 3, we corrected for the continuum contributions in the line spectral channels, revealing the photocenter displacement corresponding to the line emission only. Besides the determined photocenter oﬀsets, for comparison we also show the size of the H-band continuum-emitting disk (black ellipse, FWHM Gaussian, Sect. 4.2) and of the stellar photosphere (grey ellipse, Sect. 4.1). Gas distribution, kinematics, and excitation structure in the disks around β CMi and ζ Tau 5

Fig. 4.— Upper panel: VLTI/AMBER spectra (1st row), visibilities (2nd row), and DPs (3rd row) measured on ζ Tau for the epoch 2010-01-01 in the Brγ (left) and Pf14-22 transitions (right). Middle/bottom panel: Photocenter displacement vectors derived from the measured DPs (including line and continuum emission, middle panel) and the continuum-corrected DPs (tracing the line emission only, bottom panel). For comparison, we also show the size of the H-band continuum-emitting disk, as determined by Schaefer et al. (2010) for epoch 2009-11-10 (black ellipse, FWHM Gaussian) and of the stellar photosphere (grey ellipse). 6 over a more extended region. Accordingly, the displace- (Tycner et al. 2005). Ground-based photometric moni- ment would be strongest at zero velocities, and then sym- toring has provided no clear indications for significant metrically decrease towards higher velocities. The super- variability (Pavlovski et al. 1997), while there is some position of this weak displacement (perpendicular to the marginal evidence for long-term variations in the Hα- disk plane) with the displacement due to Keplerian ro- profile (Pollmann 2002; Hesselbach 2009). tation (parallel to the disk plane) might result in the In the following, we present a refined model for the observed arc-shaped structure, which would also provide photospheric emission of β CMi taking the near-critical a unique tool to determine the orientation of the disk in stellar rotation into account (Sect. 4.1), followed by our space and the disk rotation sense. In the case of β CMi, modeling of the interferometric data in continuum emis- this implies that the north-eastern part of the disk is fac- sion (Sect. 4.2). Given that the disk rotation signatures ing towards the observer (based on the displacement of in the Brγ emission line are overlayed on the rotation the low-velocity channels in this direction) and that the signatures of the star (in Brγ absorption), we investigate disk is in clockwise rotation (based on the location of the influence of the stellar rotation on our measurements the red-shifted photocenter displacements in the north- (Sect. 4.3). In the following, we construct a position- western quadrant). velocity diagram (Sect. 4.4) and present a full kinemati- Another intriguing feature in the β CMi data is the cal modeling in Sect. 4.5. W -shaped profile, which we observe in the wavelength- dependent visibilities and DP on our longest interfero- 4.1. Constraining the stellar parameters metric baselines (Fig. 3, 2nd and 3rd row). Likely, this In order to obtain a model for the photospheric emis- profile indicates that the visibility function of the line- sion of β CMi, we searched for a consistent set of stellar emitting region passes through a visibility null and tran- parameters, taking the evidence for near-critical rota- sits from the first to the second visibility lobe. Since tion into account (Saio et al. 2007). For this purpose, we this effect would reverse the direction of the photocen- employed our rapid rotator code (Monnier et al. 2007; ter vector, we manually correct the phase sign in these Che et al. 2011), which simulates the stellar oblateness corresponding spectral channels close to the line center. and surface temperature distribution using the modi- For ζ Tau, we can determine the rotation axis for the fied von Zeipel theorem (gravity darkening coefficient = Brγ (123.9 ± 2.0◦) and the nine Pf14-Pf22 transitions ◦ 0.188), and computes model images in the continuum (127.6 ± 1.5 ) separately and find that the line-emitting and in the photospheric Brγ absorption line. gas rotates in the same disk plane within the observa- ◦ In order to constrain the stellar parameters, we com- tional uncertainties of ∼ 2 . puted a small parameter grid, in which we systemati- Using the aforementioned procedure, it is possible to cally varied the inclination angle i and the fractional reliably measure the direction of the photocenter dis- angular velocity ω/ωcrit (where ωcrit is the critical an- placement, while the length of the displacement vector gular velocity). For a given inclination angle and frac- is biased by the contributions from the underlying con- tional angular velocity, we start our iterative process tinuum emission. In order to remove these contributions by assuming a stellar mass M , from which we com- from the measured observables (F , V , φ), we apply the ⋆ pute the polar radius Rpole using the rotation velocity method outlined by Weigelt et al. (2007) and interpolate − v sin i = 244 ± 6 kms 1 (Yudin 2001). The polar tem- the continuum flux (F ) and continuum visibility (F ) c c perature T is derived from the V -band magnitude from the adjacent continuum. The visibility and DP of pole (V = 2.89). Any disk continuum flux in the optical the pure line emitting-region (V , φ ) are then given by l l band will tend to dilute the photospheric spectrum caus- 2 2 2 |FlVl| = |F V | + |FcVc| − 2 · F V · FcVc · cos φ (2) ing absorption lines to appear weaker than expected. We checked for such line reduction in blue spectra of |F V | sin φl = sin φ , (3) |FlVl| where Fl = F − Fc denotes the flux contribution from the spectral line. The continuum-corrected DPs are then used to derive the photocenter displacement of the pure line-emitting region (Figs. 3 and 4, bottom panel). Ap- plying this correction will provide the real centroid offset of the line emission, but also introduce noise from the visibility and flux measurements, resulting in an increased scatter in the position angle distribution. There- fore, we decided to measure the position angle of the rotation axis from the uncorrected line+continuum photocenter displacements (Figs. 3 and 4, middle panel), while the continuum-corrected photocenter displacements (bottom panel) will be used later on to construct a position- velocity diagram and to compare the spatial origin in different line tracers.

4. DISCUSSION ON β CANIS MINORIS Fig. 5.— Spectrum of the He I 447.1 nm and Mg II 448.1 nm line of β CMi (solid line). The dashed line shows a photospheric model β CMi (HR 2845) is a relatively quiet B8V-type clas- without disk contributions, while the dotted line corresponds to a +2.4 model where 20% disk continuum contributions have been added. sical Be star located at a distance of 52.2−2.2 pc Gas distribution, kinematics, and excitation structure in the disks around β CMi and ζ Tau 7

TABLE 2 Grid of stellar parameters, obtained for β CMi using our fast rotator model (Sect. 4.1).

app app Model i ω/ωcrit M⋆ Rpole Req Tpole P L T ◦ eﬀ [ ][M⊙][R⊙][R⊙] [K] [days] [L⊙] [K] A 45 0.99 3.70 3.32 4.61 13200 0.676 267 11558 B 45 0.96 4.11 3.03 3.94 14600 0.577 343 13085 C 45 0.93 4.47 2.87 3.59 15650 0.527 413 14235 D 40 0.99 3.98 2.95 4.10 14550 0.546 337 12871 E 40 0.96 4.49 2.74 3.56 16050 0.475 435 14505 F 40 0.93 4.91 2.60 3.25 17300 0.433 532 15847 G 35 0.99 4.43 2.62 3.64 16200 0.433 437 14447 H 35 0.96 5.03 2.44 3.17 18100 0.377 590 16474

β CMi obtained by Grundstrom (2007). Fig. 5 shows with the observed spectrum. This comparison suggests the observed profiles (solid line) of He I 447.1 nm and that the disk continuum emission in the blue range is Mg II 448.1nm based upon the average of 11 spectra negligible, and thus we will ignore any disk contribution made between 2005 and 2009 with the Kitt Peak Na- to the optical flux in the following discussion. tional Observatory Coude Feed Telescope (spectral reso- With the derived polar radius and bolometric lumi- lution R = 12, 500). The observed spectrum is compared nosity, we are able to locate the stellar position on the with a theoretical flux spectrum (dashed line) from the HR diagram, after correcting for the rotational effect grid of LTE models by Rodr´ıguez-Merino et al. (2005). (Che et al. 2011). The stellar mass from the HR diagram We adopted average parameters for the visible hemi- is used as initial value for the next iteration step, until sphere of the star of Teff = 11, 800 K, log g = 3.8, and convergence between the assumed mass and the mass es- V sin i = 230km s−1 (Frémat et al. 2005), and the model timated from the HR diagram is reached. spectrum was convolved with a simple rotational broad- For each parameter combination, our model provides ening function for a linear limb darkening coefficient of the rotation period, apparent luminosity Lapp, and ap- app ǫ = 0.42 (Wade & Rucinski 1985). Note that this apparent effective temperature Teff , as listed in Tab. 2. proach assumes a spherical star, ignores gravity darken- Comparing the model effective temperature and luminos- ing, and neglects changes in the local intensity spectrum ity with the observational constraints for β CMi (Teff = with orientation, but these simplifications are acceptable 12, 050 K; L = 195 ± 60 L⊙; Saio et al. 2007), we find for our purpose of checking for systematic line depth dif- good agreement for models A, B, and D. Due to its con- ferences. We see that the match of the Mg II λ4481 pro- sistency with the inclination and effective temperature file is satisfactory, but the wings of the He I λ4471 line (Sects. 4.2 and 4.5), we favor model D, suggesting i = ◦ −1 appear too high, probably due to low level emission in 40 , ω/ωcrit =0.99, Tpole = 14, 550, v sin i = 244 kms , this transition. Adding a 20% disk flux contribution and app Teff = 12871 K, and a polar and equatorial radius of renormalizing the spectrum yields a diluted version (dot- Rpole = 2.95 R⊙ = 0.26 mas and Req = 4.10 R⊙ = ted line), which appears to be much too weak compared 0.36 mas, respectively. We plot the values corresponding to our model grid in the HR diagram shown in Fig. 6 and compare it to the evolutionary tracks from Yi et al. (2003) and Demarque et al. (2004), yielding an evolutionary age of ∼ 0.1 − 0.15 Gyr. Obviously, an important input parameter for our modeling procedure is the rotation velocity v sin i, for which we use an average value from the literature. Townsend et al. (2004) argued that most v sin i measurements on Be stars might systematically underestimate the true projected rotation value due to the effect of gravity darkening. In order to test this scenario, we have artificially increased the measured v sin i-value by 10% and repeated our modeling procedure. We find that the increase in v sin i results in a significant increase in the stellar mass, apparent effective temperature, and apparent luminosity, making the model prediction much less consistent with the observed values. Therefore, we suggest that in the case of β CMi, the v sin i-measurements Fig. 6.— Hertzsprung-Russell diagram, including the solutions are not significantly biased by gravity darkening, which from our small rapid rotator grid for β CMi (Tab. 2). The num- is likely a result of the intermediate inclination angle bers on the left side give the stellar mass of the corresponding evolution track, while the numbers at the bottom gives the evolu- of β CMi (the effect would be stronger for an equator- tionary age corresponding to the shown isochrones Yi et al. (2003); on viewing angle) and the use of our moderate gravity 0.188 Demarque et al. (2004). As outlined in Sect. 4.1, only models A, B, darkening law (Teff ∝ g , Che et al. 2011), which is and D satisfy the observational constraints on the stellar effective eff temperature. Furthermore, model D satisfies best our inclination likely more relatistic than the classical von Zeipel darken- ◦ 0.25 obtained with interferometry (38.5 ± 1 ). ing coefficient (Teff ∝ geff ) adopted by Townsend et al. 8 oes egv h igmjrai imtr ( diameter. axis major ring the give we models, of Note (2004). Clearly, the accurate determination of v sin i re- i mains a fundamental problem for constraining the stellar ( orientation. disk face-on a denotes 0 = parameters of near-critical rotating stars. Therefore, we ( — . are currently working on incorporating the computation of model spectra in our modelling procedure, which will a then be fitted to observed spectra together with the pho- maximum half half-width the give we model, Gaussian the For ) ailitniyindex intensity Radial velocity Rotation radius disk Outer oiinangle Position size axis Major ratio Flux Inclination tometric and interferometric constraints (Che et al., in index law Rotation mass Stellar prep.). The modeling provides spectra as well as model images in the continuum emission and in hydrogen photospheric absorption lines, which we will use as a representation

of the stellar brightness distribution for our modeling of ( e the disk in the following sections. ) 4.2. Continuum disk geometry e f hogotti ae,alpsto nlsaemaue al measured are angles position all paper, this Throughout ) F With projected baseline lengths up to B = 314 m, kep disk c R M hsprmtrhsbe etﬁe hogottemdln p modeling the throughout ﬁxed kept been has parameter This )

our CHARA/MIRC interferometric observations allow χ 1AU) (1 β out a q θ i r 2 ⋆ us to constrain the disk geometry in eight spectral /F channels in the H-band with an eﬀective resolution ⋆ of λ/2B = 0.5 mas. A commonly applied model- [mas] [mas] [ M [ ing approach is to approximate the disk emission with [ ◦ ◦ ⊙ ] ]

a Gaussian intensity proﬁle, which is superposed on ]

a uniform disk representing the stellar surface (e.g. of disk the for results fitting Model Rn oe)(asinmodel) (Gaussian model) (Ring 1 . 138 0 19

Tycner et al. 2005; Meilland et al. 2009; Pott et al. 2010; 39 . 16 +0 − .71.77 1.77 . .

Schaefer et al. 2010). A caveat of this approach is that 9 2 +0 − – – – – – – – – – – 0 +7 − . . +3 − 0 41 24 a signiﬁcant part of the disk emission is concentrated 6 . . 7 . . 07 03 . 2 . 4 ( 1 8

at radii r < R , and therefore wrongly attributed to the b ⋆ Continuum ) disk instead of the stellar flux. Accordingly, these models H will systematically overestimate the disk emission with -band HH)o h afGusa nest rfie esrdalo measured profile, intensity half-Gaussian the of (HWHM) respect to the stellar emission (F /F ), resulting in in- 0 AL 3 TABLE disk ⋆ . 139 0 33 40 . consistencies with SED fitting results. In order to avoid 18 +0 − . . 1 5 0 +0 − +5 − . . +4 these biases, we employ an elliptical Gaussian model, − 18 14 0 9 . . 6 . . 05 03 . . 9 where the radial intensity profile at r > R is given by 2 (

⋆ 4 3 a a half-Gaussian. resulting in a proper estimation of the ) Fdisk/F⋆ ratio. As alternative model, we considered an

inclined ring model, where the radial intensity proﬁle is model) (Gaussian given by a Gaussian centered at radius a from the star, Continuum n h ikmjrai n ato North. of East and axis major disk the ong 0 0 139 K . . 40 33 with a ﬁxed fractional width of 25% (see Monnier et al. 20 β 1.17 -band . ( CMi. +0 − – – – – – . 1

2006a for details). Based on radiative transfer simul- a 5 ( 0 ) ( c , oeue(e et.42ad45frdtis.( details). for 4.5 and 4.2 Sects. (see rocedure . . c ( ) 15 tations of the brightness distribution in classical Be star 15 ) c disks (e.g. Waters 1986), we consider the Gaussian model ) a better representation of the expected emission proﬁle in a Be disk, while the ring model is more suited for compar- Kpeinrt)(rertto law) rotation (Free rot.) (Keplerian

ison with results from the literature (e.g. Meilland et al. − 3 5 38 136 − 1 . . 1 5 8

2006). . .61.66 1.66 0 . 6 . 5 0 ± ± . – – – – . 5 ± For all models, the stellar photosphere is represented ( 8 ± c ( ( 0 0 c ) c 0 . . ) 38 1 ) by the aforementioned rapid rotator model (Sect. 4.1) 3 5 2 2 . 2 and the position angle θ of the stellar equator and the Brackett

disk major axis are aligned. Besides the disk major axis Line a (ring major radius or Gaussian half-width half maximum), the PA θ, and the inclination i, we also treat γ − gtemjrai.( axis. major the ng 1 136 − 0 .

the flux ratio between the stellar and disk component 0 . 1 5 . . . 5 8 ± . . 5 6 ± 8 (Fdisk/F⋆)H as a free parameter. The model was fitted to ( ( ( c c ( ( 0 c c ) ) c 0 ) . ) ) the MIRC squared visibilities, closure phases, and triple 1 . amplitudes using a least-square fitting procedure, result- 1 ing in the best-fit model shown in Fig. 7. The best-fit d model parameters are displayed in Tab. 3, including 1σ- angle inclination An ) errors, which we have determined using a boot strapping b procedure. We also tested whether the agreement can ring the For ) be improved using a skewed ring model (Monnier et al. 2006a), but find that introducing any disk asymmetry only marginally improves the fit. Our model includes minor asymmetries due to the brightened polar region in our rapid rotator model, which results in small clo- Gas distribution, kinematics, and excitation structure in the disks around β CMi and ζ Tau 9

Fig. 7.— Left: Continuum visibilities measured on β CMi with CHARA/MIRC in the H-band (red points) and CHARA/CLIMB in the K-band (blue points). In order to show the position-angle dependence of the visibility function, we have binned the data by the projected baseline PA (number in lower left corner of each panel). For the H-band, the best-fit ring model (red line) and Gaussian model (green dotted line) are shown, with both models yielding nearly indistinguishable visibility profiles. For K-band, we show the best-fit Gaussian model (blue dashed line; see Sect. 4.2). Right: Intensity distribution corresponding to our best-fit ring (top) and Gaussian model (bottom). sure phases (. 1.2◦), consistent with the measurement which is tracing the stellar rotation and is underlying the (Fig. 2). Brγ-line emission. The determined disk position angles of 139.2◦ (ring In order to simulate the effect of stellar rotation on our model) and 139.3◦ (Gaussian model) are in excellent interferometric observables, we employ our rapid rotator agreement with the gas disk rotation axis determined code and compute the stellar surface brightness distribu- with our VLTI/AMBER photocenter analysis (140.0 ± tion around the Brγ-line with a similar resolution as our 1.7◦). AMBER observations (Fig. 8, top). For this, we align Using the detailed information about the H-band ge- the stellar rotation axis with the measured rotation axis ometry obtained with our extensive MIRC data set, we (θ = 227.5◦, Sect. 3). Brγ absorption on one side of then fitted our Gaussian model to the CLIMB data the photosphere will shift the photocenter towards the to determine the K-band continuum geometry. Given opposite direction. Taking this into account, we adjust the lower amount of observational constraints for the the orientation of our model photosphere, so that the K-band, we treated only the disk-to-star flux ratio Brγ absorption in the blue-shifted line wing matches the (Fdisk/F⋆)K as a free parameters and kept the remaining measured direction of the photocenter displacement at parameters fixed. The resulting best-fit visibility curves red-shifted wavelengths. In addition to the photospheric +0.09 are shown in Fig. 7. With (Fdisk/F⋆)K = 0.25−0.08, we emission, we include disk emission in our model, assum- do not find any significant deviations between the K- ing (Fdisk/F⋆)K =0.25, as determined in Sect. 4.2. band and H-band flux ratio. From the model, we compute the corresponding DPs and find signatures φ . 0.5◦ on all baselines (Fig. 8, bot- 4.3. Investigating the differential phase signatures of tom). This result reflects the fact that the fraction of the the stellar rotation in the Brγ absorption line photospheric absorption to the total K-band emission is rather small (< 8%) and that the absorption appears With the currently achievable accuracy, DP measure- as a rather diffuse structure on the extended stellar sur- ments can already reveal photocenter displacements two face. We conclude that the photocenter signatures due to three orders smaller than the formal angular resolution to stellar rotation (in photospheric absorption lines) are (λ/2B). For instance, our 2009-12-31 AMBER observa- ◦ a minor effect compared to the rotation signatures of the tions exhibit a DP accuracy of ∼ 0.8 (standard devia- disk (in line emission), but nevertheless result in DP sig- tion over all continuum spectral channels), corresponding natures comparable to the currently achievable DP accu- to a photocenter displacement of ∼ 8 micro-arcsecond racy and are therefore important, in particular for more (0.008 mas) on the employed 115 m baseline. Given that detailed future studies. this is much smaller than the equatorial stellar diameter (∼ 660 micro-arcsecond), it is important to investigate 4.4. Constraints on the disk kinematics from the whether the measured DP signatures might also contain position-velocity diagram contributions from the photospheric Brγ absorption line, 10

Fig. 8.— Top: Model images from our rapid rotator model for β CMi for some representative wavelengths (Sect. 4.3). Bottom: Model diﬀerential phases (including the stellar photosphere only) computed for the baselines employed by our β CMi VLTI measurements from 2009-12-31 (left) and 2010-04-23 (right).

Position-velocity (p − v) diagrams provide a powerful displacements) trace the light barycenter of each annu- tool for the interpretation of a disk velocity field and are lus, corresponding to the B-C and A-D curves in Fig. 9 commonly employed in radio astronomy to derive the (for more details see Weintroub et al. 2008). We fit this rotation profile of circumnuclear galactic disks or proto- simple model to the β CMi p − v diagram assuming a stellar disks using, for instance, maser (e.g. Miyoshi et al. fixed inclination of i = 40◦ (as determined in Sect. 4.2), 1995; Pestalozzi et al. 2009) or molecular emission trac- which yields best agreement for M⋆ = 3.6 ± 0.3 M⊙, ers (e.g. Sofue & Rubin 2001; Isella et al. 2007). Using Rout =2.4 ± 0.2 mas, and Rin < 1.0 mas (solid curve in our continuum-corrected photocenter displacement mea- Fig. 9). surements, we can construct an equivalent diagram from our VLTI interferometric data by measuring the length of 4.5. Detailed modeling of the disk kinematics the continuum-corrected photocenter displacement vec- By reducing the information content to purely astro- tors projected on the disk plane θ = 140.0◦. Performing metric data, the p−v diagram analysis method presented this projection on the disk plane also allows us to avoid in the last section provides a very intuitive and power- a potential bias due to opacity effects, since these effects ful method to constrain the gas velocity field in the ob- would move the photocenter only perpendicular to the served line tracer. In this section, we will use a more disk plane (Sect. 3). The resulting p−v diagram (Fig. 9) sophisticated model in order to make full use of the rich shows a symmetric rotation curve, which we interpret information contained in our spectro-interferometric ob- using the simple model of a thin Keplerian-rotating disk servations, including spectra, wavelength-differential vis- (Weintroub et al. 2008). The disk extends from an in- ibilities, DPs, and CPs. The Brγ-emitting gas is assumed ner (Rin) to an outer (Rout) radius and the line-of-sight to be located in a thin disk plane, which is justified due (LOS) velocity v of the emission element located at ra- to the expected small opening angle of Be star disks (e.g. dius r and at angle ϑ in the disk plane is given by Bjorkman & Cassinelli 1993) and the intermediate inclination angle under which β CMi is observed (Sect. 4.2). GM⋆ The gas kinematics is parameterized with a rotation v (r, ϑ)= sin ϑ · cos i, (4) β kep r r profile |~v(r)| = fkep(Rref ) · (r/Rref ) , where β = −1 for rotation with constant angular momentum, β = −0.5 where ϑ is measured against the LOS and G is the for Keplerian rotation, β = 0 for constant rotation, gravitational constant. Considering only the emission and β = +1 for solid body rotation (e.g. Stee 1996). from the outer-most disk annulus (i.e. at Rout and ϑ = fkep = |~v(Rref )|/|~vkep(Rref )| is the orbital velocity at the −π...+π) will result in a straight line in the p−v diagram reference radius Rref = 1 AU, expressed as fraction of −1/2 (line A-B in Fig. 9). Adding the emission from the re- the Keplerian velocity |~vkep(r)| = (GM⋆/r) . In the maining disk annuli (Rin to Rout) will result in a charac- model, the line emission extends from the stellar surface teristic “bowtie”-shaped filled region in the p-v diagram, to an outer radius Rout with a radial power-law inten- q such as commonly observed in CO imaging observations sity profile Il(r) ∝ r . To include thermal line broad- (Sofue & Rubin 2001; Isella et al. 2007). However, astro- ening in our kinematic model, we adopt a constant ra- metric observations (such as radio masers or photocenter dial gas temperature of 0.6 · Teff (=7860K for β CMi), Gas distribution, kinematics, and excitation structure in the disks around β CMi and ζ Tau 11

Fig. 9.— Left: Position-velocity diagram, derived from the continuum-corrected photocenter displacement vectors in Fig. 3 (bottom, left). The color of the data points matches the color-coding of the wavelength channels in Fig. 3. The model curves show the Keplerian rotation profile for three stellar masses with Rin = 0.68 mas, and Rout = 2.4 mas (Sect. 4.4). Right: In order to illustrate the interpretation of the p − v-diagram, we show the Brγ-line model images from our best-fit kinematical model (Sect. 4.5, Fig. 10) for some representative velocities and mark the corresponding photocenter displacement (corresponding to the position offset plotted on the abscissa in the diagram) in each channel map with a blue cross. The yellow vertical line marks the position of the star, which corresponds to the photocenter of the continuum emission and the zero position in the p − v diagram. as suggested by the radiative transfer modeling from Fig. 11 and the uncertainties have been derived using the Carciofi & Bjorkman (2006). bootstrapping technique. We include both line and continuum emission in our Our model fits show that the intriging phase inversion model and compute the interferometric observables for observed at the line center on our longest VLTI base- our given VLTI array configurations and the covered lines (Fig. 10, 3rd row) can be explained with the phase wavelength channels. The wavelength-dependent model jumps in the Fourier phase crossing a visibility null. visibilities and phases are then fitted to the VLTI inter- These phase jumps appear in the same spectral chan- 2 ferometric data using a reduced χr goodness-of-fit esti- nels where we measure visibility minima in the W -shaped mator (see Kraus et al. 2009 for a definition), including visibility profile, indicating that the visibility function of our flux, visibility, DP, and CP constraints. the pure line-emitting geometry transits here from the For the continuum emission, we assume the geometry first to the second visibility lobe (Fig. 10, 2nd and 3rd determined with our CHARA observations (Sect. 4.2) row). At the same time, the continuum emission is only and use this model to renormalize the AMBER contin- marginally resolved, which results in the measured com- uum visibilities. In order to incorporate the underly- posite line+continuum visibility with a rather high con- ing photospheric Brγ absorption, we include the photo- trast of ∼ 0.7. The phase jumps and visibility mini- sphere model discussed in Sect. 4.3, although, as dis- mums are a basic property of the Fourier transform of cussed above, the influence on the differential phase strongly resolved objects and are reproduced very natu- (. 0.5◦ on all baselines) is still within the measurement rally and without finetuning from our kinematical mod- uncertainities. eling. Therefore, our results do not support the idea The disk position angle is fixed to θ = 140.0◦, as de- outlined by Stefl et al. (2011) that these features might termined by our model-independent photocenter analysis indicate secondary dynamical effects or polar mass out- (Sect. 3). The remaining six parameters in our modeling flows. are the outer disk radius Rout, the inclination i, the stel- As best-fit value for the radial intensity index, we yield lar mass M⋆, the radial intensity power-law index q, the q = −1.6±0.2. Isothermal viscous decretion disk models, disk rotation index β, and the velocity at the reference such as discussed in Bjorkman & Carciofi (2005), pre- −2 radius Rref (expressed as fraction of the Keplerian veloc- dict a radial disk surface density law Σ(r) ∝ r , cor- ity, fkep). In a first step, we fix the velocity profile to responding to q = −2 in the optically thin case. More Keplerian rotation (β = −0.5, fkep = 1) and vary the re- realistic non-LTE disk models suggest a more shallow maining four parameters systematically on a grid, yield- surface density profile (Σ(r) ∝ r−2...−1) due to a steep ing the best-fit model shown in Fig. 10. The correspond- temperature drop in the inner few stellar radii of the ing parameters and uncertainties are listed in column 4 disk Carciofi & Bjorkman (2008). Therefore, we con- of Tab. 3. In a second step, we test also non-Keplerian clude that our measured intensity profile is well consis- velocity fields, yielding the best-fit values listed in col- tent with these models. umn 5 (Fig. 12). The resulting χ2-surfaces are shown in The most significant deviation of our simple kinemat- 12

Fig. 10.— Upper panel: Comparison of the VLTI/AMBER β CMi spectra (1st row), visibilities (2nd row), DPs (3rd row), and CPs (4th row) with our Keplerian disk model (solid lines). The three diﬀerent colors correspond to the three baselines, where the associated projected baseline lengths and PAs are labeled in the 2nd row. Lower panel: To illustrate our kinematical model, we show intensity channel maps for ﬁve representative wavelengths. Gas distribution, kinematics, and excitation structure in the disks around β CMi and ζ Tau 13

Fig. 11.— χ2-surface derived from our β CMi model grid (Sect. 4.5) around the best-ﬁt solution assuming a Keplerian velocity ﬁeld (Tab. 3).

Hα emission line (Rivinius et al. 2006; Steflˇ et al. 2009). These V/R variations exhibit a period of ∼1429 days and are generally attributed to a global one-armed density oscillation in a Keplerian (or nearly Keplerian) disk (Okazaki 1991). A kinematic model for the Brγ-line of some earlier AMBER MR observations on ζ Tau was presented by Steflˇ et al. (2009) and Carciofi et al. (2009) and suggested a Keplerian rotation profile. Our spectro- interferometric observations on ζ Tau provide a Brγ measurement at a new epoch in the Hα V/R-phase and cover, for the first time, also the wavelength region around the hydrogen Pfund lines. Besides our Brγ and Pfund series Fig. 12.— χ2-surface for the rotation law index β derived from our β CMi model grid including also non-Keplerian velocity fields data, we include Hα-sizes from the literature, namely the (Sect. 4.5). GI2T photocenter measurement by Vakili et al. (1998, ∼ 7 R⋆ or ∼ 1.33 mas) and the Mark III and NPOI ical model from the data is in the precise shape of the measurements by Quirrenbach et al. (1994, 1997) and Brγ-spectrum. For instance, the measured spectra show Tycner et al. (2004). It is important to note that the a weak emission component in the line center (Fig. 10), Mark III and NPOI results were based on visibility am- which is not reproduced by the model and which might be plitudes instead of differential phases. Also, a direct com- related either to the presence of a uniformly distributed parison with these earlier observations is complicated by low-velocity gas halo or, more likely, to radiative transfer the known V/R variability of ζ Tau, although we note effects. that our measurement at phase 0.484 is reasonably close Another deviation between the model and the observa- to the Hα-measurement at phase 0.577, presented by tion concerns the weak V/R-asymmetry which can be ob- Tycner et al. (2004). We convert the Gaussian FWHM served in our VLTI/AMBER Brγ spectra at both epochs derived by these studies to photocenter displacements (Fig. 3, top row) and which is not reproduced by our ax- by computing the centroid of the corresponding Gaus- ialsymmetric kinematical model. Within the 113 days sian brightness distributions. We find that Brγ emerges covered by our AMBER observations, no changes in the V/R asymmetry could be observed, which is consistent with the conclusion by Tycner et al. (2005) and Jones et al. (2011) that the Hα profile does not show significant long-term variability. It is interesting to compare the average photocenter displacement between the blue- and red-shifted line emission (2.5 mas) with the characteristic size of the Hα-emitting region (2.13 mas, Tycner et al. 2005), which suggests that Brγ emerges from a similar spatial region as Hα and a much more extended region than the H- and K-band continuum emission (RBrγ ∼ RHα > Rcont, Fig. 3), in agreement with the prediction from Carciofi (2010). However, a detailed comparison is difficult due to the non-simultaneity of the different observations.

5. DISCUSSION ON ζ TAURI Fig. 13.— ζ Tau disk position angle, plotted as function of the Hα V/R phase, including data presented in Schaefer et al. (2010) ζ Tau is a particularly well-studied classical Be star and in this study. Our new phase measurement (phase 0.484) which shows a cyclic variability in the flux ratio of is in good agreement with the disk precession law proposed by the violet- and red-shifted wing of the double-peaked Schaefer et al. 14 from a similar spatial region than Hα, but a more ex- rial is located in an equatorial disk which extends out- tended region than the near-infrared continuum emis- wards from the stellar radius with a constant vertical sion, as already found for β CMi (Sect. 4.5). The Pfund density per unit volume and a half-opening angle of emission originates from intermediate stellocentric radii Θ=5◦. Then, we integrate for each radius r the op- (RHα ∼ RBrγ > RPf > Rcont; Fig. 4). tical depth τν in vertical direction, assuming hydrogen under LTE conditions and an isothermal temperature 5.1. Signatures of the known one-armed oscillation distributions. The number density for the different ex- Both in Brγ and in the Pf14-Pf22 lines, we detect a citation levels and ionization stages is computed using double-peaked line profile and clear rotation signatures the Saha and Boltzmann equation assuming a Gaussian in the differential phases. All line transitions exhibit a line profile with thermal line broadening (Wilson et al. photocenter displacement with a stronger amplitude in 2009). The Einstein coefficients for the different hy- the south-eastern (red-shifted) than in the north-western drogen transitions are estimated using the series expan- (blue-shifted) lobe (Fig. 4, middle panel). Such an asym- sion published by Omidvar & McAllister (1995). The radial temperature and surface density profile are param- metric displacement is consistent with the presence of −0.5 ˇ eterized with T (r) = Teff (r/R⋆) (Stee & de Araujo a one-armed oscillation in the disk (Stefl et al. 2009; −2 −2 Carciofi et al. 2009). CHARA/MIRC observations by 1994) and Σ(r)=Σ0(r/R⋆) with Σ0 = 2.1 gcm Schaefer et al. (2010) showed that the one-armed oscilla- (Carciofi et al. 2009). For the stellar temperature, equation pattern can also be observed as an asymmetry in the torial stellar radius, and distance, we assume Teff = H-band continuum emission. Using multi-epoch data, 19 370 K, R⋆ = 7.7 R⊙, and d = 126 pc, respectively they find that the position angle of the asymmetry is cor- (Carciofi et al. 2009). related with the spectroscopic variability and precesses Using the radiative transfer equation for LTE condi- −τν(r) around the star with the Hα V/R period. Our observations (Iν (r)= Bν (T )(1−e )), we compute the emit- tion (2010-01-01) adds an additional epoch for this analy- ted intensity per unit area as function of radius r in sis at the V/R phase 0.484 (assuming maximum phase at the disk, where Bν is the Planck spectrum for temper- JD=2454505.0 and a period of 1429 days, Schaefer et al. ature T at frequency ν. From the radial intensity pro- 2010). We find that the spiral density maximum is lo- files (Fig. 14) we compute the centroid of the brightness cated towards the south-east of the central star and the distribution, which is then compared to the stellocentric measured PA is consistent with the sinusoidal PA mod- emission radius measured in the different line transitions. ulation proposed by Schaefer et al. (2010), as shown in Given that our model assumes a simplified vertical den- Fig. 13. sity structure and does not include inclination effects, The reported polarization angle for ζ Tau (32.9−33.1◦, we do not aim to match the absolute sizes of the emit- Quirrenbach et al. 1997; 35 ± 4◦, Ghosh et al. 1999; 32 − ting region in all line transitions, but focus instead on 33.5◦, McDavid 1999) is in excellent agreement with the the relative sizes. For this step, we normalize both the rotation axis position angle (35.8 ± 2.0) found by our measured and the model photocenter offsets relative to spectro-interferometric observations. Brγ. The comparison between the Hα, Brγ, and Pf14-17 relative sizes and our LTE model is shown in Fig. 15. 5.2. Probing the disk excitation structure using We find that we can reproduce important observational multi-transition spectro-interferometry features, in particular that In order to test whether the measured differences in the stellocentric emitting radius of the Brγ and Pfund transitions are consistent with the expected excitation structure in the disk, we construct a simple radiative transfer model assuming local thermodynamic equilib- rium (LTE). We assume that the line-emitting mate-

Fig. 15.— Origin of the Brγ and Pfund line emission from ζ Tau, as measured by the length of the photocenter displacement vector, averaged over all velocity channels for a given line transition (the error bars represent the standard deviation in the different velocity channels). In addition, we include the Hα-line photocenter displacement by Vakili et al. (1998, black data point) and the Hα size estimates from Quirrenbach et al. (1994), Quirrenbach et al. Fig. 14.— Radial intensity profile, as computed with our LTE (1997), Tycner et al. (2004), which we converted from Gaussian model for different line transitions (Sect. 5.2). The intensity has FWHM to the corresponding photocenters (grey data point). The been weighted by the emitting area and normalized to the peak dashed line shows the prediction from our LTE excitation model intensity. (Sect. 5.2). Gas distribution, kinematics, and excitation structure in the disks around β CMi and ζ Tau 15

a) Brγ originates from a similar spatial region in the is in excellent agreement with earlier spectroscopic stud- disk as Hα (RBrγ ∼ RHα). This result is also in ies (Saio et al. 2007). The inclination of the system is agreement with the predictions from more sophis- 38.5 ± 1◦, as determined with our CHARA continuum ticated non-LTE radiative transfer computations, and VLTI line observations. As shown with our de- such as made by Carciofi (2010). In order to tailed kinematical modeling, the rotation law is Keple- better characterize the differences between the rian (β = −0.5 ± 0.1) and we do not have to include Brγ and Hα-emitting region, contemporaneous an expanding velocity component in order to explain our observations in these two wavelength bands will data. Furthermore, our kinematical model allowed us to be required. identify the origin of the phase inversion, which has now been observed in the differential phases in five out of eight b) Brγ originates from a more extended region than Be stars. These phase jumps correspond to the transition the Pfund lines (RBrγ > RPf ). Computing the from the first to the second visibility lobe, removing the weighted average of the measurements in the indi- necessity for speculations beyond the canonical star+disk vidual transitions, we yield RPf /RBrγ =0.48±0.12. paradigm (Stefl et al. 2011). For the Pfund lines, we are not aware of earlier pre- For ζ Tau, we obtained spectro-interferometric obser- dictions obtained with non-LTE radiative transfer vations covering simultaneously the Brγ and at least nine codes. transitions from the Pfund line series. For all transitions, we detect a significantly stronger photocenter dis- Our observational results confirm the finding from placement in the red-shifted line wing than in the blue- Pott et al. (2010), which found R > R on the classi- Brγ Pf shifted line wing, tracing the one-armed oscillation which cal Be star 48 Lib, and consolidates their suggestion that has been deduced for the ζ Tau disk before. Comparing the measured size differences can already be explained the photocenter displacement in the different line tran- with the expected optical depth differences between these sitions, we find that the Pfund, Brackett, and Balmer line transitions. We encourage theoreticans and model- lines originate from different stellocentric emitting re- ers to employ their more sophisticated non-LTE radiative gions (R < R < R ∼ R ), which we can re- transfer codes in order to test the influence of inclina- cont Pf Brγ Hα produce qualitatively with a simple LTE line radiative tion and non-LTE effects, although, based on the results transfer model. More work, including non-LTE radia- from Iwamatsu & Hirata (2008, e.g. Fig. 3), we expect tive transfer modeling, will be required in order to derive no significant departure from LTE for the inner disk re- quantitative constraints. gions and the high transitions traced by our observations. By detecting a purely Keplerian velocity field, our ob- Future multi-transition spectro-interferometric observa- servations are inconsistent with disk-formation mecha- tions with improved uv-coverage might also measure the nisms incorporating a strong outflowing velocity com- precise radial intensity profile in a model-independent ponent, such as the wind compression scenario, which fashion. Together with sophisticated radiative transfer predicts a strong radial velocity component compara- simulations, these observations will reveal the excitation ble to the escape velocity (Bjorkman & Cassinelli 1993). structure of the disk and constrain parameters such as On the other hand, our kinematical constraints, as well the temperature profile and the vertical disk structure, as the measured hydrogen line intensity profiles (with a which are currently difficult to access. radial power law index q = −1.6 ± 0.2) are consistent 6. with the predictions from Keplerian viscous decretion CONCLUSIONS disk models (Lee et al. 1991). As shown by Kato (1983) Using CHARA and VLTI interferometry, our study and discussed in various reviews (e.g. Carciofi 2010), vis- combined high angular resolution (with baseline lengths cous Keplerian disks are also able to produce one-armed up to 330 m) with kinematical information obtained at density oscillations, such as detected for ζ Tau. high spectral dispersion, yielding direct observational Considering that our kinematic constraints have been constraints on the gas distribution, excitation structure, obtained using a very limited number of individual mea- and kinematics of the disks around two classical Be stars. surements (2 pointings on β CMi, 1 pointing on ζ Tau), Using a model-independent photocenter analysis method our study also illustrates the high effectiveness achievable we derived the disk rotation axis for the prototypical ob- with spectro-interferometry, in particular if a very high jects β CMi and ζ Tau and spatially and spectrally re- spectral resolution is employed or several line transitions solved the disk rotation profile on scales of a few stellar are observed. radii. For both objects, we find that the determined gas rotation plane agrees well with the orientation of the continuum-emitting disk as resolved by CHARA, al- We thank A. Carciofi for helpful discussions on polar- though there is also clear evidence for substructure in ization effects in the infrared, A. Mérand for validating the disk around ζ Tau, revealing a one-armed oscillation, our spectral line imaging code, N. Morrison for providing as indicated by different displacement amplitudes in the Ritter spectra, and C. Jones for providing us information blue- and red-shifted line wing. about the photometric variability of β CMi. This work Using our data set on β CMi we constructed a position- was done in part under contract with the California Insti- velocity diagram, which can be interpreted using the tute of Technology (Caltech), funded by NASA through well-established procedures from radio interferometry, the Sagan Fellowship Program (SK is a Sagan fellow). but probes the milli-arcsecond scale position displace- JDM and GHS acknowledge support for this work pro- ments resulting from the rotating disk around this Be vided by the National Science Foundation under grants star. From our kinematical constraints, we derive the AST-0707927 and AST-1009080. The MIRC beam com- dynamical mass of the central star to 3.5 ± 0.2, which biner was developed with funding from the University of 16

Michigan. The CHARA Array is funded by the Geor- tion, by the NASA Exoplanet Science Institute, and the gia State University, by the National Science Foundation David and Lucile Packard Institute. through grant AST-0908253, by the W.M. Keck Founda- Facilities: CHARA, VLTI.

REFERENCES Barnes, T. G., Evans, D. S., & Moffett, T. J. 1978, MNRAS, 183, Monnier, J. D., et al. 2006b, in Presented at the Society of 285 Photo-Optical Instrumentation Engineers (SPIE) Conference, Bjorkman, J. E., & Carciofi, A. C. 2005, in Astronomical Society Vol. 6268, Society of Photo-Optical Instrumentation Engineers of the Pacific Conference Series, Vol. 337, The Nature and (SPIE) Conference Series Evolution of Disks Around Hot Stars, ed. R. Ignace & —. 2007, Science, 317, 342 K. G. Gayley, 75–+ Okazaki, A. T. 1991, PASJ, 43, 75 Bjorkman, J. E., & Cassinelli, J. P. 1993, ApJ, 409, 429 Omidvar, K., & McAllister, A. M. 1995, Phys. Rev. A, 51, 1063 Bonneau, D., et al. 2006, A&A, 456, 789 Pavlovski, K., Harmanec, P., Bozic, H., Koubsky, P., Hadrava, P., Carciofi, A. C. 2010, ArXiv e-prints Kriiz, S., Ruzic, Z., & Stefl, S. 1997, A&AS, 125, 75 Carciofi, A. C., & Bjorkman, J. E. 2006, ApJ, 639, 1081 Pestalozzi, M. R., Elitzur, M., & Conway, J. E. 2009, A&A, 501, —. 2008, ApJ, 684, 1374 999 Carciofi, A. C., Okazaki, A. T., Le Bouquin, J., Stefl,ˇ S., Rivinius, Petrov, R. G., et al. 2007, A&A, 464, 1 T., Baade, D., Bjorkman, J. E., & Hummel, C. A. 2009, A&A, Pollmann, E. 2002, H-alpha Monitoring of Be Stars 504, 915 Porter, J. M., & Rivinius, T. 2003, PASP, 115, 1153 Che, X., et al. 2011, ApJ, 732, 68 Pott, J., et al. 2010, ApJ, 721, 802 Chelli, A., Utrera, O. H., & Duvert, G. 2009, A&A, 502, 705 Quirrenbach, A., Buscher, D. F., Mozurkewich, D., Hummel, Delaa, O., et al. 2011, A&A, 529, A87+ C. A., & Armstrong, J. T. 1994, A&A, 283, L13 Demarque, P., Woo, J.-H., Kim, Y.-C., & Yi, S. K. 2004, ApJS, Quirrenbach, A., et al. 1997, ApJ, 479, 477 155, 667 Rivinius, T., Baade, D., Stefl, S., Stahl, O., Wolf, B., & Kaufer, Draper, Z. H., Wisniewski, J. P., Bjorkman, K. S., Haubois, X., A. 1998, A&A, 333, 125 Carciofi, A. C., Bjorkman, J. E., Meade, M. R., & Okazaki, A. Rivinius, T., Stefl,ˇ S., & Baade, D. 2006, A&A, 459, 137 2011, ApJ, 728, L40+ Rodr´ıguez-Merino, L. H., Chavez, M., Bertone, E., & Buzzoni, A. Duflot, M., Figon, P., & Meyssonnier, N. 1995, A&AS, 114, 269 2005, ApJ, 626, 411 Frémat, Y., Zorec, J., Hubert, A.-M., & Floquet, M. 2005, A&A, Saio, H., et al. 2007, ApJ, 654, 544 440, 305 Schaefer, G. H., et al. 2010, AJ, 140, 1838 Ghosh, K., Iyengar, K. V. K., Ramsey, B. D., & Austin, R. A. Sofue, Y., & Rubin, V. 2001, ARA&A, 39, 137 1999, AJ, 118, 1061 Stee, P. 1996, A&A, 311, 945 Grundstrom, E. D. 2007, PhD thesis, Georgia State University Stee, P., & de Araujo, F. X. 1994, A&A, 292, 221 Hesselbach, E. N. 2009, PhD thesis, The University of Toledo Stefl, S., Rivinius, T., Baade, D., & Carciofi, A. 2011, ArXiv Isella, A., Testi, L., Natta, A., Neri, R., Wilner, D., & Qi, C. e-prints 2007, A&A, 469, 213 Sturmann, J., Ten Brummelaar, T., Sturmann, L., & McAlister, Iwamatsu, H., & Hirata, R. 2008, PASJ, 60, 749 H. A. 2010, in Presented at the Society of Photo-Optical Jones, C. E., Tycner, C., & Smith, A. D. 2011, AJ, 141, 150 Instrumentation Engineers (SPIE) Conference, Vol. 7734, Kato, S. 1983, PASJ, 35, 249 Kervella, P., & Fouqué, P. 2008, A&A, 491, 855 Society of Photo-Optical Instrumentation Engineers (SPIE) Kraus, S., Hofmann, K., Malbet, F., Meilland, A., Natta, A., Conference Series Schertl, D., Stee, P., & Weigelt, G. 2009, A&A, 508, 787 Tatulli, E., et al. 2007, A&A, 464, 29 Le Bouquin, J., Bauvir, B., Haguenauer, P., Schöller, M., ten Brummelaar, T. A., et al. 2005, ApJ, 628, 453 Rantakyrö, F., & Menardi, S. 2008, A&A, 481, 553 Townsend, R. H. D., Owocki, S. P., & Howarth, I. D. 2004, Le Bouquin, J.-B., Absil, O., Benisty, M., Massi, F., Mérand, A., MNRAS, 350, 189 & Stefl, S. 2009, A&A, 498, L41 Tycner, C., et al. 2004, AJ, 127, 1194 Lee, U., Osaki, Y., & Saio, H. 1991, MNRAS, 250, 432 —. 2005, ApJ, 624, 359 McDavid, D. 1999, PASP, 111, 494 Stefl,ˇ S., et al. 2009, A&A, 504, 929 Meilland, A., Stee, P., Chesneau, O., & Jones, C. 2009, A&A, Vakili, F., et al. 1998, A&A, 335, 261 505, 687 Wade, R. A., & Rucinski, S. M. 1985, A&AS, 60, 471 Meilland, A., Stee, P., Zorec, J., & Kanaan, S. 2006, A&A, 455, Waters, L. B. F. M. 1986, A&A, 162, 121 953 Weigelt, G., et al. 2007, A&A, 464, 87 Meilland, A., et al. 2007, A&A, 464, 59 Weintroub, J., Moran, J. M., Wilner, D. J., Young, K., Rao, R., —. 2011, ArXiv e-prints & Shinnaga, H. 2008, ApJ, 677, 1140 Miyoshi, M., Moran, J., Herrnstein, J., Greenhill, L., Nakai, N., Wilson, T. L., Rohlfs, K., & Hüttemeister, S. 2009, Tools of Diamond, P., & Inoue, M. 1995, Nature, 373, 127 Radio Astronomy (Springer-Verlag) Monnier, J. D., et al. 2006a, ApJ, 647, 444 Yi, S. K., Kim, Y.-C., & Demarque, P. 2003, ApJS, 144, 259 Yudin, R. V. 2001, A&A, 368, 912