Issue |
A&A
Volume 668, December 2022
|
|
---|---|---|
Article Number | A119 | |
Number of page(s) | 6 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/202244921 | |
Published online | 13 December 2022 |
Spatially resolved mid-infrared observations of the circumstellar environment of the born-again object FG Sge★
Instituto de Astrofísica, Universidad Andrés Bello,
Fernández Concha 700,
Las Condes, Santiago, Chile
e-mail: k1.ohnaka@gmail.com
Received:
8
September
2022
Accepted:
31
October
2022
Context. FG Sge has evolved from the hot central star of the young planetary nebula Hen 1–5 to a G–K supergiant in the last 100 yr. It is one of the three born-again objects that have been identified as of yet, and they are considered to have undergone a thermal pulse in the post-asymptotic giant branch evolution.
Aims. We present mid-infrared spectro-interferometric observations of FG Sge and probe its dusty environment.
Methods. FG Sge was observed with MIDI at the Very Large Telescope Interferometer at baselines of 43 and 46 m between 8 and 13 µm.
Results. The circumstellar dust environment of FG Sge was spatially resolved, and the Gaussian fit to the observed visibilities results in a full width at half maximum of ~10.5 mas. The observed mid-infrared visibilities and the spectral energy distribution can be fairly reproduced by optically thick (τV ≈ 8) spherical dust shell models consisting of amorphous carbon with an inner radius rin of ~30 R★ (corresponding to a dust temperature of 1100 ± 100 K). The dust shell is characterized with a steep density profile proportional to r−3.5±0.5 from the inner radius rin to (5–10) × rin, beyond which it changes to r−2. The dust mass is estimated to be ~ 7 × 10−7 M⊙, which translates into an average total mass-loss rate of ~ 9 × 10−6 M⊙ yr−1 as of 2008 with a gas-to-dust ratio of 200 being adopted. In addition, the 8–13 µm spectrum obtained with MIDI with a field of view of 200 mas does not show a signature of the polycyclic aromatic hydrocarbon (PAH) emission, which is in marked contrast to the spectra taken with the Spitzer Space Telescope six and 20 months before the MIDI observations with wide slit widths of 3″.6–10″. This implies that the PAH emission originates from an extended region of the optically thick dust envelope.
Conclusions. The dust envelope of FG Sge is much more compact than that of the other born-again stars’ Sakurai’s object and V605 Aql, which might reflect the difference in the evolutionary status. The PAH emission from the extended region of the optically thick dust envelope likely originates from the material ejected before the central star became H-deficient, and it may be excited by the UV radiation from the central star escaping through gaps among dust clumps and/or the bipolar cavity of a disk-like structure.
Key words: stars: AGB and post-AGB / infrared: stars / circumstellar matter / stars: carbon / techniques: interferometric / stars: individual: FG Sge
© K. Ohnaka and B. A. Jara Bravo 2022
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.
1 Introduction
FG Sge is the central star of the young planetary nebula Hen 1–5 and a very rare object whose evolution can be observed in a lifetime of a human being. In the last 100 yr, it has rapidly evolved redward on the Hertzsprung-Russell (H-R) diagram, from the hot central star of the planetary nebula (effective temperature Teff ≈ 40000 K) to a G supergiant (Teff ≈ 5500 K, Jeffery & Schönberner 2006). The rapid evolution of FG Sge is considered to be caused by a late thermal pulse (LTP) that occurs after the asymptotic giant branch (AGB), when the star is evolving blue-ward and almost horizontally on the H-R diagram. The nuclear energy produced by the thermal pulse results in the expansion of the thin outer layers on a timescale of 50 to a few 100 yr, and the star evolves rapidly redward on the H-R diagram, turning into a so-called born-again object (e.g., Langer et al. 1974; Iben 1976; Blocker 2001). When the temperature of the expanding outer layers becomes low enough for convection to penetrate to deep layers, the nuclear-processed material, such as C, He, and s-process elements, is dredged up to the surface. Since the mass of the H-rich outer layers is already small, this dredge-up process makes the atmosphere H-deficient and C-rich. After the surface temperature reaches its minimum of several 1000 K, the star evolves again blueward on a timescale of 103–104 yr finally to a white dwarf.
The final flash episode is so brief that objects such as FG Sge are very rare. Only three objects are known to date: FG Sge, Sakurai’s object (V4334 Sgr), and V605 Aql. The evolution of FG Sge back to the AGB over ~ 100 yr is consistent with the LTP scenario, while the return of Sakurai’s object and V605 Aql to the AGB is even more rapid on the order of several years. They are considered to have experienced the so-called very late thermal pulse (VLTP) – a thermal pulse that occurs when the star is evolving downward to a white dwarf on the H-R diagram (Clayton et al. 2006). In addition, Reindl et al. (2017) recently found that SAO244567, the central star of the Stingray nebula, has been rapidly returning to the AGB since 2002, and they conclude that it has undergone an LTP. The final flash episode is considered to be a fairly common phenomenon, which occurs in 10-20% of stars evolving off the AGB (Blocker 2001; Herwig 2005). Therefore, studies of final flash objects provide a unique opportunity to obtain insight into the evolution of post-AGB stars.
Although FG Sge is discussed as an LTP object in the literature (e.g., Schönberner 2008), the evolution of its chemical composition is not yet fully understood. The detailed spectroscopic analysis of Jeffery & Schönberner (2006) shows that the atmosphere of FG Sge became H-deficient some time between 1960 and 1995, and the detection of C2 (Iijima & Strafella 1993; Kipper & Kipper 1993) reveals that the atmosphere is C-rich. However, this photospheric chemical composition of FG Sge and its evolution on the H-R diagram cannot be consistently explained by the present stellar evolutionary theory as discussed by Jeffery & Schönberner (2006). The calculations of born-again giants (Blöcker 2001) predict that the atmosphere becomes H-deficient only after the star begins to evolve blueward again (i.e., when the temperature is increasing again). However, FG Sge became H-deficient between 1960 and 1995, when the surface temperature was still decreasing. In the case of a VLTP, hydrogen in the outer layers becomes significantly depleted due to the hydrogen burning during the helium flash (Herwig et al. 1999). This means that the star is already H-deficient when it begins to evolve redward. However, FG Sge was H-normal in 1960 when it was evolving redward. The alternative double-loop scenarios proposed by Lawlor & MacDonald (2003) and Miller Bertolami et al. (2006), in which FG Sge, Sakurai’s object, and V605 Aql can be considered to represent different evolutionary stages of the same evolutionary path, cannot entirely explain the observed chemical composition of FG Sge either, as Jeffery & Schönberner (2006) have discussed.
Another remarkable change occurred in FG Sge in 1992, when its visual brightness decreased suddenly by ∆V ≈ 5 magnitude and recovered gradually over the course of ~1 yr. Since then, FG Sge has shown sudden, unpredictable deep declines at an interval of 1–2 yr; although, such events have become so frequent that the central star has been seldom visible since 2008 (Arkhipova et al. 2022). At the same time, it started to show a noticeable infrared excess characterized by carbonaceous dust with ~1000 K (Gehrz et al. 2005). The light curve of FG Sge is similar to those of R CrB stars – a class of F–G supergiants in which optically thick dust clouds are ejected in random directions – and an occasional, deep decline is observed when a cloud happens to form in the line of sight (Clayton 1996). While some authors propose that FG Sge is undergoing such a random dust cloud formation (Gonzalez et al. 1998), Gehrz et al. (2005) have argued for a continuous outflow with abrupt, episodic changes of the mass-loss rate. In either case, however, it is unknown where and how dust can form around FG Sge, whose present temperature (4500–5500 K) is deemed to be too high for dust to condensate.
To obtain a clearer picture of the dusty environment of FG Sge, high spatial resolution mid-infrared observations are effective. In this paper, we present the first spatially resolved mid-infrared (8–13 µm) observations of FG Sge and radiative transfer modeling of its circumstellar envelope.
2 Observations
FG Sge was observed with the MIDI instrument (Leinert et al. 2003) at ESO’s Very Large Telescope Interferometer (VLTI) on 2008 June 22 (MJD 54639) using the unit telescope (UT) configuration UT2–UT3 with projected baseline lengths of 43–46 m (Program ID: 081-D0224, P.I.: K. Ohnaka). A prism with a spectral resolution of λ/∆λ ≃ 30 (at 10 µm) was used to obtain spectrally dispersed fringes in the N band between 8 and 13 µm. The data were taken in SCI-PHOT mode, in which photometric data of each telescope are recorded simultaneously with interferometric data. A detailed description of the observing procedure is given in Przygodda et al. (2003), Leinert et al. (2004), and Chesneau et al. (2005). Table 1 summarizes the observations’ FG Sge and the calibrator.
We used the MIA+EWS package ver.2.01 to reduce the MIDI data (Leinert et al. 2004; Jaffe 2004). We adopted a uniform-disk diameter of 3.1 mas for the calibrator HD 168723 (Bourges et al. 2017). It should be noted that there are significant differences in airmass among the data sets. To avoid systematic errors, we calibrated the data sets of FG Sge with those of the calibrator taken at a similar airmass: the data set SCI1 was calibrated with CAL1, while the data sets SCI2 and SCI3 were calibrated with CAL2.
We assessed the data quality and derived the calibrated visibilities as described by Ohnaka et al. (2008). Because it is difficult to estimate the errors in the calibrated visibilities properly from two calibrator data sets taken at very different airmasses, we assumed a total relative error of 10% as adopted by Zhao-Geisler et al. (2011) for SCI-PHOT mode data. Because the projected baseline length and position angle of SCI2 and SCI3 are very close, we averaged the calibrated visibilities of these two data sets. The visibilities obtained with MIA and EWS agree well, and we only provide the results derived with EWS in the discussion below.
The absolutely calibrated N-band spectra of FG Sge were obtained from the same three SCI-CAL pairs taken at similar airmasses (SCI1-CAL1, SCI2-CAL2, and SCI3-CAL2), using the method described in Ohnaka et al. (2007). The average and the standard deviation of the three spectra were taken as the final spectrum and its error, respectively.
We also obtained the 19.5 µm flux of FG Sge with the mid-infrared imager and spectrograph VLT/VISIR (Lagage et al. 2004). The VISIR observations of FG Sge and a photometric standard star, HD 173780, were carried out on 2008 May 30, about three weeks before the MIDI observations, with the Q3 filter (central wavelength = 19.5 µm, filter full width at half maximum = 0.4 µm) and a pixel scale of 07075. The data with chopping and nodding were first reduced using the VISIR pipeline ver.3.0.0 to obtain the images with the sky background almost subtracted. From these images, we further subtracted the residual sky background measured in an annular region with an inner and outer radius of 1″.5 and 2″.25, respectively. The flux was obtained within an aperture radius of 175 both for FG Sge and HD 173780, and the data of FG Sge were flux-calibrated with a Q3-band flux of 1.7276 Jy for HD 1737802. The error in the calibrated flux of FG Sge was estimated by the quadratic sum of the photon noise of the source and that of the background. The observed 19.5 µm flux of FG Sge is 6.2 ± 0.49 Jy.
MIDI observations of FG Sge and the calibrator: time of observation (coordinated universal time=UTC), projected baseline length Bp, position angle of the projected baseline on the sky (P.A.), seeing in the visible, and the airmass at the times of the observations of FG Sge and the calibrator.
![]() |
Fig. 1 VLTI/MIDI observations of FG Sge. a) Observed visibilities of the data set SCI1 (red dots) and the average of SCI2 and SCI3 (blue dots). The error bars are only shown for the data set SCI1 for visual clarity. b) FWHM obtained by a Gaussian fit to the observed visibilities. |
3 Results
Figure la shows the measured N-band visibilities, which slightly increase from 0.8 at 8 µm to 0.9 at 13 µm. Figure 1b shows the full width at half maximum (FWHM) obtained by the Gaussian fitting, which suggests that the observed visibilities correspond to an approximately constant FWHM of ~10.5 mas. The observed visibilities suggest that the dust envelope of FG Sge is compact compared to the other born-again objects’ Sakurai’s object and V605 Aql. Chesneau et al. (2009) showed that the MIDI visibilities of Sakurai’s object obtained at baseline lengths of 42–46 m – the same as in our observations of FG Sge – range from 0.05 to 0.2, which is much lower (thus much more extended) than the 0.8–0.9 obtained for FG Sge; although, the distance to two objects is comparable (FG Sge: 3.0 kpc, Sakurai’s object: 3.8 kpc, Evans et al. 2020). Clayton et al. (2013) report that V605 Aql was so large that it was resolved out in the MIDI observations with UT2 and UT3 (the same as used for FG Sge) in spite of its distance of 4.6 kpc. This is primarily because the effective temperature of FG Sge (4500–5500 K, see Sect. 4) is noticeably lower than those of Sakurai’s object (12 000 K, van Hoof et al. 2007) and V605 Aql (95 000 K, Clayton et al. 2006). In addition, the luminosity of these latter objects (~104 L⊙, Schönberner 2008; Chesneau et al. 2009) is higher than that of FG Sge ( L⊙, Sect. 4), which can also make them appear larger than FG Sge.
Figure 2 shows the MJDI spectrum (red solid line with the error bars) and the 19.5 µm flux measured with VISIR (green filled diamond). To complement the MIDI N-band spectrum and the VISIR Q-band photometric data, we downloaded the mid-infrared spectrum taken with the InfraRed Spectrograph3 (IRS, Houck et al. 2004) onboard the Spitzer Space Telescope (Werner et al. 2004)4. The data were taken on 2007 December 11 (MJD 54446) in the short and long low mode with spectral resolutions of 64–128, covering from 5 to 40 µm (Program ID: 40061, P.I.: A. Evans).
The 19.5 µm flux measured with VISIR is in good agreement with the Spitzer data, which were taken approximately six months before the MIDI observations. The MIDI spectrum is also consistent with the Spitzer/IRS spectrum longward of ~10 µm. However, the Spitzer/IRS spectrum shows a bump at ~8 µm, which is not seen in the MIDI spectrum. Evans et al. (2015) reported this feature and the one at ~6.2 µm as the unidentified infrared emission, which is commonly attributed to polycyclic aromatic hydrocarbons or PAH (e.g., Allamandola et al. 1999). Evans et al. (2015) show that the Spitzer spectrum taken on 2006 October 24 (MJD 54032) also exhibit the PAH features at 6.2 and 8 µm, which suggests that the PAH emission was stably present over a little more than a year from October 2006 to December 2007. Therefore, it is not very likely – if not entirely excluded – that the PAH emission disappeared entirely within six months between the Spitzer observation in December 2007 and our MIDI observations in June 2008. The absence of the PAH emission in the MIDI spectrum can be due to a smaller field of view of 200 mas of MIDI compared to the slit widths of 3″.6-10″.7 of Spitzer/IRS, which suggests that the PAH emission may originate from an extended region. Given that the central star is now H-deficient, the recently rejected inner circumstellar matter should also be H-deficient, from which PAH emission is not expected. Therefore, the PAH emission likely originates from the older, outer circumstellar material ejected before the star became H-deficient. This is consistent with the above observations. However, the UV radiation required to give rise to the PAH emission cannot penetrate the optically thick dust envelope suggested by our modeling presented below. We return to this point in Sect. 4.
![]() |
Fig. 2 Observed mid-infrared spectrum of FG Sge. The solid red line with the error bars shows the MIDI spectrum. Black dots correspond to the Spitzer/IRS spectrum. The green diamond is the VISIR Q3-band (19.5 µm) flux. |
4 Modeling of the dust envelope
We carried out radiative transfer modeling with DUSTY (Ivezic et al. 1999) for the observed N-band visibilities and the spectral energy distribution (SED). The circumstellar envelope of FG Sge is considered to be complex, given the R CrB-like light curve that suggests the formation of clumpy dust clouds or episodic mass loss. However, the present MIDI data are insufficient for constraining such complex structures. Therefore, we assumed spherical shell models to constrain the global properties of the dusty environment of FG Sge.
In the present work, we adopted a distance of kpc derived by Chornay & Walton (2021)5 based on the Gaia Early Data Release 3 (Gaia Collaboration 2021). The bolometric flux of the central star was obtained from the SED obtained before the onset of the mass ejection in 1992. Jurcsik & Montesinos (1999) presented the photometric data from the V to N band obtained from June to August in 1983, which they fit with a photospheric model with Teff = 5500 K and log g = 1.5. We obtained a bolometric flux of 1. 72 × 10−11 W m−2 by integrating this SED, which translates into a luminosity of
L⊙ with the distance of
kpc. It should be noted that the large uncertainty in the distance (hence in the luminosity) does not affect the fit to the SED and visibility, because the luminosity of the model star increases as ∝ d2, but the flux predicted at the distance of the Earth decreases as ∝ d−2, resulting in the same predicted flux. Likewise, the inner radius of the dust shell increases as
, but its angular size remains unaffected because it follows ∝ d−1.
To obtain the SED observed as contemporaneously as possible with our MIDI observations, we used the visible (V, R, I, and Rc bands) photometric data taken on MJD 54642 (just three days after our MIDI observations) by Arkhipova et al. (2009) and the infrared (JHKLM) photometric data taken on MJD 54659.5 (20 days after the MIDI observations) by Shenavrin et al. (2011). As presented above, the MIDI spectrum and VISIR 19.5 µm flux are in agreement with the Spitzer spectrum obtained in December 2007 except for the PAH emission. Therefore, we used these data to cover the mid-infrared region of the SED. We adopted E(g − r) = 0.23 based on the 3D extinction map of Green et al. (2019)6 for the distance of 3.0 kpc and converted it to E(B − V) = 0.202 using the relations given in Schlafly & Finkbeiner (2011). This translates into AV = 0.63 with Rv = 3.1 assumed. The observed SED was dereddened by applying the wavelength dependence of the interstellar extinction from Cardelli et al. (1989).
The effective temperature of the central star is controversial. Jeffery & Schönberner (2006) concluded that FG Sge had been maintaining Teff ≈ 5500 K since 1980 (as of 1995) and, therefore, proposed that this Teff might be the lowest limit in its redward evolution. On the other hand, Fadeyev (2019) derived Teff = 4445 K based on theoretical modeling of FG Sge’s pulsation. We treated Teff as a free parameter and calculated models with Teff = 4500, 5000, and 5500 K. Combined with the luminosity of 4850 L⊙, these values of Teff correspond to a stellar radius of 115, 93, and 77 R⊙, respectively.
Given the C-rich nature of the photosphere and the circum-stellar envelope, we assumed amorphous carbon for the dust grains around FG Sge. We adopted a grain size of 0.01 µm as concluded by Arkhipova et al. (2022) to explain the SED before 2019. The density distribution of dust grains was assumed to be ∝ r−p, where p is a free parameter, and characterized by the dust condensation temperature (Tc), which determines the inner radius of the envelope (rin), and the optical depth in the radial direction at 0.55 µm (τV). The outer radius (rout) of the circumstellar envelope was estimated in the following manner. Gonzalez et al. (1998) measured an outflow velocity of ~200 km s−1 from the NaI D lines. Because the mass ejection started in August 1992, the material should have traveled 1.0 × 1016 cm (670 AU) at the time of our MIDI observations in June 2008. On the other hand, DUSTY models show that the inner radius of the envelope ranges from 25 to 40 R★ for Teff = 4500 to 5500 K for a condensation temperature of ~1000 K. Using the above values of R★, the inner radius is estimated to be (1 – 3) × 1014 cm. Therefore, the maximum distance that the ejected material traveled (1.0 × 1016 cm) corresponds to 33–100 rin. We adopted a fixed value of 70 rin.
We computed models with Teff = 4500, 5000, and 5500 K, p = 2.0 … 4.0 (∆p = 0.5), Tc = 800 … 1600 K (∆Tc = 200 K), and τV = 5 … 30 (∆τV = 5). However, we found no model that can reasonably reproduce the observed SED and N-band visibilities simultaneously. For example, the models with the density profile ∝ r−2, which corresponds to stationary mass loss with a constant outflow velocity, can fit the observed SED very well, but they systematically predict the N-band visibility to be much lower (≲0.5) than the observed data, because the dust envelope appears to be much more extended than the observed data suggest. The models with steeper density profiles can explain the observed N-band visibilities. However, these models predict the flux longward of ~10 µm to be too low compared to the observed data, because there is too little cold material in the outer region of the envelope.
It is possible that the density profile is steeper than the r−2 in the inner region of the dust envelope before the velocity reaches its terminal velocity. We attempted to explained the observed data with models in which the density profile falls as r−p from the inner radius rin to a certain radius rtr and then changes to r−2. Because the above model grid shows that τV should be below ~10 to explain the visible part of the observed SED, we calculated models with the following parameter ranges: Teff (K) = 4500, 5000, and 5000; τV = 4.0 … 10.0 (∆τV = 2.0); Tin (K) = 800 … 1600 (∆Tin = 100 K); p = 2.5 … 4.5 (∆p = 0.5); and rtr/rin = 2.5, 5, 10, 15, and 20. Figure 3 shows one of the best-fit models, which reproduces the observed SED and N-band visibilities fairly well. The model is characterized by Teff = 5000 K, τV = 8, Tin = 1100 K (rin = 33 R★), p = 3.5, and rtr = 5 rin. The model still predicts the flux to be somewhat lower than the observed data longward of 20 µm, and the model visibilities are slightly lower than the MIDI data. Decreasing p and/or rtr makes the fit to the SED better, but such models predict the N-band visibilities to be even lower, because their density profiles become similar to r−2. Our model grid constrains the parameter ranges as follows: τV = 8 ± 2, Tin = 1100 ± 100 K, p = 3.5 ± 0.5, and rtr = 10 ± 5. We cannot constrain Teff of the central star within the adopted range between 4500 and 5500 K, because it is heavily obscured.
Recently, Arkhipova et al. (2022) modeled the SED observed from 0.4 to 5.0 µm in 2019 with a density profile of ∝ r−2 and obtained an optical depth of 4.5 at 0.55 µm, which is significantly lower than the 8 ± 2 derived in the present work. This is because they modeled the photometric data taken in 2019, when the dust clouds temporarily cleared. The central star became visible with V ≈ 13, which is ~4 magnitude fainter than the V ≈ 9 before the start of the dust formation in 1992. With τV = 0.921 × ∆V (if scattered light and dust thermal emission are ignored), τV is estimated to be 3.7, which is in broad agreement with 4.5 as derived from their modeling. As mentioned above, the models with the density profile r−2 can explain the observed SED, but not the N-band visibilities. Given that Arkhipova et al. (2022) only modeled the SED taken at a different epoch (May 2019) from ours, we cannot draw the conclusion that the difference in the density profile reflects the lack of spatially resolved data in the modeling or a temporal variation in the density profile.
The remaining discrepancy between the best-fit models and the observed data is likely due to the complex, clumpy nature of the circumstellar environment of FG Sge as inferred from the irregular, R CrB-like declines in the light curve. Furthermore, a disk-like or ring-like structure may also be present, as observed in the other two born-again objects’ Sakurai’s object and V605 Aql (Chesneau et al. 2009; Hinkle & Joyce 2014; Hinkle et al. 2020; Clayton et al. 2013). As described in Sect. 2, if the PAH emission indeed originates from the outer region of the circumstellar envelope, the UV radiation should escape the optically thick envelope. The clumpy and/or disk-like structures can provide a natural explanation for the UV radiation to escape and excite the carrier of the PAH emission even if it is located farther away from the star.
![]() |
Fig. 3 Comparison of a spherical shell model with the observed SED and N-band visibilities of FG Sge. The model is characterized with Tin = 1100 K, a piecewise power-law density profile with r−3.5 (rin ≤ r ≤ rtr with rin = 33 R★ and rtr = 5 rin) and r−2 (r > rtr), τV = 8, and Teff = 5000 K. a) Observed SED. Orange filled triangles correspond to visible photometric data from Arkhipova et al. (2009). Orange filled diamonds are the infrared photometric data from Shenavrin et al. (2011). The photometric data were dereddened with AV = 0.63. The thick solid red line is the Spitzer/IRS spectrum. The thin solid blue line corresponds to the MIDI spectrum. The blue cross corresponds to VISIR photometry. The DUSTY model is shown with the solid black line. b) Observed visibilities are shown in the same manner as in Fig. 1a. The solid red and blue lines represent the model visibilities predicted for the 43.0 and 46.0 m baselines, respectively. |
5 Discussion and concluding remarks
The dust mass obtained from the model shown in Fig. 3 is 7.4 × 10−7 M⊙, adopting a bulk density of 2.25 g cm−3 for amorphous carbon (Gilman 1974). If we assume a gas-to-dust ratio of 200, the total envelope mass is 1.5 × 10−4 M⊙. Given that the mass ejection started in 1992, the average mass-loss rate at the time of the MIDI observations in 2008 is 9.3 × 10−6 M⊙ yr−1. The dust mass of FG Sge seems to be lower than 6 × 10−5 M⊙ as found in Sakurai’s object (Chesneau et al. 2009) and ~ 2 × 10−3 M⊙ as derived for V605 Aql (Clayton et al. 2013). The difference in the ejected mass may be related to the difference in the nature between FG Sge and the other two final-flash objects – Sakurai’s object and V605 Aql are VLTP objects, while FG Sge is likely to be an LTP object. However, follow-up observations of FG Sge would be necessary to study the long-term variation of the mass of the ejected dust.
As the light curve published by Arkhipova et al. (2022) shows, FG Sge is still experiencing the continuous ejection of dusty material. If we adopt the outflow velocity of 200 km s−1 as inferred from the Na I D lines, the material at the outer radius of the dust envelope in 2008 should have traveled 560 au (~1200 R★) by now, which translates into an angular displacement of 190 mas at the distance of 3.0 kpc. While the farthest (thus the coldest) material at the outer boundary may not be detectable in the 10 µm region, temporal variations in the circumstellar envelope can be expected. Furthermore, given the presence of disk-like structures in Sakurai’s object and V605 Aql, it is of great interest to study whether FG Sge has a disk or bipolar structures. To better constrain the structure and evolution of the circumstellar environment of FG Sge, interferometric measurements at more baselines and at multiple epochs are indispensable, which can be carried out with the VLTI/MATISSE instrument (Lopez et al. 2022).
Acknowledgement
We thank the VLTI team at ESO and the MIDI team for carrying out the observations and making the data reduction software publicly available. K.O. acknowledges the support of the Agencia Nacional de Investigación Científica y Desarrollo (ANID) through the FONDECYT Regular grant 1210652. This research made use of the SIMBAD database, operated at the CDS, Strasbourg, France.
References
- Allamandola, L. J., Hudgins, D. M., & Sandford, S. A. 1999, ApJ, 511, L115 [Google Scholar]
- Arenou, F., Grenon, M., & Gómez, A. 1992, A&A, 258, 104 [NASA ADS] [Google Scholar]
- Arkhipova, V. P., Esipov, V. F., Ikonnikova, N. P., Komissarova, G. V., & Shugarov, S. Yu. 2009, Astron. Lett., 35, 534 [NASA ADS] [CrossRef] [Google Scholar]
- Arkhipova, V. P., Ikonnikova, N. P., Shenavrin, V. I., et al. 2022, Astron. Lett., 48, 329 [NASA ADS] [CrossRef] [Google Scholar]
- Blöcker, T. 2001, Ap&SS, 275, 1 [CrossRef] [Google Scholar]
- Bourges, L., Mella, G., Lafrasse, S., et al. 2017, VizieR On-line Data Catalog: II/346 [Google Scholar]
- Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
- Chesneau, O., Meilland, A., Rivinius, T., et al. 2005, A&A, 435, 275 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Chesneau, O., Clayton, G. C., Lykou, F., et al. 2009, A&A, 493, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Clayton, G. C. 1996, PASP, 108, 225 [NASA ADS] [CrossRef] [Google Scholar]
- Clayton, G. C., Kerber, F., Pirzkal, N., et al. 2006, ApJ, 646, L69 [NASA ADS] [CrossRef] [Google Scholar]
- Clayton, G. C., Bond, H. E., Long, L. A., et al. 2013, ApJ, 771, 130 [NASA ADS] [CrossRef] [Google Scholar]
- Chornay, N., & Walton, N. A. 2021, A&A, 656, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Evans, A., Gehrz, R. D., Helton, L. A., & Woodward, C. E. 2015, EAS Pub. Ser., 71, 281 [CrossRef] [EDP Sciences] [Google Scholar]
- Evans, A., Gehrz, R. D., Woodward, C. E., et al. 2020, MNRAS, 493, 1277 [Google Scholar]
- Fadeyev, Yu. A. 2019, Astron. Lett., 45, 655 [NASA ADS] [CrossRef] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gehrz, R. D., Woodward, C. E., Temim, T., Lyke, J. E., & Mason, C. G. 2005, ApJ, 623, 1105 [NASA ADS] [CrossRef] [Google Scholar]
- Gilman, R. C. 1974, ApJS, 28, 397 [NASA ADS] [CrossRef] [Google Scholar]
- Gonzalez, G., Lambert, D. L., Wallerstein, G., et al. 1998, ApJS, 114, 133 [NASA ADS] [CrossRef] [Google Scholar]
- Green, G. M., Schlafly, E., Zucker, G., Speagle, J. S., & Finkbeiner, D. 2019, ApJ, 887, 93 [NASA ADS] [CrossRef] [Google Scholar]
- Herwig, F. 2005, ARA&A, 43, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Herwig, F., Blöcker, T., Langer, N., & Driebe, T. 1999, A&A, 349, L5 [NASA ADS] [Google Scholar]
- Hinkle, K. H., & Joyce, R. R. 2014, ApJ, 785, 146 [NASA ADS] [CrossRef] [Google Scholar]
- Hinkle, K. H., Joyce, R. R., Matheson, T., Lacy, J. H., & Richter, M. J. 2020, ApJ, 904, 34 [NASA ADS] [CrossRef] [Google Scholar]
- Houck, J., Roellig, T., van Cleve, J., et al. 2004, ApJS, 154, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Iben, I., Jr. 1976, ApJ, 208, 165 [NASA ADS] [CrossRef] [Google Scholar]
- Iijima, T., & Strafella, F. 1993, Inf. Bull. Variab. Stars, 3959, 1 [Google Scholar]
- Ivezic, Ž., Nenkova, M., & Elitzur, M. 1999, User Manual for DUSTY, University of Kentucky Internal Report, accessible at http://www.pa.uky.edu/~moshe/dusty [Google Scholar]
- Jaffe, W. 2004, SPIE Proc., 5491, 715 [Google Scholar]
- Jeffery, C. S., & Schönberner, D. 2006, A&A, 459, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jurcsik, J., & Montesinos, B. 1999, New Astron. Rev., 43, 415 [Google Scholar]
- Kipper, T., & Kipper, M. 1993, A&A, 276, 389 [NASA ADS] [Google Scholar]
- Lagage, P. O., Pel, J. W., Authier, M., et al. 2004, The Messenger, 117, 12 [NASA ADS] [Google Scholar]
- Langer, G. E., Kraft, R. P., & Anderson, K. S. 1974, ApJ, 189, 509 [NASA ADS] [CrossRef] [Google Scholar]
- Lawlor, T. M., & MacDonald, J. 2003, ApJ, 583, 913 [NASA ADS] [CrossRef] [Google Scholar]
- Leinert, Ch., Graser, U., & Przygodda, F. 2003, Ap&SS 286, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Leinert, Ch., van Boekel, R., Waters, L. B. F. M., et al. 2004, A&A, 423, 537 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lopez, B., Lagarde, S., Petrov, R. G., et al. 2022, A&A, 659, A192 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miller Bertolami, M. M., Althaus, L. G., Serenelli, A. M., & Panei, J. A. 2006, A&A, 449, 313 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ohnaka, K., Driebe, T., Weigelt, G., & Wittkowski, M. 2007, A&A, 466, 1099 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ohnaka, K., Driebe, T., Hofmann, K.-H., Weigelt, G., & Wittkowski, M. 2008, A&A, 484, 371 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Przygodda, F., Chesneau, O., Graser, U., Leinert, Ch., & Morel, S. 2003, Ap&SS, 286, 85 [NASA ADS] [CrossRef] [Google Scholar]
- Reindl, N., Rauch, T., Miller Bertolami, M. M., Todt, H., & Werner, K. 2017, MNRAS, 464, L51 [Google Scholar]
- Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
- Schönberner, D. 2008, ASP Conf. Ser. 391, 139 [Google Scholar]
- Shenavrin, V. I., Taranova, O. G., & Nadzhip, A. E. 2011, Astron. Rep., 55, 31 [NASA ADS] [CrossRef] [Google Scholar]
- van Hoof, P. A. M., Hajduk, M., Zijlstra, A. A., et al. 2007, A&A, 471, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Zhao-Geisler, R., Quirrenbach, A., Köhler, R., Lopez, B., & Leinert, C. 2011, A&A, 530, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Available at http://www.strw.leidenuniv.nl/~nevec/MIDI
Downloaded from https://sha.ipac.caltech.edu/applications/Spitzer/SHA/
All Tables
MIDI observations of FG Sge and the calibrator: time of observation (coordinated universal time=UTC), projected baseline length Bp, position angle of the projected baseline on the sky (P.A.), seeing in the visible, and the airmass at the times of the observations of FG Sge and the calibrator.
All Figures
![]() |
Fig. 1 VLTI/MIDI observations of FG Sge. a) Observed visibilities of the data set SCI1 (red dots) and the average of SCI2 and SCI3 (blue dots). The error bars are only shown for the data set SCI1 for visual clarity. b) FWHM obtained by a Gaussian fit to the observed visibilities. |
In the text |
![]() |
Fig. 2 Observed mid-infrared spectrum of FG Sge. The solid red line with the error bars shows the MIDI spectrum. Black dots correspond to the Spitzer/IRS spectrum. The green diamond is the VISIR Q3-band (19.5 µm) flux. |
In the text |
![]() |
Fig. 3 Comparison of a spherical shell model with the observed SED and N-band visibilities of FG Sge. The model is characterized with Tin = 1100 K, a piecewise power-law density profile with r−3.5 (rin ≤ r ≤ rtr with rin = 33 R★ and rtr = 5 rin) and r−2 (r > rtr), τV = 8, and Teff = 5000 K. a) Observed SED. Orange filled triangles correspond to visible photometric data from Arkhipova et al. (2009). Orange filled diamonds are the infrared photometric data from Shenavrin et al. (2011). The photometric data were dereddened with AV = 0.63. The thick solid red line is the Spitzer/IRS spectrum. The thin solid blue line corresponds to the MIDI spectrum. The blue cross corresponds to VISIR photometry. The DUSTY model is shown with the solid black line. b) Observed visibilities are shown in the same manner as in Fig. 1a. The solid red and blue lines represent the model visibilities predicted for the 43.0 and 46.0 m baselines, respectively. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.