# Large-scale environments of binary AGB stars probed by Herschel

A. Mayer, A. Jorissen, C. Paladini, F. Kerschbaum, D. Pourbaix, C. Siopis, R. Ottensamer, M. Mečina, N. Cox, M. Groenewegen, et al.

► To cite this version:

A. Mayer, A. Jorissen, C. Paladini, F. Kerschbaum, D. Pourbaix, et al. Large-scale environments of binary AGB stars probed by Herschel. Astronomy and Astrophysics - A&A, 2014, 570, pp.A113. 10.1051/0004-6361/201424465. hal-03526443

HAL Id: hal-03526443

https://hal.science/hal-03526443

Submitted on 15 Jan 2022

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
# Large-scale environments of binary AGB stars
## probed by Herschel

## II. Two companions interacting with the wind of $\pi^1$ Gruis ⋆, ⋈

A. Mayer<fnref n="1" />, A. Jorissen<fnref n="2" />, C. Paladini<fnref n="2" />, F. Kerschbaum<fnref n="1" />, D. Pourbaix<fnref n="2" />, ⋈, C. Siopis<fnref n="2" />, R. Ottensamer<fnref n="1" />, M. Mečina<fnref n="1" />, N. L. J. Cox<fnref n="3" />, M. A. T. Groenewegen<fnref n="4" />, D. Klotz<fnref n="1" />, G. Sadowski<fnref n="2" />, A. Spang<fnref n="5" />, P. Cruzalèbes<fnref n="5" />, and C. Waelkens<fnref n="3" />

Received 24 June 2014 / Accepted 14 August 2014

### ABSTRACT

**Context.** The Mass loss of Evolved StarS (MESS) sample observed with PACS on board the Herschel Space Observatory revealed that several asymptotic giant branch (AGB) stars are surrounded by an asymmetric circumstellar envelope (CSE) whose morphology is most likely caused by the interaction with a stellar companion. The evolution of AGB stars in binary systems plays a crucial role in understanding the formation of asymmetries in planetary nebulae (PNe), but at present, only a handful of cases are known where the interaction of a companion with the stellar AGB wind is observed.

**Aims.** We probe the environment of the very evolved AGB star $\pi^1$ Gruis on large and small scales to identify the triggers of the observed asymmetries.

**Methods.** Observations made with Herschel/PACS at 70 µm and 160 µm picture the large-scale environment of $\pi^1$ Gru. The close surroundings of the star are probed by interferometric observations from the VLTI/AMBER archive. An analysis of the proper motion data of HIPPARCOS and Tycho-2 together with the HIPPARCOS Intermediate Astrometric Data help identify the possible cause for the observed asymmetry.

**Results.** The Herschel/PACS images of $\pi^1$ Gru show an elliptical CSE whose properties agree with those derived from a CO map published in the literature. In addition, an arc east of the star is visible at a distance of 38" from the primary. This arc is most likely part of an Archimedean spiral caused by an already known GOV companion that is orbiting the primary at a projected distance of 460 au with a period of more than 6200 yr. However, the presence of the elliptical CSE, proper motion variations, and geometric modelling of the VLTI/AMBER observations point towards a third component in the system, with an orbital period shorter than 10 yr, orbiting much closer to the primary than the GOV star.

**Key words.** stars: AGB and post-AGB – binaries: general – circumstellar matter – stars: winds, outflows – stars: individual: $\pi^1$ Gruis – infrared: stars

## 1. Introduction

The evolution of low- and intermediate-mass stars ends with an ascent of the asymptotic giant branch (AGB) in the Hertzsprung-Russell diagram. This phase involves an increase of mass loss that strips off the envelope through a slow and dust-enriched wind ($v_w = 5–20 \text{ km s}^{-1}$) blown into the interstellar medium (ISM), where it ranks among the dominant contributors of heavy elements in the Galaxy. Finally, at the end of

the AGB phase, the remaining envelope as a whole is ejected. The hot remnant stellar core ionizes the ejecta, forming what is known as a planetary nebula (PN). PNe show a manifold of morphological diversity, including highly asymmetric and bipolar forms that can only be adequately described by a binary star model (e.g. Nordhaus et al. 2007; De Marco et al. 2008; Miszalski et al. 2009a,b).

Paczyński (1971), Livio & Soker (1988), and Theuns & Jorissen (1993) have theoretically shown that already in the AGB phase the stellar winds must be heavily distorted in binary systems depending on the size of the system and the evolutionary type of the companion. For binary systems with small separations, the primary AGB star fills the Roche lobe and transfers mass onto the companion (Paczyński 1965, 1971).

In detached systems, the stellar AGB wind fills the Roche lobe, and up to half of the material is accreted by the companion (named wind Roche lobe overflow – WRLOF, Mohamed & Podsiadlowski 2011; Abate et al. 2013). Jets and bipolar

<footnote n="1">¹ University of Vienna, Department of Astrophysics, Sternwartestraße 77, 1180 Wien, Austria
e-mail: a.mayer@univie.ac.at</footnote>

<footnote n="2">² Institut d'Astronomie et d'Astrophysique, Université Libre de Bruxelles, CP. 226, Boulevard du Triomphe, 1050 Brussels, Belgium</footnote>

<footnote n="3">³ Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan, 200D, 3001 Leuven, Belgium</footnote>

<footnote n="4">⁴ Koninklijke Sterrenwacht van België, Ringlaan 3, 1180 Brussels, Belgium</footnote>

<footnote n="5">⁵ Laboratoire Lagrange, UMR 7293, Université de Nice-Sophia Antipolis, CNRS, Observatoire de la Côte d'Azur, Bd de l'Observatoire, CS 34229, 06304 Nice Cedex 4, France</footnote>

<footnote>* Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA. This paper makes use of data from ESO programmes 076.D-0624, 077.D-0620, 078.D-0122, 080.D-0076, 187.D-0924.</footnote>

<footnote>** Maps of Fig. 4 are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/570/A113</footnote>

<footnote>*** Senior Research Associate, F.R.S.-FNRS, Belgium.</footnote>
outflows have also been observed for some of these systems, e.g., o Cet (Meaburn et al. 2009), R Aqr (Wallerstein & Greenstein 1980; Kafatos et al. 1989), and V Hya (Hirano et al. 2004), where the accretion disc of the companion is fed by a strong wind from the AGB primary (e.g. Morris 1987; Soker & Rappaport 2000; Huggins 2007). But in general, the companion affects the circumstellar envelopes (CSEs) of the AGB star in two ways. First, the material that is transferred via the WRLOF is focused by the gravitational potential of the companion and forms a density wake that trails the orbital motion of the companion. The result is a wind pattern shaped as an Archimedean spiral, as predicted by hydrodynamic simulations in Theuns & Jorissen (1993), Mastroemos & Morris (1998, 1999), and Kim & Taam (2012b). Second, the presence of the companion also manifests itself via the gravitational force that it exerts on the primary, causing it to move around the centre of mass of the binary system (Soker 1994; Kim & Taam 2012c,a). Furthermore, Kim & Taam (2012c) and Kim et al. (2013) recently demonstrated that the combination of the two effects leads to a spiral wind pattern exhibiting knots where the two structures intersect.

Observationally, spiral patterns were found around a small number of AGB or proto-PN objects: AFGL 3068 (Mauron & Huggins 2006), CIT 6 (Dinh-V.-Trung & Lim 2009; Kim et al. 2013), o Cet (Mayer et al. 2011), R Scl (Maercker et al. 2012), and W Aql (Mayer et al. 2013), all of which are wide binary systems with an orbital separation in the range of ≈50–160 au. Recently, Mauron et al. (2013) found that 50% of a sample of 22 AGB stars have elliptical emission, which the authors attributed to binaries whose envelopes are flattened by a companion.

This work continues our study of large-scale environments of binary AGB stars from the Herschel Mass loss of Evolved StarS sample (MESS; Groenewegen et al. 2011). Contrary to Mayer et al. (2013, hereafter Paper I), which concentrated exclusively on Hersche/PACS observations of the large-scale structures (around R Aqr and W Aql), here we also explore the close surroundings of the star using HIPPARCOS Intermediate Astrometry Data (IAD; van Leeuwen & Evans 1998) as well as unpublished archive observations obtained with the Very Large Telescope Interferometer Astronomical Multi-BEam combineR (VLTI/AMBER, Petrov et al. 2007). We do this by analysing the structures around the binary AGB star π<fnref n="1" /> Gru at angular scales from 0.′02 to 60″. Section 2 discusses the fundamental properties of π<fnref n="1" /> Gru. In Sect. 3, the observations of Herschel/PACS and VLTI/AMBER are described, with their results presented in Sect. 4. The different interaction scenarios that can produce asymmetries in the extended environment of the star on both small and large scales are discussed in Sect. 5.

## 2. General properties of π<fnref n="1" /> Gruis

The S5.7 star π<fnref n="1" /> Gru (HIP 110478) is an SRb variable, and because of its proximity (d = 163 pc; van Leeuwen 2007) one of the brightest and best studied intrinsic<fnref n="1" /> S stars (Keenan 1954). The intrinsic nature of the S star π<fnref n="1" /> Gru is defined from the presence of spectral lines of the element Tc (Jorissen et al. 1993). The pulsation period of the star was initially derived to be ≈150 days by Eggen (1975) and used in various publications since then. A new derivation of the period based on the light curve provided by the All Sky Automated Survey (ASAS) Photometric V-band Catalogue (Pojmanski & Maciejewski 2005) revealed, however, π<fnref n="1" /> Gru varies with a pulsation period of 195 days (see Fig. 1).

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** ASAS-3 light curve of π<fnref n="1" /> Gru with an adopted pulsation period of 195 days covering 14 cycles.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Evolutionary tracks from STAREVOL for stars with initial masses 1.5 M☉ (black line), 1.7 M☉ (red dotted line) and 2.0 M☉ (green dashed line) from the pre-main-sequence to the end of the AGB. π<fnref n="1" /> Gru is represented by the large blue circle.</figcaption></figure>

The Hertzsprung-Russell diagram for HIPPARCOS S stars from Van Eck et al. (1998) point to the very evolved nature of this star, close to the tip of the AGB. With values of T<sub>eff</sub> = 3100 K and log L/L<sub>☉</sub> = 3.86 (Van Eck et al. 1998), the location of π<fnref n="1" /> Gru in the Hertzsprung-Russell diagram can be compared with evolutionary tracks (Fig. 2) computed from the STAREVOL code (Siess 2006; Siess & Arnould 2008) with a metallicity Z = 0.02. It appears that π<fnref n="1" /> Gru falls on the track of a star of initial mass 2.0 M☉, but by the time that star has reached log L/L<sub>☉</sub> = 3.86, its mass has dropped to about 1.5 M☉.

From CO(1–0) observations, a present-day mass loss rate of $\dot{M} = 2.73 \times 10^{-6} M_{\odot} \text{yr}^{-1}$ (Winters et al. 2003), an expansion velocity of 14.5 km s<sup>-1</sup> (Guandalini & Busso 2008), and a gas-to-dust ratio of 380 (Groenewegen & de Jong 1998) is derived.

<footnote n="1">¹ About intrinsic S stars, see Van Eck & Jorissen (1999).</footnote>
<table><caption>**Table 1.** Journal of the AMBER observations of $\pi^1$ Gru.</caption><thead><tr><th>Date and UT time</th><th>Config.</th><th>Baselines [m]</th><th>PA [deg]</th><th>Seeing ["]</th><th>Airmass</th></tr></thead><tbody><tr><td>09 Oct. 2007 T23:51:44.10</td><td>E0-G0-H0</td><td>16–31–47</td><td>230</td><td>0.68–0.71</td><td>1.18</td></tr><tr><td>10 Oct. 2007 T00:53:41.50</td><td>"</td><td>16–32–48</td><td>242</td><td>0.70–0.78</td><td>1.09</td></tr><tr><td>10 Oct. 2007 T01:05:39.23</td><td>"</td><td>16–32–48</td><td>244</td><td>0.67–0.63</td><td>1.08</td></tr><tr><td>10 Oct. 2007 T01:50:37.99</td><td>"</td><td>16–32–48</td><td>252</td><td>1.05–1.10</td><td>1.07</td></tr><tr><td>10 Oct. 2007 T03:05:37.54</td><td>"</td><td>15–30–45</td><td>265</td><td>0.67–0.68</td><td>1.11</td></tr><tr><td>10 Oct. 2007 T03:54:47.47</td><td>"</td><td>14–28–42</td><td>274</td><td>0.56–0.47</td><td>1.18</td></tr><tr><td>10 Oct. 2007 T04:39:09.68</td><td>"</td><td>13–26–39</td><td>283</td><td>0.57–0.56</td><td>1.30</td></tr></tbody></table>

$\pi^1$ Gru is known to have a faint G0V companion with an apparent visual magnitude of 10.4 (Feast 1953; Ake & Johnson 1992). From the HIPPARCOS parallax (6.13 ± 0.76 mas or 163 pc; van Leeuwen 2007), the distance modulus is 6.1, yielding an absolute visual magnitude of 4.3 for the G0V star, in accordance with its spectral classification. The companion must thus be physically associated with the S star, but the orbital period is quite long, since the relative position did not change significantly over the past century, according to the list of relative positions collected by the Washington Double Star Observations catalogue, and kindly communicated to us by B. Mason (see Appendix A.2 and Table A.2). Assuming that the observed angular separation (≈2.′′8) corresponds to the semi-major axis and adopting 2.5 $M_\odot$ as the total mass of the system ($\pi^1$ Gru+G0V companion), the system parallax implies an orbital separation of the order of 460 au and an orbital period of about 6200 yr.

## 3. Observations

### 3.1. Herschel/PACS

The observations presented here are part of the MESS guaranteed time key programme (Groenewegen et al. 2011) for the Herschel Space Observatory (Pilbratt et al. 2010) using the Photodetector Array Camera and Spectrometer (PACS, Poglitsch et al. 2010) and the Spectral and Photometric Image Receiver (SPIRE, Griffin et al. 2010) on board the spacecraft. $\pi^1$ Gru was observed on May 21 2010. For the following analysis, we exclusively used the PACS data since the instrument offers a resolution that best suits our purpose. Observations of $\pi^1$ Gru were obtained at 70 and 160 $\mu$m at a FWHM of 5.′′6 and 12.′′, respectively. The adopted data processing and image reconstruction for $\pi^1$ Gru was made in the same way as for the data presented in Paper I, following Groenewegen et al. (2011) and Roussel (2013). We oversampled the reconstructed images by a factor 3.2 to achieve a pixel size of 1′′ in the blue and 2′′ in the red band.

An overview of the MESS objects is given by Cox et al. (2012) and detailed studies of individual objects are presented in Ladjal et al. (2010a), Kerschbaum et al. (2010), Jorissen et al. (2011), Mayer et al. (2011, 2013), Decin et al. (2011, 2012), van Hoof et al. (2013), and Mečina et al. (2014).

### 3.2. VLTI/AMBER

$\pi^1$ Gru was observed with VLTI/AMBER and VLTI/MIDI in the framework of ESO programmes 076.D-0624, 077.D-0620, 078.D-0122, 080.D-0076, and 187.D-0924. The detailed description of the MIDI data reduction and modelling is given in Sacuto et al. (2008) and Paladini et al. (in prep.). The MIDI observations do not deviate from spherical symmetry. These

data, however, sample only the low spatial frequencies, and it is known that for AGBs, asymmetric structures are usually detected at high spatial frequencies. The only information that we can extract from the MIDI observations is the overall size of the envelope. For this reason the MIDI data are not discussed here.

We retrieved seven VLTI/AMBER observations from the ESO archive. The data were recorded in low-resolution mode ($R = 30$) on the night of October 9, 2007 and cover the J, H, and K bands. The log of these observations is presented in Table 1, while Fig. 3 shows the corresponding uv-plane coverage.

We reduced the data with amdlib v.3.0.8 (Tatulli et al. 2007; Chelli et al. 2009) using the K3III star $\lambda$ Gru as calibrator (Bordé et al. 2002; Cruzalèbes et al. 2010, with a stellar diameter of 2.62 ± 0.03 mas). The data analysis is limited to the H and K bands because reliability for the wavelengths shorter than 1.46 $\mu$m is not guaranteed by the current pipeline version.

The medium- and long-baseline visibilities sample the second and even third lobe, meaning that the star is fully resolved. There is evidence for deviation from centro-symmetry (or rather left/right symmetry for this triplet of aligned baselines), as judged from the non-zero closure phase. We return to this in Sect. 5.3.2.

## 4. Results

Figure 4 depicts the images obtained with Herschel/PACS at 70 $\mu$m (blue band) and 160 $\mu$m (red band)$^2$. Since the blue band

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** uv-plane coverage of the VLTI/AMBER data of $\pi^1$ Gru.</figcaption></figure>

<footnote>$^2$ Both maps are available as FITS files from CDS/VizieR at http://cdsweb.u-strasbg.fr/</footnote>
<figure><img src="image_5.png" /><figcaption>Fig. 4. Deconvolved Herschel/PACS images of $\pi^1$ Gru at 70 µm (left panel) and 160 µm (right panel).</figcaption></figure>

offers a better spatial resolution, all following discussions are based on this image unless stated otherwise. The PACS 70 µm image is dominated by two features, an elliptical CSE and an arc east of the star, which are described in the remainder of this section.

The elliptical CSE is oriented east-west with its major axis at PA ≈ 105°. The total size of the emission is ≈72" × 60" [11750 × 9790 au]. A confirmation of the CSE size is provided by the 870 µm image of $\pi^1$ Gru obtained with the APEX bolometer LABoCa (Ladjal et al. 2010b). The authors found an elongation of the CSE in east-west direction with a total size of the structure of about 60" × 40".

The size of the CSE obtained from the PACS 70 µm image is much smaller than that inferred by Young et al. (1993) from the IRAS 60 µm dust emission (4.9'). A similar discrepancy for the CSE sizes obtained from IRAS and Herschel data was found for the targets X Her and TX Psc analysed previously (Jorissen et al. 2011), most likely owing to the size of the IRAS PSF (FWHM = 1.6 at 60 µm), which causes difficulties in probing asymmetries of the order of 1' as found in the Herschel data.

The second main feature visible in the far-IR emission of $\pi^1$ Gru is an arc east of the central star. It emerges 38" away from the stellar system, in the direction of the major axis (≈105°). This finding diminishes the probability that the elliptical emission is caused by the interaction with the ISM. The observed arc is curved towards the north and extends in that direction for almost 25" and 16" to the east. On the PACS 160 µm image, the arc suffers from the low resolution but remains recognisable as a clump. We further note that in Fig. A.2 of Ladjal et al. (2010b), a spike is visible on the eastern side, that might reflect the arc seen on the 70 µm PACS image.

# 5. Discussion

## 5.1. Origin of the far-infrared dust arc

Arcs or arms around AGB stars were found recently for CIT 6 (Dinh-V.-Trung & Lim 2009; Kim et al. 2013), TX Cam (Castro-Carrizo et al. 2010), o Cet (Mayer et al. 2011), and W Aql (Mayer et al. 2013). They are most likely part of Archimedean spirals caused by a combination of (i) the accretion wake of the companion when it orbits the mass-losing primary and the wind pushes the wake radially outwards (Theuns & Jorissen 1993; Mastrodemos & Morris 1998, 1999); and (ii) the motion of the primary around the centre of mass caused by the

gravitational drag of the companion (Soker 1994; Kim & Taam 2012c,a). Spectacular examples of complete Archimedean spirals have been observed around AFGL 3068 (Mauron & Huggins 2006) and R Scl (Maercker et al. 2012).

It is not clear, however, whether the far-infrared arc seen around $\pi^1$ Gru has anything to do with a spiral pattern produced by the orbital motion of a companion. If it does, the properties of the arc have to be consistent with the properties of the $\pi^1$ Gru+G0V binary system. From the orbital period $P \approx 6200$ yr, the parallax $\sigma = 6.13$ mas, and the wind velocity $v_w = 14.5$ km s$^{-1}$, one derives an arm separation of

$$ \rho = v_w P \sigma = 116'' \quad (1) $$

This is more than twice as large as the separation of 52" between the central part of the arc and $\pi^1$ Gru as seen on the PACS 70 µm image. If the arc is part of a spiral, it only represents a part of the first spiral twist.

In the upper panel of Fig. 5 we plot a spiral with $\rho = 116''$ over the PACS image that seems to match the observed arc well. However, this assumption has to be made with care since the curvature of the arm allows many solutions. A unique solution is obtained when the start of the spiral coincides with the current position of the G0V star. The lower panel of Fig. 5 displays the same scenario in a polar-radial diagram, which facilitates the illustration of the Archimedean spiral (dashed line). The current position of an anticlockwise orbiting companion is given by

$$ \Phi_{\text{comp}} = \Phi_{\text{arc}} + \frac{2\pi R}{\rho}, \quad (2) $$

where $R$ is the distance from a given part of the spiral to its origin, and $\Phi_{\text{arc}}$ the PA of that given part<fnref n="3" />. Assuming $R = 52'' - 2''8$, $\Phi_{\text{arc}} = 85^\circ$, and $\rho = 116''$, the PA of the companion is ≈240°, which has to be compared with the observed PA $\Phi_{\text{comp}} = 203^\circ$ (see Table A.2).

A spiral that matches both the slope of the arc and the position of the G0V star has an arm spacing $\rho = 168''$. This value is about 45% higher than that derived from the wind velocity of the primary, however, and the orbital period and this difference may originate from the uncertainties on these values.

<footnote n="3">³ This is a simplification since the origin of the Archimedean spiral is the primary. For a binary system this wide, it is assumed, however, that the accretion wake of the companion causes the spiral, not the reflex motion of the primary. The spiral thus follows the involute of the circular orbit. Since the outcome is almost identical (Kim & Taam 2012c), we keep the description of an Archimedean spiral.</footnote>
<figure><img src="image_6.png" /></figure>

**Fig. 5.** *Upper panel:* contour plot of the Herschel/PACS 70 µm image of π¹ Gru over-plotted with an Archimedean spiral. The spiral spacing of 116'' was derived from the wind velocity and the orbital period of the known GOV companion. *Lower panel:* polar-radial intensity profile of the same image. The position angle is measured north over east and the colour code is given in Jy arcsec⁻². The dashed line shows the same Archimedean spiral as in the upper panel. The (yellow) star illustrates the current position of the GOV companion.

Maercker et al. (2012) remarked that the arm separation of the spiral around R Scl changed significantly during the past 1800 yr. The outer (older) part of the spiral shows a larger separation than the inner (younger) part. For the authors, this indicates a modulation of the mass-loss rate by a factor of 30 caused by a thermal pulse. At the beginning of that phase, the wind velocity increased by 40% and subsequently declined to the present-day value within 1200 yr. This is measurable in the spiral-arm separation, which varied by this value.

Such a thermal pulse presumably also occurred in π¹ Gru given its location at the tip of the AGB. Furthermore, Knapp et al. (1999) suggested that the mass-loss rate of π¹ Gru has increased in the past 1000 years to explain the presence of the CO disc.

The second uncertainty on the arm separation ρ stems from the assumption that the observed orbital separation of the π¹ Gru+GOV system (≈2''8) is de-projected, meaning that the binary orbit is seen face-on. For instance, to obtain an arm spacing of 168'' with a constant wind velocity of 14.5 km s⁻¹, an orbital period of 9000 yr is required, which implies an inclination i = 46° of the orbit w.r.t. the plane of the sky.

The appearance of inclined spirals was studied in the hydrodynamic simulations of Mastrodemos & Morris (1999), Mohamed & Podsiadlowski (2011), and Kim & Taam (2012c).

The authors found that the spiral shape is preserved up to an inclination angle of ≈70° and then changes its appearance to broken concentric shells. Hence, a spiral pattern inclined by 46° would still be recognisable as such.

## 5.2. Origin of the elliptical emission

A puzzling fact is the presence of the elliptical emission together with the spiral arc. While the arc can be explained by the GOV companion interacting with the primary's wind (see Sect. 5.1), this cannot be the case for the elliptical emission.

An inclined disc structure around π¹ Gru was proposed by several CO and SiO line studies in the past 20 years (e.g. Sahai 1992; Winters et al. 2003; Knapp et al. 1999; Chiu et al. 2006). In all of the observations, the CO(2–1) profiles show an asymmetric, double-peaked structure and extended emission wings. On a corresponding map (Chiu et al. 2006), the envelope is elongated in the east-west direction with a size of ≈40'' [6530 au] and a velocity gradient in the north-south direction.

Knapp et al. (1999) interpreted these observations in terms of an expanding disc with a radius of 3340 au [20''] inclined by 35° to the plane of the sky with the northern part of the disc tilted away from the observer. The same inclination was found from the axis ratio of the PACS far-IR emission along with the same orientation of the projected major axis, which assumes that both structures are identical, but visible on different scales. A similar conclusion was drawn by Ladjal et al. (2010b) from the APEX/LABoCa observations where the authors mention an inclination of the structure of about 70° (a probable misprint for 48° = arccos[40/60]).

According to Knapp et al. (1999), the disc is produced by a constant mass-loss rate of 1.2 × 10⁻⁶ M☉ yr⁻¹ expanding in the plane of the disc with a velocity of 13 ± 2 km s⁻¹, in agreement with the velocity of the SiO(6–5) line. The expansion velocity derived from the CO lines increases steadily to 18 km s⁻¹ towards the pole. In addition, a fast molecular wind (with velocities of at least 70 km s⁻¹) is observed, which is most likely a continuation of the velocity increase towards the poles. This is in conflict with Sahai (1992), who interpreted the spatially-separated horn features in the CO lines as arising from a bipolar flow perpendicular to the disc. In the most recent study of the CO emission, Chiu et al. (2006) adhered to the Knapp et al. (1999) model and found, moreover, that the high-velocity outflow is oriented along PA = 30–210°.

A natural explanation for the elliptical emission is that the wind of the AGB star does not propagate in a spherically symmetric fashion but is focused towards a plane. Three known mechanisms can account for this: (i) a fast differential internal rotation creates a gradient in the wind velocity between the equator of the star and its poles. However, only one AGB star is known to show a fast rotation (V Hya; Barnbaum et al. 1995), and generally slow rotation rates are found among white dwarf stars (Kawaler et al. 1999); (ii) a bipolar magnetic field causes the wind to become denser at the magnetic equator (Matt et al. 2000). Surface magnetic fields have been recently reported for some AGB stars (Vlemmings et al. 2011; Leal-Ferreira et al. 2013; Lèbre et al. 2014), but there is no observational evidence for stellar magnetic fields to shape stellar AGB winds; (iii) the most common explanation is gravitational focusing of the AGB wind on the orbital plane of a companion. This scenario is supported by hydrodynamic simulations (Mastrodemos & Morris 1999; Kim & Taam 2012a) and by observations (e.g. van Winckel et al. 2009).
If the elliptical emission is part of the primary's wind that is focused towards the orbital plane of the GOV companion, that star is orbiting within the disc. Since the high-velocity receding lobe is almost centred on the position of the secondary (see Fig. 4 of Chiu et al. 2006), it is hard to imagine that the GOV companion does not produce any disturbance in the low-velocity CO disc. Therefore, the orientation of the CO disc might not be aligned with the orbital plane of the GOV companion.

The enormous separation of the system ($d > 460$ au) additionally makes it unlikely that the main-sequence companion is able to focus the AGB wind towards the orbital plane. Mastrodemos & Morris (1999) used models with binary separations of 3.6 au to 50.4 au and $M_{\text{prim}}/M_{\text{comp}} = 0.75-6$ in their hydrodynamic simulations, but only the models up to 12 au were able to form bipolar or elliptical CSEs (see Table 2 of Mastrodemos & Morris 1999). The density contrast between the mass-accretion rate of the secondary and the mass-loss rate of the primary thereby defines the degree of focusing of the stellar wind and can be estimated after Morris (1990) as

$$ \alpha_{\text{foc}} \equiv \frac{\dot{M}_{\text{acc}}}{\dot{M}_{\text{prim}}} = \left( \frac{GM_{\text{comp}}}{d} \right)^2 \frac{1}{v_{\text{w}}} \left( v_{\text{w}}^2 + \frac{G(M_{\text{comp}} + M_{\text{prim}})}{d} \right)^{-3/2} . \quad (3) $$

Assuming $M_{\text{prim}} = 1.5 M_{\odot}$, $M_{\text{comp}} = 1.0 M_{\odot}$, $v_{\text{wind}} = 14.5 \text{ km s}^{-1}$, and $d = 460$ au, the focusing ratio is $\alpha_{\text{foc}} = 8 \times 10^{-5}$. According to Han et al. (1995), strong and mild focusing is expressed by $\alpha_{\text{foc}} > 0.1$, which is more than three orders of magnitude higher than the value found for the $\pi^1$ Gru+GOV system. The shaping agent of the elliptical CSE observed by CO and dust emission might therefore be another object that is located much closer to the AGB star than the GOV companion. This hypothesis was first expressed by Chiu et al. (2006) and is discussed here in Sect. 5.3.

If the elliptical emission is indeed an inclined disc that is not located in the orbital plane of the GOV companion, the question arises whether the mass accreted by the star is large enough to form the arc. Given the large system separation of more than 460 au, the accretion rate on the companion is very low, even if the stellar AGB wind expands isotropically. The focusing ratio $\alpha_{\text{foc}} = 8 \times 10^{-5}$ is a factor of four lower than the lowest value simulated by Mastrodemos & Morris (1999) in their model M9. But even at this low rate, a spiral pattern is forming. Nevertheless, the accretion rate of the GOV star is surely enhanced when the star moves through the disc, that is, in the region where the orbital plane and the disc intersect. Unfortunately, the orbital parameters are unknown because of the long period of the GOV companion.

Another interpretation of the CSE ellipticity is that it represents a deformed asterosphere caused by the stellar wind interacting with the ISM. Ueta et al. (2006), Jorissen et al. (2011), and Cox et al. (2012) showed that fast-moving AGB stars can alter the wind bubble and produce a bow shock in the direction of the space motion at the interface of the wind and the ISM. For stars with a low space velocity, the CSE appears elliptical (Weaver et al. 1977). Adopting the long time-scale proper motion from the Tycho-2 catalogue (see Table 2), the direction of the space motion is at PA = 103.1° ± 4.0°, which means that it is aligned with the major axis at PA ≈ 105°. However, the velocity of the space motion is only 15.0 ± 2.8 km s⁻¹, which is comparatively low to cause the elongation. Cox et al. (2012) nevertheless showed that even stars with a space velocity similar to $\pi^1$ Gru are able to form bow shocks (R Leo and UU Aur). The presence of the arc in the direction of the space motion, however,

<table><caption>**Table 2.** Kinematic data of $\pi^1$ Gru from the HIPPARCOS (van Leeuwen 2007) and Tycho-2 (Høg et al. 2000) catalogues.</caption><thead><tr><th></th><th>HIPPARCOS</th><th>Tycho-2</th></tr></thead><tbody><tr><td>μ<sub>α*</sub> (mas/yr)</td><td>28.48 ± 0.94</td><td>33.4 ± 1.1</td></tr><tr><td>μ<sub>δ</sub> (mas/yr)</td><td>-12.14 ± 0.60</td><td>-17.7 ± 1.2</td></tr><tr><td>v<sub>LSR</sub> (km s<sup>-1</sup>)</td><td>11.7 ± 2.5</td><td>15.0 ± 2.8</td></tr><tr><td>i (°)</td><td>-40.4 ± 15.0</td><td>-30.4 ± 13.2</td></tr><tr><td>PA (°)</td><td>81.3 ± 2.9</td><td>103.1 ± 4.0</td></tr></tbody></table>

Notes. $v_{\text{LSR}}$ is the velocity of the star corrected for the solar motion, $i$ and PA the inclination to the sky plane and position angle of the space motion. The parallax $\bar{\sigma} = 6.13 \pm 0.76$ mas from van Leeuwen (2007), the radial velocity $V_r = -5.7 \pm 1.5$ km s⁻¹ from Van Eck et al. (2000), and the solar motion vector $(U, V, W)_\odot = (8.50 \pm 0.29, 13.38 \pm 0.43, 6.49 \pm 0.26)$ from Coşkunoğlu et al. (2011) were used to derive $v_{\text{LSR}}$, $i$, and PA.

diminishes the possibility that the elliptical emission is shaped by the oncoming ISM.

## 5.3. A hidden companion in the $\pi^1$ Gruis system?

As shown above, the elliptical emission around $\pi^1$ Gru cannot be shaped by the GOV companion given its distance and mass. Therefore, an object closer to the star might focus the primary wind on the orbital plane, making $\pi^1$ Gru a hierarchical triple system. Interestingly, Chiu et al. (2006) found a central cavity with a radius of 200 au (1′′2) in their CO map. The authors note that it is large enough to host a putative close companion, but not the GOV companion orbiting the S star at an angular distance of 2′′8. In the following sections, we discuss further indications for a second companion from astrometric and interferometric observations.

### 5.3.1. $\Delta\mu$ behaviour and the HIPPARCOS Intermediate Astrometric Data

$\pi^1$ Gru is found to be a $\Delta\mu$ binary (Makarov & Kaplan 2005; Frankowski et al. 2007), meaning that its long-term proper motion (Tycho-2; Høg et al. 2000) is different from its short-term proper motion (HIPPARCOS; van Leeuwen 2007), because the latter is altered by the orbital motion while the orbital motion averages out on the long-term proper motion (see Table 2 and also Appendix A.1).

In the following we investigate whether the $\Delta\mu$ binary nature of $\pi^1$ Gru can be caused by the GOV companion separated by more than 460 au. The observed proper motion discrepancy is derived as

$$ \begin{aligned} \Delta\mu &= \left[ (\mu_{\alpha*}(\text{HIP}) - \mu_{\alpha*}(\text{TYC2}))^2 + (\mu_{\delta}(\text{HIP}) - \mu_{\delta}(\text{TYC2}))^2 \right]^{1/2} \\ &= 7.4 \pm 1.4 \text{ mas yr}^{-1}. \end{aligned} \quad (4) $$

Frankowski et al. (2007) have shown, thanks to a comparison with known spectroscopic binaries from the *Ninth Catalogue of Spectroscopic Binary Orbits* (S$_B^*$; Pourbaix et al. 2004), that the binaries detectable by the $\Delta\mu$ approach must have orbital periods in the range of 1500 to 30 000 days. This already implies that the GOV companion with an orbital period of ≈2.2 × 10⁶ days can hardly account for the proper motion variation. For a more detailed analysis, Makarov & Kaplan (2005) showed that $\Delta\mu$ is related to the orbital parameters in the following way:

$$ \Delta\mu \le \frac{2\pi\sigma R_0 M_2}{(M_1 + M_2)^{2/3} P^{1/3}}, \quad (5) $$
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Constraints on the orbital period and the companion mass $M_2$ set by the $\Delta\mu$ criterion (Eq. (7); above the blue line) and by the Roche lobe criterion (Eq. (10); below the red line) adopting a stellar radius of 420 $R_\odot$ and a primary mass of 1.5 $M_\odot$. The latter condition imposes the companion to be less massive than 1.5 $M_\odot$ (below the dashed horizontal line). The ability to detect a $\Delta\mu$ binary moreover imposes 1500 < $P$(days) < 30 000 (between both vertical dashed lines). The admissible region is enclosed within these boundaries (green area).</figcaption></figure>

where $M_1$ and $M_2$ are the primary and secondary masses, $\bar{\omega}$ is the parallax, $P$ is the orbital period, and $R_0$ is a time-dependent orbital phase term,

$$R_0 = \left( \frac{1 + e \cos E}{1 - e \cos E} \right)^{1/2}, \qquad (6)$$

where $e$ is the orbital eccentricity and $E$ the eccentricity anomaly. In the following, $R_0 = 1$ is assumed, equivalent to a circular orbit.

The equality in Eq. (5) corresponds to an orbit seen face-on, thus all our estimations in the remainder of this section correspond to lower bounds. With $\bar{\omega} = 6.13$ mas (van Leeuwen 2007) and $\Delta\mu = 7.4$ mas yr$^{-1}$, Eq. (5) yields

$$0.192 P^{1/3} \le \frac{M_2}{(M_1 + M_2)^{2/3}}. \qquad (7)$$

Another constraint on the orbital period comes from the Roche lobe radius $R_1$ (Paczyński 1971). Since $\pi^1$ Gru appears to be close to the tip of the AGB (Van Eck et al. 1998, and Sect. 2), its large radius limits the admissible orbital separation and period

$$R_1 \le a \left(0.38 + 0.2 \log \frac{M_1}{M_2}\right). \qquad (8)$$

By substituting $a$ from the third Kepler law

$$P^2 = \frac{4\pi^2}{G(M_1 + M_2)} a^3, \qquad (9)$$

one obtains

$$R_1 \le P^{2/3} \left( \frac{G(M_1 + M_2)}{4\pi^2} \right)^{1/3} \left( 0.38 + 0.2 \log \frac{M_1}{M_2} \right), \qquad (10)$$

where $G$ is the gravitational constant. Equations (7) and (10) allow us to restrict the range of possible values for $P$ and $M_2$ (Fig. 6) if we adopt 420 $R_\odot$ for the radius of $\pi^1$ Gru (12 mas

<table><caption>**Table 3.** Possible orbital solutions obtained from the analysis of the 55 HIPPARCOS Intermediate Astrometric Data of $\pi^1$ Gru (HIP 110478), by scanning a $P-e$ grid.</caption><thead><tr><th>e</th><th>P<br>(yr)</th><th>χ<sup>2</sup></th><th>F2</th><th>ω<br>(mas)</th><th>μ<sub>a*</sub><br>(mas yr<sup>-1</sup>)</th><th>μ<sub>δ</sub><br>(mas yr<sup>-1</sup>)</th></tr></thead><tbody><tr><td>0.5</td><td>8.3</td><td>60.45</td><td>1.74</td><td>6.55</td><td>26.1</td><td>-19.1</td></tr><tr><td>0.5</td><td>4.6</td><td>59.28</td><td>1.64</td><td>6.53</td><td>26.4</td><td>-16.4</td></tr><tr><td>0.7</td><td>4.6</td><td>56.74</td><td>1.42</td><td>6.55</td><td>27.3</td><td>-17.8</td></tr><tr><td>0.9</td><td>9.7</td><td>55.66</td><td>1.32</td><td>6.55</td><td>34.3</td><td>-21.9</td></tr><tr><td>0.9</td><td>4.6</td><td>54.01</td><td>1.17</td><td>6.66</td><td>28.9</td><td>-17.9</td></tr><tr><td>0.9</td><td>6.3</td><td>53.18</td><td>1.09</td><td>6.68</td><td>31.2</td><td>-18.8</td></tr></tbody></table>

Notes. F2 is the goodness-of-fit, as defined by Eq. (11).

at 163 pc) as derived by Cruzalèbes & Sacuto (2006) from VLTI/MIDI observations, and $M_1 = 1.5 M_\odot$ for the primary mass from its location in the Hertzsprung-Russell diagram (Van Eck et al. 1998). Figure 6 presents these constraints in graphical form.

We thus conclude from this simple analysis that the GOV companion separated by at least 460 au with an orbital period of several thousand years cannot be the cause of the $\Delta\mu$ binary. A close second companion is required instead. An analysis of the HIPPARCOS IAD was performed along the method outlined by Pourbaix & Jorissen (2000)<fnref n="4" /> and especially by Jorissen (2004), to search for the possible presence of a close binary companion in the HIPPARCOS IAD without any a priori knowledge of the orbital elements. Basically, the routine scans a grid in eccentricity – period and searches for the best possible solution (in terms of $\chi^2$ value) by including orbital motion at each grid point. Satisfactory solutions (i.e., with $\chi^2$ in the range of 53 to 60, because 55 data points are available, or goodness-of-fit<fnref n="5" /> values between 1.0 and 1.8) are obtained for eccentricities higher than 0.5 and orbital periods in the range of 5 to 11 yr. Although the available data do not allow to fully constrain the orbit, confidence in the orbital solutions obtained from the HIPPARCOS IAD is bolstered because the proper motion derived from the analysis of the HIPPARCOS data now agrees with the long-term Tycho-2 proper motion. Possible solutions are listed in Table 3. The orbital arc corresponding to the best fitting among these possible solutions (last entry in Table 3) is presented in Fig. 7. Moreover, the favoured orbital periods (4.6–9.7 yr, or 1680–3540 d) are in the range considered to be likely from the analysis of Eqs. (7) and (10) (4.1–52.3 yr). According to Fig. 6, periods as short as 4.6 yr are only marginally possible, but solutions with orbital periods around 6 yr or longer are perfectly admissible and imply masses for the companion in the range of 0.5 to 1.5 $M_\odot$ that would correspond to spectral types K9 to F3 on the main sequence.

A system similar to $\pi^1$ Gru + close companion was modelled by Mastrodemos & Morris (1999) in their models M10 (1 $M_\odot$ companion) and M17 (0.5 $M_\odot$ companion). The radius of the primary ($R_p = 452.4 R_\odot$) and the separation of the system

$$F2 = \sqrt{\frac{9\nu}{2} \left[ \sqrt[3]{\frac{\chi^2}{\nu}} + \frac{2}{9\nu} - 1 \right]}. \qquad (11)$$

<footnote n="4">⁴ Except for the fact that the condition imposing a positive parallax has since been found to be inappropriate and has been lifted in recent applications of the Pourbaix & Jorissen (2000) method.</footnote>

<footnote n="5">⁵ If $\chi^2$ does follow a chi-square distribution with $\nu$ degrees of freedom, the goodness-of-fit follows a $N(0, 1)$-distribution irrespective of $\nu$ (see e.g., Pasquato et al. 2011):</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Orbital arc derived from the analysis of the HIPPARCOS Intermediate Astrometric Data of π¹ Gru (see last entry in Table 3). The HIPPARCOS observations are 1D data, the bars therefore indicate the possible position of the photocentre perpendicular to the orbital segment.</figcaption></figure>

$(a = 6.3 \text{ au})$ are the same for both models and resemble the results obtained from interferometric observations (Sacuto et al. 2008) and the $\Delta\mu$ estimate above.

In model M10, the wind morphology is indeed collimated and oblate due to the gravitational force of the companion on the spherical wind of the AGB star. It is thus conceivable that a close 1 $M_\odot$ companion is the cause for the disc observed in the CO and dust emission. A lower-mass companion, as in M17, however, prevents this behaviour and a well-defined spiral pattern occurs instead. The size of the spiral pattern can be evaluated with Eq. (1), and results in $\rho \approx 0.2 (\approx 33 \text{ au})$. This is much smaller than the 1" pixel size of the 70 µm PACS image, which does not allow us to favour either of the models.

### 5.3.2. Interferometric observations

The angular resolution of VLTI/AMBER is perfect to investigate the deformation of the envelope induced by a close secondary companion.

The usual first step for interpreting interferometric observations with limited uv-coverage, like the AMBER data presented here, is the comparison with geometric models. For this purpose we made use of the software GEM-FIND described by Klotz et al. (2012). It is sufficient here to report that none of the geometric toy models could fit the data in a satisfactory way. Given the presence of closure-phase signatures (i.e. asymmetric structures), 1D model atmospheres cannot reproduce the observations either. Therefore, we decided to switch to a more realistic physical model. The presence of a binary companion very close to the primary would trigger a tidal deformation of the primary star. If the primary is close to filling its Roche lobe (and this possibility is not excluded by the results presented in Sect. 5.3.1 since Roche lobe fitting giants are located along the left curve

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Intensity map of the Roche lobe model for one of the AMBER low-resolution spectral channels. The lines represent the direction of the projected baselines used to simulate the interferometric data. Since the location of the close companion and the orientation of the orbital plane are unknown, the depicted Roche lobe orientation illustrates only one possible solution.</figcaption></figure>

of Fig. 6), the shape will resemble that of a pear, as shown in Fig. 8. This kind of geometry will produce a signature in the closure phase.

Siopis & Sadowski (2012) developed the software package Gaia Eclipsing System Simulator and Solver (GESSS) with the primary aim of modelling the light curve of eclipsing binaries for the Gaia survey. This tool is very flexible and can also be used to model other binary configurations. The code was recently adapted to extract interferometric observables (Paladini et al. 2014).

For our experiment, we used as starting point a MARCS model atmosphere (Gustafsson et al. 2008) with stellar parameters $T_{\text{eff}} = 3200$ K, $\log g = 0.3$, and solar metallicity. The Gaia eclipsing binary software identifies the stellar surface with Roche equipotentials, which are numerically computed for each component of the binary system using a dense mesh of points. This mesh defines a scalar field of intensities (calculated from the MARCS model that incorporates the limb darkening), which is then linearly interpolated to produce a synthetic 1020-by-1020-pixel image of the system. We produced a set of 27 intensity maps spread across the H and K bands, with the spectral resolution of AMBER ($R = 30$). This set of images was produced for two Roche lobe models, one with mass ratio 1/3, and one with mass ratio 1. As we do not know the orientation of the orbital plane, we assumed for simplification that for both models the companion is currently located at one of the orbital nodes, that is, at the intersection of the orbital plane with the plane of the sky. Thus, the Roche lobe is seen face-on. An example of the intensity maps in one of the AMBER low-resolution spectral channels is shown in Fig. 8.

The orientation of the π¹ Gru system on the sky with respect to the projected baselines is unknown. In principle, the
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Comparison between AMBER data (crosses) and the synthetic Roche lobe observations (full line for the *upper three panels*, triangles for *the others*; see Fig. 8). The *first three panels* depict the visibilities of the E0-G0-H0 configuration (see Table 1). The *other panels* show the closure phases produced by the triplet of baselines B at position angle PA.</figcaption></figure>

full parameter space of azimuthal *and* polar angles describing
the orientation of the system on the plane of the sky has to be
explored. Together with the position angle of the baselines, sim-
ilar images to that of Fig. 8 have to be generated for all these
possible orientations to compute the corresponding visibilities
and closure phases. This effort is beyond the scope of this paper,
and we here restrict ourselves to showing that the Roche lobe
model may yield closure-phase variations qualitatively similar
to those observed. We chose the system orientation that yielded

the highest asymmetry (shown in Fig. 8), and searched for the
baseline position angle that yielded closure phases that varied in
the same way as the observations. Starting from a first baseline
with a position angle of either –60, –40, 0, 10, 20, or 30 degrees,
we extracted visibilities and closure phases for the same baseline
pattern as displayed in Fig. 3.

In Fig. 8, we show the orientation of the baseline pattern of
Fig. 3 needed to obtain the model curves displayed in Fig. 9.
These model curves should by no means be considered as best
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Illustration of the surroundings of π¹ Gru including the close companion, the disc, the GOV companion, and the spiral arc. The coloured areas represent the field of view or detection range of the used instrument: olive represents the VLTI/AMBER range, horizontal grey lines the HIPPARCOS Δμ constraint, light blue the SMA CO(2-1) observations, and dark orange the partly overlapping Herschel/PACS field of view.</figcaption></figure>

fits, but they illustrate the good prospects offered by the Roche model. Nevertheless, it has to be kept in mind that not only the Roche lobe geometry can cause non-zero signatures in the closure phase, but so can stellar spots or flares (e.g.; Chiavassa et al. 2010; Wittkowski et al. 2011).

The Roche lobe model suggests that a solution with the smallest orbital separation (red line in Fig. 6, corresponding to Roche lobe filling) is consistent with both the AMBER and HIPPARCOS data. The AMBER data do not allow us to unambiguously select one among the possible HIPPARCOS solutions (Table 3), but they offer confidence for the hypothesis of a close companion in the π¹ Gru system.

As a summary, the whole scenario including the close companion, the disc, the GOV companion, and the spiral arc are illustrated in Fig. 10. The coloured areas represent the field of view or detection range of the respective observing facility. The horizontal grey lines around the close companion indicate the range of acceptable orbital separations using the HIPPARCOS Δμ constraint. This is refined to the region indicated by the “error bar” which is based on the results from the HIPPARCOS IAD fitting (see Table 3).

**6. Conclusions and summary**

We have analysed the CSE of the highly evolved AGB star π¹ Gru based on Herschel observations at 70 µm and 160 µm as part of the MESS sample. The images show an asymmetric stellar wind morphology with two main features, namely an elliptical CSE and an arc east of the star. The arc emerges 38'' away from the star along the major axis of the ellipse and is curved towards the north-east before it becomes too faint. The arc is most likely a small part of an Archimedean spiral caused by the interaction of the stellar AGB wind with a companion. π¹ Gru has a physically related GOV companion that has been known for over a century and is separated by 460 au (2.'8) from the primary. We were able to fit the observations with a spiral given the properties of the GOV star. For a perfect match, however, the wind velocity

has to be adjusted to higher values, as was suggested by Knapp et al. (1999) to explain their CO observations.

The second far-IR feature, the elliptical CSE, stretches over 72'' × 60'' [11 750 × 9790 au] and represents the dusty counterpart of CO emission with the same axis ratio, but it is slightly smaller (Knapp et al. 1999; Chiu et al. 2006; Sahai 1992). All of the authors interpreted their findings as a disc structure inclined by 35° to the plane of the sky, which is also supported by the axis ratio of the far-IR emission. In the CO map by Chiu et al. (2006) the disc has an inner radius of 1.'2 and an outer radius of ≈20''. Based on the focusing ratio, it can be ruled out that the known GOV companion focuses the AGB wind towards the orbital plane given the enormous separation of the system. Furthermore, the GOV companion would be orbiting within the disc without causing any observable disturbance. Because of this, we assumed that the disc is not located in the orbital plane of the GOV companion and followed the hypothesis of Chiu et al. (2006) that π¹ Gru may have a close second companion.

We found support for this assumption in several observations. π¹ Gru is a known Δμ binary, meaning that its long-term and short-term proper motions are significantly different (Makarov & Kaplan 2005; Frankowski et al. 2007). An analysis of the HIPPARCOS (van Leeuwen 2007) and Tycho-2 (Høg et al. 2000) data eliminates the GOV companion as the source of disturbance and suggests a 0.5–1.5 M☉ companion with an orbital period in the range of 4–50 yr. This result is strengthened by the HIPPARCOS IAD, which reveal an orbital motion of the photocentre of π¹ Gru. This motion is best fitted by a highly eccentric orbit with a period of 4.6–9.7 yr. Mastrodemos & Morris (1999) used this configuration in their hydrodynamic simulations where the companion was indeed able to focus the primary’s wind towards the orbital plane.

To obtain direct indications for the close companion we used archival data from VLTI/AMBER that show closure-phase signatures. Although the interferomeric observations can be qualitatively reproduced by a Roche lobe model, we cannot exclude
<table><caption>**Table A.1.** Positional data of $\pi^1$ Gru from various catalogues over the past century.</caption><thead><tr><th></th><th>SAO<sup>(1)</sup></th><th>AC 2000.2<sup>(2)</sup></th><th>PPM<sup>(3)</sup></th><th>YZC<sup>(4)</sup></th><th>CPC-2<sup>(5)</sup></th><th>SHCB<sup>(6)</sup></th></tr></thead><tbody><tr><td>Mean epoch of RA</td><td>1896</td><td>1903.7</td><td>1950.9</td><td>1961.8</td><td>1962.1</td><td>1974.6</td></tr><tr><td>Mean epoch of Dec</td><td>1890</td><td>1903.7</td><td>1951.0</td><td>1961.8</td><td>1962.1</td><td>1974.6</td></tr><tr><td>μ<sub>α*</sub> (mas yr<sup>-1</sup>)</td><td>9 ± 11</td><td>–</td><td>38.2 ± 4</td><td>39</td><td>38.6 ± 3</td><td>–</td></tr><tr><td>μ<sub>δ</sub> (mas yr<sup>-1</sup>)</td><td>-7 ± 8</td><td>–</td><td>-17 ± 4</td><td>-22</td><td>-13 ± 3</td><td>–</td></tr><tr><td>RA at obs. epoch [°]</td><td>335.682791</td><td>335.682854</td><td>335.683458</td><td>335.683616</td><td>335.683791</td><td>335.684319</td></tr><tr><td>Error in RA [10<sup>-5</sup> °]</td><td>5.833</td><td>7.056</td><td>2.5</td><td>–</td><td>1.361</td><td>6.111</td></tr><tr><td>Dec at obs. epoch [°]</td><td>-45.947583</td><td>-45.947530</td><td>-45.947752</td><td>-45.947802</td><td>-45.947802</td><td>-45.947926</td></tr><tr><td>Error in Dec [10<sup>-5</sup> °]</td><td>5</td><td>2.694</td><td>2.5</td><td>–</td><td>1.778</td><td>5.556</td></tr></tbody></table><table><thead><tr><th></th><th>SPM4<sup>(7)</sup></th><th>FOCAT-S<sup>(8)</sup></th><th>Tycho-2<sup>(9)</sup></th><th>PPMXL<sup>(10)</sup></th><th>USNO-B<sup>(11)</sup></th></tr></thead><tbody><tr><td>Mean epoch of RA [yr]</td><td>1983.3</td><td>1985.3</td><td>1991.4</td><td>1991.2</td><td>2000.0</td></tr><tr><td>Mean epoch of Dec [yr]</td><td>1983.3</td><td>1985.3</td><td>1991.4</td><td>1991.2</td><td>2000.0</td></tr><tr><td>μ<sub>α*</sub> [mas yr<sup>-1</sup>]</td><td>26.9 ± 4.8</td><td>32 ± 1</td><td>33.4 ± 1.1</td><td>35.3 ± 2</td><td>32</td></tr><tr><td>μ<sub>δ</sub> [mas yr<sup>-1</sup>]</td><td>-20.2 ± 4.3</td><td>-12 ± 4</td><td>-17.7 ± 1.2</td><td>-16.3 ± 2</td><td>-18</td></tr><tr><td>RA at obs. epoch [°]</td><td>335.684025</td><td>335.684000</td><td>335.684094</td><td>335.684177</td><td>335.684203</td></tr><tr><td>Error in RA [10<sup>-5</sup> °]</td><td>1.183</td><td>5</td><td>0.111</td><td>0.056</td><td>–</td></tr><tr><td>Dec at obs. epoch [°]</td><td>-45.947895</td><td>-45.947872</td><td>-45.947916</td><td>-45.947915</td><td>-45.947959</td></tr><tr><td>Error in Dec [10<sup>-5</sup> °]</td><td>1.025</td><td>6.667</td><td>0.111</td><td>0.056</td><td>–</td></tr></tbody></table>

**Notes.** All positions are given in the J2000.0 system.

**References.** (1): Smithsonian Astrophysical Observatory Star Catalog (SAO Staff 1966); (2): The Astropiche Catalogue on the HIPPARCOS System (Urban et al. 2001); (3): Positions and Proper Motions – South (Bastian & Röser 1993); (4): Yale Zone Catalogues Integrated (Fallon 1983); (5): Cape Photographic Catalogue 2 (Nicholson et al. 1984); (6): Southern Hemisphere Catalogue of Bordeaux (Rousseau et al. 1996); (7): Yale/San Juan Southern Proper Motion Catalog 4 (Girard et al. 2011); (8): Pulkovo photographic Catalogue of Southern Hemisphere (Bystrov et al. 1994); (9): The Tycho-2 Catalogue of the 2.5 Million Brightest Stars (Høg et al. 2000); (10): The PPMXL catalog of positions and proper motions on the ICRS (Roeser et al. 2010); (11): The USNO-B Catalog (Monet et al. 2003).

that a more complex model including the presence of spots or flares will also be able to reproduce these observations. The main restriction comes from the uv-plane coverage. A VLTI/PIONIER imaging programme would probably help to break the degeneracy in the currently available data. Based on the current observational status, $\pi^1$ Gru is most likely a hierarchical triple system in which the close companion shapes the disc observed in CO and dust emission, while the (previously known) G0V companion is located farther outside and causes the spiral arc visible in the Herschel/PACS images.

**Acknowledgements.** We thank B. Mason (USNO) for providing the individual position measurements of the companion of $\pi^1$ Gru from the WDS catalogue and the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. We also thank I. Platais for advice with the proper motion evaluation. This work was supported in part by the Belgian Federal Science Policy Office via the PRODEX Programme of ESA. A.M. and F.K. acknowledge funding by the Austrian Research Promotion Agency FFG under project number FA 538019, RO under project number I163-N16. Furthermore, A.M. acknowledges funding by the Abschlusstipendium 2014 of the University of Vienna. This work used data from the ESO archive (Programmes 076.D-0624, 077.D-0620, 078.D-0122, 079.D-0138, 080.D-0076, 187.D-0924) and from the HST archive (Programs ID 3603 and 10185). This research has made use of the AMBER data reduction package of the Jean-Marie Mariotti Center$^6$.

## Appendix A: Positional data

### A.1. $\pi^1$ Gru

The difference between the long- and short-term proper motions is illustrated in Fig. A.1, which shows the positions of $\pi^1$ Gru in the last 100 years from various catalogues (see Table A.1 for a list of these catalogue positions). The position at the observation epoch of a given catalogue, if not directly given by the Vizier database at the Centre de Données Stellaires (Strasbourg),

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

**Fig. A.1.** Positions of $\pi^1$ Gru from catalogues spanning more than 100 years. The solid red line corresponds to a linear fit through the positions, weighted by the inverse square of the position uncertainties. The dashed blue line shows the positions extrapolated according to the HIPPARCOS proper motion, while the dashed orange line represents the same for the Tycho-2 proper motion.

<footnote>$^6$ Available at http://www.jmmc.fr/amberdrds</footnote>
<table><caption>**Table A.2.** Separations and position angles of the GOV companion of π¹ Gru, from the Washington Double Star Observations catalogue (courtesy of B. Mason).</caption><thead><tr><th>Obs. Epoch<br>(yr)</th><th>PA<br>(°)</th><th>d<br>(")</th><th>Reference</th></tr></thead><tbody><tr><td>1896.8</td><td>190</td><td>2.5</td><td>Innes (1897)</td></tr><tr><td>1900.75</td><td>201.6</td><td>3.04</td><td>Innes (1905)</td></tr><tr><td>1900.76</td><td>201.4</td><td>2.70</td><td>Lunt (1908)</td></tr><tr><td>1912.63</td><td>201.4</td><td>2.27</td><td>Innes (1914)</td></tr><tr><td>1926.02</td><td>200.6</td><td>2.74</td><td>van den Bos (1928)</td></tr><tr><td>1929.02</td><td>200.4</td><td>2.78</td><td>Rossiter (1955)</td></tr><tr><td>1936.74</td><td>201.2</td><td>2.75</td><td>van den Bos (1938)</td></tr><tr><td>1943.49</td><td>202.1</td><td>2.71</td><td>Voûte (1955)</td></tr><tr><td>1956.446</td><td>200.83</td><td>2.831</td><td>The (1975)</td></tr><tr><td>1960.80</td><td>202.8</td><td>2.63</td><td>van den Bos (1961)</td></tr><tr><td>1966.81</td><td>201.0</td><td>2.82</td><td>Knipe (1969)</td></tr><tr><td>1975.722</td><td>201.0</td><td>2.79</td><td>Worley (1978)</td></tr><tr><td>1989.86</td><td>200.4</td><td>2.45</td><td>Sahai (1992)</td></tr><tr><td>2003.47</td><td>203.1</td><td>2.82</td><td>HST (PI: Sahai)</td></tr></tbody></table>

has been derived by applying the catalogue proper motion to the listed J2000 epoch position. Care has been exercised to ensure consistency between the equinox of the proper motion and the position. Our own estimate for the proper motion is obtained from a linear fit on all these positions, weighted by the inverse square of the uncertainty on the position.

The resulting space velocity $v_{LSR}$ (corrected for the solar motion vector $(U, V, W)_\odot = (8.50 \pm 0.29, 13.38 \pm 0.43, 6.49 \pm 0.26)$; Coşkunoğlu et al. 2011) based on these proper motions and the radial velocity derived by Van Eck et al. (2000) is listed in Table 2. Since the long-term Tycho-2 proper motion appears slightly more precise than the one we derived from the positions displayed in Fig. A.1, the Tycho-2 proper motion is adopted in the remainder of this paper.

## A.2. GOV companion

Table A.2 shows the separations and position angles of the GOV companion. The observations do not show significant changes of the companion’s position over a timespan of more than 100 years. However, in 1989.86, the system was observed by Sahai (1992), who reported the position of the companion with a projected separation of 2.45 at PA = 200.4°, clearly different from previous observations. The author remarks that such a fast motion for an object located at least 460 au away from the primary would imply unrealistically large stellar masses. Moreover, given the large uncertainties induced by the 0.99 seeing that prevailed during the observations, no reliable conclusions regarding the orbital period could be drawn by Sahai (1992).

## References

* Abate, C., Pols, O. R., Izzard, R. G., Mohamed, S. S., & de Mink, S. E. 2013, A&A, 552, A26

* Ake, T. B., & Johnson, H. R. 1992, in BAAS, 24, Meeting Abstracts #180, 788

* Barnbaum, C., Morris, M., & Kahane, C. 1995, ApJ, 450, 862

* Bastian, U., & Röser, S. 1993, PPM Star Catalogue. Positions and proper motions of 197179 stars south of -2.5 degrees declination for equinox and epoch J2000.0. Vol. III: Zones -00deg to -20deg. Vol. IV: Zones -30deg to -80deg

* Bordé, P., Coudé du Foresto, V., Chagnon, G., & Perrin, G. 2002, A&A, 393, 183

* Bystrov, N. F., Polojentsev, D. D., Potter, H. I., et al. 1994, Bulletin d'Information du Centre de Données Stellaires, 44, 3

* Castro-Carrizo, A., Quintana-Lacaci, G., Neri, R., et al. 2010, A&A, 523, A59

* Chelli, A., Utrera, O. H., & Duvert, G. 2009, A&A, 502, 705

* Chiavassa, A., Lacour, S., Millour, F., et al. 2010, A&A, 511, A51

* Chiu, P.-J., Hoang, C.-T., Dinh-V-Trung, et al. 2006, ApJ, 645, 605

* Coşkunoğlu, B., Ak, S., Bilir, S., et al. 2011, MNRAS, 412, 1237

* Cox, N. L. J., Kerschbaum, F., van Marle, A.-J., et al. 2012, A&A, 537, A35

* Cruzalèbes, P., & Sacuto, S. 2006, A&A, 452, 885

* Cruzalèbes, P., Jorissen, A., Sacuto, S., & Bonneau, D. 2010, A&A, 515, A6

* De Marco, O., Hillwig, T. C., & Smith, A. J. 2008, AJ, 136, 323

* Decin, L., Royer, P., Cox, N. L. J., et al. 2011, A&A, 534, A1

* Decin, L., Cox, N. L. J., Royer, P., et al. 2012, A&A, 548, A113

* Dinh-V.-Trung, & Lim, J. 2009, ApJ, 701, 292

* Eggen, O. J. 1975, ApJ, 195, 661

* Fallon, F. W. 1983, Transactions of the Astronomical Observatory of Yale University, 32, 1

* Feast, M. W. 1953, MNRAS, 113, 510

* Frankowski, A., Jancart, S., & Jorissen, A. 2007, A&A, 464, 377

* Girard, T. M., van Altena, W. F., Zacharias, N., et al. 2011, AJ, 142, 15

* Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3

* Groenewegen, M. A. T., & de Jong, T. 1998, A&A, 337, 797

* Groenewegen, M. A. T., Waelkens, C., Barlow, M. J., et al. 2011, A&A, 526, A162

* Guandalini, R., & Busso, M. 2008, A&A, 488, 675

* Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

* Han, Z., Podsiadlowski, P., & Eggleton, P. P. 1995, in Asymmetrical Planetary Nebulae, eds. A. Harpaz, & N. Soker, Annals of the Israel Physical Society, 11, 66

* Hirano, N., Shinnaga, H., Dinh-V-Trung, et al. 2004, ApJ, 616, L43

* Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27

* Huggins, P. J. 2007, ApJ, 663, 342

* Innes, R. T. A. 1897, MNRAS, 57, 533

* Innes, R. T. A. 1905, Ann. Cape Obs., 2, 4

* Innes, R. T. A. 1914, Circular of the Union Observatory Johannesburg, 14, 97

* Jorissen, A. 2004, in Asymptotic Giant Branch Stars, eds. H. Habing, & H. Olofsson (New York: Springer Verlag), 461

* Jorissen, A., Frayer, D. T., Johnson, H. R., Mayor, M., & Smith, V. V. 1993, A&A, 271, 463

* Jorissen, A., Mayer, A., van Eck, S., et al. 2011, A&A, 532, A135

* Kafatos, M., Hollis, J. M., Yusef-Zadeh, F., Michalitsianos, A. G., & Elitzur, M. 1989, ApJ, 346, 991

* Kawaler, S. D., Sekii, T., & Gough, D. 1999, ApJ, 516, 349

* Keenan, P. C. 1954, ApJ, 120, 484

* Kerschbaum, F., Ladjal, D., Ottensamer, R., et al. 2010, A&A, 518, L140

* Kim, H., & Taam, R. E. 2012a, ApJ, 759, L22

* Kim, H., & Taam, R. E. 2012b, ApJ, 744, 136

* Kim, H., & Taam, R. E. 2012c, ApJ, 759, 59

* Kim, H., Hsieh, I.-T., Liu, S.-Y., & Taam, R. E. 2013, ApJ, 776, 86

* Klotz, D., Sacuto, S., Paladini, C., Hron, J., & Wachter, G. 2012, in SPIE Conf. Ser., 8445, 14

* Knapp, G. R., Young, K., & Crosas, M. 1999, A&A, 346, 175

* Knipe, G. F. G. 1969, Republic Observatory Johannesburg Circular, 128, 177

* Ladjal, D., Barlow, M. J., Groenewegen, M. A. T., et al. 2010a, A&A, 518, L141

* Ladjal, D., Justtanont, K., Groenewegen, M. A. T., et al. 2010b, A&A, 513, A53

* Leal-Ferreira, M. L., Vlemmings, W. H. T., Kemball, A., & Amiri, N. 2013, A&A, 554, A134

* Lèbre, A., Aurière, M., Fabas, N., et al. 2014, A&A, 561, A85

* Livio, M., & Soker, N. 1988, ApJ, 329, 764

* Lunt, J. 1908, Ann. Cape Obs., 2, 4

* Maercker, M., Mohamed, S., Vlemmings, W. H. T., et al. 2012, Nature, 490, 232

* Makarov, V. V., & Kaplan, G. H. 2005, AJ, 129, 2420

* Mastrodemos, N., & Morris, M. 1998, ApJ, 497, 303

* Mastrodemos, N., & Morris, M. 1999, ApJ, 523, 357

* Matt, S., Balick, B., Winglee, R., & Goodson, A. 2000, ApJ, 545, 965

* Mauron, N., & Huggins, P. J. 2006, A&A, 452, 257

* Mauron, N., Huggins, P. J., & Cheung, C.-L. 2013, A&A, 551, A110

* Mayer, A., Jorissen, A., Kerschbaum, F., et al. 2011, A&A, 531, L4

* Mayer, A., Jorissen, A., Kerschbaum, F., et al. 2013, A&A, 549, A69 (Paper I)

* Meaburn, J., López, J. A., Boumis, P., Lloyd, M., & Redman, M. P. 2009, A&A, 500, 827

* Mečina, M., Kerschbaum, F., Groenewegen, M. A. T., et al. 2014, A&A, 566, A69

* Miszalski, B., Acker, A., Moffat, A. F. J., Parker, Q. A., & Udalski, A. 2009a, A&A, 496, 813

* Miszalski, B., Acker, A., Parker, Q. A., & Moffat, A. F. J. 2009b, A&A, 505, 249

* Mohamed, S., & Podsiadlowski, P. 2011, in Why Galaxies Care about AGB Stars II: Shining Examples and Common Inhabitants, eds. F. Kerschbaum, T. Lebzelter, & R. F. Wing, ASP Conf. Ser., 445, 355

* Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984
Morris, M. 1987, PASP, 99, 1115

Morris, M. 1990, in From Miras to Planetary Nebulae: Which Path for Stellar Evolution, eds. M. O. Mennessier, & A. Omont, 520

Nicholson, W., Penston, M. J., Murray, C. A., & de Vegt, C. 1984, MNRAS, 208, 911

Nordhaus, J., Blackman, E. G., & Frank, A. 2007, MNRAS, 376, 599

Paczyński, B. 1965, Acta Astron., 15, 89

Paczyński, B. 1971, ARA&A, 9, 183

Paladini, C., Jorissen, A., Siopis, C., et al. 2014, in SPIE Conf. Ser. 9146, eds. J. K. Rajagopal, M. J. Creech-Eakman, & F. Malbet

Pasquato, E., Pourbaix, D., & Jorissen, A. 2011, A&A, 532, A13

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

Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1

Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2

Pojmanski, G., & Maciejewski, G. 2005, Acta Astron., 55, 97

Pourbaix, D., & Jorissen, A. 2000, A&AS, 145, 161

Pourbaix, D., Tokovinin, A. A., Batten, A. H., et al. 2004, A&A, 424, 727

Roeser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440

Rossiter, R. A. 1955, Publications of Michigan Observatory, 11, 1

Rousseau, J. M., Perie, J. P., & Gachard, M. T. 1996, VizieR Online Data Catalog: I/230

Roussel, H. 2013, PASP, 125, 1126

Sacuto, S., Jorissen, A., Cruzalèbes, P., et al. 2008, A&A, 482, 561

Sahai, R. 1992, A&A, 253, L33

SAO Staff. 1966, Star Catalog (Washington DC: Smithsonian Institution)

Siess, L. 2006, A&A, 448, 717

Siess, L., & Arnould, M. 2008, A&A, 489, 395

Siopis, C., & Sadowski, G. 2012, in Orbital Couples: Pas de Deux in the Solar System and the Milky Way, eds. F. Arenou & D. Hestroffer (Observatoire de Paris), 59

Soker, N. 1994, MNRAS, 270, 774

Soker, N., & Rappaport, S. 2000, ApJ, 538, 241

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

The, P. S. 1975, A&AS, 21, 51

Theuns, T., & Jorissen, A. 1993, MNRAS, 265, 946

Ueta, T., Speck, A. K., Stencel, R. E., et al. 2006, ApJ, 648, L39

Urban, S. E., Corbin, T. E., Wycoff, G. L., et al. 2001, in Meeting Abstracts, BAAS, 33, 129.04

van den Bos, W. H. 1928, Annalen van de Sterrewacht te Leiden, 14, D1

van den Bos, W. H. 1938, Circular of the Union Observatory Johannesburg, 100, 481

van den Bos, W. H. 1961, Circular of the Union Observatory Johannesburg, 120, 353

Van Eck, S., & Jorissen, A. 1999, A&A, 345, 127

Van Eck, S., Jorissen, A., Udry, S., Mayor, M., & Pernier, B. 1998, A&A, 329, 971

Van Eck, S., Jorissen, A., Udry, S., et al. 2000, A&AS, 145, 51

van Hoof, P. A. M., Van de Steene, G. C., Exter, K. M., et al. 2013, A&A, 560, A7

van Leeuwen, F. 2007, A&A, 474, 653

van Leeuwen, F., & Evans, D. W. 1998, A&AS, 130, 157

van Winckel, H., Lloyd Evans, T., Briquet, M., et al. 2009, A&A, 505, 1221

Vlemmings, W. H. T., Humphreys, E. M. L., & Franco-Hernández, R. 2011, ApJ, 728, 149

Voûte, J. 1955, Journal des Observateurs, 38, 109

Wallerstein, G., & Greenstein, J. L. 1980, PASP, 92, 275

Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, R. 1977, ApJ, 218, 377

Winters, J. M., Le Bertre, T., Jeong, K. S., Nyman, L.-Å., & Epchtein, N. 2003, A&A, 409, 715

Wittkowski, M., Boboltz, D. A., Ireland, M., et al. 2011, A&A, 532, L7

Worley, C. E. 1978, Publications of the US Naval Observatory Second Series, 24, 1

Young, K., Phillips, T. G., & Knapp, G. R. 1993, ApJ, 409, 725