Near-infrared interferometry of eta Carinae with
spectral resolutions of 1 500 and 12 000 using
AMBER/VLTI

G. Weigelt, S. Kraus, T. Driebe, R. G. Petrov, K.-H. Hofmann, F. Millour, O. Chesneau, D. Schertl, F. Malbet, J. D. Hillier, et al.

► To cite this version:

G. Weigelt, S. Kraus, T. Driebe, R. G. Petrov, K.-H. Hofmann, et al. Near-infrared interferometry of eta Carinae with spectral resolutions of 1 500 and 12 000 using AMBER/VLTI. Astronomy & Astrophysics - A&A, 2007, 464, pp.87-106. 10.1051/0004-6361:20065577. hal-00199354

HAL Id: hal-00199354

https://hal.science/hal-00199354v1

Submitted on 5 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Near-infrared interferometry of η Carinae with spectral resolutions of 1 500 and 12 000 using AMBER/VLTI*

G. Weigelt¹, S. Kraus¹, T. Driebe¹, R. G. Petrov², K.-H. Hofmann¹, F. Millour²,³, O. Chesneau⁴, D. Schertl¹,
F. Malbet³, D. J. Hillier⁶, T. Gull⁷, K. Davidson⁸, A. Domiciano de Souza²,⁴, P. Antonelli⁴, U. Beckmann¹, Y. Bresson⁴,
A. Chelli³, M. Dugué⁴, G. Duvert³, S. Gennari⁵, L. Glück³, P. Kern³, S. Lagarde⁴, E. Le Coarer³, F. Lisi⁵, K. Perraut³,
P. Puget³, F. Rantakyrö⁹, S. Robbe-Dubois², A. Roussel⁴, E. Tatulli³,⁵, G. Zins³, M. Accardo⁵, B. Acke³,¹⁶, K. Agabi²,
E. Altariba³, B. Arezki³, E. Aristidi², C. Baffa⁵, J. Behrend¹, T. Blöcker¹, S. Bonhomme⁴, S. Busoni⁵, F. Cassaing¹⁰,
J.-M. Clausse⁴, J. Colin⁴, C. Connot¹, A. Delboulbé³, P. Feautrier³, D. Ferruzzi⁵, T. Forveille³, E. Fossat², R. Foy¹¹,
D. Fraix-Burnet³, A. Gallardo³, E. Giani⁵, C. Gil³,¹⁷, A. Glentzlin⁴, M. Heiden¹, M. Heininger¹, O. Hernandez Utrera³,
D. Kamm⁴, M. Kiekebusch⁹, D. Le Contel⁴, J.-M. Le Contel⁴, T. Lesourd¹², B. Lopez⁴, M. Lopez¹², Y. Magnard³,
A. Marconi⁵, G. Mars⁴, G. Martinot-Lagarde¹²,⁴, P. Mathias⁴, P. Mège³, J.-L. Monin³, D. Mouillet³,¹⁸, D. Mourard⁴,
E. Nussbaum¹, K. Ohnaka¹, J. Pacheco⁴, C. Perrier³, Y. Rabbia⁴, S. Rebattu⁴, F. Reynaud¹³, A. Richichi¹⁴, A. Robini²,
M. Sacchettini³, M. Schöller⁹, W. Solscheid¹, A. Spang⁴, P. Stee⁴, P. Stefanini⁵, M. Tallon¹¹, I. Tallon-Bosc¹¹,
D. Tasso⁴, L. Testi⁵, F. Vakili², O. von der Lühe¹⁵, J.-C. Valtier⁴, M. Vannier²,⁹,¹⁹, N. Ventura³,
K. Weis²⁰, and M. Wittkowski¹⁴

(Affiliations can be found after the references)

Received 9 May 2006 / Accepted 1 August 2006

ABSTRACT

**Aims.** We present the first NIR spectro-interferometry of the LBV η Carinae. The observations were performed with the AMBER instrument of the ESO Very Large Telescope Interferometer (VLTI) using baselines from 42 to 89 m. The aim of this work is to study the wavelength dependence of η Car's optically thick wind region with a high spatial resolution of 5 mas (11 AU) and high spectral resolution.

**Methods.** The observations were carried out with three 8.2 m Unit Telescopes in the K-band. The raw data are spectrally dispersed interferograms obtained with spectral resolutions of 1500 (MR-K mode) and 12 000 (HR-K mode). The MR-K observations were performed in the wavelength range around both the He I 2.059 µm and the Bry 2.166 µm emission lines, the HR-K observations only in the Bry line region.

**Results.** The spectrally dispersed AMBER interferograms allow the investigation of the wavelength dependence of the visibility, differential phase, and closure phase of η Car. In the K-band continuum, a diameter of 4.0 ± 0.2 mas (Gaussian FWHM, fit range 28–89 m baseline length) was measured for η Car's optically thick wind region. If we fit Hillier et al. (2001, ApJ, 553, 837) model visibilities to the observed AMBER visibilities, we obtain 50% encircled-energy diameters of 4.2, 6.5 and 9.6 mas in the 2.17 µm continuum, the He I, and the Bry emission lines, respectively. In the continuum near the Bry line, an elongation along a position angle of 120° ± 15° was found, consistent with previous VINCI/VLTI measurements by van Boekel et al. (2003, A&A, 410, L37). We compare the measured visibilities with predictions of the radiative transfer model of Hillier et al. (2001), finding good agreement. Furthermore, we discuss the detectability of the hypothetical hot binary companion. For the interpretation of the non-zero differential and closure phases measured within the Bry line, we present a simple geometric model of an inclined, latitude-dependent wind zone. Our observations support theoretical models of anisotropic winds from fast-rotating, luminous hot stars with enhanced high-velocity mass loss near the polar regions.

**Key words.** stars: individual; η Carinae – stars: mass-loss, emission-line Be, circumstellar matter, winds, outflows – infrared: stars – techniques: interferometric

**1. Introduction**

The enigmatic object η Car is one of the most luminous and most massive (M ~ 100 M⊙) unstable Luminous Blue Variables suffering from an extremely high mass loss (Davidson & Humphreys 1997). Its distance is approximately 2300 ± 100 pc (Davidson & Humphreys 1997; Davidson et al. 2001; Smith 2006). η Car, which has been subject to a variety of studies over the last few decades, is surrounded by the expanding bipolar Homunculus nebula ejected during the Great Eruption in 1843. The

inclination of the polar axis of the Homunculus nebular with the line-of-sight is ~41° with the southern pole pointing towards us (Davidson et al. 2001; Smith 2006). The first measurements of structures in the innermost sub-arcsecond region of the Homunculus were obtained by speckle-interferometric observations (Weigelt & Ebersberger 1986; Hofmann & Weigelt 1988). These observations revealed a central object (component A) plus three compact and surprisingly bright objects (components B, C, and D) at distances ranging from approximately 0.1" to 0.2". HST observations of the inner 1" region (Weigelt et al. 1995) provided estimates of the proper motion of the speckle objects B, C, and D (velocity ~50 km s⁻¹; the low velocity suggests that the speckle objects are located within the equatorial plane), and follow-up HST spectroscopy unveiled their unusual spectrum

<footnote>* Based on observations collected at the European Southern Observatory, Paranal, Chile, within the AMBER guaranteed time programme 074.A-9025 and the VLTI science demonstration programme 074.A-9024.</footnote>
(Davidson et al. 1995). The central object (speckle object A) showed broad emission lines, while the narrow emission lines came from the speckle objects B, C, and D. Therefore, A is certainly the central object while B, C, and D are ejecta. Recent observations of η Car by Chesneau et al. (2005) using NACO and MIDI/VLTI revealed a butterfly-shaped dust environment at 3.74 and 4.05 µm and resolved the dusty emission from the individual speckle objects with unprecedented angular resolution in the NIR. Chesneau et al. also found a large amount of corundum dust peaked ~1" south-east of the central object.

Spectroscopic studies of the Homunculus nebula showed that the stellar wind of η Car is aspherical and latitude-dependent, and the polar axes of the wind and the Homunculus appear to be aligned (bipolar wind model; Smith et al. 2003). Using Balmer line observations obtained with HST/STIS, Smith et al. (2003) found a considerable increase of the wind velocity from the equator to the pole and that the wind density is higher in polar direction (parallel to the Homunculus; PA of the axis ~132°; Davidson et al. 2001) than in equatorial direction by a factor of ~2. van Boekel et al. (2003) resolved the optically thick, aspheric wind region with NIR interferometry using the VINCI/VLTI instrument. They measured a size of 5 mas (50% encircled-energy diameter), an axis ratio of 1.25 ± 0.05, and a position angle (PA) of the major axis of 134° ± 7°, and derived a mass-loss rate of $1.6 \times 10^{-3} M_{\odot} \text{yr}^{-1}$. The aspheric wind can be explained by models for line-driven winds from luminous hot stars rotating near their critical speed (e.g., Owocki et al. 1996, 1998). The models predict a higher wind speed and density along the polar axis than in the equatorial plane. In addition, van Boekel et al. showed that the K broad-band observations obtained with VINCI are in agreement with the predictions from the detailed spectroscopic model by Hillier et al. (2001).

The Hillier et al. (2001, 2006) model was developed to explain STIS HST spectra. The luminosity of the primary ($5 \times 10^6 L_{\odot}$) was set by observed IR fluxes (see discussion by Davidson & Humphreys 1997) and the known distance of 2.3 kpc to η Car. Any contribution to the IR fluxes by a binary companion was neglected. Modeling of the spectra was undertaken using CMFGEN, a non-LTE line blanketed radiative transfer developed to model stars with extended outflowing atmospheres (Hillier & Miller 1998). For the modeling of η Carinae, ions of H, He, C, N, O, Na, Mg, Al, Si, S, Ca, Ti, Cr, Mn, Fe, Ni, and Co were included. The mass loss was derived from the strength of the hydrogen lines and their associated electron scattering wings. Due to a degeneracy between the mass-loss rate and the He abundance, the H/He helium abundance ratio could not be derived, but was set at 5:1 (by number), which is similar to that found by Davidson et al. (1986) from nebula studies. CNO abundances were found to be consistent with those expected for full CNO processing. With the exception of Na (which was found to be enhanced by at least a factor of 2), the adoption of solar abundances for other metal species was found to yield satisfactory fits to the STIS spectra. A more recent discussion of the basic model, with particular reference to the UV and outer wind, is given by Hillier et al. (2006).

Because the wind is optically thick, the models are fairly insensitive to the radius adopted for the hydrostatic core (i.e., the radius at which the velocity becomes subsonic). One exception was the He I lines, which decreased in strength as the radius increased and, in general, were very sensitive to model details. Additional HST STIS observations show that the He I lines are strongly variable and blue-shifted throughout most of the 5.54-year variability period. These observations cannot be explained in the context of a spherical wind model. It now

appears likely that a large fraction of the He I line emission originates in the bow shock and an ionization zone, associated with the wind-wind interaction zone in a binary system (Davidson et al. 1999; Davidson 2001; Hillier et al. 2006; Nielsen et al. 2006). Consequently, the hydrostatic radius derived by Hillier et al. (2001) is likely to be a factor of 2 to 4 too small. Because the wind is so thick, a change in radius will not affect the Bry formation region, and it will only have a minor influence on the Bry continuum emitting region. If this model is correct, the He I emission will be strongly asymmetrical and offset from the primary star.

A variety of observations suggest that the central source of η Car is a binary. Damineli (1996) first noticed the 5.5-year periodicity in the spectroscopic changes of this object (see Damineli et al. 1997, 2000; Duncan et al. 1999; Ishibashi et al. 1999; Davidson et al. 1999, 2000; van Genderen et al. 2003; Steiner & Damineli 2004; Whitelock et al. 2004; Corcoran 2005; Weis et al. 2005). On the other hand, to date, the binary nature of the central object in η Car and its orbital parameters are still a matter of debate (see, e.g., Zanella et al. 1984; Davidson 1999, 2001; Davidson et al. 1999, 2000, 2005; Ishibashi et al. 1999; Smith et al. 2000; Feast et al. 2001; Ishibashi 2001; Pittard & Corcoran 2002; Smith et al. 2003; Martin et al. 2006).

The 1997.9 X-ray peak with the subsequent rapid drop to a few-month-long minimum was detected by RXTE (see Corcoran 2005). Then the first spectra with HST/STIS were obtained at 1998.0, demonstrating changes in both the central star and the aforementioned speckle objects (Davidson et al. 1999; Gull et al. 1999). Pittard & Corcoran (2002) demonstrated that the CHANDRA X-ray spectrum can be explained by the wind-wind collisions of the primary star ($\dot{M} = 2 \times 10^{-4} M_{\odot} \text{yr}^{-1}$ at 500 km s$^{-1}$) and a hot companion ($\dot{M} = 10^{-5} M_{\odot} \text{yr}^{-1}$ at 3000 km s$^{-1}$). Verner et al. (2005) used models calculated with the CLOUDY code to demonstrate that during the spectroscopic minimum, the excitation of the speckle objects is supported by the primary stellar flux, but that the UV flux of a hot companion consistent with an O7.5V, O9I, or early WN star was probably necessary to excite the speckle objects during the broad spectroscopic maximum.

In this paper we present the first spectro-interferometric K-band observations of η Car obtained with the VLTI beam-combiner instrument AMBER with medium and high spectral resolution and in the projected baseline range from 28 to 89 m.

The paper is organized as follows: in Sect. 2 we give an overview of the AMBER observations of η Car and describe the data reduction procedure in detail, and in Sect. 3, the analyses of the continuum data and the measurements within the Bry and He I lines are discussed individually.

## 2. AMBER observations and data processing

AMBER (Petrov et al. 2003, 2007) is the near-infrared (J, H, K band) beam-combiner instrument of ESO's Very Large Telescope Interferometer, which allows the measurement of visibilities, differential visibilities, differential phases, and closure phases (Millour et al. 2006; Petrov et al. 2005; Tatulli et al. 2007). AMBER offers three spectroscopic modes: low (LR mode; $R=\lambda/\Delta\lambda = 75$), medium (MR mode; $R=1500$), and high (HR mode; $R=12000$) spectral resolutions. The fibers in AMBER limit the field-of-view to the diameter of the fibers on the sky (~60 mas). In AMBER the light is spectrally dispersed using a prism or grating. The AMBER detector is a Hawaii array detector with $512 \times 512$ pixels.
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Spectrally dispersed AMBER/VLTI Michelson interferograms of η Car. The two panels show the spectrally dispersed fringe signal (IF) as well as the photometric calibration signals from the three telescopes (P1–P3) in high (HR, *upper panel*) and medium spectral resolution mode (MR, *lower panel*). In both panels, the bright regions are associated with the Doppler-broadened Bry emission line.</figcaption></figure>

Figure 1 shows two AMBER raw interferograms taken in the wavelength range around the Bry line in HR (top) and MR (bottom) mode. In the MR data sets, the Doppler-broadened Bry line covers ~8 spectral channels, whereas in HR mode, the line is resolved by ~50 spectral channels.

η Car was observed with AMBER on 2004 December 26, 2005 February 25, and 2005 February 26 with the three 8.2 m Unit Telescopes UT2, UT3, and UT4. With projected baseline lengths up to 89 m, an angular resolution of ~5 mas was achieved in the K band. As listed in Table 1, the MR-K observations were performed in the wavelength range around both the He I 2.059 µm and the Brγ 2.166 µm emission lines. The HR-K observations were only performed in a wavelength range around the Bry line. The widths of the wavelength windows of the obtained MR-K and HR-K observations are approximately 0.05 µm and 0.02 µm, respectively.

For the reduction of the AMBER data, we used version 2.4 of the *amdlib*<fnref n="1" /> software package. This software uses the P2VM (pixel-to-visibilities matrix) algorithm (Tatulli et al. 2007) in order to extract complex visibilities for each baseline and each spectral channel of an AMBER interferogram. From these three complex visibilities, the amplitude and the closure phase are

derived. While the closure phase is self-calibrating, the visibilities have to be corrected for atmospheric and instrumental effects. This is done by dividing the η Car visibility through the visibility of a calibrator star measured on the same night. In order to take the finite size of the calibrator star into account, the calibrator visibility is corrected beforehand through division by the expected calibrator star visibility (see Table 1). In the case of the MR measurement performed on 2005-02-25, the interferograms recorded on the calibrator contain only fringes corresponding to the shortest baseline (UT2-UT3). Thus, the η Car visibility for this night could only be calibrated for this shortest baseline.

Besides the calibrated visibility and the closure phase, the spectral dispersion of AMBER also allows us to compute differential observables; namely the differential visibility and the differential phase (Millour et al. 2006; Petrov et al. 2005; Tatulli et al. 2007). These quantities are particularly valuable, as they provide a measure of the spatial extent and spatial offset of the line-emitting region with respect to the continuum emission. Since the measured complex visibilities are affected by wavelength-dependent atmospheric piston (optical path difference), the piston has to be estimated and subtracted. This was done using the *ammyorick*<fnref n="1" /> tool (version 0.56).

Since a large fraction of the interferograms is of low contrast (probably due to vibration; see Malbet et al. 2007), we removed a measurement from the data sets if (a) the intensity ratio of two of the photometric channel signals is larger than 4 (a large ratio means that the interferograms are very noisy since the signal is very weak in one channel) or (b) it belongs to the 70 percent of the interferograms with the lowest fringe contrast SNR (with the SNR defined as in Tatulli et al. 2007). In order to optimize the selection for each baseline of the telescope triplet, both of these criteria are applied for each telescope pair individually. Furthermore, the first 10 frames in each new sequence of recorded interferograms are removed since they are degraded by electronic noise.

Figures 2 and 4 show the spectra as well as the wavelength dependence of the visibilities, differential visibilities, differential phases, and closure phases derived from the AMBER interferograms for the observations around the Brγ and He I emission lines. The uv coverage of the observations is displayed in Fig. 3.

The η Car spectra were corrected for instrumental effects and atmospheric absorption through division by the calibrator spectrum. For the HR 2005-02-26 measurement, we found that the calibrator itself (L Car) shows prominent Bry line absorption (see Fig. A.1). Therefore, we had to remove this stellar line by linear interpolation before the spectrum could be used for the calibration. The wavelength calibration was done using atmospheric features, as described in more detail in Appendix A.

In order to test the reliability of our results, we split each of the raw data sets into 5 subsets, each containing the same number of interferograms. The results obtained with these individual subsets allowed us to test that the major features detected in the visibility, differential visibility, differential phase, and closure phase are stable, even without any frame selection applied. As an exception, we found that for a small wavelength range of the HR 2005-02-26 data set (hatched areas in the two lower right panels of Fig. 2), the differential phase corresponding to the middle and longest baselines and the closure phase vary strongly within the subsets and are, therefore, unreliable. This is likely due to the very low visibility value on these two baselines, resulting in a low fringe SNR within this wavelength range. Furthermore, with this method we found that the differential visibility, differential phase, and closure phase extracted

<footnote n="1">¹ This software package is available from
http://amber.obs.ujf-grenoble.fr</footnote>
<table><caption>**Table 1.** Summary of the AMBER η Car observations using the UT2, UT3, and UT4 telescopes.</caption><thead><tr><th>Date [UT]</th><th>Time [UT] start</th><th>Time [UT] end</th><th>Orbital phase<sup>e</sup></th><th>Spectral mode</th><th>Line within spectral range</th><th>DIT<sup>d</sup> [ms]</th><th>N<sub>η Car</sub><sup>f</sup></th><th>Calibrator</th><th>N<sub>calib.</sub><sup>g</sup></th><th>Calibrator uniform disk diameter [mas]</th></tr></thead><tbody><tr><td rowspan="2">2004 Dec. 26</td><td>07:52</td><td>08:16</td><td>0.267</td><td>MR-K</td><td>Bry</td><td>40</td><td>7500</td><td>HD 93030</td><td>5000</td><td>0.39<sup>a</sup></td></tr><tr><td>08:19</td><td>08:32</td><td>0.267</td><td>MR-K</td><td>He I 2.059 µm</td><td>40</td><td>5000</td><td>HD 93030</td><td>5000</td><td>0.39<sup>a</sup></td></tr><tr><td rowspan="2">2005 Feb. 25</td><td>04:33</td><td>04:43</td><td>0.298</td><td>MR-K</td><td>Bry</td><td>50</td><td>5000</td><td>HD 89682</td><td>2500</td><td>3.08<sup>b</sup></td></tr><tr><td>04:55</td><td>05:05</td><td>0.298</td><td>MR-K</td><td>He I 2.059 µm</td><td>50</td><td>5000</td><td>HD 89682</td><td>2500</td><td>3.08<sup>b</sup></td></tr><tr><td>2005 Feb. 26</td><td>08:16</td><td>08:57</td><td>0.299</td><td>HR-K</td><td>Bry</td><td>82</td><td>7500</td><td>L Car</td><td>2500</td><td>2.70<sup>c</sup></td></tr></tbody></table>

<sup>a</sup> Uniform disk (UD) diameter estimated using the method described by Dyck et al. (1996). <sup>b</sup> UD diameter taken from the CHARM2 catalog (Richichi et al. 2005). <sup>c</sup> UD diameter of L Car at the time of the AMBER high-resolution observations derived from the limb-darkened diameter $d_{LD} = 2.80$ mas at the L Car pulsation phase $\phi = 0.0$ (Kervella et al. 2004b, 2006) and $d_{UD}/d_{LD} = 0.966$ (Kervella et al. 2004a).

<sup>d</sup> Detector integration time per interferogram. <sup>e</sup> The orbital phase was computed assuming a zero point at JD 2 450 800.0 and a period of 2024 days (Corcoran 2005). <sup>f</sup> Number of η Car interferograms. <sup>g</sup> Number of calibrator interferograms.

from the MR 2005-02-25 He I data set are very noisy and not reliable. Therefore, these differential quantities and closure phases were dropped from our further analysis.

The subsets were also used to compute statistical errors. We estimated the variance for each spectral channel and derived formal statistical errors for both the continuum and line wavelength ranges. In each panel of Figs. 2 and 4, we show two types of error bars corresponding to these regions, which not only take these statistical errors but also a systematic error (e.g. resulting from an imperfect calibration) into account.

## 3. Observational results and interpretation

### 3.1. Comparison of the observed wavelength dependence of the visibility with the NLTE radiative transfer model of Hillier et al. (2001)

For the analysis presented in this chapter, we used the AMBER data sets from 2004 Dec. 26 and 2005 Feb. 25 and 26, presented in Figs. 2 and 4, and compared the AMBER visibilities and spectra with the NLTE radiative transfer model of Hillier et al. (2001). To directly compare the AMBER measurements with this model, we derived monochromatic model visibilities for all wavelengths between 2.03 and 2.18 µm (with $\Delta\lambda = 10^{-4}$ µm) from the model intensity profiles, assuming a distance of 2.3 kpc for η Car. The comparison is visualized in Fig. 5 for the individual AMBER HR and MR measurements. The first row displays the AMBER and model spectra, while all other panels show the AMBER and model visibilities for the different projected baselines. We note that for the comparison shown in Fig. 5, we used the original model of Hillier et al. (2001) without any additional size scaling or addition of a background component.

As the figure reveals, the NLTE model of Hillier et al. (2001) can approximately reproduce the AMBER continuum observations for all wavelengths (i.e. 2.03–2.18 µm) and all baselines. Moreover, the wavelength dependence of the model visibilities inside the Bry line is also similar to the AMBER data. There is a slight tendency for the model visibilities in the Bry line to be systematically lower, which can be attributed to the overestimated model flux in the line. On the other hand, there is an obvious difference in the wavelength dependence of the visibility across the He I line between the observations and the model predictions. This difference probably indicates that the primary wind model does not completely describe the physical origin and, hence, the spatial scale of the He I line-forming region. The discrepancy is possibly caused by additional He I emission from the wind-wind interaction zone between the binary components and by

the primary's ionized wind zone caused by the secondary's UV light illuminating the primary's wind (e.g., Davidson et al. 1999; Davidson 2001; Pittard & Corcoran 2002; Steiner & Damineli 2004; Hillier et al. 2006; Nielsen et al. 2006; Martin et al. 2006), as discussed in Sects. 1, 3.4.2, and 3.7.2 in more detail.

Figure 6 shows the AMBER and model visibilities as a function of spatial frequency and the corresponding model center-to-limb intensity variations (CLVs) for seven selected wavelengths (2 continuum wavelengths; center, blue-shifted, and red-shifted wings of Bry emission; center of both He I emission and absorption). As Fig. 6 reveals, at several wavelengths we find a very good agreement between the visibilities measured with AMBER and the visibilities predicted by the model of Hillier et al. (2001). This is especially true for the continuum data (upper two panels).

From the model CLVs, FWHM model continuum diameters of 2.24 mas and 2.33 mas can be derived for $\lambda = 2.040$ and $2.174$ µm, respectively. If we allow for a moderate rescaling of the size of the model, we find that the best $\chi^2$ fit at both continuum wavelengths can be obtained with scaling factors of 1.015 and 1.00, respectively. This means that the model size has to be increased by only 1.5% at $\lambda = 2.040$ µm and that the best fit at $2.174$ µm is indeed obtained with the original Hillier model with a scaling factor of 1.0. Thus, taking the slight rescaling for the best $\chi^2$ fit into account, we can conclude that, based on the NLTE model from Hillier et al. (2001, 2006) and the AMBER measurements, the apparent FWHM diameters of η Car in the K-band continuum at $\lambda = 2.040$ µm and $2.174$ µm are 2.27 mas and 2.33 mas, respectively (see Table 2), corresponding to a physical size of approximately 5 AU.

Since the deviations between the model and the measurements are larger in the case of the Bry and He I line data (lower 5 panels in Fig. 6), the scaling factors corresponding to the best $\chi^2$ fit in the lines show stronger deviations from unity. For the Bry emission line, we find scaling factors of 0.74, 0.76, and 0.78 for $\lambda = 2.1661, 2.1635$, and $2.1669$ µm, corresponding to FWHM diameters of 1.83, 9.52 and 2.02 mas (see Table 2).

For the He I emission line, rescaled models with scaling factors of 1.24 and 1.11 provide the best $\chi^2$ fit for the peaks of the emission and absorption within the He I line ($\lambda = 2.058$ and $2.056$ µm), resulting in FWHM diameters of 4.24 and 4.19 mas, respectively.

In addition to the inner CLV core, at several wavelengths, the CLVs show a very extended wing corresponding to the extended Bry and He I line emission regions. Since the intensity in the wing is much lower than 50% of the peak intensity, the FWHM diameter is not very sensitive to this part of the CLV.
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** AMBER observables derived from our η Car data around the Bry line for three independent measurements (Left: MR, 2004 December 26, Middle: MR, 2005 February 25, Right: HR, 2005 February 26). The first row shows the continuum-normalized spectra as extracted from the interferometric channels, followed by the derived calibrated visibilities and the differential visibilities. In the fourth and fifth row, the differential phase and the closure phase are presented. In the spectra we mark the wavelength regimes, which we defined as continuum for our analysis (shaded regions). The vertical grey line marks the rest-wavelength of Bry ($\lambda_{vac}$ = 2.1661 µm; the small correction due to the system velocity of –8 km s⁻¹ (Smith 2004) has been neglected). We show different error bars within each panel: the left error bars correspond to the total (including statistical and systematic) error estimated for the continuum wavelength range, and the error bar towards the right visualizes the total error for the wavelength range within the line. For the HR 2005-02-26 measurement, data splitting showed that for a small wavelength range (hatched areas in the two lower right panels), the differential phase for the longest and middle baseline as well as the closure phase become very noisy and are therefore not reliable. Furthermore, the HR differential phase of the longest baseline is noisy at all wavelengths. See Sect. 2 for further details.</figcaption></figure>

In other words, in the case of CLVs with multiple or very extended components, a *FWHM* diameter can be quite misleading. In such a case, it seems to be more appropriate to use, for instance, the diameter measured at 10% of the peak intensity

($d_{10%}$) or the 50% encircled-energy diameter ($d_{50\% EED}$). For example, at $\lambda = 2.1661$ µm we obtain $d_{10\%} = 9.39$ mas and $d_{50\% EED} = 9.58$ mas, while for the continuum at 2.174 µm we find $d_{10\%} = 5.15$ mas and $d_{50\% EED} = 4.23$ mas. Thus, based
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** uv coverage of the AMBER η Car observations. The data obtained with medium spectral resolution in Dec. 2004 and Feb. 2005 are indicated by red squares and green bullets, respectively, while the high-resolution measurements are shown as blue triangles. Filled symbols denote observations around the Bry line, and open symbols denote observations around the He I line.</figcaption></figure>

on $d_{50\%\text{EED}}$, η Car appears ~2.2 times larger at $\lambda = 2.1661$ µm compared to the continuum at $\lambda = 2.174$ µm. The best-fit model diameters at the other wavelengths are listed in Table 2. The errors of the diameter measurements are ±4% for the two continuum diameters and ±10% for the line diameters, derived from the visibility errors and the uncertainty of the fitting procedure.

## 3.2. Continuum visibilities

### 3.2.1. Comparison of the continuum visibilities with the Hillier et al. (2001) model predictions

The comparison of the AMBER continuum visibilities with the NLTE model from Hillier et al. (2001, 2006) is shown in the two upper left panels of Fig. 6 for the continuum near the He I 2.059 µm and Bry 2.166 µm emission lines (the exact wavelengths are described in Fig. 6). Taking a slight rescaling into account, we concluded in the previous section that, based on the NLTE model from Hillier et al. (2001) and the AMBER measurements, the apparent 50% encircled-energy diameters $d_{50\%\text{EED}}$ of η Car in the K-band continuum at $\lambda = 2.040$ µm and 2.174 µm are 3.74 mas and 4.23 mas, respectively (see Table 2). These diameters are in good agreement with the 50% encircled-energy K-band diameter of 5 mas reported by van Boekel et al. (2003).

For comparison, we also fitted the AMBER visibilities with simple analytical models such as Gaussian profiles, as described in more detail in Sect. B in the Appendix. From a Gaussian fit of the AMBER visibilities, we obtain a FWHM diameter of $d_{\text{Gauss}} \sim 4.0 \pm 0.2$ mas in the K-band continuum. As outlined in Sect. B, the diameter value strongly depends on the range of projected baselines used for the fit, since a Gaussian without an additional background component is not a good representation of the visibility measured with AMBER. As discussed in Appendix B, using a Gaussian fit with a fully resolved background component as a free parameter results in a best fit with a 30% background flux contribution (see also Petrov et al. 2005).

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Similar to Fig. 2, but showing the MR measurement from 2004 December 26 covering the region around the He I line. The vertical grey line marks the He I rest-wavelength ($\lambda_{\text{vac}} = 2.0586$ µm).</figcaption></figure>

### 3.2.2. Comparison of the VINCI and AMBER continuum visibilities

In Fig. 6 (left, second row) displaying the averaged Bry continuum data, the visibilities of η Car obtained with VINCI/VLTI are
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Comparison of the AMBER spectra and visibilities with the NLTE model predictions of Hillier et al. (2001). The figure displays the spectra (upper row) and visibilities (lower three rows, see labels for projected baselines) of the four AMBER measurements (green lines) and the corresponding data of the Hillier et al. NLTE model (red lines). The errors of the AMBER continuum and line visibility measurements are indicated by the two vertical error bars (see Figs. 2 and 4; the left bar is the continuum error bar), and the uncertainty of the AMBER wavelength calibration is indicated by the horizontal error bar. As the figure shows, we find good agreement between the AMBER data and the model predictions for the continuum visibilities as well as the shape and depth of the visibility inside the Bry line. In the case of the He I line, the wavelength dependence of the model visibility inside the line differs considerably from the AMBER measurements, indicating a different physical process involved in the line formation (see Sect. 3.1). The He I wavelength shift can be attributed to a combination of both Doppler shift (e.g. Nielsen et al. 2006; Hillier et al. 2006) and uncertainties in the wavelength calibration of the AMBER data. Note that no additional scaling has been applied to the Hillier et al. model. The model spectra and visibilities have a spectral resolution ($\lambda/\Delta\lambda \sim 20\,000$) comparable to the HR measurements.</figcaption></figure>

shown in addition to the AMBER data. These VINCI measure-
ments were carried out in 2002 and 2003 using the 35 cm test
siderostats at the VLTI with baselines ranging from 8 to 62 m
(for details, see van Boekel et al. 2003). Like AMBER, VINCI
is a single-mode fiber instrument. Therefore, its field-of-view is
approximately equal to the Airy disk of the telescope aperture
on the sky, which is ~1.4" in the case of the siderostats. From
the VINCI measurements and using only the 24 m baseline data,
van Boekel et al. (2003) derived a *FWHM* Gaussian diameter of
$d_{\text{Gauss}} \sim 7$ mas for the wind region of $\eta$ Car. At first glance, this
diameter measurement seems to contradict the $d_{\text{Gauss}} \sim 4.0$ mas
*FWHM* diameter derived from the AMBER data. This is not
the case, however, since the diameter fit is very sensitive to the
baseline (or spatial frequency) fit range, because a Gaussian is
not a good representation of the visibility curve at all, as can
be seen in Fig. B.1. If only the VINCI data points are fitted,
which have spatial frequencies <60 cycles/arcsec (correspond-
ing to projected baselines <28 m), $d_{\text{Gauss}} \sim 7$ mas provides the
best fit. On the other hand, if the data point at 136 cycles/arcsec

(corresponding to a projected baseline of ~62 m) is included
in the fit, we obtain $d_{\text{Gauss}} \sim 4.3$ mas (see also the discussion
Sect. B). Thus, when using comparable baseline ranges for the
Gaussian fits, there is good agreement between the AMBER and
VINCI measurements.

To account for the background contamination of the VINCI
data caused by nebulosity within VINCI's large 1.4" field-of-
view (in which, for instance, all speckle objects B, C, and D
are located), van Boekel et al. introduced a background com-
ponent (derived from NACO data) providing 55% of the total
flux. Adding this background component to the model of Hillier
et al. (2001), they found a good match between the model and
the observations. Since our AMBER observations were carried
out with the 8.2 m Unit Telescopes of the VLTI, the field-of-
view of the AMBER observations was only ~60 mas. Thus,
the background contamination of the AMBER data can be ex-
pected to be much weaker, if not negligible, compared to the
VINCI measurements. To check this, we first performed a fit
of the Hillier et al. (2001) model, which not only contains the
<figure><img src="image_7.png" /></figure>

**Fig. 6.** Left: comparison of the AMBER visibilities (filled green squares; baseline range 28–89 m) as a function of spatial frequency with the NLTE model predictions of Hillier et al. (2001) (solid red lines) for two continuum wavelengths (upper two panels; see labels for the exact wavelengths), the central wavelength of the Bry emission line (third), wavelengths in the blue and red wing of the Bry emission line (fourth and fifth row; at these wavelengths, the AMBER data show the strongest differential and closure phase signals), and the central wavelengths of the HeI emission and P Cygni absorption (lower two panels). The blue triangles are the background-corrected VINCI K-band measurements from van Boekel et al. (2003). Right: center-to-limb variation (CLV; i.e. intensity as a function of angular radius) of the monochromatic Hillier et al. (2001) NLTE models for the wavelengths indicated by the labels. As the figures show, the typical FWHM diameter of the models is of the order of ~2–4 mas for an assumed distance of 2.3 kpc for η Car. See text for details.

size scaling as a free parameter, but also a fully resolved background component. As we expected, we found the best fit (smallest $\chi^2$) with no background contamination<fnref n="2" />. Therefore, when we finally compared the AMBER observations with the model from

Hillier et al. (2001), we did not introduce a background component. In Fig. 6 (second row, left) we plot both the AMBER visibilities (no background correction required) plus the background-corrected VINCI data (assuming a 55% background contribution; blue triangles). As can be seen from the figure, these VINCI points nicely match the AMBER data and the corresponding fit of the NLTE model from Hillier et al. (2001).

<footnote n="2">² An additional argument in favor of only a very faint background contribution in the AMBER UT observation can be found in the shape of the high spectral resolution line: the light from the speckle objects B, C, and D is produced in areas with velocities smaller than 50 km s⁻¹. Therefore, it produces a narrow emission line which should appear in</footnote>

<footnote>the center of the broad Bry line. Just looking at the shape of the line, it can be concluded that such an effect is negligible.</footnote>
<table><caption>**Table 2.** Diameters obtained by fitting Hillier et al. (2001) model visibilities to the measured AMBER visibilities. Errors are ±4% for the diameters in the continuum and ±10% in the lines (see text).</caption><thead><tr><th>Spectral region</th><th>Wavelength [µm]</th><th>d<sub>FWHM</sub> [mas]</th><th>d<sub>10%</sub> [mas]</th><th>d<sub>50% EED</sub> [mas]</th></tr></thead><tbody><tr><td>continuum</td><td>2.0400</td><td>2.27</td><td>4.85</td><td>3.74</td></tr><tr><td>continuum</td><td>2.1740</td><td>2.33</td><td>5.15</td><td>4.23</td></tr><tr><td>Bry (center)</td><td>2.1661</td><td>1.83</td><td>9.39</td><td>9.58</td></tr><tr><td>Bry (blue wing)</td><td>2.1635</td><td>9.52</td><td>16.46</td><td>9.60</td></tr><tr><td>Bry (red wing)</td><td>2.1669</td><td>2.02</td><td>9.61</td><td>9.78</td></tr><tr><td>He I (absorption)</td><td>2.0560</td><td>4.24</td><td>8.22</td><td>5.36</td></tr><tr><td>He I (emission)</td><td>2.0580</td><td>4.19</td><td>4.30</td><td>6.53</td></tr></tbody></table>

$d_{FWHM} = FWHM$ diameter; $d_{10\%} = \text{diameter measured at 10\% peak intensity}; d_{50\% EED} = 50\% \text{ encircled-energy diameter.}$

Therefore, from the analysis of the continuum data, we can conclude that the background contamination in the AMBER measurements is negligible and that the AMBER measurements are in good agreement with both the previous VINCI measurements and the model predictions from Hillier et al. (2001).

### 3.3. Elongated shape of the continuum intensity distribution

To look for detectable elongations of the continuum intensity distribution, we fitted an elliptically stretched 2-D version of the radiative transfer model visibilities from Hillier et al. (2001) to the measured visibilities. Our best $\chi^2$ fit reveals a projected axis ratio of $\xi = 1.18 \pm 0.10$ and PA = $120 \pm 15^\circ$. Comparison with the results found by van Boekel et al. (2003) shows that the projected axis ratio $\xi$ derived from the AMBER data is in basic agreement with the K-broad-band values of $\xi = 1.25 \pm 0.05$ and PA = $138 \pm 7^\circ$ from van Boekel et al. (2003).

We also studied the elongation inside the Bry emission line at $\lambda = 2.166 \mu m$, following the same procedure as in the continuum; i.e., we fitted an elliptically stretched 2-D version of the Hillier et al. model shown in Fig. 6 to the AMBER data. However, since the global shape of the model function at $\lambda = 2.166 \mu m$ shows stronger deviations from the measurements than in the continuum, the elongation determination suffers from larger uncertainties, resulting in large error bars of the fit parameters. For instance, for $\lambda = 2.166 \mu m$ we obtained $\xi = 1.66 \pm 0.60$ and PA = $81 \pm 40^\circ$ from the best ellipse fit.

The 2-D ellipse fitting was also performed for the continuum near the He I emission line and in the center of the He I line ($\lambda = 2.057 \mu m$), where our model fits give an axis ratio of $\xi = 1.35 \pm 0.30$ and a PA of the major axis of $98 \pm 40^\circ$ in the continuum, and $\xi = 1.74 \pm 0.60$ and PA = $159 \pm 40^\circ$ in the center of the He I emission line. It should be noted that for the He I line region, only four visibility points are available, covering the small PA range of only $60^\circ$. Because of this limited number of data points and the small PA coverage, we conclude that the He I elongation measurements in the continuum as well as the line region are not reliable and abandoned in the further elongation analysis of the He I data.

From the K-band VINCI data, van Boekel et al. (2003) derived a PA of $138 \pm 7^\circ$ for the major axis, very well aligned with the Homunculus (132°, Davidson et al. 2001) and in agreement with our results (PA = $120 \pm 15^\circ$). Van Boekel's and our continuum elongation measurements favor the physical model according to which $\eta$ Car exhibits an enhanced mass loss in polar direction as proposed, for instance, by Owocki et al. (1996, 1998) or Maeder & Desjacques (2001) for stars rotating close to

their critical rotation speed. Axis ratios of the order of 1.2 appear reasonable in the context of such polar-wind models. Suppose, for example, that the wind's polar/equatorial density ratio is 2 at any given radius $r$, as reported by Smith et al. (2003) to explain latitude-dependent changes in the Balmer line profiles. Relevant absorption and scattering coefficients have radial dependencies between $n_e \sim r^{-2}$ (Thomson scattering) and $n_e^2 \sim r^{-4}$ (most forms of thermal absorption and emission). A meridional map of projected optical thickness through the wind would show cross-sections of prolate spheroids, correlated with the appearance of the configuration. With the radial dependencies and polar/equatorial density ratio mentioned above, these spheroids have axial ratios between about 1.2 and 1.4; i.e., appreciably less than 2. Viewed from an inclination angle $i \approx 45^\circ$ (Davidson et al. 2001), the apparent (projected) axis ratios are between 1.1 and 1.2. This is merely one example, and we have omitted many details, but it illustrates that the polar/equatorial density ratio is around 2, in agreement with Smith et al. (2003).

Finally, Smith et al. (2003) suggested that the stellar wind should become basically spherical during an event at periastron. This prediction can be tested if AMBER/VLTI data are obtained at the next periastron passage.

### 3.4. Continuum-corrected visibilities

#### 3.4.1. Continuum-corrected visibility in the Bry emission line

To investigate the brightness distribution in the Bry line in more detail, we tried to disentangle the continuum and pure line emission from both the AMBER data as well as the model data to derive the size of the pure Bry line-emitting region. Since the visibility measured inside an emission or absorption line is the composite of a pure line component and an underlying continuum, the measured line visibility, $V_{tot}$ (see Fig. 7 top), has to be corrected for the continuum contribution to obtain the visibility $V_{line}$ of the line emitting (absorbing) region. As discussed in Malbet et al. (2007), $V_{line}$ can be calculated if the continuum level within the line is known. If the continuum within the line is equal to the continuum level outside the line for an optically thin environment, as assumed in Fig. 7, we obtain:

$$V_{\text{line}} = \frac{V_{\text{tot}} \cdot F_{\text{tot}} - V_{\text{cont}} \cdot F_{\text{cont}}}{F_{\text{line}}} \quad (1)$$

with $F_{tot} = F_{cont} + F_{line}$ being the total measured flux and $V_{tot}$ being the measured visibility (also see Fig. 7 top for illustration). As outlined in the Appendix, taking a non-zero differential phase $\Phi'$ into account introduces an additional term in Eq. (1), leading to

$$V_{\text{line}} = \frac{\sqrt{|V_{\text{tot}}F_{\text{tot}}|^2 + |V_{\text{cont}}F_{\text{cont}}|^2 - 2V_{\text{tot}}F_{\text{tot}}V_{\text{cont}}F_{\text{cont}} \cdot \cos \Phi'|}}{F_{\text{line}}} \quad (2)$$

We applied Eq. (2) to the line-dominated AMBER data in the 2.155–2.175 $\mu m$ Bry wavelength range to derive the continuum-corrected visibility of the region emitting the Bry line radiation. One big uncertainty in this correction is the unknown continuum flux within the line. Due to intrinsic absorption, the continuum flux might be considerably lower within the line than measured at wavelengths outside the line. Especially the blue-shifted wing of the Bry emission line might be affected by the P Cygni-like absorption, as discussed below. In case such an absorption component is present, our continuum-correction would overestimate
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Top: the solid and dashed red lines show the continuum-corrected visibility inside the Bry line of the HR measurement with the shortest projected baseline (29 m). The continuum correction follows Eq. (2). Except for the phase, all quantities entering Eq. (2) are indicated in the figure. We assume $F_{\text{tot}} = F_{\text{cont}} + F_{\text{line}}$ and that $F_{\text{cont}}$ inside the line is approximately equal to the level outside the line. The continuum-corrected visibility in the blue-shifted region of the emission line is shown with a dashed line since it is highly uncertain due to the presence of the P Cygni-like absorption component (see text for details). Bottom: continuum-corrected AMBER visibilities (filled blue squares) in the red region of the Bry line as a function of spatial frequency. To derive the visibilities according to Eq. (2), the data in the wavelength range 2.1661–2.1670 µm (red-shifted line region) were averaged before the continuum correction. The solid red line shows the continuum-corrected model visibility. Obviously, observations at shorter baselines are needed to further test the model predictions.</figcaption></figure>

the size of the line-emitting region. The influence of this effect on the red-shifted wing of the line is likely to be much smaller, as P Cygni-like absorption mainly affects the blue-shifted emission.

The visibility across the Bry line is shown in Fig. 7 (top) for the HR data corresponding to the shortest projected baseline. The results for the other data sets are similar. From Fig. 7 (top), one can see that after the subtraction of a continuum contribution equal to the continuum outside the line, the visibility reaches very small values in the center of the emission line. This means that the pure line-emitting region is much larger than the region providing the continuum flux.

Figure 7 (top) shows a strong asymmetry between the blue- and red-shifted part of the visibility in the line with respect to the spectrum. While the visibility ($V_{\text{tot}}$ as well as $V_{\text{line}}$) rises concomitantly with the drop of the line flux on the red side, the

situation is very different on the blue side line center. In agreement with the model predictions from Hillier et al. (2001), this indicates the existence of a P Cygni-like absorption component in this wavelength region. In fact, at $\lambda = 2.1625$ µm, we see a small dip in the Bry spectrum in both the model spectrum from Hillier et al. (2001) and the HR AMBER observations. If such an absorption component is present, it can explain the asymmetric behaviour of the line visibility with respect to the spectrum. The P Cygni absorption in the blue wing of the Bry line makes the continuum correction of the visibility uncertain for wavelengths shorter than the central wavelength $\lambda_c$ of the emission line. Because of this uncertainty, in Fig. 7 (top) the continuum corrected visibility is shown with a dashed line for $\lambda < \lambda_c$, and the following discussion is restricted to the red-shifted region of the Bry line emission.

The continuum-corrected AMBER visibilities in the red-shifted region of the Bry line are displayed in Fig. 7 (bottom) for all data sets. To derive the visibilities in the red region, the data in the wavelength range 2.1661–2.1670 µm were averaged before the continuum correction. To now compare the continuum-corrected AMBER visibilities with the model predictions (2.1661–2.1670 µm), we constructed a model intensity profile of the pure Bry emission line region by subtracting the Hillier et al. intensity profile of the nearby continuum from the combined line + continuum profile.

As Fig. 7 illustrates, the model prediction is in agreement with the low visibilities found for spatial frequencies beyond 60 cycles/arcsec. On the other hand, the figure also clearly indicates that measurements at smaller projected baselines are needed to further constrain the Hillier et al. model in the line-emitting region. With the baseline coverage provided by the current AMBER measurements, we obtain a *FWHM* diameter of ≥15.4 mas (lower limit) for the (continuum-corrected) line-emitting region in the red line wing.

### 3.4.2. Continuum-corrected visibility in the He I emission line

As can be seen in Fig. 4, the AMBER spectrum of the He I line shows a P Cygni-like profile with a prominent absorption and emission component. This is in agreement with earlier findings by Smith (2002) from long-slit spectroscopy using OSIRIS on the CTIO 4m telescope.

To estimate the spatial scale of the region emitting the He I emission line, we followed the same approach as outlined in the previous section for the Bry line; i.e., we first applied Eq. (2) and then compared the continuum-corrected visibility with the continuum-corrected radiative transfer model of Hillier et al. (2001). Figure 8 (top) shows the measured flux and visibility for the MR-2004-12-26 measurement with the shortest projected baseline (43 m) as well as the continuum-corrected visibility across the He I emission component (solid red line). Because of the P Cygni-like absorption component, the continuum subtraction is highly uncertain in the blue region of the emission line (dashed red line in Fig. 8), as already discussed in the context of the Bry line in Sect. 3.4.1.

In Fig. 8 (bottom), the continuum-corrected visibility of all AMBER data in the red region of the He I emission line (averaged over the wavelength range 2.057–2.058 µm) is shown as a function of spatial frequency. As the figure reveals, similar to the Bry emission, the visibilities inside the He I emission line region reach rather low values. As the comparison shows, the line visibilities predicted by the model are much higher than the line visibilities measured with AMBER, indicating that the size of the line-emitting region in the model is too small. Rescaling of
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Visibility in the He I line. Top: the figure is similar to Fig. 7 (top) but displays the MR-2004-12-26 data with the shortest projected baseline (43 m). As in the case of the Bry emission, the figure reveals that the region dominated by line emission is fully resolved by the AMBER measurements. As in Fig. 7 (top), the continuum-corrected visibility is shown with a dashed line in the blue-shifted region of the emission line to indicate that in this region, the line visibility is highly uncertain due to the presence of the P Cygni absorption. Bottom: continuum-corrected AMBER visibilities (filled blue squares) in the red region of the He I emission line as a function of spatial frequency and continuum-corrected visibility curve (solid red line) according to the NLTE model of Hillier et al. (2001). To derive the visibilities according to Eq. (2), the data in the wavelength range 2.057–2.058 µm were averaged before the continuum correction. The model curve predicts much larger visibility values in the line, which indicates that the size of the line-emitting region in the model is probably underestimated by the model. A rescaling of the size of the model by a factor of ~2.3 (green line) reveals better agreement with the AMBER measurements.</figcaption></figure>

the model size by a factor of 2.4 results in a much better agreement between the model and observations (green curve in Fig. 8, bottom) and a *FWHM* diameter of ≥8.2 mas, which is 3.6 times larger than the *FWHM* diameter of 2.3 mas in the continuum. Due to the lack of interferometric data at small projected baselines, this value can only give a rough lower limit of the size.

The results for the visibility inside the He I line can possibly be explained in a qualitative way in the framework of the binary model for the central object in η Car (e.g. Davidson et al. 1999; Davidson 2001; Pittard & Corcoran 2002; Hillier et al. 2006; Nielsen et al. 2006). In a model of this type, He I emission should arise near the wind-wind interaction zone between the

binary components. The hot secondary star is expected to ionize helium in a zone in the dense primary wind, adjoining the wind-wind interaction region. Such a region can produce He I recombination emission<fnref n="3" />. The wind-wind shocked gas, by contrast, is too hot for this purpose, while the density of the fast secondary wind is too low. Since the AMBER measurements (Dec. 2004 and Feb. 2005, at orbital phases φ = 0.268 and φ = 0.299, see Table 1) were obtained at an intermediate phase between periastron in July 2003 and apoastron in April 2006, the extension of the He I emission zone is expected to be rather diffuse and larger than the continuum size. In other words, the He I emission zone should be fairly extended and larger than the Hillier et al. model prediction, which is in agreement with the AMBER data.

### 3.5. Differential phases and closure phases

The measurement of phase information is essential for the reconstruction of images from interferometric data, but such an image reconstruction is only possible with an appropriate coverage of the uv plane. Nevertheless, even single phase measurements, in particular of the closure phase and differential phase, provide important information.

The closure phase (CP) is an excellent measure for asymmetries in the object brightness distribution. In our AMBER measurements, as illustrated in Figs. 2 and 4, we find that the CP in the continuum is zero within the errors for all the various projected baselines of the UT2-UT3-UT4 baseline triplet, indicating a point-symmetric continuum object. However, in the line emission, we detect a non-zero CP signal in all data sets. In both MR measurements covering the Bry line, we find the strongest CP signal in the blue wing of the emission line at λ = 2.164 µm (−34° and −20°) and a slightly weaker CP signal in the red wing of the emission line at λ = 2.167 µm (+12° and +18°). We also detected non-zero CP signals in the HR measurement around Bry taken at a different epoch. In the case of the He I line, a non-zero CP could only be detected at λ = 2.055 µm, just in the middle between the emission and absorption part of the P Cygni line profile.

The differential phase (DP) at a certain wavelength bin is measured relative to the phase at all wavelength bins. Therefore, the DP measured within a wavelength bin containing line emission yields approximately the Fourier phase of the combined object (continuum plus line emission) measured relative to the continuum. This Fourier phase might contain contributions from both the object phase of the combined object and a shift phase, which corresponds to the shift of the photocenter of the combined object relative to the photocenter of the continuum object. Significant non-zero DPs were detected in the Doppler-broadened line wings of the Bry line. Particularly within the blue-shifted wings, we found a strong signal (up to ~−60°), whereas the signals are much weaker within the red-shifted line wings. These DPs might correspond to small photocenter shifts, possibly arising if the outer Bry wind region consists of many clumps which are distributed asymmetrically. The small differential phases of up to ~−15° for the different baselines of the blue-shifted light in the He I line can perhaps also be explained by the above-mentioned asymmetries or within the

<footnote n="3">³ For a qualitative sketch of the geometry, see “zone 4” in Fig. 8 of Martin et al. (2006), even though this figure was drawn to represent He++ in a different context. For reasonable densities, the predicted He⁺ zone has a quasi-paraboloidal morphology. In addition, some extremely dense cooled gas, labeled “zone 6” in the same figure, may also produce He I emission.</footnote>
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Left: illustration of the components of our geometric model for an optically thick, latitude-dependent wind (see text for details). For the weak aspherical wind component, we draw the lines of latitudes to illustrate the 3D-orientation of the ellipsoid. Right (a), b): the upper row shows the brightness distribution of the modeled aspherical wind component (item (3) in the text) for two representative wavelengths. The figures below show the total brightness distribution after adding the contributions from the two spherical constituents of our model.</figcaption></figure>

framework of the binary model discussed in previous sections. In the binary model, a large fraction of the He I is possibly emitted from the wind-wind collision zone, which is located between the primary and the secondary (Davidson et al. 1999; Davidson 2001; Pittard & Corcoran 2002; Hillier et al. 2006; Nielsen et al. 2006).

### 3.6. Modeling with an inclined aspherical wind geometry

The goal of the modeling presented in this section is to find a model which is able to explain several remarkable features in our data; in particular, (a) the asymmetry in the Bry line profile (showing less emission in the blue-shifted wing than in the red-shifted wing) and the P Cygni-like absorption dip in the blue-shifted Bry wing, (b) the strong DP in the blue-shifted wing and a weaker DP signal in the red-shifted wing, and (c) the structure of the CP, showing a change in the sign between the blue- and red-shifted line wing. We aimed for a geometrical but physically motivated model which would reproduce these features at all wavelength channels simultaneously. For this, we concentrate on the Bry line, as this line shows a stronger phase signal than the He I line and was measured with a better uv coverage.

As Smith et al. (2003) convincingly showed, the stellar wind from η Car seems to be strongly latitude-dependent, with the highest mass flux and velocities at the poles. This anisotropy can be understood in the context of theoretical models (see, e.g., Maeder & Desjacques 2001), which take the higher temperatures at the poles ($g_{eff}$-effect) and the equatorial gravity darkening on a rapidly rotating star into account (von Zeipel effect, Zeipel 1925). As these models are quite successful in explaining the bipolar structure of the Homunculus nebula, we investigated

whether such bipolar geometries with a latitude-dependent velocity distribution might also be suited to explain our interferometric data.

Due to its success in reproducing both the spectrum and the measured visibilities, we based our wind model on the spherical Hillier et al. (2001) model and superposed a weak aspherical stellar wind geometry, which is inclined with respect to the line-of-sight. Our model includes three components (see Fig. 9); namely,

(1) a continuum component (using the Hillier et al. continuum CLV, see Fig. 6 top) with a blue-shifted absorption component;

(2) a spherical stellar wind (using the Hillier et al. continuum-subtracted Bry CLV); and

(3) an aspherical wind geometry, represented by a 41° inclined ellipsoid.

The relative contribution of these different constituents to the total flux is given by the input spectra shown in the upper row of Fig. 10. For the spherical and aspherical wind component, we assume Gaussian-shaped spectra. The original Hillier et al. CLVs slightly underestimate the size of the observed structures (see Fig. 5). Therefore, we rescaled them by 10% to obtain a better agreement for the visibilities at continuum wavelengths.

The aspherical wind of η Car is simulated as an ellipsoid with an inclination similar to the inclination angle of the Homunculus (41°, Smith 2006). While the south-eastern pole (which is inclined towards the observer) is in sight, the north-western pole is obscured. The latitude-dependent velocity distribution expected for the η Car wind was included in our model by coupling the latitude-dependent brightness distribution of the ellipsoid to
the wind velocity. At the highest *blue-shifted* velocities, mainly the *south-eastern polar region* contributes to the emission (see Fig. 9a). In the *red-shifted* line wing, mainly the (obscured) *north-western pole* radiates (see Fig. 9b). The axis of the ellipsoid was assumed to be oriented along the Homunculus polar axis (PA 132°, Smith 2006) and its axis ratio was fixed to 1.5.

As our simulations show, such an asymmetric geometry can already explain the measured DPs and CPs with a rather small contribution of the asymmetric structure to the total flux (see black line in Fig. 10, upper row). Although the large number of free parameters prevented us from scanning the whole parameter space, we found reasonable agreement with a size of the ellipsoid major axis of 8 mas. Figure 10 shows the spectrum, visibilities, DPs, and CPs computed from the model.

As our model was inspired by physical models, but does not take the complicated radiation transport and hydrodynamics involved in reality into account, we would like to note that our model allows us to check for consistency between the considered geometry and the AMBER spectro-interferometric data, but can neither constrain the precise parameters of a possible aspherical latitude-dependent stellar wind around η Car, nor can it rule out other geometries. We summarize some qualitative properties of our wind model as follows:

- The strong underlying spherical component mainly accounts for the very low visibilities measured within the line.

- The aspherical wind component introduces the asymmetry required to roughly explain the measured phase signals. In particular, it reproduces the larger DPs and CPs within the blue-shifted line wing compared to the red-shifted wing, as the red-shifted emission region is considerably obscured. It also accounts for the flip in the CP sign, as the photocenter of the line emission shifts its location between the blue- and red-shifted wing relative to the continuum photocenter.

- The absorption component which we introduced in the blue-shifted wing of the Bry line allows us to reproduce the asymmetry measured in the shape of the Bry emission line profile (showing an increase of flux towards red-shifted wavelengths) and the weak dip observed at far-blue-shifted wavelengths. Furthermore, with the decrease of the continuum contribution, the absorption component helps to lower the visibilities in the blue-shifted line wing, simultaneously increasing the asymmetry in the brightness distribution (increasing the phase signals). Finally, with the interplay between the absorption and emission component, our simulation reproduced a “bump” in the visibility similar to the one observed on the shortest baseline of our HR measurement ($\lambda \approx 2.163$ µm).

## 3.7. Feasibility of the detection of the hypothetical hot companion and the wind-wind interaction zone

One of the most intriguing questions regarding η Car is whether or not its central object is a binary, as suggested to explain cycle (e.g., Damineli 1996).

### 3.7.1. A simple binary continuum model

To investigate whether the AMBER measurements presented here can shed more light on the binarity hypothesis, we used the following approach: we constructed a *simple binary model* consisting of a primary wind component with a CLV according to the continuum model of Hillier et al. (2006; FWHM diameter ~2.3 mas; see upper panels in Fig. 6) and an unresolved binary companion represented by a point-like source (uniform disk

with ~0.1 mas FWHM diameter). The secondary component is predicted to be approximately located at PA -36° with a separation of 8 mas from the primary for the time of the AMBER observations (Nielsen et al. 2006). The continuum flux ratio *q* was treated as a free parameter. We would like to note that in our model, we assumed that all *K*-band light from the secondary is reaching us unprocessed; i.e. we ignored a possible dilution or re-distribution of the secondary’s radiation.

We calculated the 2D visibility function (see Fig. 11c) of this model intensity distribution for different values of *q*, as well as the closure phases for the baselines and PAs corresponding to our AMBER measurements (Fig. 11d). Finally, we compared the results with those obtained from a single component model where only the primary wind is present (Fig. 11b). The differences of the visibilities and closure phases between the single star and the binary model (at the baselines and PAs corresponding to our AMBER measurements) are displayed in Fig. 11d as a function of the *K*-band flux ratio of the binary components.

Figure 11d shows two interesting results: First, the closure phase is more sensitive to the binary signature than the visibilities and, thus, a more suitable observable to constrain the binary hypothesis. And second, given the accuracies of our first AMBER visibility and closure phase measurements (indicated by the horizontal dashed-dotted lines), we can conclude, for the particular model shown in Fig. 11, that the AMBER closure phases put an intensity ratio limit $q_{min} \approx 110$ on the binary *K*-band flux ratio. This limit is in line with the estimate $q \approx 200$ given by Hillier et al. (2006). Thus, based on the model shown in Fig. 11, the AMBER measurements are not in conflict with recent model predictions for the binary.

To investigate whether we can put similar constraints on the minimum *K*-band flux ratio $q_{min}$ for arbitrary separations and PAs, we calculated a larger grid of binary models and compared the residuals of visibilities and closure phases analogue to the example shown in Fig. 11.

For these grid calculations, we used values in the range from 4 to 14 mas for the binary separation with increments of 1 mas, and PAs of the secondary in the whole range from 0° to 360° in steps of 10°. The *K*-band flux ratio *q* of the binary components was varied in the range from 1 to 250 with $\Delta q = 4$. As a result of the grid calculation, we obtained the minimum *K*-band flux ratio $q_{min}$ as a function of binary separation and orientation.

Whereas the study with a fixed companion position presented above allowed us to put rather stringent constraints on $q_{min}$ (see Fig. 11), this systematic study revealed that due to the rather poor uv coverage, a few very specific binary parameter sets exist where we are only sensitive to $q_{min} \approx 10$. Nevertheless, for the above-mentioned separation interval (4 to 14 mas), we found that we are able to detect companions up to $q_{min} = 50$ at more than 90% of all PAs. In order to push this sensitivity limit in future observations, a better uv coverage will be required. Together with the expected higher closure phase accuracy, AMBER will be sensitive up to $q > 200$ and, therefore, have the potential to probe the currently favored binary models.

### 3.7.2. Can AMBER detect a He I wind-wind interaction zone shifted a few mas from the primary wind?

In the context of the binary hypothesis, it is also important to discuss the implications for the interpretation of the AMBER He I measurements. According to the binary model, a large fraction of the He I line emission should arise from the wind-wind collision zone expected between the primary and the secondary (Davidson et al. 1999; Davidson 2001; Hillier et al. 2006;
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Observables computed from our optically thick, latitude-dependent wind model (see Fig. 9 for a model illustration). The points (crosses) represent the measurements (as also shown in Fig. 2), and the solid lines give the observables computed from our model. The upper row shows the contributions from the various model components to the total flux. Besides the continuum emission (purple line), we introduced a spherical (blue line) and an aspherical (black line) wind component.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Simulation illustrating the signatures of a binary companion at the predicted position (for the orbital phase at the time of our continuum observations around the Bry line; separation ~8 mas, PA ~−36°; see Nielsen et al. 2006). For the primary component, we assume a continuum emission CLV from Hillier et al. (2006), to which we add a uniform disk of 0.1 mas diameter to account for the binary companion (see a). b) Shows the 2D-visibility function of a single component model containing only the continuum-emitting object with Hillier-type CLV. c) Shows the 2D-visibility corresponding to the binary model shown in panel a) for a small K-band flux ratio $q = 10$ between the continuum emitting region (Hillier-type CLV) and the companion. In d) the residuals of the visibilities (upper panel) and closure phases (lower panel) between the models with and without a companion are shown as a function of flux ratio $q$ for baselines and PAs corresponding to the seven AMBER measurements. The dashed horizontal lines indicate the uncertainty of the measurements. Given these uncertainties, the figure illustrates that a companion signature should be visible from the AMBER closure phase measurements, if the binary is less than ~110 times fainter compared to the primary. See text for details.</figcaption></figure>

Nielsen et al. 2006). The exact intensity ratio of primary He I wind and He I emission from the wind-wind interaction zone is not known. Figure 5 suggests that during the AMBER observations, the total He I flux was roughly two times larger than the model prediction of Hillier et al. (2001) for the primary He I wind.

At the orbital phases of the AMBER measurements, the wind-wind collision zone should be at resolvable distances from η Car's primary (resolution ~5 mas; companion separation ~8 mas, PA ~−36°, Nielsen et al. 2006). Looking at the AMBER He I data, we see that the differential as well as the closure phases are zero everywhere except for the transition region between the absorption and emission part of the He I line, where

we find differential phases of ~10–20° and a closure phase of ~−30°; i.e., the phases measured across the He I line are significantly weaker compared to the Bry line. The question is now, why AMBER measured weaker phase signals within the He I line and if this result is in line with the predictions of the wind-wind collision model.

One possible explanation for the small measured phases could be the orientation of the binary orbit. If the orbit's major axis is nearly aligned with the line-of-sight, the photocenter shift inside the He I line will be very small. In addition, the deviations from point symmetry would be rather small. Therefore, in the case of this special geometry, both differential phases and the closure phase would be small, in qualitative agreement with
the AMBER data. Another explanation could be that the contri-
bution of the wind-wind collision zone to the He I line emission
is much weaker than that of the primary wind. However, this is
not very likely (see Hillier et al. 2006).

A different explanation for the weak phases can be found
from a modeling approach similar to the one for the Bry line
region outlined in Sect. 3.6. Based on the results presented in
Sects. 3.1 and 3.4.2, we constructed a simple He I model consist-
ing of a spherical primary wind component with a Hillier-type
CLV (2.5 mas *FWHM* diameter) and an extended spherical He I
line-emitting region with Gaussian CLV and a 7 mas *FWHM* di-
ameter (i.e., for simplicity, we assumed that all He I flux is emit-
ted from the wind-wind interaction region; however, some frac-
tion of He I is also emitted from the primary wind; see Hillier
et al. 2006; and Fig. 5 of the present paper). The center of the
line-emitting component of this model is located 3 mas away
from the primary wind component towards PA 132°; i.e., in the
direction of the Homunculus axis. The spectra of the continuum
and line-emitting components were chosen in such a way that the
combined spectrum resembles the observed He I line spectrum.

The modeling results show that this simple model is approx-
imately able to simultaneously reproduce the observed spec-
trum and the wavelength dependence of visibilities, differential
phases (10–20°), and closure phases (~30°). Thus, our simple
model example illustrates that the AMBER measurements can
be understood in the context of a binary model for η Car and
the predicted He I wind-wind collision scenario (e.g. Davidson
et al. 1999; Davidson 2001; Hillier et al. 2006; Nielsen et al.
2006). We note that the model parameter values given above
are of preliminary nature. A more detailed, quantitative model-
ing is in preparation and will be subject of a forthcoming paper.
Furthermore, we would like to emphasize, as already discussed
in previous sections, that there are likely to be *three* sources of
He I emission – the primary wind, a wind-wind interaction zone
(bow shock), and the ionized wind zone caused by the ionization
of the secondary. For both the bow shock and the ionized wind
zone, the ionizing UV radiation field of the secondary is of cru-
cial importance. On the basis of the observed blue-shift and the
weakness of the He I during the event, we believe that the pri-
mary wind contribution is small. It is not yet possible to decide
on the relative contributions of the bow shock and the ionized
wind region.

4. Conclusions

In this paper we present the first near-infrared spectro-
interferometry of the enigmatic Luminous Blue Variable η Car
obtained with AMBER, the 3-telescope beam combiner of
ESO's VLTI. In total, three measurements with spectral reso-
lutions of R = 1500 and R = 12 000 were carried out in Dec.
2004 (φ = 0.268) and Feb. 2005 (φ = 0.299), covering two
spectral windows around the He I and Bry emission lines at
λ = 2.059 and 2.166 µm, respectively. From the measurements,
we obtained spectra, visibilities, differential visibilities, differen-
tial phases, and closure phases. From the analysis of the data, we
derived the following conclusions:

- In the *K*-band continuum, we resolved η Car's optically thick wind. From a Gaussian fit of the *K*-band continuum visibilities in the projected baseline range from 28–89 m, we obtained a *FWHM* diameter of 4.0±0.2 mas. Taking the different fields-of-view into account, we found good agreement between the AMBER measurements and previous VINCI/VLTI observations of η Car presented by van Boekel et al. (2003).

- When comparing the AMBER continuum visibilities with the NLTE radiative transfer model from Hillier et al. (2001), we find very good agreement between the model and observations. The best fit was obtained with a slightly rescaled version of the original Hillier et al. model (rescaling by 1–2%), corresponding to FWHM diameters of 2.27 mas at λ = 2.040 µm and 2.33 mas at λ = 2.174 µm.

- If we fit Hillier et al. (2001) model visibilities to the observed AMBER visibilities, we obtain, for example, 50% encircled-energy diameters of 4.2, 6.5, and 9.6 mas in the 2.17 µm continuum, the He I, and the Bry emission lines, respectively.

- In the continuum around the Bry line, we found an asymme-
  try towards position angle PA = 120° ± 15° with a projected
  axis ratio of ξ = 1.18 ± 0.10. This result confirms the earlier
  finding of van Boekel et al. (2003) using VINCI/VLTI and
  supports theoretical studies which predict an enhanced mass
  loss in polar direction for massive stars rotating close to their
  critical rotation rate (e.g. Owocki et al. 1996, 1998).

- For both the Bry and the He I emission lines, we measured non-zero differential phases and non-zero closure phases within the emission lines, indicating a complex, asymmetric object structure.

- We presented a physically motivated model which shows that the asymmetries (DPs and CPs) measured within the wings of the *Bry* line are consistent with the geometry expected for an aspherical, latitude-dependent stellar wind. Additional AMBER/VLTI measurements and radiative transfer modeling will be required to determine the precise parameters of such an inclined aspherical wind.

- Using a simple binary model, we finally looked for a possible binary signature in the AMBER closures phases. For separations in the range from 4 to 14 mas and arbitrary PAs, our simple model reveals a minimum *K*-band flux ratio of ~50 with a 90% likelihood.

Our observations demonstrate the potential of AMBER/VLTI
observations to unveil new structures of η Car on the scales
of milliarcseconds. Repeated observations will allow us to
trace changes in observed morphology over η Car's spectro-
scopic 5.5 yr period, possibly revealing the motion of the wind-
wind collision zone as predicted by the η Car binary model.
Furthermore, future AMBER observations with higher accuracy
might be sensitive enough to directly detect the hypothetical hot
companion.

*Acknowledgements.* We thank the referee Dr. N. Smith and Dr. A. Damineli for very valuable comments and suggestions which helped to considerably improve the manuscript. The NSO/Kitt Peak FTS data used here to identify the telluric lines in the AMBER data were produced by NSF/NOAO.

The AMBER project<fnref n="4" /> was founded by the French Centre National de la Recherche Scientifique (CNRS), the Max Planck Institute für Radioastronomie (MPIfR) in Bonn, the Osservatorio Astrofisico di Arcetri (OAA) in Firenze, the French Region "Provence Alpes Côte d'Azur" and the European Southern Observatory (ESO). The CNRS funding has been made through the Institut National des Sciences de l'Univers (INSU) and its Programmes Nationaux (ASHRA, PNPS, PNP).

The OAA co-authors acknowledge partial support from MIUR grants to the Arcetri Observatory: *A LBT interferometric arm, and analysis of VLTI interferometric data* and *From Stars to Planets: accretion, disk evolution and planet formation* and from INAF grants to the Arcetri Observatory *Stellar and Extragalactic Astrophysics with Optical Interferometry*. C. Gil work was supported in part by the Fundação para a Ciência e a Tecnologia through project POCTI/CTE-AST/55691/2004 from POCTI, with funds from the European program FEDER.

The preparation and interpretation of AMBER observations benefit from the tools developed by the Jean-Marie Mariotti Center for optical interferometry

<footnote n="4">⁴ The structure and members of the AMBER Consortium can be found in the website: http://amber.obs.ujf-grenoble.fr</footnote>
<figure><img src="image_13.png" /></figure>

**Fig. A.1.** *Top:* spectral calibration of the AMBER observations. The figure shows the spectrum of η Car (blue) and the calibrator star L Car (green) around the Bry line obtained from the HR measurements in Feb. 2005. In addition, a telluric spectrum from the Kitt Peak Observatory is shown in red. The original telluric spectrum with a spectral resolution of *R* = 40 000 was spectrally convolved to match the resolution of the AMBER measurements. For the flux calibration of the η Car spectrum, the Bry line in the L Car spectrum was interpolated before the division. See text for further details. *Bottom:* same as top panel, but for both MR measurements around the He I emission line. Together with the η Car spectra, the spectra of the calibrators HD 93030 and HD 89682 are displayed. In addition, the telluric spectrum obtained at Kitt Peak with *R* = 40 000 and a spectrally convolved telluric spectrum is shown, which matches the spectral resolution of the MR measurements (*R* = 1500). The figure illustrates that the forest of telluric lines forms a quasi-continuum which modulates the AMBER spectra.

JMMC<fnref n="5" /> and from the databases of the Centre de Données Stellaires (CDS) and of the Smithsonian/NASA Astrophysics Data System (ADS).
The data reduction software amdlib is freely available on the AMBER site http://amber.obs.ujf-grenoble.fr. It has been linked to the public do-main software Yorick<fnref n="6" /> to provide the user-friendly interface ammyorick.

**Appendix A: Wavelength calibration**

To obtain both an accurate wavelength calibration of the AMBER raw data and properly calibrated spectra of η Car, we compared the AMBER raw spectra of η Car as well as the calibrator stars L Car, HD 93030, and HD 89682 with a K-band telluric spectrum recorded at the Kitt Peak Observatory with

a spectral resolution of 40 000. For the comparison with the AMBER spectrum, this telluric spectrum was spectrally convolved to match the spectral resolution of the AMBER measurements with high (*R* ~ 12 000) and medium (*R* ~ 1500) spectral resolution.

The result of the comparison is shown in Fig. A.1. In the upper panel, the high spectral resolution AMBER spectra of η Car and the calibrator L Car are shown together with the telluric spectrum with *R* = 10 000. From the comparison with the telluric spectrum, we identified 7 prominent telluric absorption features in the L Car spectrum, which are indicated by the dashed vertical lines. The strongest absorption line seen in the L Car spectrum is not telluric, but can be identified as intrinsic Bry absorption in L Car. Therefore, to properly calibrate the η Car spectrum with the L Car spectrum, we had to interpolate the Bry line region in the L Car spectrum before dividing the two spectra. From the spectral calibration shown in Fig. A.1, we estimated a wavelength calibration error of the AMBER data Δλ = 3 × 10⁻⁴ μm

The lower panel in Fig. A.1 shows the wavelength calibration of the medium spectral resolution data in the wavelength region around the He I line. The figure contains the two η Car MR spectra and the spectra of the two corresponding calibrator stars, HD 93030 and HD 89682, as well as the telluric spectra with spectral resolutions of *R* = 40 000 and *R* = 1500. As the telluric spectra reveal, there is a forest of telluric lines in the spectral region around the He I line. As can be seen in Fig. 4, the modulation of the continuum flux introduced by the telluric quasi-continuum cancels out completely when the η Car spectra are divided by the corresponding calibrator spectra, which show no prominent intrinsic line features. Since there are no sharp spectral features in the 2.03–2.08 μm region of either the calibrator or telluric spectra which could be used for the spectral calibration, we estimated a wavelength calibration error Δλ = 6 × 10⁻⁴ μm for the MR He I data. On the other hand, for the MR data around the Bry line, we found Δλ = 4 × 10⁻⁴ μm.

**Appendix B: Continuum uniform disk and Gauss diameter fits**

For each spectral channel as well as for an averaged continuum, we performed 1-D fits to the visibility data using simple uniform disk (UD) and Gaussian models. In this step of the analysis, possible asymmetries were ignored and all visibility points at a given wavelength were fitted together, regardless of the position angle of the observations. The results of these 1-D fits are illustrated in the two upper panels of Fig. B.1 for the averaged continuum data in the wavelength ranges 2.03–2.08 μm and 2.155–2.175 μm, respectively. As the figure reveals, neither a uniform disk nor a single Gaussian provides a good fit to the continuum data. At least, this is true as long as no contamination by a fully resolved background component is taken into account.

The wavelength dependence of the apparent size obtained from the UD and GAUSS model fits for the individual spectral channels is shown in the lower left panel of Fig. B.1. This panel illustrates that the equivalent UD and GAUSS K-band diameters of η Car derived from the AMBER data are ~4 and ~6.5 mas, respectively.

It should be added here that a good fit of the AMBER data using, for instance, a Gaussian can indeed be obtained when a certain amount of contamination due to a fully resolved background component is taken into account (see also Petrov et al. 2005). To illustrate that, we performed Gaussian fits to the AMBER data, where we introduced such a fully resolved component as a free fitting parameter. We found that the

<footnote n="5">⁵ The JMMC is a center providing software tools for optical interferometry described at the website: http://www.jmmc.fr</footnote>

<footnote n="6">⁶ http://yorick.sourceforge.net</footnote>
<figure><img src="image_14.png" /><figcaption>**Fig. B.1.** 1-D visibility fits of the AMBER continuum measurements. The upper two panels show uniform-disk and Gaussian fits for the averaged continuum around the He I and Brγ lines, respectively. As both panels illustrate, neither a uniform disk nor a Gaussian is a good representation of the AMBER measurements. The panel on the left summarizes the fitting attempt for all continuum wavelength channels by displaying the fitted diameters of the UD and GAUSS models as a function of wavelength.</figcaption></figure>

best Gaussian fit is obtained with a FWHM diameter $d_{\text{Gauss}}$ = 3.01 mas and a ~29% background contamination for the He I continuum region and $d_{\text{Gauss}}$ = 3.32 mas and a ~30% background contamination for the Brγ continuum region. Thus, from this fit we would derive a background contamination which is only ~50% smaller than in the VINCI data. We think that such a large amount of background contamination is not very likely given the small AMBER fiber aperture (60 mas) of the 8.2 m telescopes. We think that the large amount of background contamination needed to find a reasonable Gauss fit just reflects the fact that a Gaussian is not appropriate to describe the observations. This is confirmed by the fact that for the fit of the radiative transfer model of Hillier et al. (2001), no background component has to be taken into account to reproduce the AMBER measurements.

We would like to note here that the Gaussian FWHM diameters of typically ~4 mas found from the AMBER measurements are not in contrast to the value $d_{\text{Gauss}}$ ~ 7 mas found by van Boekel et al. (2003) from VINCI/VLTI observations for the following reason: Since a Gaussian is not a good representation of both the VINCI and the AMBER visibilities, the diameter resulting from a Gaussian fit strongly depends on the fit range. This is illustrated in Fig. B.2 for the four VINCI measurements given in Fig. 1 of van Boekel et al. (2003). As the figure shows, from a Gaussian fit of all four data points, $d_{\text{Gauss}}$ = 6.0 mas is obtained. If only the data point with $q$ = 45 cycles/arcsec is fitted (corresponding to a projected baseline length of ~24 m), we get $d_{\text{Gauss}}$ = 7.65 mas. This is in agreement with the values given in van Boekel et al. (2003) for the elliptical Gaussian fit of the large number of VINCI measurements with a projected baseline of 24 m (see their Fig. 2). On the other hand, if we fit only the VINCI data point corresponding to the longest projected baseline ($q$ = 136 cycles/arcsec), a Gaussian fit provides $d_{\text{Gauss}}$ = 4.39 mas (see Fig. B.2), which is very close

<figure><img src="image_15.png" /><figcaption>**Fig. B.2.** Dependence of the Gaussian FWHM diameter on the fit range. The figure shows the background-corrected visibilities obtained with VINCI/VLTI (see van Boekel et al. 2003) as well as Gaussian fits of (a) all four data points (long-dashed green line), (b) only the point with $q$ = 45 cycles/arcsec (short-dashed blue), (c) only the point corresponding to the longest baseline (dotted purple), and (d) only the point corresponding to the shortest baseline (dashed-dotted light blue). See the labels for the Gaussian FWHM diameters resulting from the different fits. The figure illustrates that the fitted diameter strongly depends on the spatial frequency range which is used to fit the data. The strong diameter variation (in this case, the diameter changes by a factor of ~3) occurs since a Gaussian is not a good representation of the measured visibility function. See text for further discussion.</figcaption></figure>

to the diameter we obtain from the AMBER measurements for $\lambda$ ~ 2.174 µm ($d_{\text{Gauss}}$ = 4.35 mas). This is not surprising since the spatial frequency of this VINCI data point agrees with the average spatial frequency of our AMBER observations ($q$ ~ 50–200 cycles/arcsec). Thus, it can be concluded that good agreement between the Gaussian FWHM diameters derived from the AMBER and VINCI measurements is found if a comparable spatial frequency range is used for the fit.
## Appendix C: Visibility and differential phase of an emission line object

We assume that the target's intensity distribution can be described by two components: the continuum spectrum $o_{\text{cont}}(x, y; \lambda)$ and the emission line spectrum $o_{\text{line}}(x, y; \lambda)$. In the part of the spectrum containing the emission line, both $o_{\text{cont}}$ and $o_{\text{line}}$ contribute to the total intensity distribution $o_{\text{tot}}$. According to the van-Zittert-Zernike theorem, the Fourier transforms $O_{\text{cont}}(B/\lambda)$ and $O_{\text{line}}(B/\lambda)$ of $o_{\text{cont}}(x, y; \lambda)$ and $o_{\text{line}}(x, y; \lambda)$ are measured with an optical long baseline interferometer at wavelength $\lambda$ and projected baseline vector $B$. In the following, we assume that all Fourier spectra are normalized to 1 at frequency zero. The complex Fourier spectrum $O_{\text{tot}}(B/\lambda)$ of the intensity distribution $o_{\text{tot}}$ measured at the emission line $\lambda_{\text{line}}$ is given by

$$O_{\text{tot}} = \frac{1}{F_{\text{cont}} + F_{\text{line}}} (F_{\text{cont}} \cdot O_{\text{cont}} + F_{\text{line}} \cdot O_{\text{line}}), \quad (C.1)$$

where $F_{\text{cont}}$ and $F_{\text{line}}$ are the fluxes of the continuum component $o_{\text{cont}}$ and the line component $o_{\text{line}}$, respectively. In the emission line, the total flux measured is $F_{\text{tot}} := F_{\text{cont}} + F_{\text{line}}$.

From the spectrally dispersed interferometric data, we can derive the differential phase, which is the difference of the Fourier phases of the continuum component $o_{\text{cont}}$ and the total intensity $o_{\text{tot}}$ in the emission line. The differential phase $\Phi'(B/\lambda_{\text{line}})$ in the emission line at $\lambda_{\text{line}}$ is given by

$$O_{\text{cont}} \cdot O_{\text{tot}}^{*} = V_{\text{cont}} \cdot V_{\text{tot}} \cdot e^{i\Phi'}, \quad (C.2)$$

where $\Phi'(B/\lambda_{\text{line}}) := \Phi_{\text{cont}}(B/\lambda_{\text{line}}) - \Phi_{\text{tot}}(B/\lambda_{\text{line}})$. $\Phi_{\text{cont}}(B/\lambda_{\text{line}})$ is the Fourier phase of the continuum component, and $\Phi_{\text{tot}}(B/\lambda_{\text{line}})$ denotes the Fourier phase of $o_{\text{tot}}$ at the position of the emission line $\lambda_{\text{line}}$. The asterisk * in this equation denotes conjugate complex operation. $V_{\text{cont}}(B/\lambda_{\text{line}})$ describes the visibility of the continuum component at the position of the emission line $\lambda_{\text{line}}$, and $V_{\text{tot}}(B/\lambda_{\text{line}})$ is the visibility measured at the position of the emission line $\lambda_{\text{line}}$. Inserting Eq. (C.1) into Eq. (C.2) yields

$$V_{\text{cont}} \cdot V_{\text{tot}} \cdot e^{i\Phi'} = \frac{V_{\text{cont}}}{F_{\text{cont}} + F_{\text{line}}} \cdot (F_{\text{cont}} V_{\text{cont}} + F_{\text{line}} V_{\text{line}} e^{i\Delta\Phi}), \quad (C.3)$$

where $\Delta\Phi(B/\lambda_{\text{line}})$ denotes the difference of the Fourier phases of the continuum and line components; i.e., $\Delta\Phi(B/\lambda_{\text{line}}) := \Phi_{\text{cont}}(B/\lambda_{\text{line}}) - \Phi_{\text{line}}(B/\lambda_{\text{line}})$. $\Phi_{\text{cont}}(B/\lambda_{\text{line}})$ and $\Phi_{\text{line}}(B/\lambda_{\text{line}})$ are the Fourier phases of the continuum and line components, respectively.

In the vector representation of complex numbers, the three quantities $F_{\text{cont}}V_{\text{cont}}$, $F_{\text{line}}V_{\text{line}}$, and $F_{\text{tot}}V_{\text{tot}}$ form a triangle with one corner placed at the center of the coordinate system. According to the law of cosines, the correlated flux of the line component is given by Eq. (2) (see Sect. 3.4.1):

$$|F_{\text{line}} V_{\text{line}}|^2 = |F_{\text{tot}} V_{\text{tot}}|^2 + |F_{\text{cont}} V_{\text{cont}}|^2 - 2 \cdot F_{\text{tot}} V_{\text{tot}} \cdot F_{\text{cont}} V_{\text{cont}} \cdot \cos(\Phi'). \quad (C.4)$$

Since the flux $F_{\text{line}}$ can be calculated from the measured fluxes $F_{\text{cont}}$ and $F_{\text{tot}}$, the visibility $V_{\text{line}}$ of the line component can be derived using Eq. (C.4).

Applying the law of sines to this triangle in the complex plane yields the differential phase $\Delta\Phi(B/\lambda_{\text{line}})$, which is the

difference between the Fourier phase $\Phi_{\text{cont}}(B/\lambda_{\text{line}})$ of the continuum component and the Fourier phase $\Phi_{\text{line}}(B/\lambda_{\text{line}})$ of the line component:

$$\sin(\Delta\Phi) = \sin(\Phi') \cdot \frac{|F_{\text{tot}}V_{\text{tot}}|}{|F_{\text{line}}V_{\text{line}}|}, \quad (C.5)$$

where $\Phi'(B/\lambda_{\text{line}})$ is the differential phase measured at the position of the emission line $\lambda_{\text{line}}$.

## References

* Chesneau, O., Min, M., Herbst, T., et al. 2005, A&A, 435, 1043

* Corcoran, M. F. 2005, AJ, 129, 2018

* Damineli, A. 1996, ApJ, 460, L49

* Damineli, A., Conti, P. S., & Lopes, D. F. 1997, New Astron., 2, 107

* Damineli, A., Kaufer, A., Wolf, B., et al. 2000, ApJ, 528, L101

* Davidson, K. 1999, in ASP Conf. Ser., 179, Eta Carinae at The Millennium, ed. J. A. Morse, R. M. Humphreys, & A. Damineli, 304

* Davidson, K. 2001, in Eta Carinae and Other Mysterious Stars: The Hidden Opportunities of Emission Spectroscopy, ed. T. R. Gull, S. Johannson, & K. Davidson, ASP Conf. Ser., 242, 3

* Davidson, K., Dufour, R. J., Walborn, N. R., & Gull, T. R. 1986, ApJ, 305, 867

* Davidson, K., Ebbets, D., Weigelt, G., et al. 1995, AJ, 109, 1784

* Davidson, K., & Humphreys, R. M. 1997, ARA&A, 35, 1

* Davidson, K., Ishibashi, K., Gull, T. R., & Humphreys, R. M. 1999, in Eta Carinae at The Millennium, ed. J. A. Morse, R. M. Humphreys, & A. Damineli, ASP Conf. Ser., 179, 227

* Davidson, K., Ishibashi, K., Gull, T. R., Humphreys, R. M., & Smith, N. 2000, ApJ, 530, L107

* Davidson, K., Martin, J., Humphreys, R. M., et al. 2005, AJ, 129, 900

* Davidson, K., Smith, N., Gull, T. R., Ishibashi, K., & Hillier, D. J. 2001, AJ, 121, 1569

* Duncan, R. A., White, S. M., Reynolds, J. E., & Lim, J. 1999, in ASP Conf. Ser., 179, Eta Carinae at The Millennium, ed. J. A. Morse, R. M. Humphreys, & A. Damineli, 54

* Dyck, H. M., Benson, J. A., van Belle, G. T., & Ridgway, S. T. 1996, AJ, 111, 1705

* Feast, M., Whitelock, P., & Marang, F. 2001, MNRAS, 322, 741

* Gull, T. R., Ishibashi, K., & The Cycle 7 STIS Go Team. 1999, in Eta Carinae at The Millennium, ed. J. A. Morse, R. M. Humphreys, & A. Damineli, ASP Conf. Ser., 179, 144

* Hillier, D. J., & Miller, D. L. 1998, ApJ, 496, 407

* Hillier, D. J., Davidson, K., Ishibashi, K., & Gull, T. 2001, ApJ, 553, 837

* Hillier, D. J., Gull, T., Nielsen, K., et al. 2006, ApJ, 642, 1098

* Hofmann, K.-H., & Weigelt, G. 1988, A&A, 203, L21

* Ishibashi, K. 2001, in Eta Carinae and Other Mysterious Stars: The Hidden Opportunities of Emission Spectroscopy, ed. T. R. Gull, S. Johannson, & K. Davidson, ASP Conf. Ser., 242, 53

* Ishibashi, K., Corcoran, M. F., Davidson, K., et al. 1999, ApJ, 524, 983

* Kervella, P., Fouqué, P., Storm, J., et al. 2004a, ApJ, 604, L113

* Kervella, P., Nardetto, N., Bersier, D., Mourard, D., & Coudé du Foresto, V. 2004b, A&A, 416, 941

* Kervella, P., Mérand, A., Perrin, G., & Coudé de Foresto, V. 2006, A&A, 448, 623

* Maeder, A., & Desjacques, V. 2001, A&A, 372, L9

* Malbet, F., Benisty, M., De Wit, W.-J., et al. 2007, A&A, 464, 43

* Martin, J. C., Davidson, K., Humphreys, R. M., Hillier, D. J., & Ishibashi, K. 2006, ApJ, 640, 474

* Millour, F., Vannier, M., Petrov, R. G., et al. 2006, in EAS Publ. Ser., in press

* Nielsen, K. E., Corcoran, M. F., Gull, T. R., et al. 2006, ApJ submitted

* Owocki, S. P., Cranmer, S. R., & Gayley, K. G. 1996, ApJ, 472, L115

* Owocki, S. P., Cranmer, S. R., & Gayley, K. G. 1998, Ap&SS, 260, 149

* Petrov, R. G., Malbet, F., Weigelt, G., et al. 2003, in Interferometry for Optical Astronomy II, ed. W. A. Traub, Proc. SPIE, 4838, 924

* Petrov, R., Millour, F., Chesneau, O., et al. 2005, in The power of optical/IR interferometry: recent scientific results and 2nd generation VLTI instrumentation, ed. F. Delplancke, A. Richichi, [astro-ph/0509208]

* Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1

* Pittard, J. M., & Corcoran, M. F. 2002, A&A, 383, 636

* Richichi, A., Percheron, I., & Khristoforova, M. 2005, A&A, 431, 773

* Smith, N. 2002, MNRAS, 337, 1252

* Smith, N. 2004, MNRAS, 351, L15

* Smith, N. 2006, ApJ, 644, 1151

* Smith, N., Morse, J. A., Davidson, K., & Humphreys, R. M. 2000, AJ, 120, 920
Smith, N., Davidson, K., Gull, T. R., Ishibashi, K., & Hillier, D. J. 2003, ApJ, 586, 432

Steiner, J. E., & Damineli, A. 2004, ApJ, 612, L133

Tatulli, E., Millour, F., Chelli, A., et al. 2007, 464, 29

van Boekel, R., Kervella, P., Schöller, M., et al. 2003, A&A, 410, L37

van Genderen, A. M., Sterken, C., Allen, W. H., & Liller, W. 2003, A&A, 412, L25

Verner, E., Bruhweiler, F., & Gull, T. 2005, ApJ, 624, 973

Weigelt, G., & Ebersberger, J. 1986, A&A, 163, L5

Weigelt, G., Albrecht, R., Barbieri, C., et al. 1995, in Rev. Mex. Astron. Astrofis. Conf. Ser., ed. V. Niemela, N. Morrell, & A. Feinstein, 11

Weis, K., Stahl, O., Bomans, D. J., et al. 2005, AJ, 129, 1694

Whitelock, P. A., Feast, M. W., Marang, F., & Breedt, E. 2004, MNRAS, 352, 447

Zanella, R., Wolf, B., & Stahl, O. 1984, A&A, 137, 79

Zeipel, H. V. 1925, MNRAS, 85, 678

¹ Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany

² Laboratoire Universitaire d'Astrophysique de Nice, UMR 6525 Université de Nice - Sophia Antipolis/CNRS, Parc Valrose, 06108 Nice Cedex 2, France

³ Laboratoire d'Astrophysique de Grenoble, UMR 5571 Université Joseph Fourier/CNRS, BP 53, 38041 Grenoble Cedex 9, France

⁴ Laboratoire Gemini, UMR 6203 Observatoire de la Côte d'Azur/CNRS, BP 4229, 06304 Nice Cedex 4, France

⁵ INAF-Osservatorio Astrofisico di Arcetri, Istituto Nazionale di Astrofisica, Largo E. Fermi 5, 50125 Firenze, Italy

⁶ Department of Physics and Astronomy, University of Pittsburgh, 3941 O Hara Street, Pittsburgh, PA 15260, USA

⁷ Laboratory for Extraterrestrial Planets and Stellar Astrophysics, Goddard Space Flight Center, 20771 Greenbelt, Maryland, USA

⁸ School of Physics and Astronomy, University of Minnesota, 116 Church Street SE, Minneapolis, MN 55455, USA

⁹ European Southern Observatory, Casilla 19001, Santiago 19, Chile

¹⁰ ONERA/DOTA, 29 avenue de la Division Leclerc, BP 72, 92322 Chatillon Cedex, France

¹¹ Centre de Recherche Astronomique de Lyon, UMR 5574 Université Claude Bernard/CNRS, 9 avenue Charles André, 69561 Saint Genis Laval Cedex, France

¹² Division Technique INSU/CNRS UPS 855, 1 place Aristide Briand, 92195 Meudon Cedex, France

¹³ IRCOM, UMR 6615 Université de Limoges/CNRS, 123 avenue Albert Thomas, 87060 Limoges Cedex, France

¹⁴ European Southern Observatory, Karl Schwarzschild Strasse 2, 85748 Garching, Germany

¹⁵ Kiepenheuer Institut für Sonnenphysik, Schöneckstr. 6, 79104 Freiburg, Germany

¹⁶ Instituut voor Sterrenkunde, KU-Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium

¹⁷ Centro de Astrofísica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal

¹⁸ Laboratoire Astrophysique de Toulouse, UMR 5572 Université Paul Sabatier/CNRS, BP 826, 65008 Tarbes Cedex, France

¹⁹ Departamento de Astronomia, Universidad de Chile, Chile

²⁰ Astronomisches Institut, Ruhr-Universität Bochum, Universitätsstr. 150, 44780 Bochum, Germany; & Lise-Meitner fellowship