Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.sao.ru/hq/balega/PUBL_BAL/PUB_2009/2009A&A...497..195K.pdf Äàòà èçìåíåíèÿ: Wed Feb 24 16:01:54 2010 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 12:59:42 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: ï ï ð ð ð ð ð ð ð ï ï ï ï ï ï ï

A&A 497, 195­207 (2009) DOI: 10.1051/0004-6361/200810368
c ESO 2009

Astronomy & Astrophysics

Tracing the young massive high-eccentricity binary system 1 Orionis C through periastron passage
S. Kraus1 ,G.Weigelt1 , Y.Y.Balega2 , J. A. Docobo3 , K.-H. Hofmann1 , T. Preibisch4 , D. Schertl1 , V. S. Tamazian3 , T. Driebe1 , K. Ohnaka1 ,R.Petrov5 ,M.SchÆller6 , and M. Smith7
1 2 3 4 5 6 7

Max-Planck-Institut fÝr Radioastronomie, Auf dem HÝgel 69, 53121 Bonn, Germany e-mail: skraus@mpifr-bonn.mpg.de Special Astrophysical Observatory, Russian Academy of Sciences, Nizhnij Arkhyz, Zelenchuk region, Karachai-Cherkesia 357147, Russia Astronomical Observatory R. M. Aller, University of Santiago de Compostela, Galicia, Spain UniversitÄts-Sternwarte MÝnchen, Scheinerstr. 1, 81679 MÝnchen, Germany Laboratoire Universitaire d'Astrophysique de Nice, UMR 6525 UniversitÈ de Nice/CNRS, Parc Valrose, 06108 Nice Cedex 2, France European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany Centre for Astrophysics & Planetary Science, University of Kent, Canterbury CT2 7NH, UK

Received 11 June 2008 / Accepted 27 January 2009
ABSTRACT

Context. The nearby high-mass star binary system 1 Ori C is the brightest and most massive of the Trapezium OB stars at the core of the Orion Nebula Cluster, and it represents a perfect laboratory to determine the fundamental parameters of young hot stars and to constrain the distance of the Orion Trapezium Cluster. Aims. By tracing the orbital motion of the 1 Ori C components, we aim to refine the dynamical orbit of this important binary system. Methods. Between January 2007 and March 2008, we observed 1 Ori C with VLTI/AMBER near-infrared ( H - and K -band) longbaseline interferometry, as well as with bispectrum speckle interferometry with the ESO 3.6 m and the BTA 6 m telescopes ( B and V -band). Combining AMBER data taken with three different 3-telescope array configurations, we reconstructed the first VLTI/AMBER closure-phase aperture synthesis image, showing the 1 Ori C system with a resolution of 2 mas. To extract the astrometric data from our spectrally dispersed AMBER data, we employed a new algorithm, which fits the wavelength-differential visibility and closure phase modulations along the H - and K -band and is insensitive to calibration errors induced, for instance, by changing atmospheric conditions. Results. Our new astrometric measurements show that the companion has nearly completed one orbital revolution since its discovery in 1997. The derived orbital elements imply a short-period (P 11.3 yr) and high-eccentricity orbit (e 0.6) with periastron passage around 2002.6. The new orbit is consistent with recently published radial velocity measurements, from which we can also derive the first direct constraints on the mass ratio of the binary components. We employ various methods to derive the system mass ( Msystem = 44 ± 7 M ) and the dynamical distance (d = 410 ± 20 pc), which is in remarkably good agreement with recently published trigonometric parallax measurements obtained with radio interferometry.
Key words. stars: formation ­ stars: fundamental parameters ­ stars: individual: 1 Orionis C ­ binaries: close ­ techniques: interferometric ­ stars: imaging

1. Introduction
The Orion Nebula (M42) is one of the closest (d 400-450 pc) and most prominent star-forming regions (see O'Dell 2001 for a review). It contains a massive cluster of very young (1 â 106 yr) stars (cf. Herbig & Terndrup 1986; McCaughrean & Stauffer 1994; Hillenbrand 1997), which is known as the Orion Nebula Cluster (ONC). Due to its relatively close distance and its favorable celestial location (which makes it accessible to observatories on the northern and southern hemispheres), the ONC is probably the best investigated young cluster in the whole sky and has been observed at virtually every wavelength. It is a perfect laboratory for observations of young stellar objects over the
Based on observations made with ESO telescopes at the La Silla Paranal Observatory under the OT and VISA-MPG GTO programme IDs 078.C-0360(A), 080.C-0541(A,B,C,D), 080.D-0225(B), and 080.C-0388(A).

full mass range, from very low-mass brown dwarfs to massive O-type stars. The brightest star in the cluster is the massive O7­O5.5 type1 star 1 Ori C, which is known to be a close visual binary system. After the initial discovery of the companion at a separation of 0. 033 (33 mas, corresponding to about 15 AU) with nearinfrared (NIR) bispectrum speckle interferometry by Weigelt et al. (1999), Schertl et al. (2003) have presented further observations and reported the first detection of orbital motion. Kraus et al. (2007) presented the first speckle observations at visual wavelengths, the first NIR long-baseline interferometric observations of 1 Ori C using the IOTA interferometer, and produced
1 SimÑn-DÌaz et al. (2006) determined the stellar effective temperature of 1 Ori C to be T eff = 39 000 ± 1000 K with log g = 4.1 dex. In recent stellar atmosphere models (Martins et al. 2005), this corresponds to a spectral type of O5.5­O6, while the same parameters correspond to a later spectral type of O7.5­O8 using earlier models (e.g. Vacca et al. 1996).

Article published by EDP Sciences


196

S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

an aperture-synthesis image of the system. They also performed a joint analysis of all existing interferometric measurements that covered a period of more than 9 years and clearly revealed orbital motion. After reaching a maximum value of 42 mas in 1999, the separation of the system steadily decreased to 13 mas in 2005. Detailed modeling of these data yielded a preliminary orbit solution with a high eccentricity (e 0.91) and a period of 10.9 yrs. According to this solution, the periastron passage should have occurred around July 2007 with a closest separation of less than 2 AU. Patience et al. (2008) recently presented additional interferometric observations of 1 Ori C obtained with NPOI at visual wavelengths. Extending the orbital coverage by about 1.2 yrs, they measured a companion position which deviates 4 mas from the position predicted by the orbital solution of Kraus et al. (2007) and concluded that the orbit has a considerably lower eccentricity (e 0.16) and a longer period (26 yr). This shows a clear need for further interferometric observations that will clarify this apparent discrepancy. A reliable and accurate orbit solution will also provide unique constraints to the basic stellar parameters and, most importantly, will yield the masses of the stars and the distance. The orbital motion of the companion might also be responsible for the radial velocity variations reported by Stahl et al. (1996, 2008) and others. Due to the incomplete coverage and the significant scatter in the derived velocities, it is not yet possible to derive the orbital elements of the spectroscopic orbit (Stahl et al. 2008). In spite of its importance and the huge number of studies, the distance to the ONC is, even after decades of investigation, still not well known and an issue of ongoing discussion. The "canonical" value of 450 pc (Herbig & Terndrup 1986; Genzel & Stutzki 1989) that was widely used during the last two decades was recently challenged by some studies finding significantly smaller values. Numerous new distance determinations (e.g. Stassun et al. 2006; Jeffries 2007; Mayne & Naylor 2008), also including the first direct trigonometric parallax measurements on radio sources in the ONC by Hirota et al. (2007, 437 ± 19 pc), Sandstrom et al. (2007, 389+24 pc), and Menten -21 et al. (2007, 414 ± 7 pc), yielded distances mainly in the range between 390 pc and 440 pc. The 13% difference in these distance estimates, which to some extent might also include an intrinsic distance spread of the studied stars, leads to an 30% uncertainty in the derived stellar luminosities and correspondingly affects any age and mass estimates for the stars. This is a serious limitation to the usefulness of the exceptionally well-studied young stellar population; e.g., for the calibration of pre-mainsequence evolutionary models or in the investigation of the spatial relationship between the young stellar groups in the wider region of the Orion association. Therefore, trigonometric parallax measurements of stars in the core of the Trapezium OB star cluster (e.g. on GMR 12 = 1 Ori A2, for which Menten et al. 2007, measured 418 ± 9 pc) are highly desirable. An alternative way to obtain such distance estimates is through astrometric measurements on close binary systems such as 1 Ori C, yielding the dynamical parallax of the system.

2. Observations and data reduction
2.1. Bispectrum speckle interferometry

between 1997 and 2004 at wavelengths between 2.1 m and 545 nm ( = 43 ... 24 mas; Weigelt et al. 1999; Schertl et al. 2003; Kraus et al. 2007) using the Russian Big Telescope Altazimuthal (BTA) 6.0 m telescope, 1 Ori C became unresolvable with 6 m-class telescopes in 2006. In November 2007 ( = 20 mas) and January 2008, the system could again be resolved with the BTA 6 m telescope using a V medium-band filter (c = 550 nm, = 20 nm). For these observations, a 1280 â 1024 pixel CCD with a multi-alkali S25 intensifier photocathode was used. In January 2008, additional speckle measurements using V (c = 550 nm, = 30 nm) and B (c = 440 nm, = 16 nm) medium-band filter were obtained with the ESO 3.6 m telescope. The speckle observations were carried out in the course of ESO open-time programme 080.C-0388(A) using our visitor speckle camera, which employs a Marconi/EEV electronmultiplying CCD. For the speckle observations, we recorded interferograms of 1 Ori C and of nearby unresolved stars in order to compensate for the atmospheric speckle transfer function. The calibrator stars, the number of recorded interferograms and the detector integration times (DITs) are listed in Table 1. The modulus of the Fourier transform of the object (visibility) was obtained with the speckle interferometry method (Labeyrie 1970). For image reconstruction we used the bispectrum speckle interferometry method (Weigelt 1977; Weigelt & Wirnitzer 1983; Lohmann et al. 1983; Hofmann & Weigelt 1986). Besides providing an independent astrometric measurement, our bispectrum speckle interferometric measurement is of special importance as it allows us to solve the 180-ambiguity, which is inherent to long-baseline interferometric investigations which do not include the closure phase (e.g. Patience et al. 2008) or for which the instrumental closure phase sign has not yet been calibrated (as for VLTI/AMBER). Solving this ambiguity is essential for deriving the orbit of the system. Therefore, we paid special attention while deriving the field orientation of our speckle images, using calibration measurements which were taken with the same instrument setup during the same night as the 1 Ori C observations on 1 Ori A, 1 Ori B and the well-studied multi-component object Carinae. For the January 2008 observations, a position angle calibration with an accuracy of 0.3 was done using K -band observations covering the Trapezium stars 1 Ori A,B,and E, followed by the V and B observations of 1 Ori B, which is a binary system with a separation of 1 , where the fainter component is a close binary with a separtion of 0. 15. This allows us to unambiguously determine that in January 2008, the fainter component (C2) was located to the southwest of the primary star (Figs. 1 and A.1). Performing a detailed re-analysis on all Speckle data taken by our group between 1997 and 2005, we found a 180 -calibration problem which affected the position angle of the Speckle measurements at epochs 2003.8, 2003.9254, 2003.928, and 2004.8216, as published in Kraus et al. (2007). Since these Speckle measurements were used for the calibration of the IOTA and NPOI long-baseline interferometric observations by Kraus et al. (2007) and Patience et al. (2008), a revision of the preliminary orbital solutions presented in these papers is required (as presented in Sect. 5.1). In Table 3, we list all available astrometric data, taking the quadrant correction into account.
2.2. VLTI/AMBER spectro-interferometry

Bispectrum speckle interferometry is a powerful technique to overcome atmospheric perturbations and to reach the diffractionlimited resolution of ground-based telescopes at visual wavelengths. After monitoring the orbital evolution of the system

AMBER is the NIR interferometric beam-combiner instrument (Petrov et al. 2007) of the Very Large Telescope Interferometer (VLTI), which is located on Cerro Paranal/Chile and operated


S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage Table 1. Observation log for our new bispectrum speckle and long-baseline interferometric observations. Instrument BTA 6 m/Speckle ESO 3.6 m/Speckle ESO 3.6 m/Speckle BTA 6 m/Speckle VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER VLTI/AMBER Date [UT] 2007 2008 2008 2008 2007 2007 2007 2007 2007 2007 2007 2008 2008 2008 2008 2008 2008 UT Nov. 25 Jan. 10 Jan. 10 Jan. 26 Jan. 08 Dec. 03 Dec. 03 Dec. 03 Dec. 05 Dec. 05 Dec. 05 Feb. 22 Feb. 22 Feb. 24 Feb. 24 Mar. 03 Mar. 03 Spectral Mode 550 nm/20 nm 440 nm/16 nm 550 nm/30 nm 550 nm/20 nm UT1-UT3-UT4 L R- K A0-D0-H0 LR-HK A0-D0-H0 LR-HK A0-D0-H0 LR-HK A0-K0-G1 LR-HK A0-K0-G1 LR-HK A0-K0-G1 LR-HK A0-D0-H0 LR-HK A0-D0-H0 LR-HK A0-K0-G1 LR-HK A0-K0-G1 LR-HK D0-H0-G1 LR-HK D0-H0-G1 LR-HK Telescope Triplet DIT 20 ms 10 ms 10 ms 20 ms 26 ms 26 ms 50 ms 100 ms 26 ms 50 ms 100 ms 50 ms 100 ms 50 ms 100 ms 50 ms 100 ms No. Interferograms Target/Calibrator 1940/1940 18 000/10 000 10 000/6000 1500/2000 12 000/10 000 5000/10 000 7500/7500 2000/750 5000/5000 10 000/15 000 3500/2500 5000/10 000 2500/3500 10 000/5000 6000/4000 5000/5000 1500/1500 Calibrator(s)

197

04:36, 07:19 07:34, 07:48, 06:05 05:49, 06:17, 03:44 04:21 03:11 02:37 02:31 02:52

06:32 08:44 08:54 07:28 07:46

1 Ori D 36 Ori 36 Ori 1 Ori D HD 41547 HD 33833 HD 33833 HD 33833 HD 33833 HD 33833 HD 33833 HD 37128, HD 50281 HD 37128, HD 50281 HD 37128 HD 37128 HD 43023 HD 43023

by ESO. For the beam combination, either three 8.2 m unit telescopes (UTs) or three of the allocatable 1.8 m auxiliary telescopes (ATs) can be used. One outstanding feature of AMBER is its spectral capability, allowing one to observe several spectral bands with low spectral resolution simultaneously (LR- JHK mode, R = / = 35, covering the J -band around 1.2 m, the H -band around 1.6 m, and the K -band around 2.2 m). The wide wavelength coverage offered by this spectral mode allows us to fit the binary separation vector with a high accuracy and observing efficiency (as discussed in Sect. 3). The VLTI/AMBER data was recorded in the course of three ESO programmes, yielding the uv-plane coverage shown in Fig. 2. For each science observation, we also recorded interferograms for at least one interferometric calibrator star (see Table 2 for the characteristics of our calibrator stars), which allows us to calibrate the visibilities and the closure phases for instrumental as well as atmospheric effects. For AMBER data reduction, we employed the amdlib2-data reduction software2 (release 2.1), which is based on the P2VM algorithm (Tatulli et al. 2007). For the LR-HK data, we applied the data selection criteria outlined in Kraus et al. (2008)in order to remove frames which were either degraded by atmospheric effects or were recorded significantly offset from zero optical path delay. For the accuracy achievable with our binary model fits, the absolute calibration of the wavelength scale is of special importance. Therefore, we carefully refined the wavelength calibration using the telluric gaps between the observed spectral bands. Using this procedure (which is described in Appendix A), we reach a calibration accuracy of about 0.03 m, which is still the dominating limiting factor on the total achievable astrometric accuracy (2%). While the target/calibrator observations from January 2007 were taken under good and stable atmospheric conditions, some of the observations in December 2007 and February 2008 suffer from strongly variable seeing conditions and short atmospheric coherence times. It is known that short coherence times can result in a decrease of the measured fringe contrast, which might not be completely calibrated out using calibrator measurements. This effect can impose errors on the absolute calibration, which
2 The amdlib2 software package is available from the website http://www.jmmc.fr/data_processing_amber.htm

are expected to increase with longer DIT and towards shorter wavelengths. To illustrate this effect, in Fig. 4 we plot the calibrated visibilities measured during our December 2007 observation campaign and compare them to the cosine visibility modulation of a binary source. In particular, in the H -band the resulting calibration errors can be on the order of 20%. However, since all spectral channels of an AMBER interferogram are recorded at the same time, the wavelength-differential observables (in particular, the differential visibility V ()) are practically insensitive to this degradation. Therefore, in Sect. 3 we employ a fitting algorithm in which we fit only differential visibilities and closure phases in order to determine the 1 Ori C binary separation vector. A fundamental problem one encounters when interpreting VLTI/AMBER 3-telescope data concerns the lack of a calibration measurement for the closure phase sign, resulting in a 180 uncertainty on the position angle of binary star observations. Since our 1 Ori C VLTI/AMBER observations from December 2007 and February 2008 bracket the bispectrum speckle measurement from January 2008, our data set allows us to unambigously define the closure phase sign for these observations, providing a direct calibration of the 180 uncertainty for VLTI/AMBER for the first time.

3. Modeling In order to derive the binary separation and position angle3 from the speckle data, we used the same algorithm as in our earlier studies on this system (Weigelt et al. 1999; Schertl et al. 2003; Kraus et al. 2007), fitting a cosine function directly to the 2-D speckle visibilities. For illustration, in the Appendix (Fig. A.1) we show power spectra and Fourier spectra determined from three independent data sets obtained during our observing run with the ESO 3.6 m telescope and compare them with the model power and Fourier spectra corresponding to a binary star. For the modeling of the AMBER data, we employed an optimized algorithm which fits closure phases (CP) and differential visibilities V . Each AMBER LR-HK measurement records 16 K -band plus 11 H -band spectral channels in the wavelength range from 1.51 to 2.55 m, covering about 40% of the object Fourier spectrum in radial direction (see uv-plane tracks in
3 Following convention, we measure the position angle (PA) East of North.


198

S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage
1997.8, BTA 6 m, H -band 1999.7, BTA 6 m, J -band 2000.8, BTA 6 m, V -band

2003.9, BTA 6 m, V -band

2004.8, BTA 6 m, V -band

2005.9, IOTA/IONIC3, H -band

2007.9, VLTI/AMBER, K -band

2007.9, ESO 3.6 m, B -band

2008.1, BTA 6 m, V -band

Fig. 1. Selection of interferometric images of the 1 Ori C system obtained by our group between 1997 and 2008, revealing the orbital motion of the companion. The images were reconstructed either from H -, J -, V -, or B -band speckle interferograms recorded at the BTA 6 m or the ESO 3.6 m telescopes, or from H -band Michelson interferograms recorded with the IOTA 3-telescope interferometer (see Kraus et al. 2007, for details). For epoch 2007.9, we show the aperture synthesis image which we reconstructed from VLTI/AMBER data using the procedure described in Sect. 4. For each image, 10% intensity level contours are shown. In the speckle images, the fitted component positions are marked with an asterisk.

Fig. 2). Strongly resolved objects (such as binaries with a separation of / B, where is the recorded spectral window and B is the projected baseline length) can already show significant visibility modulation over this range of spatial frequencies (see illustration in Fig. 3). For instance, our 1 Ori C AMBER measurements from December 2007 probe spatial frequencies out to the eighth lobe of the cosine binary visibility modulation (Fig. 4) and the wavelength-differential visibilities recorded in a single AMBER H -/K -band measurement sample up to 3 visibility cycles. This wavelength-differential visibility modulation already provides all information required to determine the binary separation and orientation. In this way, the absolute calibration

of the visibility, which is subject to many adverse atmospheric effects, becomes dispensable. To determine the differential visibilities V as used in our fitting algorithm, we first follow the standard data reduction and calibration procedure in order to correct for wavelengthdependent instrumental effects (using a calibrator measurement taken during the same night). Then, we remove the absolute calibration by subtracting the average visibility of the considered spectral window V () = V () - V ()
band

,

(1)


S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

199

Table 2. AMBER calibrator stars and their characteristics, including uniform disk (UD) diameters. Star H H H H H D D D D D 33833 37128 41547 43023 50281 V 5. 1. 5. 5. 6. 9 7 9 8 6 H 3. 2. 5. 3. 4. 9 4 1 7 3 K 3. 2. 5. 3. 4. 8 3 0 5 1 Spectral Type G7III B0I F4V G8III K3V Adopted UD diameter [mas] 0.83 ± 0.06a 0.86 ± 0.16b 0.41 ± 0.03a 0.98 ± 0.07a 0.77 ± 0.10c

Notes: a UD diameter computed with ASPRO (http://www.jmmc. fr/aspro_page.htm). b UD diameter taken from Mozurkewich et al. (1991). c UD diameter taken from getCal tool (http://mscweb. ipac.caltech.edu/gcWeb/).

Fig. 2. uv-plane coverage obtained at four epochs on VLTI/AMBER. The radial extension of the uv-tracks spectral coverage of our AMBER interferometric data, ther the K -band (LR- K mode, 2007 Jan. 08) or the H (LR-HK mode, 2007 Dec. 03­05, 2008 Feb. 22­24, 2008

1 Ori C with reflects the covering eiand K -band Mar. 03).

yielding the differential visibility V , which we compute separately for each spectral band ( H and K -band). A similar approach is applied to the model visibilities Vmodel () (which we compute using Eqs. (7), (11), and (12) from Kraus et al. 2005): V
model

source structure, it might, in some cases, not be applicable or might lead to biased results. Therefore, it is highly desirable that optical interferometers such as VLTI/AMBER have the capability of recovering the source brightness distribution free of any assumptions. The aim of this section is to apply state-ofthe-art aperture synthesis imaging techniques for the first time to real VLTI/AMBER data in order to independently confirm the scientific results obtained for 1 Ori C in the last section, and, simultaneously, to demonstrate the imaging capabilities of VLTI/AMBER on a relatively well-studied astrophysical target with limited intrinsic complexity. In order to obtain the uv-coverage required for aperture synthesis imaging, we combined the 1 Ori C AMBER data sets taken between December 2007 and March 2008 on three different 3-telescope array configurations and at several hour angles. Assuming that the source morphology does not change significantly over the K -band (which seems well justified based on our earlier measurements of the wavelength-dependent binary flux ratio; Kraus et al. 2007), we make use of AMBER's spectral coverage, yielding radial tracks in the uv-plane (Fig. 7, left). Since the K -band visibilities provide a more reliable absolute calibration (see Sect. 2.2), we did not use the H -band data for image reconstruction and rejected also one measurement taken under particularly poor and variable conditions (2007 Dec. 05, UT 07:46). In order to correct for the binary orbital motion over the 3-month period ( = 2.2mas, = 6.6 ), we apply a rotation-compensating coordinate transformation of the uv-plane (see Kraus et al. 2005, for a description of this procedure) using the astrometric data given in Table 3. Then, we employed our Building Block Mapping software, which is based on the algorithm described by Hofmann & Weigelt (1993). Starting from an initial single -function, this algorithm adds components to a model image in order to minimize the deviations between the measured bispectrum and the bispectrum of the model image. Finally, the image is convolved with a clean beam of 1.5 â 3 mas, reflecting the elongation of the sampled uv-plane. The resulting K -band aperture synthesis image (Figs. 1 and 7, right) yields a direct, model-independent representation of our VLTI/AMBER interferometric data. Measuring the companion position in the reconstructed image ( = 19.3mas, = 241, epoch 2007.9), we find good agreement with the results obtained with our model fitting approach using wavelength-differential observables (Sect. 3). The noise features within the image show an intensity amplitude below 2% of the peak brightness in the image.

() = a (V

model

() - b),

(2)

where a and b are adjusted to minimize the residuals between the model visibilities and the measured wavelength-differential visibilities V () before a Levenberg-Marquardt least-square fitting algorithm is used to determine the best-fit model parameters. The measured closure phases are fit simultaneously with the differential visibilities, weighting each data point according to the determined error bars (see Eqs. (8)­(10) in Kraus et al. 2005). Parameters in our binary star model are the binary separation , the position angle , the intensity ratio FC2 /FC1 , and the uniform disk diameter of the components. Since our earlier speckle measurements indicated that FC2 /FC1 is practically constant in the NIR wavelength range (see Kraus et al. 2007), we assume that the intensity ratio does not vary over the H - and K -bands. For the model fitting, we fix the apparent diameters of the two stars to 0.22 mas for the primary star (corresponding to 10.6 R at a distance of 450 pc) and 0.15 mas for the companion (7.2 R , using the spectral type determined by Kraus et al. 2007;and the stellar evolutionary models of Martins et al. 2005). Our AMBER LR measurements and best-fit models are shown in Figs. 5 and 6. The derived astrometric data are given in Table 3.

4. VLTI/AMBER aperture synthesis imaging
Currently, the most commonly applied procedure to extract scientific information from optical interferometric data is model fitting, which requires the selection of an appropriate (geometrical or astrophysically motivated) model, whose parameters are then adjusted to fit the interferometric observables. Since this model fitting approach requires a priori knowledge about the


200

S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

Table 3. Astrometric measurements for the 1 Ori C binary system. Telescope BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle BTA 6 m/Speckle IOTA NP OI VLTI/AMBER NP OI NP OI NP OI NP OI NP OI BTA 6 m/Speckle VLTI/AMBER ESO 3.6 m/Speckle ESO 3.6 m/Speckle BTA 6 m/Speckle VLTI/AMBER VLTI/AMBER Date 1997. 1998. 1999. 1999. 2000. 2001. 2003. 2003. 2003. 2004. 2005. 2006. 2007. 2007. 2007. 2007. 2007. 2007. 2007. 2007. 2008. 2008. 2008. 2008. 2008. 784 838 737 8189 8734 184 8 9254 928 8216 92055 1486 0192 1425 1507 1753 2055 2137 9014 9233 0274 0274 0712 1479 1726 Filter H K J G V J J V V V H V K V V V V V V H+K B V V H+K H+K Flux ratio FC2 /FC1 0.26 ± 0.02 0.32 ± 0.03 0.31 ± 0.02 0.35 ± 0.04 0.35 ± 0.03 0.29 ± 0.02 0.30 ± 0.02 ­ ­ 0.34 ± 0.04 0.28 ± 0.03 ­ 0.31 ± 0.03 ­ ­ ­ ­ ­ 0.29 ± 0.02 0.24 ± 0.07 0.29 ± 0.02 0.31 ± 0.02 0.31 ± 0.02 0.23 ± 0.09 0.26 ± 0.06 [ ] 226.0 ± 3 222.0 ± 5 214.0 ± 2 213.5 ± 2 210.0 ± 2 208.0 ± 2 19.3 ± 2 19.0 ± 2 19.1 ± 2 10.5 ± 4 342.74 ± 2 332.3 ± 3.5 274.9 ± 1 268.1 ± 5.2 272.9 ± 8.8 266.6 ± 2.1 265.6 ± 1.9 263.0 ± 2.3 238.0 ± 2 241.2 ± 1 237.0 ± 3 236.5 ± 3 236.2 ± 2 234.6 ± 1 236.4 ± 1 [mas] 33 ± 2 37 ± 4 43 ± 1 42 ± 1 40 ± 1 38 ± 1 29 ± 2 29 ± 2 29 ± 2 24 ± 4 13.55 ± 0. 11.80 ± 1. 11.04 ± 0. 11.94 ± 0. 12.13 ± 1. 12.17 ± 0. 12.28 ± 0. 12.14 ± 0. 19.8 ± 2 19.07 ± 0. 19.7 ± 3 19.6 ± 3 20.1 ± 2 21.17 ± 0. 21.27 ± 0. Ref. a a b c c b c c c c c d ­ d d d d d ­ ­ ­ ­ ­ ­ ­

5 11 5 31 58 37 41 43 5

5 5

Notes ­ The position angles given in Col. 5 are measured East of North and were corrected for the 180 -calibration problem described in Sect. 2.1. References ­ (a) Weigelt et al. (1999), (b) Schertl et al. (2003), (c) Kraus et al. (2007), (d) Patience et al. (2008).

5. Results
5.1. Dynamical orbit of the 1 Ori C binary system

Since our new astrometric data extend the orbital coverage for the 1 Ori C system by about 12 months and solves the calibration problem described in Sect. 2.1, we can significantly improve the orbital solution of the system. To derive a refined orbital solution, we applied two independent orbit fitting approaches. First, we used the method presented by Docobo (1985), which generates a family of Keplerian orbits, whose apparent orbits pass through three base points. These base points might be selected from the actual astrometric measurements or represent observationally favored areas in the (, , t) parameter space. From the generated class of possible solutions, the orbit which best agrees with the measured separations and PAs is selected. We use the error bars of the individual measurements as weight. The orbital elements for the determined best-fit orbit solution are given in Table 4. In addition, we implemented a grid search algorithm which scans the parameter space of the dynamical elements eccentricity e, period P, and time of periastron passage T 0 . As described by Hilditch (2001), at each grid point, the geometrical orbital elements a, i, , and can be determined by fitting the Thiele-Innes constants to the observational data. We scanned the parameter space between P = 9 ... 16 yrs (in increments of 0.001 yrs), e = 0.0 ... 1.0 (in increments of 0.005), and T 0 = 2000.0 ... 2010.0 (in increments of 0.001 yrs) and determined the least-square distance between the N measured positions (i , i ) and the corresponding orbit positions (i , i ) for each orbit solution: i - i 2 i - i 2 1 · 2 r = + (3) N i=1...N i i

Using the 2 value determined at each grid point, we built a r 2 (P, e, T 0) data cube, which we projected to determine the lor cal minimum and the associated uncertainty for each parameter. Figure 9 shows the determined 2 minima curves, from which r we determine the best-fit orbital elements given in Table 4.Given that the orbital coverage has been substantially improved since our 2007 study, we do not have to impose a priori constraints on the distance/system mass in order to eliminate unphysical orbit solutions. In Fig. 8c we compare the orbital solutions determined with the Docobo (1985) algorithm and with the grid search algorithm to the available astrometric data.
5.2. Constraining the binary mass ratio

Besides the 15.424-day period which is associated with the wind from the primary, long-term radial velocity variations were also found (Vitrichenko 2002; Stahl et al. 2008). Using a large data set covering more than 15 yrs of spectroscopic observations (plus three archival measurements, which extend the coverage to more than 64 yrs), Stahletal. (2008) showed that these variations are consistent with the orbital motion of a high-eccentricity binary system. Although the strong scatter within the radial velocity measurements prevents us from solving for the precise spectroscopic orbit, the combination of these data with our new orbital solution can be used to provide a first direct constraint on the mass ratio of the components in the 1 Ori C system. In Fig. 10 we plot the available radial velocity data as a function of orbital phase using the orbital period P and time for periastron passage T 0 determined independently from our interferometric measurements (Sect. 5.1). Using the method from Pourbaix (1998), we compute the radial velocity variations corresponding to our full set of orbital elements (Table 4) and perform a least-square fit between the measured (vi ) and the


S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

201

Fig. 3. Illustration of the basic principle of fitting binary parameters using wavelength-differential interferometric observables. Left: the VLTI/AMBER 3-telescope interferometer measures the interferometric observables in various spectral channels from 1.4 to 2.4 m (as indicated by the color of the dots) and towards three different position angles, probing different regions in the two-dimensional Fourier spectrum of the source brightness distribution. The figure shows the visibility spectrum for a binary with = 19.07 mas and = 241.2 (as inferred for 1 Ori C on 2007 Dec. 03) and the uv-sampling obtained with one of our AMBER observations on the same date. Middle and right: the visibilities and phases show a wavelength-differential modulation which is independent of the absolute calibration (see Sect. 3) and which can be fitted to analytical models. Table 4. Orbital elements as determined with the algorithm from Docobo (1985) and our grid search algorithm (described in Sect. 5.1). Docobo algorithm 11.05 ± 0.03 2002.87 ± 0.40 0.534 ± 0.050 40.00 ± 3.00 100.7 ± 1.0 25.3 ± 1.5 290.9 ± 2.5 1.84 524 ± 130 0.21 ± 0.05 49 ± 4 456 ± 13 Grid search algorithm 11.26 ± 0.5 2002.57 ± 0.5 0.592 ± 0.07 43.61 ± 3 99.0 ± 2.6 26.5 ± 1.7 285.8 ± 8.5 0.56 645 ± 200 0.23 ± 0.05 47 ± 4 416 ± 12

Fig. 4. AMBER visibilities measured on 2007 Dec. 03 and 2007 Dec. 05 plotted versus projected distance in the uv-plane (where the projection was perpendicular to the fitted binary separation vector = 241.2 ,i.e. x = u cos(+ 90 ) - v sin(+ 90 )). The solid red line shows the theoretical cosine visibility profile for a binary star with separation 19.07 mas and intensity ratio 0.30. As indicated by the strong visibility offsets, which particularly occur at high spatial frequencies, the absolute calibration is sometimes rather poor, reflecting the changing atmospheric conditions during these nights. It can also be seen that this calibration bias is particularly important for long DITs (50 ms, 100 ms), while it is nearly negligible for short DITs (26 ms, grey & black points). As expected, the spectral dependence of the visibility is not affected by these calibration uncertainties.

Parameter P T0 e a i 2 r a3 /P2 MC1 /MC2 MC1 + MC2 ddyn

[yrs] [mas] [ ] [ ] [ ] [mas3 /yrs2 ] [M ] [pc]

Notes. Besides the orbital elements, we give the mass ratio (Sect. 5.2), dynamical distance, and system mass (Sect. 5.3), derived from both set of orbit elements. The dynamical distance and system mass was determined using the method from Baize & Romani (1946, method c in Sect. 5.3) and three different MLRs. When assuming another distance d , the dynamical system mass MC1 + MC2 must be scaled by a factor (d /ddyn )3 . The mass ratio MC1 /MC2 was also computed for the distance ddyn , but can be converted to any other distance using Eq. (5).

predicted (vi ) radial velocities in order to determine the massratio between the binary components. As least-square measure for N measurements, we use 2 = r 1 N vi - vi - V vi
0 2

,

(4)

i=1...N

where V0 is the velocity of the center of mass of the system, which is given by V0 := vi - vi . Since the radial velocities were extracted from various spectral lines (C iv,He ii,and O iii) and it is known that these lines can show, with respect to each other, systematic velocity offsets on the order of 2­3 km s-1 (Stahl et al. 2008), we used 3 km s-1 as minimum velocity error vi in order to avoid overweighting individual measurements.


202

S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

Fig. 5. Differential visibilities and closure phases measured with VLTI/AMBER in January 2007 (left), February 2008 (middle), and March 2008 (right) on 1 Ori C. The solid lines show the best-fit model corresponding to the binary parameters given in Table 3.

Fig. 6. Similar to Fig. 5, showing the AMBER data and best-fit model for our observations from December 2007.


S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

203

Fig. 7. Combining AMBER data obtained on three telescope configurations (left: uv-coverage), we reconstructed an aperture synthesis image of the 1 Ori C system with an effective resolution of 2mas (right). For a detailed description, we refer to Sect. 4.

40

New orbit solution (Docobo Algorithm) New orbit solution (Grid Search Algorithm) Data 2003.9 (Sp6)

1 Ori C system with respect to the parental cloud is smaller than previously assumed (O'Dell 2001; Stahl et al. 2008).
5.3. Dynamical masses and parallaxes

30 2004.8 (Sp6) 20

10

2005.9 (IOTA) 2006.1 (NPOI)

dDEC [mas]

0



2007.0 (VLTI) 2007.1 (NPOI) 2007.2 (NPOI) 2007.9 (Sp6) 2007.9 (VLTI) 2008.0 (Sp3.6) 2008.1 (Sp6) 2008.2 (VLTI)

-10

-20 1997.8 (Sp6) 1998.8 (Sp6) -30 (Sp6) 2001.2 (Sp6) 2000.8 -40 10 0 -10 dRA [mas] -20 -30 -40 1999.8 (Sp6)

Fig. 8. Comparison of our new orbit solutions with the available astrometric data. Each position measurement is connected to the orbit prediction with an O­C vector (arrows). Furthermore, we show the lines of nodes in grey. The plots are centered on the primary star. North is up and east is to the left.

Kepler's third law ( Msystem · 3 = a3 /P2 ) relates a and P with the product of the system mass Msystem := MC1 + MC2 and the cube of the parallax . Therefore, our astrometric measurement of a3 /P2 directly constrains certain areas in the (, Msystem ) parameter space, as shown by the shaded area in Fig. 13. Since the elements of the spectroscopic orbit are currently only weakly constrained, it is not possible to directly separate the system mass and the dynamical parallax. However, several independent methods exist to disentangle these important parameters using additional information, such as the flux ratio of the components, their mass ratio, and/or their stellar parameters. The stellar parameters of 1 Ori C1 have already been studied extensively with spectroscopy, placing the effective temperature between 37 000 and 40 000 K (e.g. Rubin et al. 1991; Baldwin et al. 1991; SimÑn-DÌaz et al. 2006; Polyakov & Tsivilev 2007). In particular, SimÑn-DÌaz et al. (2006) included non-LTE and line-blanketing effects (T eff ,C1 = 39 000 ± 1000 K) for their detailed spectroscopic modeling, which makes their results compatible with the O-star stellar parameter calibration by Martins et al. (2005). Since each of the above-mentioned stellar or observational parameters is associated with certain assumptions and uncertainties, it seems advisable to take a number of alternative approaches for deriving the underlying physical parameters in order to yield some insight into the associated uncertainties. Therefore, we follow three alternative approaches: a) based on the effective temperature determination of the primary and the binary flux ratio, one can give a reasonable mass range for each component and, thus, the system mass (green-shaded area in Fig. 13); b) one can estimate the mass of the primary from the stellar temperature and then derive the companion mass using the mass ratio constraints obtained from the radial velocities (Sect. 5.2), yielding the orange curve in Fig. 13; c) Baize & Romani (1946) presented a method which solves for the system mass and the dynamical parallax using a mass-luminosity relation (MLR), the bolometric corrections of the components and their extinction-corrected apparent magnitudes. To evaluate the influence of the MLRs on the

By varying the mass-ratio between the components, we find q(d ) := MC1 = 5.27 MC2
+1.17 -0.75

d -1 414 pc

-1

,

(5)

i.e. q(450 pc) = 0.21 ± 0.04 or q(414 pc) = 0.23 ± 0.05, which is slightly lower than the value we derived earlier by modeling the wavelength-dependent binary flux ratio of the 1 Ori C system (q = 0.45 ± 0.15, Kraus et al. 2007). For the radial velocity of the center of mass, we determine 23.6 km s-1 , which is in good agreement with the heliocentric velocity of the Orion Molecular Cloud (28 km s-1 , O'Dell 2001). This might indicate that the relative motion of the


204

S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

a) Orbital period P

b) Eccentricity e

c) Periastron passage T

0

Fig. 9. Minimum 2 curves as function of the dynamical orbital parameters P, e, T 0 . These curves were used to determine the best-fit orbit solution r and to evaluate the uncertainties on the individual parameters (see Table 4).

Fig. 10. The radial velocities measured towards 1 Ori C plotted as function of orbital phase (assuming the orbital elements determined with our grid search algorithm; see Table 4). Besides the extensive data set by Stahl et al. (2008), we included the radial velocities measured by Struve & Titus (1944), Conti (1972), and Morrell & Levato (1991). The red line shows the radial velocities corresponding to our orbital solution with a system velocity V0 = 23.6 km s-1 (blue dashed line) and the mass ratio q(414 pc) = 0.23 as determined from our least-square fit (Sect. 5.2).

Fig. 11. 2 curve for our fit of the binary mass ratio to the available r radial velocity data (assuming a distance of 414 pc). We assume the orbital elements (Table 4) derived from our interferometric measurements.

result, we used three different MLRs; namely, from Baize & Romani (1946), Heintz (1978), and Demircan & Kahraman (1991) (corresponding to the three blue data points in Fig. 13). For methods a) and b), we employ the stellar calibration by Martins et al. (2005), while for method c), three earlier calibrations are used. In method a) and b), the stellar temperature of the primary component is fixed to the value by SimÑn-DÌaz et al. (2006, T eff ,C1 = 39 000 ± 1000 K), while in c) we scan a slightly wider range of temperature values (T eff ,C1 = 37 000 ... 40 000 K, T eff ,C2 = 30 000 ... 33 000 K, using the bolometric correction by Bessell et al. 1998). Method a) requires the V -band flux ratio (FC2 /FC1 = 0.31 ± 0.02) and c) the extinction-corrected magnitude of the total system (V = 5.12 ± 0.1, AV = 1.74 ± 0.1, Hillenbrand 1997). Evidently, each method is associated with considerable uncertainties, which makes it very desirable for future observations to improve not only the astrometric orbit, but to derive the accurate spectroscopic orbit of the system as well. Nevertheless, within their large uncertainties, the methods employed cover a common area in parameter space, corresponding to a system mass of Msystem = 44 ± 7 M and a dynamical distance of

Fig. 12. Minimum 2 curve as a function of a3 /P2 = ( Msystem ) · 3 , r through which the dynamical orbital parameters can be related to the mass sum and the dynamical parallax. For illustration, we give the dynamical distance in the upper axis assuming a total system mass of 45 M .

d = 410 ± 20 pc, as marked with the red ellipse in Fig. 13. Methods a) and b) yield systematically lower system masses than method c), reflecting the recent correction in the mass calibration scale in stellar evolutionary models (Martins et al. 2005). Assuming a distance of 414 pc, as determined by Menten et al. (2007), would yield a system mass of 46 M .


S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

205

Fig. 13. The derived a3 /P2 value puts some direct constraints on the system mass Msystem := MC1 + MC2 and the dynamical distance. As described in Sect. 5.3, we employ three alternative methods to obtain constraints, which help to disentangle the system mass and the distance (method a: dashed green-shaded area, method b: orange curve, method c: blue data points). The black-shaded area marks the a3 /P2 constraints determined with our orbit grid search algorithm. With the red ellipse, we visually mark the area of best agreement between all constraints.

5.4. Possible implications on the dynamical history of the 1 Ori C system

From our long-baseline interferometric data, we have reconstructed the first model-independent VLTI/AMBER aperture synthesis image, depicting the 1 Ori C system at a resolution of 2 mas and demonstrating the imaging capabilities of this unique facility. In order to extract accurate astrometric information for all epochs, we have followed a new modeling approach which is based on wavelength-differential observables and which demonstrates the benefits of spectro-interferometry in terms of observing efficiency and robustness to poor observing conditions. Furthermore, our ESO 3.6 m and the BTA 6 m speckle observations allow us to solve 180-ambiguity and to calibrate the closure phase sign of our VLTI/AMBER observations, providing a potential reference for other AMBER observations using closure phase information4. Our new 1 Ori C astrometric data shows that since its discovery in 1997, the 1 Ori C companion has nearly completed one orbital revolution and that the system has a high eccentricity (e 0.6). Solving for the orbital elements, we determine a period of 11.3 yrs, a semi-major axis of 44 mas, and a periastron passage around 2002.6. According to our orbital solutions, the physical separation between the components decreases to 7 AU (angular separation 2.8 mas) during periastron passage, which might be too large to result in detectable signatures of wind-wind interaction between the stellar winds. Using additional information about the stellar parameters and various stellar models, we estimate the total system mass to be 44 ± 7 M and the dynamical distance to be d = 410 ± 20 pc. In the coming years, the uncertainties on these parameters could be considerably reduced with new spectroscopic and astrometric observations of this important binary system.
Acknowledgements. We would like to thank the referee, M. McCaughrean, for helpful comments which improved this paper. Furthermore, we acknowledge helpful discussions with E. Vitrichenko, F. Millour and members of the AMBER consortium.

With an eccentricity of 0.6, the orbit of 1 Ori C is located on the upper end of the eccentricity distribution of low- as well as high-mass binary stars (Mathieu 1994; Mason et al. 1998), perhaps providing important information about the dynamical history of the system. Tan (2004, 2008) proposed that the BecklinNeugebauer (BN) object, which is located 45 northwest of the Trapezium stars, might be a runaway B star ejected from the 1 Ori C multiple system approximately 4000 yrs ago. This scenario is based on proper motion measurements, which show that BN and 1 Ori C recoil roughly in opposite directions. Threebody interaction is a crucial part of this interpretation, and the high eccentricity of the 1 Ori C orbit which we derive for this system might be a direct consequence of this dynamical interaction event. However, another study (RodrÌguez et al. 2005) also aimed to identify the multiple system from which BN was ejected and identified Source I as the likely progenitor system. Later, GÑmez et al. (2005, 2008) added further evidence to this interpretation by identifying Source n as a potential third member of the decayed system. Therefore, it is still unclear whether this scenario can explain the measured properties of the 1 Ori C orbit. As an alternative explanation for the high eccentricity of the 1 Ori C system, Zinnecker & Yorke (2007) pointed out that such systems are predicted by star formation scenarios which include sub-Keplerian rotating disks or filament fragmentation (Kratter & Matzner 2006).

Appendix A: Spectral calibration of AMBER-LR data
Employing wavelength-differential visibilities and phases for model fitting requires a precise knowledge of the central wavelength of the recorded spectral channels. Therefore, we performed a re-calibration of the wavelength axes of our data sets using atmospheric absorption features and found deviations of up to 0.1 m compared to the standard wavelength table applied to the data on Paranal. In the recorded J -/H -/K -band object raw spectra, the telluric features are not pronounced enough to be used for the calibration, which is mainly due to absorption by internal optical components. Therefore, besides the raw spectrum of the object exposures, we also extracted the spectrum from the sky and the P2VM calibration exposures. The P2VM calibration files (P2VM = Pixel-to-Visibility-Matrix; see Tatulli et al. 2007) are recorded at the beginning of each observation block using a lamp located in AMBER's Calibration and Alignment Unit and provide a measure of the instrumental transmission. By subtracting the sky spectrum from the object spectrum and then dividing by the P2VM spectrum, we yield a spectrum which is corrected for most instrumental effects and shows the telluric absorption features much more clearly (see Fig. A.2, red curve). These corrected spectra were compared with standard atmospheric
4 Since our data set might be useful for VLTI/AMBER users as reference data for the calibration of the closure phase sign, we provide our data on the following website: http://www.mpifr.de/staff/skraus/files/amber.htm

6. Conclusions
We have presented new bispectrum speckle (V -/ B -band) and VLTI/AMBER (H -/K -band) interferometric observations of the Orion Trapezium star 1 Ori C covering several epochs over a time period of about 14 months.


206

S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage

Fig. A.1. Left: three two-dimensional average power spectra (top row) and Fourier phase spectra (middle row) derived from B -band speckle data recorded on 2008 Jan. 10 with the ESO 3.6 m telescope. Right: comparing this data to model power spectra (top row) and Fourier spectra (middle row) for a binary with separation = 19.7 mas and PA = 237 shows that the fainter component is located southwest of the primary star (the bispectrum speckle image reconstructed from the model visibilities and phases are shown to the bottom right). Together with the calibration data sets mentioned in Sect. 2.1, this data set unambigously defines the current orientation of the 1 Ori C binary system and shows that the system has nearly completed one orbital revolution since its discovery in 1997.

pronounced gaps between the J -/H -and H -/K -band and find best agreement assuming a linear dispersion law (i) 0.0328 · i, where i is the number of the spectral channel on the detector. The remaining wavelength calibration uncertainty of about 1 spectral channel (corresponding to about 0.03 m) is the dominant limiting factor on the achievable astrometric accuracy (2%). Therefore, for future observations, it seems highly desirable to implement an AMBER on-site spectral calibration device (e.g. NIR lasers with well-known frequencies), fascilitating an absolute spectral calibration in the standard instrument calibration procedure, pushing AMBER to its full astrometric accuracy.

Fig. A.2. To obtain a spectrum which is corrected for instrumental effects for the spectral calibration, we subtract from the RAW spectrum the SKY spectrum and then divide by the P2VM spectrum (blue dashed curve). The resulting instrument-corrected spectrum (red curve) is compared with an atmospheric transmission spectrum (green curve: R = 4000; black curve: R = 35).

References
Baize, P., & Romani, L. 1946, Ann. Astrophys., 9, 13 Baldwin, J. A., Ferland, G. J., Martin, P. G., et al. 1991, ApJ, 374, 580 Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 Conti, P. S. 1972, ApJ, 174, L79 Demircan, O., & Kahraman, G. 1991, Ap&SS, 181, 313 Docobo, J. A. 1985, Celest. Mech., 36, 143 Genzel, R., & Stutzki, J. 1989, ARA&A, 27, 41 GÑmez, L., Rodriguez, L. F., Loinard, L., et al. 2008, ArXiv e-prints, 805 GÑmez, L., RodrÌguez, L. F., Loinard, L., et al. 2005, ApJ, 635, 1166 Heintz, W. D. 1978, Geophys. Astrophys. Monogr., 15 Herbig, G. H., & Terndrup, D. M. 1986, ApJ, 307, 609 Hilditch, R. W. 2001, An Introduction to Close Binary Stars, ed. R. W. Hilditch (Cambridge, UK: Cambridge University Press), ISBN 0521241065

transmission spectra provided by the Gemini observatory5 (see Fig. A.2, green curve). We manually align the spectra using the
These ATRAN transmission spectra can be found on the website http://www.gemini.edu/sciops/telescopes- and- sites/ observing- condition- constraints/transmission- spectra.
5


S. Kraus et al.: Tracing the high-eccentricity binary 1 Ori C through periastron passage
Hillenbrand, L. A. 1997, AJ, 113, 1733 Hirota, T., Bushimata, T., Choi, Y. K., et al. 2007, PASJ, 59, 897 Hofmann, K.-H., & Weigelt, G. 1986, A&A, 167, L15 Hofmann, K.-H., & Weigelt, G. 1993, A&A, 278, 328 Jeffries, R. D. 2007, MNRAS, 376, 1109 Kratter, K. M., & Matzner, C. D. 2006, MNRAS, 373, 1563 Kraus, S., Schloerb, F. P., Traub, W. A., et al. 2005, AJ, 130, 246 Kraus, S., Balega, Y. Y., Berger, J.-P., et al. 2007, A&A, 466, 649 Kraus, S., Preibisch, T., & Ohnaka, K. 2008, ApJ, 676, 490 Labeyrie, A. 1970, A&A, 6, 85 Lohmann, A. W., Weigelt, G., & Wirnitzer, B. 1983, Appl. Opt., 22, 4028 Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 Mason, B. D., Gies, D. R., Hartkopf, W. I., et al. 1998, AJ, 115, 821 Mathieu, R. D. 1994, ARA&A, 32, 465 Mayne, N. J., & Naylor, T. 2008, MNRAS, 326 McCaughrean, M. J., & Stauffer, J. R. 1994, AJ, 108, 1382 Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 Morrell, N., & Levato, H. 1991, ApJS, 75, 965 Mozurkewich, D., Johnston, K. J., Simon, R. S., et al. 1991, AJ, 101, 2207 O'Dell, C. R. 2001, ARA&A, 39, 99 Patience, J., Zavala, R. T., Prato, L., et al. 2008, ApJ, 674, L97 Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1

207

Polyakov, A. M., & Tsivilev, A. P. 2007, Astron. Lett., 33, 34 Pourbaix, D. 1998, A&AS, 131, 377 RodrÌguez, L. F., Poveda, A., Lizano, S., & Allen, C. 2005, ApJ, 627, L65 Rubin, R. H., Simpson, J. P., Haas, M. R., & Erickson, E. F. 1991, ApJ, 374, 564 Sandstrom, K. M., Peek, J. E. G., Bower, G. C., Bolatto, A. D., & Plambeck, R. L. 2007, ApJ, 667, 1161 Schertl, D., Balega, Y. Y., Preibisch, T., & Weigelt, G. 2003, A&A, 402, 267 SimÑn-DÌaz, S., Herrero, A., Esteban, C., & Najarro, F. 2006, A&A, 448, 351 Stahl, O., Kaufer, A., Rivinius, T., et al. 1996, A&A, 312, 539 Stahl, O., Wade, G., Petit, V., Stober, B., & Schanne, L. 2008, A&A, 487, 323 Stassun, K. G., van den Berg, M., Feigelson, E., & Flaccomio, E. 2006, ApJ, 649, 914 Struve, O., & Titus, J. 1944, ApJ, 99, 84 Tan, J. C. 2004, ApJ, 607, L47 Tan, J. C. 2008, ArXiv e-prints, 807 Tatulli, E., Millour, F., Chelli, A., et al. 2007, A&A, 464, 29 Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 914 Vitrichenko, è. A. 2002, Astron. Lett., 28, 324 Weigelt, G. P. 1977, Opt. Commun., 21, 55 Weigelt, G., & Wirnitzer, B. 1983, Opt. Lett., 8, 389 Weigelt, G., Balega, Y., Preibisch, T., et al. 1999, A&A, 347, L15 Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481