author: T. Driebe; J. H. Groh; K.-H. Hofmann; K. Ohnaka; S. Kraus; F. Millour; K. Murakawa; D. Schertl; G. Weigelt; R. Petrov; M. Wittkowski; C.A. Hummel; J.B. Le Bouquin; A. Merand; M. Schöller; F. Massi; P. Stee; E. Tatulli
bibliography: 11723.bib
date: Version: 2024-10-02
title: Resolving the asymmetric inner wind region of the yellow hypergiant with VLTI/AMBER in low and high spectral resolution mode[^1]

is a massive evolved star belonging to the group of yellow hypergiants. Currently, this star is rapidly evolving through the Hertzprung-Russell diagram, crossing the so-called yellow void. is suffering from intensive mass loss which led to the formation of an extended dust shell. Moreover, the dense stellar wind of is subject to strong line emission. Our goal was to probe the photosphere and the innermost circumstellar environment of , to measure the size of its continuum- as well as the Br$\gamma\,$line-emitting region on milliarcsecond scales, and to search for evidence of an asymmetric distribution of 's dense, circumstellar gas. We obtained near-infrared long-baseline interferometry of with the AMBER instrument of ESO's Very Large Telescope Interferometer (VLTI). The measurements were carried out in May/June 2007 and May 2008 in low-spectral resolution mode in the $JHK$ bands using three Auxillary Telescopes (ATs) at projected baselines ranging from 30 to 96 m, and in October 2008 in high-spectral resolution mode in the $K$ band around the Br$\gamma\,$emission line using three Unit Telescopes (UTs) with projected baselines between 54 and 129 m. The high-spectral resolution mode observations were analyzed by means of radiative transfer modeling using CMFGEN and the 2-D Busche & Hillier codes. For the first time, we have been able to absolutely calibrate the $H$- and $K$-band data and, thus, to determine the angular size of IRC+10420's continuum- and Br$\gamma$ line-emitting regions. We found that both the low resolution differential and closure phases are zero within the uncertainty limits across all three bands. In the high-spectral resolution observations, the visibilities show a noticeable drop across the Br$\gamma$ line on all three baselines. We found differential phases up to -25  in the redshifted part of the Br$\gamma$ line and a non-zero closure phase close to the line center. The calibrated visibilities were corrected for AMBER's limited field-of-view to appropriately account for the flux contribution of 's extended dust shell. From our low-spectral resolution AMBER data we derived FWHM Gaussian sizes of $1.05\pm 0.07$ and $0.98\pm 0.10$ mas for 's continuum-emitting region in the $H$ and $K$ bands, respectively. From the high-spectral resolution data, we obtained a FWHM Gaussian size of $1.014\pm0.010$ mas in the $K$-band continuum. The Br$\gamma\,$-emitting region can be fitted with a geometric ring model with a diameter of $4.18^{+0.19}_{-0.09}~$mas, which is approximately 4 times the stellar size. The geometric model also provides some evidence that the Br$\gamma\,$line-emitting region is elongated towards a position angle of 36, well aligned with the symmetry axis of the outer reflection nebula. Assuming an unclumped wind and a luminosity of $6\times10^5\ensuremath{\mathrm{L}_{\odot}}$, the spherical radiative transfer modeling with CMGFEN yields a current mass-loss rate of $1.5 - 2.0\times 10^{-5}\ensuremath{\mathrm{M}_{\odot}{\rm yr}^{-1}}$ based on the Br$\gamma\,$equivalent width. However, the spherical CMFGEN model poorly reproduces the observed line shape, blueshift, and extension, definitively showing that the outflow is asymmetric. Our 2-D radiative transfer modeling shows that the blueshifted Br$\gamma\,$emission and the shape of the visibility across the emission line can be explained with an asymmetric bipolar outflow with a high density contrast from pole to equator (8–16), where the redshifted light is substantially diminished.

# Introduction

Due to its distance , the relatively high wind velocity (40 ${\rm kms}^{-1}$), and the remarkable photometric history, (= = ) is most likely not a post-AGB star evolving through the proto-planetary nebula stage as suggested earlier , but a yellow hypergiant . Such extremly luminous stars are extremely rare, and only a dozen are known in the Galaxy . Yellow hypergiants (YHG) have high mass-loss rates ($10^{-5}-10^{-3} \ensuremath{\mathrm{M}_{\odot}{\rm yr}^{-1}}$) and are in a short, transitional evolutionary stage, thereby rapidly crossing the so-called "yellow void" in the Hertzsprung-Russell diagram . Their link to other advanced evolutionary phases of massive stars such as Luminous Blue Variables and Wolf-Rayet stars is still an open issue in stellar evolution theory.

Because of its large number of remarkable observational features, has been subject to extensive studies over the last 30 years. The spectral type of has changed from F8 Ia$^{+}$ in 1973 through A5Ia in the mid 90s to A2 in 2000 . Correspondingly, the effective temperature has changed by more than 3000 K in the last two and a half decades. This makes a unique object for the study of stellar evolution since it is one of the very rare stars believed to be in the rapid transition from the Red-Supergiant stage to the Wolf-Rayet phase.

HST/WFPC2 images of 's surrounding nebula revealed a variety of structures including condensations or knots, ray-like features, and several arcs or loops within $2^{''}$ from the star, plus one or more distant reflection shells. All these features suggest frequent episodes of high mass loss during the past centuries. found thermal SiO emission in a huge hollow shell around , with a typical radius of $\sim10^{17}\,$cm (=6680 AU=1 at a distance of 5 kpc), a shell width smaller than half of the radius, and an expansion velocity of 35 ${\rm kms}^{-1}$. mapped 's nebula in the $^{12}{\rm CO}$ J = 2-1 and 1-0 transitions and found that the nebula shows an approximately spherical, extended halo surrounding a bright inner region, with both components clearly presenting smaller aspherical features. The CO nebula expands isotropically with an expansion rate similar to the SiO shell.

The chemical composition of 's nebula is dominated by O-rich chemistry as suggested by several observations , and is similar to O-rich AGB stars as recently found by . is also among the brightest IRAS objects and one of the warmest stellar OH maser sources known . Ammonia emission has been reported by and . From CO observations, high mass-loss rates of the order of several $10^{-4}$ M$_{\odot}{\rm yr}^{-1}$ have been derived.

As inferred by many studies, the structure of the inner circumstellar envelope (CSE; spatial scale of milliarcseconds) of appears to be as complex as that of the outer nebula, and numerous scenarios have been proposed to explain its observed features. These models include a rotating equatorial disk , bipolar outflows , infall of circumstellar material , wind blowing in a preferential direction , and even the simultaneous presence of inflowing and outflowing matter .

Despite very intensive observing campaigns, the overall geometry of 's CSE is still a matter of debate. obtained HST/STIS spatially resolved spectroscopy of and its reflection nebula. They suggested that given the stellar temperature and the high mass-loss rate, the strong stellar wind of must be optically thick and, thus, the observed variations in the apparent spectral type and temperature are due to changes in the wind and do not reflect an interior evolution on short timescales. The detection of a nearly spherical CSE by is in marked contrast to other observations that reveal a rather axis-symmetric wind geometry in , such as the earlier HST images , integral-field spectroscopy , or recent spectropolarimetry . These latter observations suggest a symmetry axis at a position angle of $45^\circ\pm15^\circ$.

Infrared interferometric and coronographic observations of were reported by, e.g., , , , , , , , , and . presented diffraction-limited 73 mas $K$-band bispectrum speckle-interferometry observations of 's dust shell. They found that the $K$-band visibility steeply drops to a plateau-like value of $\approx0.6$ at 6 m baseline and, thus, concluded that 40% of the total $K$-band flux comes from the extended dust shell. The best radiative transfer model found by to simultaneously explain the spectral energy distribution (SED) and $K$-band visibility contains a two-component shell composed of silicate dust with an inner rim at 69 stellar radii where the dust temperature is 1000 K. Moreover, found that a phase of heavy mass loss with mass-loss rates approximately 40 times higher than the current $\dot{M}= 1.4\times 10^{-5}\,M_\odot{\rm yr}^{-1}\times(d{\rm kpc}^{-1})$ must have ceased roughly 60 to 90 yrs ago.

obtained both $K$-band aperture-masking observations of using the Keck I telescope and long-baseline interferometric observations with the beam-combiner instrument FLUOR at the IOTA interferometer . While the aperture masking observations basically confirmed the results of the speckle observations of , i.e. a sharp visibility drop for baselines shorter than 2 m and a plateau-like visibility of $\sim$`<!-- -->`{=html}0.6 up to a 8 m baseline, from the IOTA measurements found a $K$-band visibility of $\sim$`<!-- -->`{=html}0.7, from which they concluded that the compact stellar component is not resolved at baselines as long as $\sim$`<!-- -->`{=html}35 m.

Recently, presented the first near-infrared long-baseline interferometric observations of obtained with the AMBER instrument at ESO's Very Large Telescope Interferometer (VLTI). The observations presented by were carried out in medium spectral resolution mode ($\lambda/\Delta\lambda=1\,500$) around the Br$\gamma\,$emission line. resolved the Br$\gamma\,$line-emitting region and derived a Gaussian FWHM size of 3.3 milliarcseconds (mas), but due to calibration problems, the size of the continuum-emitting region could not be constrained.

In this paper, we present the first VLTI/AMBER observations of in the $H$ and $K$ bands in low-spectral resolution mode and the first AMBER measurements of its Br$\gamma$-emitting region with a spectral resolution of $\lambda/\Delta\lambda=12\,000$ (high-spectral resolution mode of AMBER) and baselines up to 128 m. From these new AMBER observations, the size of the continuum-emitting region in several spectral channels across the $H$ and $K$ bands could be derived and, by comparison of the high-spectral resolution AMBER data with 2-D gas radiative transfer models, details on the size and geometry of the Br$\gamma\,$line-emitting region could be obtained.

The observations presented here are the first step of an observing campaign to better constrain the geometry of the inner wind region of at a scale of a few stellar radii.

The paper is structured as follows. In Sect.  we will present the AMBER observations of and discuss the data reduction. In Sect.  we discuss the implications of AMBER's limited field-of-view for the observations of the very extended object . In Sect.  and , the absolute calibration, radiative transfer modeling, and interpretation of the AMBER data will be discussed. The paper closes with a summary and conclusions in Sect. .

# Observations and data reduction

## General remarks on the data reduction process

In the following two subsections we describe our VLTI/AMBER observations of in low (LR) and high spectral resolution (HR) mode. AMBER is the near-infrared closure-phase beam combiner of ESO's VLTI, simultaneously operating in the $J$, $H$, and $K$ bands. This instrument is capable of interfering the light of either three 1.8 m Auxiliary Telescopes (ATs) or three of the 8.2 m Unit Telescopes (UTs) located on Paranal. A detailed description of the AMBER instrument and its optical design is given in and . The AMBER data presented here were reduced using version 2.2 of the data reduction package *amdlib* [^2]. The reduction software is based on the so-called Pixel-to-Visibility-Matrix (P2VM) algorithm, which is described in .

Following previous experiences with the reduction of low and high-spectral resolution AMBER data , the data selection was carried as follows: As primary selection criteria we used the signal-to-noise ratio (SNR) of the fringe signal. For the LR data, we kept 20% of the data with the highest fringe SNR, while for the HR data taken with the fringe tracker FINITO we kept 80% of the data with the highest fringe SNR. In both cases, the optimal fringe SNR selection value was found by increasing the fraction of discarded data frames from 0 to 80%. A stronger selection than given above only led to an increase of the noise; i.e., finally there was only a decrease of the data quality rather than a significant improvement. In addition to the fringe SNR selection, to ensure a proper deselection of all frames which are far from zero optical path difference, we also discarded in the low- as well as high-spectral resolution data all frames with a piston of more than $10\,\mu$m. For the high-spectral resolution data, this seems to be a too strong selection, since in this case the coherence length is of the order of $\lambda/\Delta\lambda\times\lambda\approx2$ cm, but due to the use of FINITO the piston variation is typically much small than 100 micron in the HR data. Since it turned out that the HR results do not show a strong variation with piston selection, we finally decided to use the same piston selection criteria for both the HR and LR data.

Concerning the error estimates of our AMBER measurements, we were facing the problem that for each observation only one suitable calibrator measurement was available. Thus, in addition to the statistical error of the observables calculated from the average over all frames of a given data set, it was problematic to address also any systematic error sources. Therefore, to account for systematic uncertainties of the absolute calibration process, we looked up previous AMBER measurements taken under similar weather conditions. From these comparisons, we concluded to add an absolute error of 0.03 to all visibility data and 3 degrees to all phases to account for the systematic uncertainties.

## The low-spectral resolution data

The low-spectral resolution (hereafter LR) AMBER observations of were obtained with the Auxiliary Telescopes (ATs) in May/June 2007 and May 2008 using the linear baseline configurations E0-G0-H0 and A0-D0-H0. A summary of the LR observations is given in Table . As this table as well as the $(u,v)$ coverage in Fig.  shows, the total baseline and position angle ranges covered by the LR observations of are 15 to 94 m and 69 to 75, respectively, i.e. the observations were essentially carried out in the same direction on the sky. All measurements were made under sub-average to average seeing conditions (see Tabs.  and  ). For the calibration of all LR AMBER data on we used the calibrator HD 190327. A record of the calibrator observations is also provided in Table .

$$\begin{array}{cccccc}
            \hline
            \noalign{\smallskip}
            {\rm date}& {\rm UTC}& {\rm AT}           & B_{\rm p}        & PA      & {\rm seeing}\\
                      &          & {\rm configuration}& [{\rm m}]        & [\degr] & [\arcsec]   \\
             \noalign{\smallskip}
            \hline
            \noalign{\smallskip}
    \multicolumn{6}{c}{{\rm \object{IRC\,+10420 }}}\\ \hline
     2007-05-17   & 08:11:53  & {\rm A0-D0-H0}    &31-93             & 75      & 1.0\\
     2007-05-17   & 09:18:51  & {\rm A0-D0-H0}    &32-96             & 72      & 1.4\\
     2007-06-18   & 08:09:22  & {\rm E0-G0-H0}    &15-48^{\mathrm{a}}& 69      & 0.7\\
     2008-05-21   & 08:12:14  & {\rm A0-D0-H0}    &31-94             & 74      & 0.8\\ \hline
    \multicolumn{6}{c}{{\rm HD~190327}}\\ \hline
         2007-05-17   & 08:48:23  & {\rm A0-D0-H0}    &31-93             & 74      & 1.2 \\
     2007-06-18   & 08:39:13  & {\rm E0-G0-H0}    &16-47             & 70      & 0.6 \\
     2008-05-21   & 08:48:41  & {\rm A0-D0-H0}    &32-95             & 74      & 0.7 \\

            \noalign{\smallskip}
            \hline
         \end{array}$$

Unfortunately, due to technical problems during the observations, from the measurements carried out in June 2007 reliable data could only be retrieved from one of the three baselines. Thus, in this case only one visibility instead of a triplet was obtained. Due to unfavorable weather conditions and limitations of the technical performance, reliable differential and closure phases from the LR observations could not be derived for any of the 2007 observations, but only for those from May 2008.

The results of the AMBER LR observations of are shown in Fig. . The panels in the left and middle column show the AMBER $K$-band visibilities as a function of wavelength for the three epochs covered by our study, the top right panel shows the the AMBER $H$-band visibilities as a function of wavelength, and the bottom right panel illustrates the differential and closure phases in the $J,H$, and $K$ bands obtained from the May 2008 observations of .

Two main results can be seen in Fig. . First, both the $H$- and $K$-band visibilities show only a weak wavelength dependence, and the visibility change with baseline for a given wavelength is only moderate. This corresponds to the fact that even with the shortest baselines of our AMBER observations, the extended dust shell of is already fully resolved, and thus, AMBER probes the compact stellar component, i.e. the continuum emission from the central star and the dense circumstellar wind. The second result from the calibrated LR observables shown in Fig.  is that both the differential and closure phases show zero phase signals across all three near-infrared bands within the uncertaintly limits. Thus, in the LR AMBER data which mainly probe the continuum emission does not show detectable deviations from point symmetry.

## The high-spectral resolution data

In addition to the low-spectral resolution $H$- and $K$-band AMBER observations of presented in the previous section, we also obtained the first AMBER measurement of in high-spectral resolution mode (spectral resolution $\lambda/\Delta\lambda = 12\,000$). The observations were carried out in October 2008 under average seeing conditions as part of an AMBER science verification run in order to test the performance of AMBER using the UTs and the fringe tracker FINITO. In the context of this science verification run, a single measurement of and the calibrator star HD 232078 was carried out in the $K$ band, centered around the Br$\gamma\,$emission line ($\lambda_{\rm c}= 2.166\,\mu$m), with an exposure time of 3 s. Due to this long exposure time, a window covering 512 pixels, i.e. all spectral channels of the AMBER detector, could be read out, corresponding to a wavelength coverage between $\sim 2.145$ and $\sim 2.19\,\mu$m.

$$\begin{array}{lccccc}
            \hline
            \noalign{\smallskip}
            {\rm object}& {\rm UTC}& N_{\rm frames}& B_{\rm p}   & PA      & {\rm seeing}\\
                      &            &               & [{\rm m}]   & [\degr] & [\arcsec]   \\
             \noalign{\smallskip}
            \hline
            \noalign{\smallskip}
        {\rm \object{IRC\,+10420 }}      & 00:26:40  & 125            & 54-129      & 36-75   & 0.9\\
        {\rm HD~232078} & 00:03:59  & 125            & 51-130      & 37-78   & 1.1\\
            \noalign{\smallskip}
            \hline
         \end{array}$$

It should be noted that a second measurement of with an exposure time of 1 s was discarded for two reasons. First of all, the data quality of this second measurement is lower than that with the longer detector integration time (DIT), simply because the SNR is lower. Second, there is no corresponding calibrator measurement with the same DIT close in time to the science target measurement. From the ESO data archive, we found two other calibrator measurements with DIT=$1\,$s from different science programs in the same night, but these measurements of the calibrators HD 902 and HD 13692 were taken 1 and 2 hours after the observations of . Since current high-spectral resolution AMBER data suffer from a time-variable high-frequency beating introduced by the VLTI Infrared Image Sensor IRIS , we decided to discard these additional measurements because a proper calibration of the science target data turned out to be impossible. This seems to be a general problem for observations in high-spectral resolution mode if the time span between the science target and calibrator measurement becomes longer than approximately 30-45 minutes. For the same reason, we also discarded two other calibrator data sets with the same exposure time (DIT=3 s) as the measurement of presented here since these calibrator observations were carried out $\sim5.5$ and 8 hours after the observations. To further reduce the effect of the IRIS fringes in our data, we spectrally binned our final visibilities and phases with a bin size of 7 and a step size of 1. To avoid a loss of spectral resolution across the Br$\gamma\,$line, this binning was only performed in the continuum region, i.e. for wavelength channels with $\lambda < 2.1660\,\mu$m or $\lambda > 2.1677\,\mu$m. A summary of the HR observations discussed in this paper is given in Table .

The calibrated observables of the AMBER high-spectral resolution measurements are presented in Fig. . The panels in the three left columns show (from left to right) the visibilities, differential phases, closure phases (top), and the spectrum (bottom) for the full wavelength range covered by the observations. All panels in the three columns on the right-hand side of Fig.  show the same quantities as the corresponding panels in the three left columns, but only for a small wavelength range close to the Br$\gamma\,$emission line. For the sake of clarity, in each panel, only a single averaged error bar is shown on the right-hand side.

The wavelength calibration of the HR data was performed by comparison of the AMBER raw spectra of both and the calibrator HD 232078 with high-spectral resolution ($\lambda/\Delta\lambda=60\,000$) telluric spectra kindly provided by the Kitt Peak Observatory[^3]. The vertical dashed line shown in all panels on the right-hand side of Fig.  marks the zero velocity with respect to earth, assuming a radial system velocity of 73 ${\rm kms}^{-1}$ as given by , and a heliocentric correction of -26 ${\rm kms}^{-1}$ for the time of the observations. The spectrum shown in Fig.  is the average over the photometric beams from the three single telescopes and the interferometric signal.

The calibration of the high-spectral resolution AMBER measurements is more uncertain compared to the low-spectral resolution data shown in Sect.  for two reasons. First of all, for the HR data, we had only one calibrator measurement instead of three as in the case of the LR data. Second, the HR data were recorded with the fringe tracker FINITO. Since the performance of FINITO depends on, e.g., a target's correlated flux and the weather conditions at the time of the observations, the final calibrated observables obtained from AMBER+FINITO observations can be biased either due to a significant brightness difference between science target and calibrator and/or due to different weather conditions (seeing, coherence time, wind) during the science target and calibrator measurements.

In the case of the HR measurements of presented in this paper, the calibration relies on only a single calibrator measurement of HD 23078 which is approximately 0.5 magnitudes fainter in the $K$ band than ($K_{\rm HD23078} = 4.2^{\rm m}; K_{\rm IRC+10420}=3.65^{\rm m}$). From the weather data[^4], we can see that at the time of the calibrator and science target measurements the seeing was varying between 0 and 1, and the wind speed on the ground was steadily declining from 9 to 7.5 ${\rm ms}^{-1}$. Thus, there is a noticeable, but not dramatic change in the weather conditions between the calibrator and observation. Together with the 0.5 magnitude brightness difference, this might explain why, for instance, the calibrated visibility at the shortest baseline (see Fig. ) is slightly larger than unity. In fact, when using data from other calibrators of the same night (HD 25680; HD 33833) which were taken several hours later and with a seeing of 1, the calibrated visibility of reaches values $> 1.2$.

While the weaker emission peaks seen in the full AMBER spectrum of in Fig.  are telluric artifacts of the calibration process and therefore do not indicate real emission features, the spectrum clearly shows a strong Br$\gamma\,$emission line arising from which contributes approximately 50% to the total flux at $\lambda=2.1666\,\mu$m. The Br$\gamma\,$line is clearly blueshifted with $v=-25\,{\rm kms}^{-1}$ with respect to the systemic velocity, in agreement with previous findings by, e.g., ($v=-30\,{\rm kms}^{-1}$) and ($v=-22\,{\rm kms}^{-1}$). From the $\lambda$-calibration process of our data, we estimate that the uncertainty of the wavelength calibration is of the order of $10\,{\rm kms}^{-1}$. The equivalent width of the Br$\gamma\,$line in the AMBER spectrum is $-6.7~\AA$, in good agreement with the value found by from medium-spectral resolution AMBER observations in June 2006. The Br$\gamma\,$line is slightly asymmetric in the sense that the redshifted tail of the line is more pronounced with velocities up to $\sim120\,{\rm kms}^{-1}$ compared to $\sim80\,{\rm kms}^{-1}$ on the blueshifted side.

At the wavelengths of the Br$\gamma\,$emission, the visibilities on all three baselines show a strong decrease, indicating that the line-emitting zone is clearly resolved by our AMBER observations and more extended than the continuum-emitting region. At the longest baseline, the visibility in the center of the Br$\gamma\,$line is as low as 0.14. It should be noted here that similar to the LR data presented in Sect.  the HR data also suffer from a strong field-of-view (FOV) effect as will be discussed in more detail in the next section.

The differential phases in the region of the Br$\gamma\,$-emission line show a clear non-zero signal within the uncertainty limits. On all three baselines, we obtained the strongest phase signal in the redshifted wing of the Br$\gamma\,$-emission line, with phases up to $\phi_{ij}=-20^\circ$ on the two shorter baselines and $\phi_{13}=-26^\circ$ on the longest baseline. The closure phase $\Phi$, on the other hand, is in good agreement with the differential closure phase $\Phi_{\rm diff} = \phi_{12} + \phi_{23} - \phi_{13}$. $\Phi$ also shows a strong signal in the redshifted wing of the Br$\gamma\,$-emission line ($\sim-18^\circ$) and the strongest non-zero signal ($\sim-27^\circ$) at the line center.

# Results and discussion

## The field-of-view effect

Since VLTI/AMBER is a single-mode fiber instrument, the field-of-view (FOV) is limited to the Airy disk of the telescope aperture. Therefore, for AMBER observations with the 1.8 m ATs of the VLTI, the FOV is 250 mas in the $K$ band, while for AMBER observations with the 8.2 m UTs, the FOV is only $\sim60$ mas. Thus, when observing an object as extended as with AMBER, it is of great importance to take AMBER's limited FOV into account, since a non-negligible fraction of the total flux will be located well outside the FOV of the observations.

To account for the FOV effect in the AMBER data of , we used the dust radiative transfer model of . This radiative transfer model was developed with the code DUSTY to simultaneously explain the spectral energy distribution (SED) and $K$-band speckle-interferometric observations obtained with the SAO 6 m telescope.

The azimuthally averaged $K$-band visibility of is shown in the top-left panel of Fig.  as a function of baseline length. As one can see, the steep visibility drop ends in a plateau at a visibility level of $\sim$`<!-- -->`{=html}60%. The plateau itself is associated with the unresolved stellar component plus the dense stellar wind of , which can be probed by interferometry with baselines longward of approximately 10 m. On the other hand, the extended dust shell is fully resolved already with a few meter baseline. According to the best-fitting model of Blöcker et al., as a result of its heavy mass loss is surrounded by an optically thick dust shell composed of silicate dust, which contributes approximately 40% to the total $K$-band flux of . The inner boundary of the dust shell which marks the dust sublimation radius is located at $69~R_{\star}\approx 69$ mas and exhibits a dust temperature of 1000 K. Moreover, Blöcker et al. found that in addition to the steady stellar wind a phase of enhanced mass loss (superwind) that ceased roughly 60 to 90 yrs ago is required to explain both the $K$-band visibility and the SED. In the intensity profiles shown in Fig. b, this temporal superwind phase results in the bump seen at an angular displacement of $\sim 170$ mas. Thus, according to Blöcker et al. 's dust shell comprises two components: An outer shell resulting from the previous superwind phase and, directly adjacent to this outer shell, an inner shell resulting from the present-day mass loss which followed the superwind phase. During the superwind phase, the mass-loss rate was approximately 40 times higher than the current mass-loss rate of $1.4\times 10^{-5}\,M_\odot{\rm yr}^{-1}$.

Since the speckle data of Blöcker et al. were taken with a comparably large FOV ($\sim14$`<!-- -->`{=html}), the infinite FOV of the DUSTY model was well suited for the speckle data analysis. In the context of our AMBER observations, we re-examined the final radiative transfer model of Blöcker et al. as follows. We used the radial $H$- and $K$-band intensity profiles of their final model (see the top-right panel in Fig. ) and calculated the visibilities from the intensity profiles, which were truncated to match a given FOV. We note that it would have been more precise to convolve the intensity profiles with the corresponding telescope aperture before the truncation, but the effect of the convolution on the final result is only weak and, therefore, has been omitted here. The $H$- and $K$-band visibilities resulting from our FOV correction are shown in the middle panels of Fig.  for the FOVs given by the plot labels. As expected, the smaller the FOV, the higher the plateau value of the visibility, since an increasing fraction of the extended dust shell lies outside the FOV and will, therefore, no longer contribute to the flux and visibility, respectively.

In the bottom right panel of Fig. , the AMBER visibility correction for is shown as a function of wavelength across the $H$ and $K$ bands. This wavelength-dependent correction has to be subtracted from the measured visibility in order to correct AMBER's finite FOV when observing with the ATs. The figure reveals that the visibility correction decreases with increasing wavelength from approximately 0.42 at $1.45\,\mu$m to 0.1 at $2.45\,\mu$m. This decrease reflects the decreasing flux contribution of the dust shell (scattering + thermal emission) with increasing wavelength (see upper right panel in Fig. ). For instance, at $1.5\,\mu$m the contribution of the attenuated stellar flux to the stellar flux is only 35%, while it is almost 60% at $2.2\,\mu$m. Correspondingly, the FOV correction is stronger in the $H$ band than in the $K$ band. At the center of the $H$ and $K$ bands the visibilities of the unresolved stellar component in the model with a 250 mas FOV are approximately higher by 0.18 and 0.35 compared to the DUSTY model with an infinite FOV. For a 60 mas FOV, the effect is even more dramatic. Here, almost the complete dust shell is truncated due to the small FOV, and the visibility approaches unity, since mainly the compact stellar component contributes to the total flux. According to Fig. , we can expect that only $\sim 6$% instead of the originally 40% of flux from the extended dust shell remains, if the FOV is as small as 60 mas.

## The low-spectral resolution data

As discussed in the previous section, the AMBER visibilities of obtained with both ATs and UTs are highly affected by the limited FOV of the observations. From the radiative transfer model of , we estimated the effect of AMBER's limited FOV on the $H$- and $K$-band visibilities for baselines which start to resolve the stellar component and its dense wind. We lowered all LR $H$- and $K$-band visibilities shown in Fig.  according to the wavelength-dependent visibility correction shown in the bottom panel of Fig. . Figure shows the rescaled $K$-band AMBER visibilities from May 2007 and May 2008 (top panels) and the rescaled $H$-band AMBER visibilities from May 2008 (bottom left panel) which would be obtained in the case of an infinite AMBER FOV.

As the bottom-right panel in Fig.  illustrates, the rescaled AMBER $K$-band data at $\lambda=2.11\,\mu$m (red bullets with error bars) are in good agreement with the visibilities from the speckle-interferometric observations by (blue bullets) within the uncertainties of both measurements. The solid green curve in Fig.  shows a Gaussian fit of the compact stellar component, assuming a fully resolved dust shell with an a priori unknown flux contribution. As indicated in the plot, a two-parameter fit with stellar diameter $d_{\star}$ and the fractional flux contribution of the dust shell as $f_{\rm dust}$ free parameters gives a Gaussian FWHM diameter of $d_\star=0.98\pm0.16$ mas and a dust shell flux contribution of $f_{\rm dust}=0.35\pm0.02$. The fitted diameter is in agreement with the stellar diameter $d=1.01$ mas derived for a bolometric flux of $8.2\times~10^{-10}\,{\rm Wm}^{-2}$, a distance of 5 kpc, and a central star effective temperature of 7000 K .

After rescaling all AMBER $H$- and $K$-band visibilities, we fitted all visibility points for a given wavelength with the simple two-component model described above (Gaussian stellar component + infinitely extended dust shell). The result of the fit procedure is shown in Fig. , where the fit parameters are displayed as a function of wavelength for the $H$- (bottom panel) and $K$-band data (top panel). As the figure reveals, within the error bars we essentially obtained both a wavelength-independent continuum diameter and a wavelength-independent flux contribution from the extended dust shell across the $K$ band. On the other hand, a slight decrease of the dust shell's flux contribution is seen across the $H$ band from $f_{\rm dust}\approx 0.64$ at the lower band edge to $f_{\rm dust}\approx 0.54$ at the upper band edge. A comparison with the lower left panel in Fig.  shows that the wavelength dependence and the absolute level of the stellar flux contribution derived from the AMBER data is in basic agreement with the predictions from the DUSTY model of . In this model, the attenuated stellar flux amount to approximately 0.58 in the $K$ band with only a moderate wavelength dependence, while the stellar flux contribution rises from approximately 0.32 to 0.5 across the $H$ band.

Averaged over all $K$-band data, for we derived a stellar continuum diameter of $d_{\rm Gauss} = 0.98\pm0.10$ mas ($d_{\rm UD} = 1.65\pm0.14$ mas assuming a uniform disk) and a fractional flux contribution of $1-f_{\rm\star} = f_{\rm dust} = 0.41\pm0.05$ from the dust shell. From the $H$-band AMBER data, we found $d_{\rm Gauss} = 0.99\pm0.07$ mas ($d_{\rm UD} = 1.63\pm0.12$ mas) and $f_{\rm dust} = 0.44\pm0.01$.

## The high-spectral resolution data

### FOV effect correction

To account for the larger uncertainties in the absolutely calibrated data, we used a slightly different strategy to correct the HR AMBER data for the FOV effect described in Sect. . Looking at Fig. e we see that for AMBER+UT observations with a FOV of $\sim60$ mas, the fractional flux contribution of the dust shell in the $K$ band is only of the order of 6%. On the other hand, from the fit of our LR AMBER data, we found that 's $K$-band continuum diameter is 0.98 mas. Therefore, we corrected our HR AMBER visibilities in the following way: For a given baseline, we first rescaled the visibility as a function of wavelength to match the value expected for a single, compact component with 0.98 mas diameter. Then, we globally lowered the rescaled visibilities by 6% to account for the diffuse, extended flux contribution of the dust shell.

The final result of this rescaling process in order to take the FOV effect into account is shown in the upper panel of Fig. . Here, we show for a small wavelength range around the Br$\gamma\,$line the $K$-band visibilities one would retrieve for in the case of an infinite AMBER FOV. As the figure reveals, the continuum visibilities are now $\sim0.9$ for the shortest baseline, 0.83 for the intermediate baseline, and 0.68 for the longest baseline.

### Simple geometrical model

Attempting to get a first idea of the size of the Br$\gamma$ line-emitting region and disregarding a possible deviation from spherical symmetry, we fitted a spherical Gaussian to all continuum data for a given wavelength with $\lambda < 2.1660\,\mu$m and $\lambda < 2.1667\,\mu$m, taking into account the 6% flux contribution from the extended dust shell. Averaged over all spectral channels (approximately 490 channels), we obtained a FWHM Gaussian size of $d_{\rm continuum} = 1.014\pm0.01$ mas for the continuum-emitting region. We note that the error given here includes not only the the pure fitting error, which is only of the order of 0.001 mas, but also reflects the uncertainties of the overall calibration process, including the FOV-related rescaling of the data.

To derive the size of the Br$\gamma\,$line-emitting zone, we fitted the visibilities in the spectral channels with noticeable Br$\gamma\,$line-emission in the spectrum, i.e. for $2.1660\,\mu{\rm m}\la\lambda\la2.1667\,\mu$m, with a three-component model. For the sake of simplicity, in this simple model, we neglect likely deviations from a spherical or even point-symmetric distribution of the Br$\gamma\,$line-emitting region. convincingly argued that the Br$\gamma\,$line-emitting region cannot have a spherical shape, but is most likely elongated towards a position angle of 20. Unfortunately, as in the case of the study of , also our single measurement with a position angle coverage of $\sim40$`<!-- -->`{=html}  is only of limited use to put reliable constraints on the shape of the emission zone. Nevertheless, we studied the emission region also in the context of a 2-D radiative transfer model as will be shown in more detail in the next section. However, independent from a more sophisticated modeling, qualitatively we can conclude that the clear non-zero detection in the differential and closure phases across the Br$\gamma\,$line in our AMBER HR measurement shows that the geometry of the line-emitting region deviates from a point-symmtric configuration. On the other hand, in the center of the Br$\gamma\,$line the differential and closure phases are of the order of $\la10$`<!-- -->`{=html}. Thus, in order to determine the typical size of the Br$\gamma\,$-emitting zone with the simple model approach discussed here we decided to neglect the non-zero phases.

Our simple three-component model consists of the following components: (a) a fully resolved, infinitely large dust shell with a fixed fractional flux contribution of 6%, (b) the stellar continuum, approximated by a Gaussian with a fixed size of 0.98 mas according to our continuum fit, and (c) the Br$\gamma\,$line-emitting region, which is represented in our simple model by a ring with infinitely thin thickness ($\Delta r = 0$). The total visibility can then be written as

$$\label{eq1}
V_{\rm total} = f_{\rm dust} \times V_{\rm dust} + f_{\rm cont} \times V_{\rm cont} + f_{\rm ring} \times V_{\rm ring}$$ with $$f_{\rm dust} + f_{\rm cont} + f_{\rm ring} = 1.0 \quad {\rm and} \quad f_{\rm dust} = 0.06 \quad {(\rm fixed)},$$ $$V_{\rm dust} = 1 \quad {\rm for} \quad B_{\rm p} = 0 \quad {\rm and} \quad V_{\rm dust} = 0 \quad {\rm elsewhere},$$ $$%V_{\rm cont} = \exp \{-(\pi d_{\rm cont}B_{\rm p}/\lambda)^2/(4.\log(2.))\}
V_{\rm cont} = \exp \left\{ -\frac{\left(\pi d_{\rm cont}B_{\rm p}/\lambda\right)^2}{4.\times\log(2.)}\right\}$$ with $d_{\rm cont} = 0.98$ mas (fixed) and $$\label{eq6}
V_{\rm ring} = J_0(\pi d_{\rm ring}B_{\rm p}/\lambda).$$

$B_{\rm p}$ is the projected baseline, $J_0$ is the Bessel function of 0th order, and $d_{\rm ring}$ is the diameter of the infinitely thin ring representing the Br$\gamma\,$line-emitting region. As can be seen from the above equations, our fit procedure contains two free parameters: the ring diameter $d_{\rm ring}$ and the fractional flux contribution from the Br$\gamma\,$line-emitting region of 's optically thick stellar wind, $f_{\rm ring}$. In principle, one could also fix $f_{\rm ring}$ and $f_{\star} = 1 - f_{\rm ring} - f_{\rm dust}$ according to the spectrum shown in Fig. , if the continuum emission is interpolated across the Br$\gamma\,$line. For instance, a linear interpolation from the spectral channels measuring pure continuum emission would give $f_{\star} = 0.415$ and $f_{\rm ring} = 0.525$. But since it is not clear a priori to which extent photospheric absorption is present , we decided to keep the fractional flux contribution from the stellar wind component as a free parameter.

Our best-fit result for the central wavelength channel in the Br$\gamma\,$line is shown in the middle panel of Fig. . Here, the visibility according to Eqs. ()-() is plotted as a function of baseline for $\lambda=2.16664\,\mu$m. The blue bullets with error bars are the three AMBER data points, and the three curves visualize the best fit model and its tolerance. We found that our AMBER data in the center of the Br$\gamma\,$emission line can be well fitted with a ring with a diameter of $d_{\rm ring} = 4.18^{+0.19}_{-0.09}~{\rm mass}\,\approx 4.12 \times d_{\star}$, whose fractional flux contribution is $f_{\rm ring} = 0.58^{+0.02}_{-0.04}$. We note that the size of the Br$\gamma\,$line-emitting region obtained from our fit is in general agreement, although slightly larger than the 3.3 mas found by from a single-component Gaussian fit of the line data. From $f_{\rm ring}$, we directly obtain $f_{\star} = 0.36^{+0.02}_{-0.04}$; i.e., from our best fit, we indeed find a photosperic absorption in the Br$\gamma\,$line of the order of $f_{\rm abs} = (f_{\star, {\rm cont}} - f_{\star, {\rm line}})/f_{\star, {\rm cont}} =
(0.415 - 0.36)/0.415 = 0.13$.

Despite the limited amount of data and the limited position angle coverage, we also briefly addressed the question of an aspherical shape of the Br$\gamma\,$line-emitting region in the context of our simple model fits. For this study, we used the three-component model described above and fixed $f_{\rm ring}$ to the value found from the 1-D fit ($f_{\rm ring} = 0.58$) of the $\lambda=2.16664\,\mu$m data. Then, we first fitted all three data points with the three-component model with only $d_{\rm ring}$ as free parameter, and finally, we fitted an ellipse to the three ring diameters as a function of position angle. The result is shown in the bottom panel of Fig. , where the fitted ring diameters and the ellipse are displayed in a polar diagram. As the figure shows, the position-angle dependence of the ring diameter can be fitted with an ellipse with major and minor axes of $4.32\pm0.44$ mas and $3.86\pm0.37$ mas, respectively, and a position angle of the major axis of $\Phi=36\pm11$`<!-- -->`{=html}. Interestingly, the major axis is nearly aligned with the elongation of the outer reflection nebula of seen in the HST images of and the symmetry axis of the H$\alpha$ emission as inferred by .

### Radiative transfer modeling of the Br$\gamma$ line emission

A simple geometrical model as presented in the previous section cannot explain the AMBER observations in detail. Thus, a more physical model is required to interpret, for instance, the blueshift of the Br$\gamma$ emission, the visibility across the emission line, and the non-zero phases. Therefore, in a second step of our analysis we interpreted the AMBER data of using the iterative, spherical symmetric, full line-blanketed, non-LTE radiative transfer code CMFGEN .

#### The radiative transfer code

In CMFGEN, a spherically-symmetric outflow in steady-state is assumed, and line and continuum formation are calculated in a non-LTE regime. Each model is specified by its effective temperature $\mathit{T}_{\rm eff}$, luminosity $\mathit{L}_{\star}$, effective gravity $\mathrm{g}_{\rm eff}$, mass-loss rate $\dot{M}$, wind terminal velocity $v_{\infty}$, and chemical abundances $Z_i$ of the included species. Moreover, a velocity law must be adopted, since the momentum equation of the wind is not solved. CMFGEN can only handle monotonically increasing velocity laws and, therefore, no inflows can be analyzed using this code. The velocity structure $v(r)$ is parameterized by a beta-type law, which is modified at depth to smoothly match a hydrostatic structure at $\ensuremath{\mathit{R}_{\star}}$. We define the hydrostatic radius as the radius where the wind velocity is equal to one third of the sonic velocity ($v=v_\mathrm{sonic}/3$), in order to avoid any effect of the dense wind on the determination of $\ensuremath{\mathit{R}_{\star}}$. Close to and below $\ensuremath{\mathit{R}_{\star}}$, the hydrostatic equation is iterated in order to produce a quasi-hydrostatic structure which extends inwards to the inner model boundary (at a Rosseland optical depth of 100). The effective temperature is defined as the temperature where the Rosseland optical depth is 2/3, i.e., $\ensuremath{\mathit{T}_{\rm eff}}=T(\tau_\mathrm{Ross}=2/3$).

While CMFGEN can handle the effects of clumping via a volume filling factor $f$, we decided to assume a unclumped wind (i.e. $f=1$) since no strong electron scattering wings are present around the Br$\gamma\,$line. The strength of such electron scattering wings is sensitive to the amount of wind clumping , and even assuming an unclumped wind model (which produces the maximum strength possible in the electron scattering wings), no detectable electron scattering wings are seen in Br$\gamma\,$.

Line blanketing affects the ionization structure of the wind and, consequently, the spectrum. CMFGEN employs the concept of super-levels to include thousands of spectral lines in non-LTE, making it feasible to solve the equations of statistical equilibrium and radiative transfer simultaneously. The atomic model used for included lines of H, He, C, N, O, and Fe.

Although CMFGEN handles only spherical-symmetric outflows, the wind asymmetry was also analyzed using a recently developed modification in CMFGEN to compute the emerging spectrum in two-dimensional geometry . We refer the reader to for further details about the code, and to that paper and for additional details about applications of the BH05 code. In the following, we briefly describe the main aspects of the code.

As input the BH05 code requires several quantities (energy-level populations, ionization structure, temperature structure, emissivities, opacities, and specific intensity $J$) calculated by the original, spherically-symmetric CMFGEN model in the co-moving frame. Latitude-dependencies of the wind density and wind terminal velocity can then be taken into account using any arbitrary latitude-dependent density/wind terminal velocity variation. Due to our very limited amount of data, we restricted the analysis to changes only in the wind density, using oblate and prolate density parameterizations as follows,

$$\mathrm{oblate}: \,\,\, \rho_{2D} \propto \rho_{1D} \{1+ [a (1-cos^2 \theta)^b]\}\,\,$$ $$\mathrm{prolate}: \,\,\, \rho_{2D} \propto \rho_{1D} (1 + a . cos^b \theta)\,\,,$$

where $\theta$ is the latitude angle (0$^{\circ}$=pole, 90$^{\circ}$=equator). A scaling law is used to ensure that the 2-D model has the same mass-loss rate as an equivalent spherically symmetric model.

The scaling laws given above are applied only for distances greater than that where the wind velocity is equal to the sonic velocity. For smaller distances, no scaling was applied. The 2-D source function, emissivities, and opacities are then calculated, assuming that these quantities depend only on the new values of the scaled density. Appropriate scaling laws are used for different physical processes (e.g. density-squared scaling for free-free and bound-free transitions, and linear-density scaling for electron scattering). Finally, for a direct comparison with the AMBER observables the spectrum in the observer's frame and intensity maps for all wavelengths of interest are computed. In our analysis, we first considered a spherically symmetric model and then a model with a prolate/oblate wind structure.

#### Spherical model

The stellar and wind parameters of our best-fitting spherically symmetric CMFGEN model are presented in Table . Our CMFGEN model has a luminosity of $\ensuremath{\mathit{L}_{\star}}=6 \times 10^5~\ensuremath{\mathrm{L}_{\odot}}$, in good agreement with previous determinations Together with an assumed mass of $M=20~\ensuremath{\mathrm{M}_{\odot}}$, we find $\ensuremath{\log \mathrm{g}_{\rm eff}}=0.8$. Since has already lost about half of its initial mass , it is likely that the He content is enhanced on the surface compared to the solar value. Thus, apart from solar metallicity our model assumes a He abundance of He/H=0.4 (by number). Finally, we assumed a distance of 3.5 kpc.

The AMBER high-resolution $K$-band spectrum covers only a short wavelength region around Br$\gamma\,$, and does not provide enough diagnostic lines to derive a precise value for the effective temperature of . However, the absence of $2.16137\mu$m and $2.16474\mu$m lines imply an upper limit for $\ensuremath{\mathit{T}_{\rm eff}}$ of about 10500 K. Since obtained $\ensuremath{\mathit{T}_{\rm eff}}\simeq9200~\mathrm{K}$ for the observations of in 2000 and an annual increase in $\ensuremath{\mathit{T}_{\rm eff}}$ of $\sim 120$ K, we assume in our models $\ensuremath{\mathit{T}_{\rm eff}}=9600~\mathrm{K}$. We note that conclusions of the present paper such as the evidence for wind asymmetry in , are only weakly affected by reasonable changes in $\ensuremath{\mathit{T}_{\rm eff}}$.

Our model requires a mass-loss rate of $\ensuremath{\dot{M}}=1.7 \times 10^{-5}~\ensuremath{\mathrm{M}_{\odot}{\rm yr}^{-1}}$ in order to reproduce the equivalent width of the observed Br$\gamma$ line in ($-6.7 {\AA}$, Fig. ). A wind terminal velocity of $\ensuremath{v_{\infty}}=70~\ensuremath{{\rm km\,s^{-1}}}$ and $\beta=1.5$ were used in order to roughly reproduce the Br$\gamma$ line width. However, the Br$\gamma\,$-emitting region of such a model is rather compact, not reproducing the spatial extension of the Br$\gamma\,$emitting region as measured by AMBER. As a consequence, this model fails to reproduce the observed visibilities (see red curves in the lower panels of Fig. ). Moreover, while the observed spectrum is blueshifted relative to the systemic velocity, the CMFGEN model spectrum is slightly redshifted ($\sim 10\ensuremath{{\rm km\,s^{-1}}}$). The redshift in the model is caused by the turbulent velocity field of the wind , which was parameterized using a microturbulent velocity of $10~\ensuremath{{\rm km\,s^{-1}}}$.

An additional model with the same parameters as the model described above, but with $\ensuremath{\dot{M}}=5 \times 10^{-5}~\ensuremath{\mathrm{M}_{\odot}{\rm yr}^{-1}}$, was calculated in order to analyze the effects of a larger $\ensuremath{\dot{M}}$ on the fit of the AMBER observables. As expected, the increase in $\ensuremath{\dot{M}}$ results in a larger Br$\gamma$-emitting region, but while such a model with an increased mass-loss rate is now able to reasonably reproduce the observed visibilities (see blue curves in the lower panels of Fig. ), it also predicts a much higher amount of Br$\gamma\,$line emission compared to the observations (lower left panel in Fig. ). On the other hand, lower turbulent velocities provide slightly lower Br$\gamma\,$redshifts. However, as a general result we conclude that spherically symmetric wind models are intrinsically unable to reproduce the observed *blueshifted* Br$\gamma\,$line. Moreover, a reasonable simultaneous fit of both the observed Br$\gamma\,$spectrum and the corresponding visibilities with a spherical model could not be obtained.

```latex
\begin{table*}
\caption{Physical parameters of the CMFGEN models for \object{IRC\,+10420 }.}
\label{irc_param}
\centering
\begin{tabular}{l c c c c c c r}
\hline
\hline
\multicolumn{7}{c}{spherical model} \\
\cline{1-7}
\hline
log \ensuremath{\mathit{L}_{\star}}/\ensuremath{\mathrm{L}_{\odot}}& \ensuremath{\mathit{T}_{\rm eff}}& \ensuremath{\log \mathrm{g}_{\rm eff}}&  $\ensuremath{\dot{M}}$                  & \ensuremath{v_{\infty}}& $\beta$ &     He/H    \\
                 &   (K)  &       &  ($\ensuremath{\mathrm{M}_{\odot}{\rm yr}^{-1}}$)              & ($\ensuremath{{\rm km\,s^{-1}}}$) &         &   (number)  \\
\hline
  5.78           &   9600 &  0.8  &  $1.7\times10^{-5}$       & 70       & 1.5     &     0.4 \\
\hline
\\
\hline
\hline
\multicolumn{7}{c}{two-dimensional model} \\
\cline{1-7}
\hline
Density distribution & $a$ & $b$ & i         & PA        & $\ensuremath{\dot{M}}$             & \\
                     &     &     & ($\degr$) & ($\degr$) & ($\ensuremath{\mathrm{M}_{\odot}{\rm yr}^{-1}}$)         & \\
\hline
Prolate              &  15 & 2   & 45        & 60        &  $3.5\times10^{-5}$   &\\
\hline
%Oblate & 15 & 2 \\
\hline
\end{tabular}
\end{table*}
```

#### 2-D modeling of the AMBER data around the Br$\gamma$ line - evidence for wind asymmetry

In the previous paragraph we showed that spherical models are unable to simultaneously reproduce the observed AMBER visibilities and the Br$\gamma$ spectrum of and, in particular, the blueshift of the Br$\gamma$ emission. Thus, one can conclude that a non-spherical outflow must be present. This is supported by the fact that non-zero differential and closure phases have been measured with AMBER. Thus, in order to obtain insights on the wind asymmetry of , we calculated 2-D axisymmetric radiative transfer models using the code described above.

Due to the limited amount of spatial frequencies and position angles sampled by our data, in our analysis we gave highest priority to the simultaneous reproduction of the observed amount of Br$\gamma$ line emission compared to the continuum, the blueshift of this emission, and the wavelength dependence of the visibility. Given the limitations of our 2-D model and the likely three-dimensional nature of 's outflow, a full analysis of the wavelength dependencies of the differential and closure phases are deferred to a future paper. Nevertheless, we believe that the modeling results presented in the following are crucial for constraining the parameter space in future studies using more sophisticated 2-D radiative transfer models such as ASTAROTH .

We analyzed both prolate and oblate density enhancements and found that, in general, both geometries require a remarkably large density contrast ($a\sim 7-15$) in order to reproduce the low visibility detected within the Br$\gamma$ line. Interestingly, the observed blueshift of the Br$\gamma\,$emission could only be reproduced by models in which the receding part of the wind, which is responsible for the major fraction of the redshifted Br$\gamma$ emission, is blocked. For a given density contrast, oblate models provided less projected separation on the sky between the redshifted and blueshifted hemispheres than prolate models, making it challenging for oblate models to produce a blueshifted Br$\gamma$ line with a similar strength and spatial extension as the observations. Therefore, the preferred geometry to fit the observations is prolate and the blocking of the redshifted hemisphere could possibly be caused either by an optically-thick disk, a putative binary companion, or by the wind being ejected by specific parts of the stellar surface.

The best fit to simultaneously reproduce the observed amount of Br$\gamma$ line emission with respect to the continuum, the line blueshift, and the wavelength dependence of the visibility was obtained by using a prolate model with a density contrast of $\sim 16$ from pole to equator (see Fig. ) oriented at $\mathrm{PA}=60 \degr$. The optimal fit to the line's blueshift was obtained by starting the blocking of the redshifted emission at a projected offset distance from the center of the star of $\Delta$RA=-0.3 mas and $\Delta$DEC=-0.3 mas (Fig. ). As Fig.  reveals, while the 2-D models can basically resemble the AMBER spectrum and visibilities across the Br$\gamma\,$line and, in the the case of the blocked model, also the observed blueshift of the emission, neither the unblocked nor the blocked model is able to reproduce the observed wavelength dependence of the differential and closure phases. Both models predict much stronger phase signals across the line, especially in the blueshifted wing of the emission line. As noted above, a consistent fit of all AMBER observables including the phases is beyond the scope of this paper and will be subject to a more detailed analysis which also includes follow-up AMBER oberservations of which will help to better constrain the 2-D model.

Compared to our spherical model presented before, the mass-loss rate of the 2-D model had to be increased to $\ensuremath{\dot{M}}=3.5 \times 10^{-5}~\ensuremath{\mathrm{M}_{\odot}{\rm yr}^{-1}}$, since a fraction of the Br$\gamma\,$emission is blocked. Lower polar density enhancements did not provide enough extension of the Br$\gamma\,$line-emitting region, while higher mass-loss rates provided a stronger Br$\gamma\,$emission line than observed with AMBER. The value of $\mathrm{PA}=60\degr$ is admittedly biased by the visibility measurement of the longest baseline, which had an orientation on the sky of $60.6\degr$. Different PAs provide a worse fit to the Br$\gamma\,$visibilities and the differential and closure phases have less amplitude than the model with $\mathrm{PA}=60 \degr$. Interferometric measurements using baselines oriented at different PAs on the sky are urgently needed to better constrain the position angle of prolate wind.

An intermediate inclination angle of $i=45 \pm 15 \degr$ is required in order to provide enough separation between the projected redshifted and blueshifted lobes and the star. In this case, the redshifted part of the line emission can be efficiently blocked and, thus, the observed blueshifted Br$\gamma\,$emission can be reproduced and yet produce sufficient spatial extension on the sky to also fit the AMBER visibilities. Models with $i > 60 \degr$ provide too much redshifted emission even if one of the hemispheres is completely blocked, while even for a polar density enhancement of 16 models with $i < 30\degr$ do not provide enough spatial extension of the Br$\gamma\,$emission region on the sky due to projection effects. If the inclination angle of the line-of-sight to is less than $30\degr$, we estimate that a very large density contrast between pole and equator of $>30$ is required, which seems unlikely.

To summarize, from the 2-D modeling of our AMBER data we conclude that the wind of has a significantly non-spherical geometry, and the viewing angle and the physical conditions must be such that a large fraction of the redshifted Br$\gamma$ emission is blocked from our view.

# Summary and conclusions

We presented the first VLTI/AMBER observations of the yellow hypergiant in low-spectral resolution (LR) mode covering the $J$, $H$, and $K$ bands and the first AMBER observation of around the Br$\gamma\,$emission line in high-spectral resolution (HR) mode with a spectral resolution of 12 000 and projected baselines between 54 and 129 m.

From the low-spectral resolution observations carried out with two linear AT arrays (baselines ranging from 15 to 96 m), we were able to extract $H$- and $K$-band visibilities at spatial scales which probe the stellar continuum emission. Differential and closure phases could be derived from the LR data for two near-infrared bands (differential phases) and all three near-infrared bands (closure phases) covered by AMBER. Within the error bars, all LR phases turn out to be zero. We corrected all visibilities for AMBER's limited field-of-view (FOV = 250 mas for AT observations) to account for the truncation of 's extended dust shell. From the visibilities we derived FWHM Gaussian diameters of 's stellar component of $1.05\pm0.07$ mas and $0.98\pm0.10$, averaged over all spectral channels in the $H$- and $K$-bands, respectively. The stellar flux contribution in the case of an infinite FOV was found to be $0.55\pm0.02$ in the $H$ band and $0.60\pm0.02$ in the $K$ band.

On the other hand, we carried out high-spectral resolution AMBER observations of centered around the Br$\gamma\,$emission line and combining three UTs. We found that the Br$\gamma\,$emission is blueshifted with $v=-25~{\rm kms}$ with respect to the systemic velocity and that the Br$\gamma\,$emitting-region of 's dense stellar wind contributes roughly 50% to the total flux at the peak wavelength of the nearly symmetric emission line. The equivalent width of the Br$\gamma\,$line is $-6.7\AA$, in agreement with recent findings . The HR visibilities of show a strong decrease at the wavelength of the Br$\gamma\,$-line emission, indicating that the region in the wind of where the Br$\gamma\,$emission arises is fully resolved by our observations and considerably more extended than the continuum-emitting region. In addition, the HR AMBER observation of revealed non-zero differential and closure phases with absolute values up to $\sim 30^\circ$ for several spectral channels covering the central and redshifted wing of the Br$\gamma\,$-emission line. This finding indicates that the shape of the Br$\gamma\,$line-emitting region is not point-symmetric.

The HR data of around the Br$\gamma\,$emission line have been analyzed by simple geometrical models as well as 1-D and 2-D gas radiative transfer models using the code CMFGEN. From simple geometrical models, we found a size of the stellar continuum-emitting region comparable to our LR observations ($d_{\star} = 1.014\pm0.01$ mas). Approximating the line-emitting region with an infinitely thin ring, we found $d_{\rm ring} = 4.18^{+0.19}_{-0.08}$ mas, i.e. $d_{\rm ring} = 4.12 \times d_{\star}$. The size of the line-emitting region is in basic agreement, although slightly larger than the value 3.3 mas found by from AMBER observations in medium-spectral resolution mode ($\lambda/\Delta\lambda=1\,500$). The difference between the two results can be explained by the differences in the details of the two modeling approaches and the different spectral resolution of the data, but probably mainly reflects the overall calibration uncertainties.

To get a first idea of the deviation from sphericity, we fitted an elliptical Gaussian ring to the 3 AMBER visibility data points at the central wavelength of the Br$\gamma\,$emission. From this fit, we find that the Br$\gamma\,$line-emitting region is elongated, with major and minor FWHM ring diameters of $4.32\pm0.44$ and $3.85\pm0.37$ mas, and an elongation towards a position angle of $36\pm11$`<!-- -->`{=html}, perfectly aligned with the elongation axis of the outer reflection seen from HST images. Follow-up observations with higher redundancy, a more robust calibration, and a better position angle coverage are indispensable to confirm this result.

Apart from simple geometrical models, we also analyzed the AMBER HR data of around the Br$\gamma\,$emission line by comparison with 1-D and 2-D gas radiative transfer models obtained with the code CMFGEN. We found that spherical models definitely fail to reproduce the observed blueshift of the Br$\gamma\,$emission while they can well reproduce the amount of line emission. Spherical wind models also fail to reproduce, at the same time, the amount of line emission and the wavelength dependence of the visibility across the emission line and, hence, the extension of the Br$\gamma\,$line-emitting region. On the other hand, we found that 2-D models with a prolate wind structure can simultaneously reproduce the amount of Br$\gamma\,$emission, the wavelength dependence of the visibility across the line and, in addition, the observed blueshift of the Br$\gamma\,$emission, but only if a substantial fraction of the redshifted Br$\gamma\,$emission is blocked. A model without any diminishing of the redshifted emission cannot explain the blueshift of the emission while *simultaneously* reproducing the amount of line emission and the spatial extension of the line-emitting region. In our best-fitting model, the outflow exhibits a strong density enhancement towards the poles with a pole-to-equator density ratio of 16, and the outflow is oriented towards a PA of 60and inclined by 45. Since the uncertainty of the position angle of the best-fitting model is approximately 20, we find that the orientation of the symmetry axis of our radiative transfer model is in rough agreement with the orientation of the elongation derived from our simple elliptical ring model ($\phi=36\pm 11^\circ$, see Sect. ). Nevertheless, it should be noted that both results might be biased towards the derived values due to the sparse position angle coverage of our measurements.

Interestingly, our best-fitting model gives a temperature of 1500 K at a radial distance of 100 mas from the central star, assuming a distance of 5 kpc for direct comparison with the results. In their best model, based on speckle-interferometric observations from 1998, the inner dust shell boundary with a dust temperature of 1000 K is located at a radial distance of $\sim70$ mas. Thus, our gas-radiative transfer model nicely fits the previously developed dust-shell model, especially if we believe that we indeed see an outward shift of the inner dust shell boundary which is driven by the noticeable increase of the stellar effective temperature within the last decades.

The signature of the differential- and closure-phase signal of detected with AMBER is totally different from what would be expected, for instance, for a rotating circumstellar disk . On the other hand, we find that the wavelength dependence of both differential and closure phases shows some qualitative similarities to measurements obtained for the Luminous Blue Variable $\eta$ Carinae . In the interpretation of the $\eta$ Car AMBER data, the phase signature was associated with enhanced mass loss of the rapidly rotating, massive primary star in polar direction, and the observer is basically looking towards the blueshifted part of the outflow. The qualitative similiarty of the phase signals between and Br$\gamma\,$motivates the assumption that also in the case of rapid rotation drives a non-spherical and polar enhanced mass loss as inferred from theoretical considerations . This would be in line with the fact that is a massive star in a highly evolved evolutionary stage and would also qualitatively fit to the observed blueshift of the Br$\gamma\,$line. However, the major difference between the differential and closure phases of $\eta$ Car and is that in the phase signal in stronger in the redshifted wing of the Br$\gamma\,$emission and not in the blueshifted wing as in the case of $\eta$ Car. This makes the modeling of the data more puzzling and probably indicates that the complexity of the overall structure of the innermost wind region goes well beyond the capabilities of a 2-D modeling as carried out in this paper.

Nevertheless, our modeling indicates that a substantial fraction of the redshifted line emission does not reach the observer. Since our modeling is restricted to an axisymmetric outflow geometry, the fact that we were not able to reproduce all AMBER observables simultaneously (in particular the phases) might indicate the following:

-   The mass outflow from is indeed intrinsically highly asymmetric. In this case, one of the hemispheres (which is pointing towards us and, thus, emitting blueshifted Br$\gamma\,$) has a much higher mass-loss rate than the other hemisphere (responsible for the redshifted light). We can only speculate on the reason for such a strong asymmetry in the outflow but, in principle, such asymmetries could be connected to large-scale asymmetries in the surface structure of the central star (e.g., hot spots) which lead to preferred directions of mass ejection . In this context, it is interesting to note that recent CO measurements detected a flux deficit in the CO shells in the south-western direction. Thus, seems to exhibit an asymmetry in the same direction as found in our study also on much larger spatial scales, indicating that this asymmetry might be of a long-term nature.

-   The wind of is axisymmetric, but the redshifted light is blocked by an optically-thick disk which would also explain the strong phase signal in the redshifted part of the Br$\gamma\,$emission. Interestingly, such a disk has to be located very close to the star ($\sim 5$ mas) in order to block the Br$\gamma\,$emission, and would have to be formed of mainly neutral+molecular gas in order to effectively absorb the redshifted light, since the presumed dust sublimation radius is located much farther out, at scales of $\sim 100$ mas. An yet undetected binary companion could also cause the presence of an equatorial disk.

-   exhibits an inflow rather than an outflow, or maybe matter infall and outflow are present at the same time , for example in terms of an equatorial infall and a polar outflow. Such an equatorial infall could lead to the formation of an optically thick disk.

Of course, the above scenarios are not mutually exclusive, and both could contribute to explain the peculiar properties seen in .

The study presented here has illustrated AMBER's great potential to probe the innermost circumstellar gas environment of evolved stars with spectral resolutions up to 12 000. New observations of covering more baselines and position angles, which are foreseen for forthcoming observing periods, will undoubtably provide a much clearer picture of the geometry of 's inner wind zone. For instance, new observations will help to support or disprove the model according to which is exhibiting a prolate wind structure with a substantial shielding of the redshifted emission. Thus, further AMBER observations probing the inner wind zone of will help to answer long-standing questions on the physical conditions and the mass-loss processes of this outstanding, rapidly evolving object and shed more light on the question of the mass-loss conditions in highly evolved massive stars.

[^1]: The low-spectral resolution data have been obtained as part of the Guaranteed Time Programme for VLTI/AMBER (program ID: 079.D-0356(B)), while the high-spectral resolution data were obtained in the context of science verification observations (program ID: 60.A-9053(D)).

[^2]: This software package is kindly provided by the Jean-Marie Mariotti Center and publicly available from http://www.jmmc.fr/data_processing_amber.htm.

[^3]: The spectra are retrievable from
    http://www.eso.org/sci/facilities/paranal/instruments/isaac/tools/spectroscopic_standards.html#Telluric.

[^4]: Apart from the information provided by ESO with the raw data itself, the ambient weather conditions during an observation can be looked up at http://archive.eso.org/asm/ambient-server.


abstract: We propose a method to study the transition to chaos in isolated quantum systems of interacting particles. It is based on the concept of delocalization of eigenstates in the energy shell, controlled by the Gaussian form of the strength function. We show that although the fluctuations of energy levels in integrable and non-integrable systems are principally different, global properties of the eigenstates may be quite similar, provided the interaction between particles exceeds some critical value. In this case the quench dynamics can be described analytically, demonstrating the universal statistical relaxation of the systems irrespectively of whether they are integrable or not.
author: L. F. Santos; F. Borgonovi; F. M. Izrailev
title: Chaos and statistical relaxation in quantum systems of interacting particles

*Introduction.* - Recent experimental progresses in the studies of various quantum systems of interacting particles (see, e.g., ) have triggered the interest in basic problems of many-body physics. One of the issues that has been widely discussed in the literature is the onset of thermalization in an isolated quantum system caused by interparticle interactions .

A prerequisite for thermalization is the statistical relaxation of the system to some kind of equilibrium and its viability has been associated with the onset of *quantum chaos*. The latter term was originally created to address specific properties of dynamical quantum systems whose classical counterparts are chaotic. Later, it was found that similar properties of spectra, eigenstates and dynamics could emerge in quantum systems without a classical limit, as well as in quantum systems with disordered potentials. Nowadays, the term quantum chaos is used in a broader context when referring to those properties, irrespectively of the existence of a classical limit.

According to studies of isolated quantum many-body systems, their eigenfunctions (EFs) in the mean-field basis spread as the interaction between particles (or quasi-particles) increases, and they may eventually become *chaotic eigenstates* . A crucial point is that due to the finite range of the interactions, only part of the *unperturbed* basis states $|n \rangle$ is directly coupled and therefore accessible to the eigenstates $|\alpha \rangle =\sum_{n} C^{\alpha}_{n} |n\rangle$. In the energy representation, this fraction constitutes the *energy shell* of the system, being either partly or fully filled by the actual eigenstates . In the first case, the eigenstates are *localized*, having a small number of non-zero elements $C^{\alpha}_{n}$. In the second case, when the energy shell is fully filled, the eigenstates can be either *sparse* or *ergodic* , both showing a very large number of principal components ($N_{pc} >> 1$) strongly fluctuating with $n$. In ergodic eigenstates the coefficients $C^{\alpha}_{n}$ become random variables following a Gaussian distribution around the "envelope" defined by the energy shell. The latter scenario is that of chaotic eigenstates and occurs when the interaction exceeds a critical value .

As argued in Ref. , the energy shell is associated with the limiting form of the *strength function* (SF) for a given interaction strength. This function is defined via the projection of unperturbed states onto the basis of perturbed (exact) eigenstates. Written in the energy representation, SF is widely used in nuclear physics and is analogous to *local density of states* in solid state physics. It has also been shown (see, e.g., ) that the shape of SF changes from Breit-Wigner (Lorentzian) to Gaussian as the interparticle interaction increases.

It should be stressed that if a quantum system has a classical limit, the shapes of both EFs and SFs in the energy representation have classical analogs. The quantum-classical correspondence of EFs and SFs has been studied in great detail for various few and many-body systems (see Ref. ). Typically, delocalization of EFs in the energy shell is directly related to the *chaotization* of the system in the classical limit, thus providing a tool to reveal the transition to quantum chaos.

The emergence of chaotic eigenstates has been related to the onset of thermalization in many-body systems, even if the latter are isolated from a heat bath . An essential point of Refs.  is that when the eigenstates become chaotic, the distribution of occupation numbers achieves a Fermi-Dirac or a Bose-Einstein form, thus allowing for the introduction of temperature. Using a two-body random matrix model, a relation between temperature and interaction strength was analytically derived . Therefore, the interparticle interaction plays the role of a heat bath when the eigenstates are chaotic. Since the components of chaotic eigenstates can be treated as random variables, the eigenstates close in energy are statistically similar. This fact has been employed to justify the agreement between the expectation values of few-body observables and the predictions from a microcanonical ensemble .

In this Letter we present a unified approach to identify the emergence of chaotic eigenstates. Two models of interacting spins-1/2 are considered: one model is completely integrable for any value of the interaction and the other is not. The different procedures used to determine the critical parameters above which the eigenstates become chaotic lead all to the same values. Our results reveal a relation between energy level statistics, the structure of EFs and SFs, and the quench dynamics. We show that statistical relaxation occurs for both models and can be described analytically with the same expression.

*The models.* We consider two models of interacting spins-1/2. One model has only nearest-neighbor couplings, which results in the complete integrability of the system. The other case has additional next-nearest-neighbor couplings, and becomes chaotic when the two coupling strengths are comparable. Assuming open boundary conditions, the Hamiltonians read as $$\begin{aligned}
&& H_1 = H_0 + \mu V_1 , \; \;\;\; H_2 = H_1 + \lambda V_2 ,
\label{model1} \\
&& H_0 = J \sum_{i=1}^{L-1}  \left(S_i^x S_{i+1}^x + S_i^y S_{i+1}^y \right) , \;\;\;\; V_1 = J \sum_{i=1}^{L-1}   S_i^z S_{i+1}^z,
\nonumber \\
&&
V_2 = \sum_{i=1}^{L-2} J\left[ \left( S_i^x S_{i+2}^x + S_i^y S_{i+2}^y
\right) + \mu S_i^z S_{i+2}^z  \right], \nonumber
\end{aligned}$$ where $\mu$ and $\lambda$ control the perturbation in Model 1 and Model 2, respectively. Here, $L$ is the number of sites, $S^{x,y,z}_i = \sigma^{x,y,z}_i/2$ are the spin operators at site $i$, with $\sigma^{x,y,z}_i$ as the Pauli matrices and $\hbar=1$. The coupling strength $J$ defines the energy scale and is set to 1. We refer to a spin "up" in the $z$ direction as an excitation, assuming the presence of a magnetic field.

In Model 1, the term $H_0$ determines the (mean field) basis in which the total Hamiltonian $H_1$ is presented. This term moves the excitations through the chain and can be mapped onto a system of noninteracting spinless fermions  or of hard-core bosons , being therefore integrable. The system remains integrable even when the Ising interaction $V_1$ is added, no matter how large the anisotropy parameter $\mu$ is. The total Hamiltonian $H_1$ is known as the XXZ Hamiltonian and can be solved exactly via the Bethe Ansatz .

The integrable XXZ Hamiltonian determines the mean field basis for Model 2. Thus, $V_2$ is treated as the "residual interaction" responsible for the onset of chaos. The parameter $\lambda$ refers to the ratio between the next-nearest-neighbor and the nearest-neighbor exchange.

Depending on the parameters of the Hamiltonians (), different symmetries are identified . For the sake of generality, we avoid them by restricting our analysis to a subspace with $L/3$ up-spins and $\mu \neq 1$. The only remaining symmetry is parity. We take it into account by studying only even states, which leads to subspaces of dimension $N \sim (1/2) L!/[(L/3)!(L-L/3)!]$. All data are given for $L=15$ and $\mu=0.5$ for Model 2.

*Spectrum statistics.* - According to the common lore, we analyze first the level spacing distribution $P(s)$ for both models, numerically obtained for different values of the control parameters $\mu$ and $\lambda$. For Model 1, $P(s)$ is close to the Poisson distribution for any value of $\mu$. For Model 2, the transition of $P(s)$ from Poisson to Wigner-Dyson as $\lambda$ increases is shown in Fig. . The standard approach of fitting $P(s)$ with the Brody distribution allows us to extract the repulsion parameter $\beta$ characterizing the transition between the two distributions. From Fig. , the transition for Model 2 with $L=15$ occurs approximately at $\lambda \approx 0.5$.

*Emergence of chaotic eigenstates.* - Much more information is contained in the structure of the eigenstates. Our data show that as the strength of the perturbations $V_1$ and $V_2$ increases, the eigenstates of *both* integrable and non-integrable models undergo a transition from localized to chaotic-like ones. Typical examples of the amplitudes $C_n^\alpha$ of such eigenstates with energy $E_{\alpha}$ from the center of the energy band are shown in Fig. . Here, the eigenstates are given as a function of the unperturbed energy $\varepsilon_n$ rather than in the basis representation, following the one-to-one correspondence between each unperturbed state $|n \rangle$ and its energy $\varepsilon_n$. The figure shows how the eigenstates spread as the interaction increases.

In order to find the critical parameters $\mu_{cr}$ and $\lambda_{cr}$ above which the perturbation is strong and the eigenstates are extended (in energy shell), different approaches may be employed. We start by analyzing the matrix elements of $H_1$ and $H_2$. It is important to take into account that in each line $n$ of the Hamiltonians, the perturbation couples directly only $M_n$ unperturbed states . We have numerically found that at the center of the energy band, $M_n \approx N/2$ and $M_n \approx N$ for Models 1 and 2, respectively. Thus, the Hamiltonian matrix of the integrable model has more vanishing elements $H_{nm}$, being more sparse than the matrix of Model 2.

To determine the critical perturbation, we compare the average value of the coupling strength, $v_n = \sum_{m\neq n} |H_{nm}|/M_{n}$, of each line $n$ with the mean level spacing $d_n$ between directly coupled states. The mean level spacing can be estimated as $d_n=[\varepsilon_n^{max}-\varepsilon_n^{min}]/M_{n}$, where $\varepsilon_n^{max}$ ($\varepsilon_n^{min}$) is the unperturbed energy corresponding to the largest (smallest) $m$ where $H_{nm}\neq0$. Our results show that for $\mu > \mu_{cr} \approx 0.5$ and $\lambda > \lambda_{cr} \approx 0.5$ the ratio $v_n / d_n$ becomes larger than 1 and the perturbation is considered to be strong. Notice that for Model 2, the obtained value of $\lambda_{cr}$ corresponds to the onset of the Wigner-Dyson statistics, as independently found from the level spacing distribution. This is remarkable if we take into account that no diagonalization was necessary to derive the above estimates.

*Strength function: From Breit-Wigner to Gaussian.* - Another way to obtain the critical values $\mu_{cr}$ and $\lambda_{cr}$ relies on the shape of SF. The latter corresponds to the dependence of $w_n^{\alpha}=|C^{\alpha}_n|^2$ on the exact energies $E_{\alpha}$ for each fixed unperturbed energy $\varepsilon_n$. It contains information about the energies $E_{\alpha}$ that become accessible to an initial state $|n \rangle$ when the perturbation is turned on. Clearly, SF is related also to the structure of EFs, since the latter is derived from the same $w_n^{\alpha}$, but now as a function of the unperturbed energies $\varepsilon_n$.

In quantum many-body systems, the form of SF typically changes from Breit-Wigner to Gaussian as the inter-particle interactions increase . This transition occurs when the half-width of the Breit-Wigner distribution becomes comparable to the width of the energy shell. In this case, as we show next, there emerge chaotic eigenstates filling the whole available energy shell.

The energy shell corresponds to the distribution of states obtained from a matrix filled only with the off-diagonal elements of the perturbation. It is associated with the maximal SF, that is the SF that arises when the diagonal part of the Hamiltonian can be neglected. We verified that the energy shell coincides with a Gaussian of variance $\sigma^2$ given by the second moment of the off-diagonal elements of the matrix Hamiltonian, $\sigma^2=\sum_{m\neq n} |H_{nm}|^2$ . Note that no diagonalization is required to derive this expression.

Our numerical data confirm that the transition from Breit-Wigner to Gaussian occurs for the same critical values obtained above, $\mu_{cr} \, , \lambda_{cr} \approx 0.5$, as indicated in Fig. . In the figure, the envelopes of the SFs were obtained by smoothing the dependence of $w_n^{\alpha}$ on $E_{\alpha}$ for fixed unperturbed energies $\varepsilon_{n}$ with $n \approx N/2$. The fit to either a Breit-Wigner or a Gaussian form was done with high accuracy, allowing us to discriminate between the two functions. It is noteworthy the excellent agreement between the Gaussian fit and the Gaussian obtained simply from the off-diagonal elements of the Hamiltonians.

*Structure of eigenstates in energy shell.* - The eigenstates may be localized, sparse or ergodically extended in the energy shell. The data in Fig.  demonstrate that for a sufficiently strong perturbation the eigenstates undergo a transition from strongly localized to extended states, somehow filling the energy shell. The transition to chaotic eigenstates occurs again at the same critical parameters $\mu_{cr} \, , \lambda_{cr} \approx 0.5$. These results confirm the predictions made on the basis of both, the estimate of $v_n/d_n$ and the Gaussian form of the strength functions.

One notices that above the critical value, the eigenstates of Model 2 appear to fill the energy shell better than those of Model 1 (cf. bottom panels of Fig. ). Indeed, for Model 1, the eigenstates are not completely extended in the energy shell even for strong interaction, $\mu =1.5$, although they do fill a large part of it. At the same token, a close inspection of the level of delocalization of individual EFs and SFs has revealed differences between the two models. Overall, delocalization measures, such as the inverse participation ratio or Shannon entropy, show larger fluctuations for Model 1 than for Model 2 . This agrees with recent results obtained for bosonic and fermionic systems .

*Statistical relaxation.* - The knowledge of the shape of SFs allows one to describe the quench dynamics in a system of interacting particles, provided the eigenstates are chaotic and delocalized in the energy shell . By quench dynamics we mean the time evolution of initial states corresponding to unperturbed basis vectors which takes place once the interaction is turned on. To see how the relaxation occurs, we study the time dependence of the Shannon entropy $S$ in the unperturbed basis. For an initial state $|k\rangle$, it is defined as $$S_k(t)=-\sum_{n=1}^N  \Omega_n^{(k)} (t) \ln \Omega_n^{(k)} (t)
\label{entropy}$$ with $\Omega_n^{(k)} (t) = |\sum_{\alpha} C^{\alpha*}_k C^{\alpha}_n e^{-iE_{\alpha}t}|^2$.

An analytical expression for $S{(t)}$ has been derived with the use of a two-body random matrix model , $$S_k(t) = -W_k(t) \ln W_k(t) -
[1-W_k(t)] \ln \left( \frac{1-W_k(t)}{N_{pc}} \right).
\label{analyt}$$ Here $W_{k}(t)$ is the probability for the system to stay in the initial state $|k \rangle$ and $N_{pc}$ is the average number of directly coupled states. We obtain $N_{pc}$ numerically according to $N_{pc}=\langle e^S \rangle$, where the average $\langle. \rangle$ is performed over a long time after the entropy saturates. As for $W_{k}(t)$, it has been shown that for a Gaussian SF, it decays as $W_{k}(t)= \exp (-\sigma^2 t^2)$, with $\sigma^2$ as previously defined.

Figure shows numerical data for the relaxation process of both models. To reduce fluctuations, we average over 5 initial basis states excited in a narrow energy range in the middle of the spectrum. Initially, the entropy grows quadratically, as given by perturbation theory. Afterwards, a clear linear growth is observed before $S(t)$ reaches relaxation. With high accuracy the linear behavior of $S(t)$ is described by the simple relation , $$S_k(t) \approx \sigma_k t \ln M_k.
\label{linear}$$ Note that Eq. () depends only on the elements of the Hamiltonian: $\sigma_k^2=\sum_{m\neq k} |H_{km}|^2$ and $M_k$ is the connectivity, that is the number of directly coupled unperturbed states in the $k$-th line of the Hamiltonian matrix. As seen in Fig. , the analytical expressions () and () give a correct description of the increase of the entropy for *both* models in the regime corresponding to the onset of chaotic-like eigenstates delocalized in the energy shell. The same relation () was found to emerge also for an integrable model of interacting bosons .

*Conclusion.* - We have studied spectrum statistics, structures of eigenstates and strength functions for two models of interacting spins, connecting the results with the onset of chaotic-like eigenstates and the statistical relaxation of the quench dynamics. The key point of our approach is the existence of an energy shell of finite range, inside which the eigenfunctions can be either localized or extended. We have shown that the critical parameters above which the eigenstates become chaotic can be equally found by simply studying the elements of the Hamiltonian or by analyzing the shape of the strength function. The latter provides us with the form of the energy shell, thus allowing one to clearly define the notion of delocalized eigenstates in the energy shell.

By studying the time dependence of the Shannon entropy, we have shown that numerical data are in full agreement with the analytical predictions of the quench dynamics, provided the eigenstates are chaotic-like. Our approach is very general and expected to apply to different systems of interacting particles, such as those currently under theoretical and experimental investigation.

*Acknowledgments.* F.M.I. acknowledges support from CONACyT grant N-161665 and thanks Yeshiva University for the hospitality during his stay in the fall 2011.