Bispectrum speckle interferometry observations and radiative transfer models of the red supergiant NML CYG

NML Cyg is a highly evolved OH/IR supergiant, one of the most prominent infrared objects due to its strong obscuration by dust, and supposed to be among the most luminous supergiants in the galaxy. We present the first diffraction-limited 2.13μm observations of NML Cyg with 73 mas resolution. The speckle interferograms were obtained with the 6 m telescope at the Special Astrophysical Observatory, and the image reconstruction is based on the bispectrum speckle-interferometry method. The visibility function declines towards the diffraction limit to ∼0.6. Radiative transfer calculations have been carried out to model the spectral energy distribution, given by ground-based photometry and ISO spectroscopy, and our 2.13 μm visibility function. Additionally, midinfrared visibility functions at 11 μm were considered. The observed dust shell properties do not appear to be in accordance with standard single-shell (uniform outflow) models but seem to require multiple components. Considering previous periods of enhanced mass-loss, various density enhancements in the dust shell were taken into account. An extensive grid of models was calculated for different locations and strenghts of such superwind regions in the dust shell. To match the observations from the optical to the sub-mm domain requires at least two superwind regions embedded in the shell. The best model includes a dust shell with a temperature of 1000 K at its inner radius of 6.2R∗, a close embedded superwind shell extending from 15.5R∗ to 21.7R∗ with an amplitude (factor of density enhancement) of 10, and a far-out density enhancement at 186R∗ with an amplitude of 5. The angular diameters of the central star and of the inner rim of the dust shell amount to 16.2 mas and 105 mas, resp. The diameter of the embedded close superwind region extends from 263 mas to 368 mas, and the inner boundary of the distant superwind region has a diameter of 3. ′′15. In the near-infrared the dust condensation zone is limb-brightened leading to a corresponding ring-like intensity distribution. The grain sizes, a, were found to be in accordance with a standard distribution function, n(a)∼ a−3.5, with a ranging between amin = 0.005 μm and amax = 0.15 μm. The bolometric flux amounts to Fbol = 3.63 10 −9 Wm−2 corresponding to a central-star luminosity of L/L = 1.13 10 5 · (d/kpc). Within the various parts of the dust shell, 1/r density distributions could be maintained differing only in their amplitude A. A slight improvement of the far-infrared properties can be obtained if a shallower density distribution of ρ ∼ 1/r is considered in the distant superwind region. The present-day mass-loss rate was determined to be Ṁ = 1.2 10−4 M /yr. The inner embedded superwind shell corresponds to a phase of enhanced mass-loss (with amplitude 10) in the immediate history of NML Cyg which began 59 yr ago and lasted for ∼18 yr. Correspondingly, the outer superwind region is due to to a high mass-loss period (amplitude 5) which terminated 529 yr ago.