We fit the SEDs of the 18 self-obscured stars using DUSTY (Elitzur & Ivezic2001; Ivezic & Elitzur1997; Ivezic et al.1999) to model radiation transfer through a spherical dusty medium surrounding a star and Figure7 shows the best fit models. We estimate the properties of a black-body source obscured by a surrounding dusty shell that would produce the best fit to the observed SED (see Figure6 for an example). We considered models with either graphitic or silicate (Draine & Lee1984) dust. We distributed the dust in a shell with a density distribution. The models are defined by the stellar luminosity (), stellar temperature (), the total (absorption plus scattering) -band optical depth (), the dust temperature at the inner edge of the dust distribution (), and the shell thickness . The exact value of has little effect on the results, and after a series of experiments with , we fixed for the final results. We embedded DUSTY inside a Markov Chain Monte Carlo (MCMC) driver to fit each SED by varying , , and . We limit to a maximum value of 30,000K to exclude unrealistic temperature regimes.
The parameters of the best fit model determine the radius of the inner edge of
the dust distribution ().
The mass of the shell is
For a comparison sample, we followed the same procedures for the SEDs of three well-studied dust obscured stars: Car (Humphreys & Davidson1994); the Galactic OH/IR star IRC+10420 (Jones et al.1993; Humphreys et al.1997; Tiffany et al.2010); and M33's VariableA, which had a brief period of high mass loss leading to dust obscuration over the last years (Humphreys et al.2006; Hubble & Sandage1953; Humphreys et al.1987). We use the same SEDs for these stars as in Khan et al. (2013). In Table3, we report , , , , , , (Equation1), and (Equation2) for the best fit models for these three sources as well as the newly identified stars. The stellar luminosities required for both dust types are mutually consistent because the optically thick dust shell acts as a calorimeter. However, because the stars are heavily obscured and we have limited optical/near-IR SEDs, the stellar temperatures generally are not well constrained. In some cases, for different dust types, equally good models can be obtained for either a hot (, such as a LBV in quiescence) or a relatively cooler (, such as a LBV in outburst) star. Indeed, for many of our 18 sources, the best fit is near the fixed upper limit of . To address this issue, we also tabulated the models on a grid of three fixed stellar temperatures, , for each dust type. The resulting best fit parameters are reported in Tables 4 and 5.
Figure9 shows the integrated luminosities of the newly identified self-obscured stars described in Section3.2 as a function of for the best fit graphitic models of each source. ObjectX, IRC, M33VarA, and Car are shown for comparison. Figure10 shows the same quantities, but for various dust models and temperature assumptions. It is apparent from Figure10 and Tables 3, 4 and 5 that the integrated luminosity and ejecta mass estimates are robust to these uncertainties. The exceptions are N2403-4 and N7793-3. Without any optical or near-IR data, many of the models of N7793-3 are unstable so we simply drop it. The only models having a luminosity in significant excess of are some of the fixed temperature models of N2403-4. These models have a poor goodness of fit and can be ignored.
One check on our selection methods is to examine the distribution of shell radii.
Crudely, we can see a shell until
it either becomes optically thin or too cold, so the probability
distribution of a shell's radius assuming a constant expansion velocity is