Dissecting the AGB star L Puppis: a torus in the making ††thanks: Based on observations from ESO programs 090.D-0541 & 090.D-0677
Key Words.:stellar evolution – AGB stars – high-angular resolution techniques – L$_2$ Pup
Aims:The circumstellar environment of L$_2$ Pup, an oxygen-rich semiregular variable, was observed to understand the evolution of mass loss and the shaping of ejecta in the late stages of stellar evolution.
Methods:High-angular resolution observations from a single 8 m telescope were obtained using aperture masking in the near-infrared (1.64, 2.30 and 3.74 ) on the NACO/VLT, both in imaging and polarimetric modes.
Results: The aperture-masking images of L Pup at 2.30 show a resolved structure that resembles a toroidal structure with a major axis of 140 milliarcseconds (mas) and an east–west orientation. Two clumps can be seen on either side of the star, 65 mas from the star, beyond the edge of the circumstellar envelope (estimated diameter is 27 mas), while a faint, hook-like structure appear toward the northeast. The patterns are visible both in the imaging and polarimetric mode, although the latter was only used to measure the total intensity (Stokes ). The overall shape of the structure is similar at the 3.74 pseudo-continuum (dust emission), where the clumps appear to be embedded within a dark, dusty lane. The faint, hook-like patterns are also seen at this wavelength, extending northeast and southwest with the central, dark lane being an apparent axis of symmetry. We interpret the structure as a circumstellar torus with inner radius of 4.2 au. With a rotation velocity of 10 km s as suggested by the SiO maser profile, we estimate a stellar mass of 0.7 M.
During the late stages of evolution of intermediate-mass stars (1–8 M) – that is, at the transition from the asymptotic giant branch (AGB) to the post-AGB phase – the ejected material is shaped into intriguing forms (e.g., bipolar nebulae, torii, disks and/or spirals) that depart from the initial spherical symmetry (van Winckel 2003, and references therein). Although there are many studies in the literature of asymmetries found in post-AGB stars, detections of such asymmetries in AGB stars are rare. It is therefore uncertain at which point the asymmetries begin to develop.
One of the brightest AGB stars in the near-infrared ( mag) near the solar neighborhood ( pc; van Leeuwen 2007) is L$_2$ Pup (HD56096, HIP34922, IRAS 07120-4433). It is a semiregular variable (M5III; days) that experienced a dimming event starting around 1995 (3 mag; Fig. 7), which was attributed to circumstellar dust (Bedding et al. 2002). The decline lasted for six years, and as of June 2014, the star has not yet begun to recover its former brightness. Dimming events are present in a few other evolved stars (e.g., V Hya; Knapp et al. 1999), which may repeat over the course of many years. L Pup had a less severe dimming event around 1960: recovery took until 1980. Unlike V Hya, there is no clear periodicity to the dimming in the data, which go back to 1930.
Only a few images of L Pup are available in the literature. Bedding et al. (2002) stated that the IRAS images are non-descriptive. The WISE satellite observed L Pup post-dimming in the near- and mid-infrared at low angular resolutions, although all images were saturated by the stellar emission. Jura et al. (2002) studied L Pup with the Keck 10 m telescope in the mid-infrared and inferred that two asymmetries are present at position angles 135° and 225° east-of-north. Past studies on the variation of polarization of L Pup are indicative of an asymmetric111In this paper, the term ‘asymmetric’ will describe any aspherical structure. dust distribution (Magalhaes et al. 1986). In addition to the stellar component, Ireland et al. (2004b) inferred from the MAPPIT/AAT (3.9 m) observations in the optical that an additional emitting component is required. They probed different spatial scales above the photosphere and found a size variation with respect to wavelength in the titanium oxide band (0.68 – 0.82 m) as well as an asymmetric brightness profile, which they attributed to scattering by dust. The indications of asymmetry in the circumstellar environment of L Pup, suggested by polarimetry (Magalhaes et al. 1986), as well as the indications of dust scattering by interferometric observations (Ireland et al. (2004b); see also Section 3.1), might indicate a disk or a companion.
Very recently, in an effort parallel to our own, Kervella et al. (2014) have presented infrared interferometric data and deconvolved adaptive-optics images, which show a disk-like structure, approximately 700 mas across, transecting the star, opaque at , translucent at and , and a thermal dust emitter at . A faint loop extends to the northeast. Kervella et al. attributed the dimming to this disk. L Pup is therefore one of the few AGB stars with evidence for asphericities in their circumstellar material, and as the closest mass-losing AGB star, is ideally placed for resolving the details of the envelope structure. We note that the technique of aperture masking used here is different and complementary to the technique of direct adaptive optics imaging used by Kervella et al. (2014), as discussed in detail in Sect. 4.5.
The circumstellar environment of L Pup was observed in the near-infrared, using aperture masking techniques on an 8.2 m telescope (NACO/VLT), to detect any asymmetries within a 05 radius. The observations are presented in the next section, while results are discussed in Sects. 3 and 4. Conclusions are discussed in Sects. 5 and 6.
|Mode||Date||Phase||Filter||() [m]||# frames||CameraCamera pixel scales are 13.22 and 27.19 mas/pixel for the S13 and L27 cameras, respectively.||[sec]||SeeingIn units of arcseconds.||Airmass|
2 Observations and data reduction
L Pup was observed on 1 December 2012 (Prog.ID: 090.D-0541, PI: D. Klotz) at NACO/VLT (Lenzen et al. 2003; Rousset et al. 2003) in sparse aperture masking mode (SAM; Tuthill et al. 2010; Lacour et al. 2011) using the 18–hole333Although the mask has 18 holes, only 17 holes are used effectively, because one is located at the edge of the pupil and is thus omitted in the data reduction. mask at 1.64, 2.30 and 3.74 m.
Additional observations were obtained on 8 February 2013 (Prog. ID: 090.D-0677, PI: A.A. Zijlstra) at 2.30 m using the SAM polarimetric mode (SAMPol). In this mode, a Wollaston prism positioned beyond the aperture mask splits the incoming beam into two polarizations (for a detailed description, see Norris et al. 2012). With the use of a half-wave plate, the plane of polarization of the beam is rotated at 0°, 22.5°, 45° , and 67.5° to obtain the polarized intensities at 0°, 45°, 90° , and 135° , that is and . The Stokes vector component, , of the polarized beam can be recovered from these intensities through the following relation of Murakawa et al. (2005):
where is the transmission of the Wollaston prism (93% at 2.30 m).
Table 1 contains details of each observing run, such as stellar phases, integration times, number of acquisition frames, and seeing conditions. The projected 136 baseline lengths were in the range of 0.5 and 7.8 meters, which within only 10 minutes of observing time provides an excellent uv-plane coverage (Fig. 1). SAM non-redundant masks allow recovering diffraction-limited images with a relative resolution of 20, 30, and 50 mas at 1.64, 2.30, and 3.74 , respectively, with minimal disturbance by atmospheric seeing.
Pupil-tracking mode was enabled to prevent the apertures from being covered by the telescope’s spiders. This leads to field rotation on the detector, indicated by the rotation of the -plane coverage as seen in Fig. 1. All speckle frames were recorded in cube mode. To decrease readout time and the volume of recorded data, two subframes were chosen on the detector: pixels for the S13 camera and pixels for the L27 camera. Different dither positions in all four quadrants of those subframes were used in a clockwise pattern to minimize readout time and detector noise (due to the variation of photon and thermal noise in each quadrant during data acquisition), as well as simultaneously recording sky background at an opposing dither position. To allow for determining the transfer function and calibration of the observables, Pup (K3Ib) was observed immediately after the science target.
The data sets were processed with a pipeline developed at the University of Sydney (Tuthill et al. 2000; Lacour et al. 2011). The interferograms in each data cube were flat-fielded, dark, bias, and sky subtracted. Each speckle frame was Fourier transformed, and the extracted Fourier components (power spectra, bispectra) were calibrated with respect to the unresolved target. Because of calibration problems, the results from the 1.64 data are inconclusive.
3 Fitting the visibilities: indications of an extended component
3.1 Past interferometric observations
Ochsenbein & Halbwachs (1982) predicted a diameter for L Pup in the visual of 13.18 milliarcseconds (mas) based on calculations from the intrinsic color and brightness. Ireland et al. (2004b) previously studied L Pup with optical interferometric techniques using aperture masking, although their projected baselines were oriented at 180, that is north–south. They measured an extended component whose size varies with wavelength (Gaussian full-width-half-maximum mas), and the diameter of a compact source (uniform disk, UD) at mas (50% of the total flux), which is almost twice the size of the theoretical prediction of Ochsenbein & Halbwachs (1982).
VINCI/VLTI data (Meisner 2008) show a similar size for the K-band emission (; mas). This VINCI diameter varies over less than a month of observations, but the visibility is dependent on the different projected baselines and on the stellar pulsation phase, that is the duration of the VINCI observations spanned over one fifth of the stellar pulsation period, as shown in Fig. 2. Kervella et al. (2014) also found evidence for an unresolved component with an upper limit of mas, which they interpreted as the central star. The resolved components (76 and 28 mas; MAPPIT) would be affected by circumstellar emission.
3.2 This work
In the SAM power spectra (Fig. 3), a dependence of the visibility on the baseline position angle (P.A.) is shown at similar baseline lengths for both wavelengths. The distinction is clearer at 3.74 , where the difference for visibilities at P.A.°and P.A.°is .
For an initial determination of the overall extent of the detected structures, a simple source morphology was assumed (centro-symmetric) and geometric models were directly fitted to the Fourier data. The initial sizes were estimated by applying a circular Gaussian distribution fit to the azimuthally averaged visibilities, and the errors were estimated from the scatter of the azimuthally averaged measurements. None of these fits were satisfactory since the dependence of the visibility on the baseline P.A. was not taken into account. A further examination of the interferometric data was applied within a 20° cone, centered on the north–south and east–west directions. Tests were performed for both non-averaged and azimuthally averaged squared visibilities. When we applied new Gaussian distribution fits in all cases, the radius of the object was larger along the east–west orientation. A similar result is seen in the SAMPol data at 2.30 . This preferential extension can also be noted in the uv-planes of the original data (Fig. 1).
Using the geometric model simulator GEM-FIND (Klotz et al. 2012b), we found that the shape of the object can be described by an ellipsoid, composed of a circular uniform disk representing the star and an elliptical Gaussian distribution representing the circumstellar environment. The best fits to the Fourier data are shown in Fig. 4. Simpler centrosymmetric functions are also shown for comparison (a uniform disk and a Gaussian distribution).
The two-dimensional fitted structures are also overplotted for comparison to the NACO images in Fig. 6. The ellipsoid extends at approximately 140 mas along the major axis and at a similar position angle (86°) in both wavelengths, while the size of the central component (uniform disk) varies from mas at 2.30 to mas at 3.74 .
4 SAM/NACO imaging
Diffraction-limited maps were reconstructed from the Fourier observables using different methods to independently validate the results. These included the BiSpectrum Maximum Entropy Method (BSMEM; Buscher 1994), the MArkov Chain IMager (MACIM; Ireland et al. 2006), and the Multi-aperture Image Reconstruction Algorithm (MiRA; Thiébaut 2008). Each of these methods uses a different algorithm and regularization techniques to produce the best-fit map from the interferometric observables (i.e., visibilities and closure phases).
4.1 Image reconstruction with BSMEM
The reconstruction algorithm BSMEM is based on the maximum entropy method (MEM; Sivia 1987; Skilling & Bryan 1984). Based on the geometric model-fitting results (Section 3), images were reconstructed from the power spectra and the closure phases of the calibrated interferograms using a Gaussian distribution as an initial model. Although the algorithm converged at , a small misfit was found at low spatial frequencies (that is, baselines shorter than one meter) between the final modeled power spectrum and the data at 2.30 , suggesting either a miscalibration of the power spectrum at these spatial frequencies or that a more complex model was required.
Nevertheless, the image reconstruction produced similar structures, both in the SAMPol and SAM data in that wavelength (Figs. 5 and 6). In the SAMPol Stokes map of February 2013, two bright clumps appear extended to 65 mas in the east–west direction at levels above 2% of the peak intensity. In the SAM maps in particular, a bridge-type structure is more prominent in the west, while a hook-like pattern emerges in the east. It is unclear whether the circumstellar environment of L Pup is polarized at 2.30 . The Stokes map shown here was produced by conventional calibration, and this method could not provide more information on the source’s polarization. However, when using a differential calibration method for the SAMPol data (Norris et al. 2012), the method is sensitive to subtle polarization features. A more detailed analysis of the polarimetric data will be presented in a future paper.
A similar approach was followed for the 3.74 data. Here, no other misfits to the power spectra were found within the errors. In the maps presented in Fig. 6, we detect an elliptical structure with sizes 136 mas and 70 mas along the major and minor axes. This structure is also aligned in the east–west direction. At levels below 5% of the peak intensity, two hook-like patterns are seen in the southwest and northeast directions.
Extended emission might be resolved out by the instrument, and/or the stellar photosphere is unresolved at a given wavelength. We emphasize that the MEM process selects the smoothest map that fits the original data, based on positivity and finite-extension of the map (Narayan & Nityananda 1986; Skilling & Bryan 1984), and it is thus possible that another solution may have fitted the data better. MEM is also known to smooth out fine structures, for instance embedded point sources (Tuthill et al. 2002). For these reasons, imaging results were reproduced with other different algorithms, as shown in the following sections.
4.2 Image reconstruction with MACIM
To verify the consistency of the previous reconstructions, we used MACIM (Ireland et al. 2006). This algorithm uses a regularization method based on simulated annealing (for a detailed description, see Ireland et al. 2006). The initial model consisted of a uniform disk of 25 mas in diameter444Based on the MAPPIT, SAM/NACO and VINCI diameters. at 10% of the total flux. The results are strikingly similar to those of BSMEM (≳1% of the peak intensity), as shown in Fig. 6.
4.3 Image reconstruction with MiRA
MiRA is an algorithm that minimizes a penalty function composed of two terms following a gradient descent method (Thiébaut 2008). The two terms are the distance to the data ( statistics) and the distance to some priors (such as ‘regularization’, which is used to fill the Fourier space between the points of the uv-plane). MiRA implements many different regularizations; here we used the total variation with edge-preserving smoothness with an adopted regularization weight (for a detailed description, see Thiébaut 2008). This regularization minimizes the total gradients in the image and subsequently smooths the flux levels in the image. The BSMEM and MACIM images are similar to MiRA, which succeeds in reproducing the overall shape of the structures, including the clumpy and hook-like patterns seen in all reconstructions (Fig. 6).
4.4 Comparing image reconstruction and GEM-FIND fitting results
In the right column of Fig. 6, the image reconstruction results are compared to the best-fit models of GEM-FIND. While improvements can be implemented, the models provide reasonable555Although the reduced statistic of each model fit exceeds unity, these were the best fits found. fits to the reconstructed images (reduced is 17.48 and 6.34 for 2.30 and 3.74 , respectively; Fig. 4):
Circular uniform disk:
the flux ratio of the central star over the circumstellar envelope is found to be 0.441 and 0.704 at 2.30 and 3.74 . The respective modeled stellar diameters coincide with contour levels at 50 and 80% of the peak intensity in the maps.
the major and minor axes, including the position angle of the ellipse, fit the structures in the maps to within 5% of the peak intensity (cf. contours in Fig. 6) for both wavelengths.
4.5 Comparison to Kervella et al. 2014
The observations presented here are taken with the same instrument as was used in Kervella et al. (2014), but with a very different technique. Kervella et al. used direct imaging with adaptive optics as well as lucky imaging with very short integrations, selecting the best images. Deconvolution was used to further improve the resolution. This gives superb resolution but depends on the accuracy of the point-spread-function (PSF) calibration and on the Strehl ratio. The current data were taken with an aperture mask. This creates a Fizeau inteferometer, giving an exactly known PSF that is independent of the seeing (the seeing does affect the signal-to-noise), but with reduced throughput because of the mask geometry.
Comparing the two techniques, aperture masking has higher angular resolution: its PSF is somewhat smaller than the diffraction limit of a single mirror and deconvolves well because it is mathematically known. Both techniques require accurate calibration (PSF for direct imaging and visibility amplitude for aperture masking). Aperture masking filters out extended emission, while adaptive optics leaves extended pedestals in the PSF. The techniques are complementary in view of their very different strengths.
Comparing our Fig. 5 with Fig. 4 of Kervella et al. (2014) illustrates this. The current data have a slightly better resolution and show smaller details, but the extended disk is better seen in Kervella et al. Allowing for this, the agreement is very good. The aperture masking -band image shows two clumps either side of the star, indicating a better resolution than Kervella et al., but the extended disk emission is suppressed. The -band image, with its shorter -spacings, shows the disk better. The overall structures agree very well, validating the images of Kervella et al. (2014) and supporting the reliability of the image reconstruction.
5 Nature of the detected structure
The NACO observations are focused at 2.30 , where the detection of carbon monoxide emission666To our knowledge, there are currently no spectroscopic observations of the CO first-overtone band available in the literature. and scattered dust emission from the circumstellar envelope is expected, and at 3.74 , where continuum emission from the dusty environment probably is more prominent. The detected ellipsoidal structure is common in both wavelengths. This was confirmed by convolving the 2.30 maps with the resolution beam of SAM/NACO at 3.74 . The overall structure may be interpreted either as an elongated outflow or a toroidal structure, depending on the inclination of the object with respect to the line of sight and the stellar rotation axis. The interpretation of the hook-like features depends on the inclination of the object as well. A few possible scenarios on the origin of the structure detected around L Pup by SAM/NACO are presented below.
5.1 Circumstellar shell hypothesis
The structure detected at 2.30 is aspherical and extends on both sides to 65 mas on average (above 3% of the peak flux east–west), which in physical units translates in 4.2 au. The structure becomes an ellipsoid at the pseudo-continuum (3.74 ), with approximate sizes 140 mas wide in the east–west direction and 70 mas wide in the north–south direction (Fig. 6).
Bedding et al. (2002) compared the IRTS mid-infrared spectrum (recorded in 1995, pre-dimming) to a model of a silicate-rich, detached shell777Nevertheless, they advised to use the best-fitting parameters of their model with caution.. The inner and outer radii of that dust shell were and cm (or 48.1 and 86.9 au). At a distance of 64 pc, these dimensions translate into 075 and 135, and are placed within the extended envelope detected at 11.7 (1818; Jura et al. 2002). With a known expansion velocity888This is the outflow velocity of gas, not dust. These are probably only loosely coupled in a star with such a low mass-loss rate (; Schöier et al. 2004). of approximately 2.5 km s(Kerschbaum & Olofsson 1999), the outer edge of the shell would have reached that distance within 165 years. Nevertheless, such a shell is much larger than the structures detected by SAM and also lies outside the region imaged by Kervella et al. (2014).
There have been previous detections of molecular shells in the extended atmospheres of AGB stars (Le Bouquin et al. 2009; Perrin et al. 2004; Eisner et al. 2007; Thompson et al. 2002; Ohnaka et al. 2012; Wittkowski et al. 2007, 2011; Ireland et al. 2004a). The sizes of these shells have been estimated through star+molecular shell models (MOLspheres) or dynamic atmosphere models and are spherically symmetric. The observations presented here showed that the case of L Pup is different because no spherically-symmetric shell was detected.
5.2 Toroidal structure hypothesis
The molecular emission of L Pup has been extensively monitored in the past fifty years; it showed the variability of the outflow velocity (Kerschbaum et al. 1996; Kerschbaum & Olofsson 1999; Winters et al. 2002; Lepine et al. 1976; Winters et al. 2002, 2003; Haikala 1990; Haikala et al. 1994; Menten & Melnick 1991; Gomez Balboa & Lepine 1986; Lepine et al. 1978; Caswell et al. 1971; Knowles & Batchelor 1978; Balister et al. 1977; Schöier et al. 2004; González Delgado et al. 2003). González Delgado et al. (2003) reported that the interior of the circumstellar envelope of L Pup is probably complex, based on their SiO observations. However, the radio interferometric measurements by Schöier et al. (2004) of “thermal SiO” were unable to probe the direction of these flows, and a circularly symmetric envelope was assumed.
Recent observations (McIntosh & Indermuehle 2013) of the variation in the positioning of the velocity centroids of SiO maser emission over a period of four years, indicate an asymmetric nature of the gaseous environment of L Pup which may be associated with a rotating edge-on disk and bipolar ejecta (see Sect. 5.3). SiO maser emission probes the material much closer to the star than the “thermal” SiO emission. Kerschbaum & Olofsson (1999) suggested that the line profile seen in the carbon monoxide spectra is caused by a gaseous disk around the star, similar to the case of SV Psc (Klotz et al. 2012a). Bujarrabal et al. (2013) showed that similar CO line profiles can be found in post-AGB stars, and they assigned these to rotation in molecular disks. In fact, three of these objects – the Red Rectangle, AC Her, and HR 4049 – have slow expansion velocities, bipolar ejecta, and circumbinary dusty disks.
The rotational axis of L Pup is unknown, which allows for a range of scenarii on the nature of the detected structures. However, the images of Kervella et al. (2014) are strongly suggestive of a rotating east-west disk, which would indicate a north-south rotation axis.
Considering that the bilateral clumps are observed farther out than the extended circumstellar emission (≳30 mas) in all 2.30 maps, but are absent from the 3.74 maps999Due to the lower resolution at the wavelength., it is possible that the edge-on structure detected with SAM/NACO is in fact a torus of gas interwoven with dust. The origin of the hook-like patterns is quite perplexing. These emerge at levels 5% of the peak intensity and extend symmetrically southwest and northeast in all the SAM maps. Assuming that the inclination of L Pup in our line of sight is perpendicular to the north–south axis, that is , the possibility of an inclined, edge-on spiral or helical structure, similar to those observed in AGB stars (e.g., AFGL3068; Mauron & Huggins 2006) and planetary nebulae (e.g., PNG356.8+03.3; Sahai et al. 2011), should not be excluded.
5.3 Stellar and torus parameters
The SiO maser spectra of McIntosh & Indermuehle (2013) show an unusual profile, with two peaks separated by 20 km s, in addition to the more usual narrow peak at the systemic velocity. A rotating gaseous disk would fit this profile, where the separated peaks show the approaching and receding side of the disk. The central velocity component may trace the stellar wind might be emission from the disk in front of the star, where the maser is amplifying the stellar emission. This model gives a rotational velocity of the disk of , where is the inclination with respect to the plane of the sky.
The two clumps seen on either side of the star are consistent with the expected structure seen at the inner radius of the disk. This gives an inner radius of 4.2 au (Kervella et al. derived 6 au from their models).
Assuming Keplerian velocities, this gives a mass of the star of . The optical obscuration events, and the images of Kervella et al. (2014), indicate an inclination , although given the size of the star, a significant range of is possible. However, using a maximum range of and still implies a stellar mass .
This mass is lower than that derived by Kervella et al. (2 M). They located the star at the beginning of the AGB (meant is the start of the TP-AGB) where the mass is still close to the initial mass. However, the short period of the star of 140 days together with a location 0.5 mag below the Mira PL relation (Bedding & Zijlstra 1998; Bedding et al. 2002) suggest a low initial mass, as found here.
5.4 Binarity hypothesis
The possible signature of a binary companion was investigated in the SAM/NACO data, although the binary nature of L Pup was previously doubted by Jorissen et al. (2009). Using the geometric model-fitting tool LITpro101010LITpro software available at http://www.jmmc.fr/litpro (Tallon-Bosc et al. 2008), an additional, unresolved component, with a flux ratio ranging from 1/100 to 1/500, was included in the model presented in Sect. 3 and was positioned randomly within the NACO field of view. All attempts failed to converge (), suggesting that no signal from a binary companion was found in the data.
The SAM and SAMPol diffraction-limited imaging revealed the complex circumstellar nature of the envelope of L Pup in both gaseous emission and dust continuum. The maps at 2.30 show an asymmetric structure extending to 65 mas in the west at levels above 5% of the peak intensity and a hook-like extension (Figs. 5 and 6). It resembles a linear structure from the star, whose extended envelope is resolved at 30 mas, although the stellar photosphere should be smaller than that. A fainter, hook-like feature emerges in the northeastern direction. The overall shape of the structure is similar at 3.74 , where it resembles an ellipsoid with sizes 140 mas and 70 mas along the major and minor axis, respectively. This warm dust component extends farther out into a fainter, helical feature.
Both structures are oriented in the east–west axis, agreeing with Kervella et al. (2014), but different from previous aperture-masking observations: Ireland et al. (2004b) positioned their mask in a north–south orientation (P.A.=180°) and failed to detect the equatorial structure.
This ellipsoidal structure, found in pseudo-continuum (3.74 ), is the less resolved counterpart of the more detailed, bilateral structure (2.30 ), since the dimensions of both are similar. The structure is contained within the extended silicate-rich dust zone (Bedding et al. 2002). As shown in the near-infrared maps, it appears that the structure detected in L Pup is seen nearly edge on. We interpret the structure as a torus with inner radius 4.2 au around L Pup. This torus is probably the cause of the long-lived dimming event (Bedding et al. 2002). It might be that a more extended dust reservoir is located between the torus and the extended envelope (Jura et al. 2002), although it is currently outside the field of view of SAM/NACO.
SiO maser profiles indicate that the torus may be rotating at 10 km s. Using the inner radius of 4.2 au, we derive a stellar mass of 0.7 M with an upper limit of 1 M. This low mass is consistent with the short period of L Pup.
Acknowledgements.We would like to thank the NACO operational team for their assistance during our observations, as well as K. Hebden, W. Nowotny and T. Lebzelter for their valuable comments. FL, DK and JH were supported by the Austrian science fund FWF project number AP23006. CP was supported by the Belgian Federal Science Policy Office via the PRODEX Programme of ESA. JK was supported by the French ANR POLCA project (Processing of pOlychromatic interferometriC data for Astrophysics, ANR-10-BLAN-0511). We acknowledge with thanks the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. This research has made use of the Jean-Marie Mariotti Center LITpro service co-developed by CRAL, LAOG and FIZEAU.
- Balister et al. (1977) Balister, M., Batchelor, R. A., Haynes, R. F., et al. 1977, MNRAS, 180, 415
- Bedding & Zijlstra (1998) Bedding, T. R. & Zijlstra, A. A. 1998, ApJ, 506, L47
- Bedding et al. (2002) Bedding, T. R., Zijlstra, A. A., Jones, A., et al. 2002, MNRAS, 337, 79
- Bujarrabal et al. (2013) Bujarrabal, V., Alcolea, J., Van Winckel, H., Santander-García, M., & Castro-Carrizo, A. 2013, A&A, 557, A104
- Buscher (1994) Buscher, D. F. 1994, in IAU Symposium, Vol. 158, Very High Angular Resolution Imaging, ed. J. G. Robertson & W. J. Tango, 91
- Caswell et al. (1971) Caswell, J. L., Robinson, B. J., & Dickel, H. R. 1971, Astrophys. Lett., 9, 61
- Eisner et al. (2007) Eisner, J. A., Graham, J. R., Akeson, R. L., et al. 2007, ApJ, 654, L77
- Gomez Balboa & Lepine (1986) Gomez Balboa, A. M. & Lepine, J. R. D. 1986, A&A, 159, 166
- González Delgado et al. (2003) González Delgado, D., Olofsson, H., Kerschbaum, F., et al. 2003, A&A, 411, 123
- Haikala (1990) Haikala, L. K. 1990, A&AS, 85, 875
- Haikala et al. (1994) Haikala, L. K., Nyman, L.-A., & Forsstroem, V. 1994, A&AS, 103, 107
- Ireland et al. (2006) Ireland, M. J., Monnier, J. D., & Thureau, N. 2006, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 6268
- Ireland et al. (2004a) Ireland, M. J., Scholz, M., & Wood, P. R. 2004a, MNRAS, 352, 318
- Ireland et al. (2004b) Ireland, M. J., Tuthill, P. G., Bedding, T. R., Robertson, J. G., & Jacob, A. P. 2004b, MNRAS, 350, 365
- Jorissen et al. (2009) Jorissen, A., Frankowski, A., Famaey, B., & van Eck, S. 2009, A&A, 498, 489
- Jura et al. (2002) Jura, M., Chen, C., & Plavchan, P. 2002, ApJ, 569, 964
- Kerschbaum & Olofsson (1999) Kerschbaum, F. & Olofsson, H. 1999, A&AS, 138, 299
- Kerschbaum et al. (1996) Kerschbaum, F., Olofsson, H., & Hron, J. 1996, A&A, 311, 273
- Kervella et al. (2014) Kervella, P., Montargès, M., Ridgway, S. T., et al. 2014, A&A, 564, A88
- Klotz et al. (2012a) Klotz, D., Sacuto, S., Kerschbaum, F., et al. 2012a, A&A, 541, A164
- Klotz et al. (2012b) Klotz, D., Sacuto, S., Paladini, C., Hron, J., & Wachter, G. 2012b, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8445
- Knapp et al. (1999) Knapp, G. R., Dobrovolsky, S. I., Ivezić , Z., et al. 1999, A&A, 351, 97
- Knowles & Batchelor (1978) Knowles, S. H. & Batchelor, R. A. 1978, MNRAS, 184, 107
- Lacour et al. (2011) Lacour, S., Tuthill, P., Ireland, M., Amico, P., & Girard, J. 2011, The Messenger, 146, 18
- Le Bouquin et al. (2009) Le Bouquin, J.-B., Lacour, S., Renard, S., et al. 2009, A&A, 496, L1
- Lenzen et al. (2003) Lenzen, R., Hartung, M., Brandner, W., et al. 2003, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. M. Iye & A. F. M. Moorwood, Vol. 4841, 944–952
- Lepine et al. (1976) Lepine, J. R. D., Paes de Barros, M. H., & Gammon, R. H. 1976, A&A, 48, 269
- Lepine et al. (1978) Lepine, J. R. D., Scalise, Jr., E., & Le Squeren, A. M. 1978, ApJ, 225, 869
- Magalhaes et al. (1986) Magalhaes, A. M., Codina-Landaberry, S. J., Gneiding, C., & Coyne, G. V. 1986, A&A, 154, 1
- Mauron & Huggins (2006) Mauron, N. & Huggins, P. J. 2006, A&A, 452, 257
- McIntosh & Indermuehle (2013) McIntosh, G. C. & Indermuehle, B. 2013, ApJ, 774, 21
- Meisner (2008) Meisner, J. 2008, in The Power of Optical/IR Interferometry: Recent Scientific Results and 2nd Generation, ed. A. Richichi, F. Delplancke, F. Paresce, & A. Chelli, 481
- Menten & Melnick (1991) Menten, K. M. & Melnick, G. J. 1991, ApJ, 377, 647
- Murakawa et al. (2005) Murakawa, K., Suto, H., Oya, S., et al. 2005, A&A, 436, 601
- Narayan & Nityananda (1986) Narayan, R. & Nityananda, R. 1986, ARA&A, 24, 127
- Norris et al. (2012) Norris, B. R. M., Tuthill, P. G., Ireland, M. J., et al. 2012, Nature, 484, 220
- Ochsenbein & Halbwachs (1982) Ochsenbein, F. & Halbwachs, J. L. 1982, A&AS, 47, 523
- Ohnaka et al. (2012) Ohnaka, K., Hofmann, K.-H., Schertl, D., et al. 2012, A&A, 537, A53
- Perrin et al. (2004) Perrin, G., Ridgway, S. T., Mennesson, B., et al. 2004, A&A, 426, 279
- Rousset et al. (2003) Rousset, G., Lacombe, F., Puget, P., et al. 2003, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. P. L. Wizinowich & D. Bonaccini, Vol. 4839, 140–149
- Sahai et al. (2011) Sahai, R., Morris, M. R., & Villar, G. G. 2011, AJ, 141, 134
- Schöier et al. (2004) Schöier, F. L., Olofsson, H., Wong, T., Lindqvist, M., & Kerschbaum, F. 2004, A&A, 422, 651
- Sivia (1987) Sivia, D. S. 1987, PhD thesis, , Cambridge Univ., (1987)
- Skilling & Bryan (1984) Skilling, J. & Bryan, R. K. 1984, MNRAS, 211, 111
- Tallon-Bosc et al. (2008) Tallon-Bosc, I., Tallon, M., Thiébaut, E., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7013
- Thiébaut (2008) Thiébaut, E. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7013
- Thompson et al. (2002) Thompson, R. R., Creech-Eakman, M. J., & van Belle, G. T. 2002, ApJ, 577, 447
- Tuthill et al. (2010) Tuthill, P., Lacour, S., Amico, P., et al. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7735
- Tuthill et al. (2002) Tuthill, P. G., Monnier, J. D., Danchi, W. C., Hale, D. D. S., & Townes, C. H. 2002, ApJ, 577, 826
- Tuthill et al. (2000) Tuthill, P. G., Monnier, J. D., Danchi, W. C., Wishnow, E. H., & Haniff, C. A. 2000, PASP, 112, 555
- van Leeuwen (2007) van Leeuwen, F. 2007, A&A, 474, 653
- van Winckel (2003) van Winckel, H. 2003, ARA&A, 41, 391
- Winters et al. (2003) Winters, J. M., Le Bertre, T., Jeong, K. S., Nyman, L.-Å., & Epchtein, N. 2003, A&A, 409, 715
- Winters et al. (2002) Winters, J. M., Le Bertre, T., Nyman, L.-Å., Omont, A., & Jeong, K. S. 2002, A&A, 388, 609
- Wittkowski et al. (2011) Wittkowski, M., Boboltz, D. A., Ireland, M., et al. 2011, A&A, 532, L7
- Wittkowski et al. (2007) Wittkowski, M., Boboltz, D. A., Ohnaka, K., Driebe, T., & Scholz, M. 2007, A&A, 470, 191
Appendix A Visual light curve
Figure 7 shows the visual light curve of L Pup from 1985 to 2013 (data from AAVSO), where the observing periods of MAPPIT, VINCI, and NACO are clearly marked. An average light curve is overplotted to guide the eye.