Free Access
Issue
A&A
Volume 655, November 2021
Article Number A98
Number of page(s) 34
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202141279
Published online 29 November 2021

© ESO 2021

1. Introduction

The origin of dust on cosmic scales is of great astrophysical interest. The cosmic importance of dust is illustrated by the fact that about 50% of starlight is reprocessed by dust (Dole et al. 2006). Thus, dust dominates how we see galaxies. Dust is very important for the thermal balance of the interstellar medium (ISM) of galaxies and accordingly for star formation (McKinnon et al. 2016, and references therein). Large amounts of dust have been observed in high redshift galaxies (z ⪆ 6, e.g. Riechers et al. 2017; Laporte et al. 2017), where low- and intermediate-mass stars did not have time to evolve and thus to inject dust into the ISM. It is widely accepted now that dust must arise from rapidly evolving (1−10 Myr) massive stars, particularly in the ejecta of core-collapse supernovae (CC-SNe) (Sugerman et al. 2006). In fact, large amounts of dust (0.1−1 M) were postulated by Todini & Ferrara (2001) and Nozawa et al. (2003), under the assumption of favourable conditions in the SN ejecta. Nevertheless, for nearby CC-SNe, the dust masses inferred from mid-infrared observations are in general of the order of 10−4 − 10−2M (e.g. Kotak et al. 2009; Fox et al. 2011). Dust masses of the order of 10−1M were reported from sub-millimetre observations of SN1987A (e.g. Indebetouw et al. 2014; Matsuura et al. 2015), suggesting that perhaps cold dust emission is not captured by mid-infrared observations. However, such estimates are difficult to compare because they are often based on different observational constraints and assumptions. It is also expected that only a fraction of the newly formed dust reaches the ISM, as it will be heavily processed and destroyed by the reverse shock (see Micelotta et al. 2018; Gall & Hjorth 2018). Dust can also grow in dense ISM from existing dust grains (e.g. Draine 2009), although this is still not sufficient, together with the contribution from SNe, to explain dust yields in galaxies (e.g. Tamura et al. 2019). The topic is also controversial with respect to the nearby Magellanic Clouds (MCs).

The MCs are two irregular dwarf galaxies with metallicities 0.5 and 0.25 times Solar for the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), respectively. Their relative proximity to the Milky Way (MW), accurate distance measurements and low line-of-sight extinction make them an ideal laboratory, resembling the physical conditions of galaxies at the peak of star formation. In the past decades, several efforts were made to estimate dust yields in these galaxies and determine the primary dust sources at low metallicities. With the advent of Herschel this became possible. Particularly the HERITAGE Herschel key project (Meixner et al. 2013), observing in five bands from 100 to 500 μm, offered a plethora of data that advanced our understanding of the dust budget problem. Recent estimates report integrated dust masses of (7.3 ± 1.7) × 105 and (8.3 ± 2.1) × 104M for the LMC and SMC, respectively (Gordon et al. 2014), although later Chastenet et al. (2017) derived dust masses about a factor of two smaller. Several teams tried to identify the main sources of dust in such galaxies. Srinivasan et al. (2009, 2016), Matsuura et al. (2009), Riebel et al. (2012), and Boyer et al. (2012) studied the dust production rate from low- and intermediate-mass stars in the MCs using mid-infrared observations. Jones et al. (2015) identified some tens of new evolved stellar objects from HERITAGE. These works generally agree on the result that the global dust production from the numerous population of asymptotic giant branch (AGB) and red supergiant (RSG) stars cannot account for the ISM dust reported from far-infrared and sub-millimetre observations. Another result is that the majority of dust input from low- and intermediate-mass stellar sources comes from “extreme” AGB (Srinivasan et al. 2009; Riebel et al. 2012; Boyer et al. 2012).

Spitzer observations from the SAGE legacy survey (Meixner et al. 2006) were used by Bonanos et al. (2009, 2010) to investigate the infrared properties of massive stars. They found an excess of emission at 24 μm due to dust from a handful of luminous blue variable (LBV) stars. Some of these, and other LBVs with known extended nebulae detected with HST (Weis 2003), were then followed up and large masses of dust were measured (Guha Niyogi et al. 2014; Agliozzo et al. 2017a,b, 2019).

In the classical scenario, LBVs are transitional objects from the main sequence of O, B stars to the Wolf-Rayet (WR) phase (Conti 1984). This represents a short phase of evolution of the most massive (M ≥ 25 M) stars, ranging in duration from 104 to 105 yr (e.g. Maeder 1989; Massey et al. 2007; Smith 2014). Observationally, LBVs are blue supergiants (BSGs) that experience instabilities, observed in the form of spectroscopic and photometric variabilities, usually in the optical and near-infrared. These instabilities are poorly understood but, on the basis of their observational properties, are phenomenologically divided into three main groups: giant eruptions, S Doradus-type variability, and micro-variabilities (Humphreys & Davidson 1994; van Genderen 2001). Although the physical mechanisms underlying these instabilities are possibly different (e.g. the vicinity to the Eddington limit, sub-photospheric instabilities, bistability jump, envelope inflation, wind-envelope interaction, fast rotation, binarity, stellar merger; Humphreys & Davidson 1994; Gallagher 1989; Vink et al. 1999; Groh et al. 2009; Gräfener et al. 2012; Portegies Zwart & van den Heuvel 2016; Owocki et al. 2017; Grassitelli et al. 2021), observations suggest increases in the mass-loss rates by about a factor of three in S Doradus-type variability (e.g. Vink & de Koter 2002) or by several orders of magnitude in giant eruptions, like that of η Carinae (e.g. Humphreys & Davidson 1994). During the quiescent state, LBV mass-loss rates range between 10−7 and 10−5M yr−1, typical of BSGs as explained by the radiatively-driven stellar-wind theory (Castor et al. 1975; Sundqvist et al. 2019).

Particularly enigmatic are the rare giant eruptions, initially often mistaken for CC-SNe, because of their brightness, and subsequently named SN impostors (e.g. Van Dyk & Matheson 2012; Mauerhan et al. 2013; Tartaglia et al. 2016; Pastorello et al. 2017; Elias-Rosa et al. 2018; Reguitti et al. 2019). Famous examples in our Galaxy are η Car and P Cygni. Observationally, mass-loss rates exceeding 10−3 − 10−2M yr−1 are estimated (e.g. Smith 2014). In the case of η Car, more than 40 M (Morris et al. 2017) were released during the Great Eruption in the XIX century, which lasted about 20 years, forming the Homunculus nebula. In eruptions of SN impostors, nebular masses of the order of 0.1 M are observed (e.g. Margutti et al. 2014). Several observational and theoretical works suggest that some LBVs could be the immediate progenitors of SNe (e.g. Kotak & Vink 2006; Trundle et al. 2008; Gal-Yam & Leonard 2009; Smith et al. 2007, 2008b; Groh et al. 2013; Groh 2014; Boian & Groh 2018).

Luminous blue variable stars are often surrounded by dusty circumstellar nebulae. In the Milky Way typical dust masses ranging between 10−3 and 10−1M are reported (Hutsemekers 1997; Umana et al. 2012; Vamvatira-Nakou et al. 2013, 2015; Agliozzo et al. 2014; Lau et al. 2014; Arneson et al. 2018). Furthermore, several tens of candidate LBVs were identified on the basis of the detection of infrared circumstellar nebulae (Egan et al. 2002; Clark et al. 2003; Gvaramadze et al. 2010) and of an infrared excess usually peaking between 24 and 70 μm. The central star, which is a BSG during quiescence, is very hot, and one does not expect dust condensation in its wind. However, the physical conditions in the optically thick and cool “pseudo-photospheres” formed during eruptions are favourable for dust condensation and growth (Kochanek 2011a). Kochanek (2011a) showed analytically that the required temperature and particle density for dust formation are met for mass-loss rates larger than 10−2.5M yr−1, that are typical of giant eruptions. Grain growth and dust evolution may be observed in S Doradus variables (which have mass-loss rates that are a couple of orders of magnitudes lower), as suggested by the disappearance of the silicate bump around 10−13 μm during the most recent S Dor outburst of the Magellanic LBV RMC 71 (Mehner et al. 2017). Another way to form dust around a massive star is in the colliding winds of close binaries, as in some Wolf Rayet stars. In these systems dust seems to form episodically, near the periastron passage, or persistently, creating pinwheel-like nebulae (e.g. Crowther 2007; Lau et al. 2020). This mechanism of dust formation is presumably the case of the Galactic LBV binary HR Car (Boffin et al. 2016), a model of which to explain the clumpy and dusty inner nebula is an Archimedean spiral centred on the binary system (alternative to the expanding bipolar lobe model, see Buemi et al. 2017, and references therein). Episodic dust formation also potentially occurs in the shocked regions of η Car’s binary colliding winds (Smith 2010).

In the MCs only up to one third of the LBV sample appears in studies of the primary dust sources (see above), and the dust masses from these objects remain unconstrained. Here we aim to expand these studies to the full population of LBVs in the MCs. Figs. 1 and 2 show the locations of the candidate and confirmed LBVs in the LMC and SMC that are analysed in this work. We have built a maximally complete multi-wavelength dataset of 31 LBVs mostly from archival space telescope images and catalogues, but also have included ground-based mid-infrared observations of three sources. In addition, we derived new photometric values for sources missing in catalogues or with low-accuracy photometry (Sects. 2 and 3). The goals of this work are:

  1. Review the status of the MC LBV sample by means of their infrared properties, combined with literature information (Sect. 4).

  2. Estimate the dust masses of LBVs in the MCs by modelling individual sources (Sect. 5.1) and modelling photometry from stacked images of the LBVs in the LMC (Sect. 5.2).

  3. Evaluate the importance of LBVs in producing dust by comparing the derived dust mass with that inferred from the stacking of lower-mass stars and with literature values of ISM dust yields (Sects. 5.3 and 6).

thumbnail Fig. 1.

Cutout of the SAGE MIPS24 image (Meixner et al. 2006) of the LMC and positions of the candidate and confirmed LBVs analysed in this work. The numbers are the IDs in Table 1, and the red circles have a radius of 30″. North is towards the top right.

thumbnail Fig. 2.

Cutout of the SAGE MIPS24 image (Meixner et al. 2006) of the SMC and positions of the candidate and confirmed LBVs analysed in this work. The numbers are the IDs in Table 1, and the red circles have a radius of 30″. North is up and east is left.

2. The multiwavelength dataset

The sample was extracted from the most recent census of LBVs in the Local Group (Richardson & Mehner 2018), which contains 31 objects in the MCs. Table 1 includes the names and celestial coordinates of all the stars included in this work. More information on individual sources is provided in Sect. 4.

Table 1.

Sources, celestial coordinates (Gaia Collaboration 2018), LBV status (previous: Richardson & Mehner 2018; new: this work), and SED-based classification proposed in this work.

2.1. VISIR observations

RMC 71, HD 38489, and HD 34664 were observed in 2016 and 2017 (ESO program IDs: 097.D-0612(A,B) and 0100.D-0469(A,B)) with the upgraded Vlt Imager and Spectrometer for the mid-InfraRed (VISIR; Lagage et al. 2004; Käufl et al. 2015; Kerber et al. 2016) on the Very Large Telescope (VLT, Enard 1991). Some of the data of RMC 71 from program 097.D-0612 were previously published in Mehner et al. (2017). The images were taken in standard imaging mode through filters J8.9 (8.7 ± 0.37 μm), PAH2_2 (11.88 ± 0.19 μm) and Q1 (17.65 ± 0.42 μm) with perpendicular nodding and a chop/nod throw of 20″ in 097.D-0612 and 10″ in 0100.D-0469, respectively. The individual observations are listed in Table 2. J8.9 was used for target acquisition only (hence the short exposure times). For flux calibration and point spread function (PSF) reference, the science observations were preceded or followed by the observation of mid-infrared standard stars (Cohen et al. 1999), which provide a systematic flux uncertainty of ≤10%. The data were reduced with the custom Python tool VISIR and ISAAC Pipeline Environment1 (VIPE). Fluxes and source sizes were measured with Gaussian-fitting photometry performed as part of the VIPE reduction.

Table 2.

Journal of VISIR observations.

Both epochs of the PAH2_2 observations of HD 34664 were affected by a chopping-induced PSF instability at the VLT, leading to an artificial elongation of the source. However, in the J8.9 and Q1 filters, the object is consistent with being unresolved (FWHM ≤ 0.4″ in J8.9 and ≤0.5″ in Q1).

The PAH2_2 observation of HD 38489 suffers from the same PSF-instability problem while, during the Q1 observation, the source was offset such that it was on the detector in only half of the chopping positions. The latter fact has been taken into account during the data reduction. Like HD 34664, HD 38489 is not obviously resolved, with upper limits in FWHM of ≤0.37″ in J8.9 and ≤0.5″ in Q1, respectively.

RMC 71 appears resolved in all images, with the following FWHM dimensions: ∼0.46″ × 0.42″ in J8.9 (PA ∼ 126 deg), ∼0.96″ × 0.83″ in PAH2_2 (PA ∼ 110 deg), and ∼0.87″ round in Q1. These are to be regarded as upper limits because none of the observations was made under diffraction-limited conditions (also because of the low declination of the object). Thus the differences between the two axes in each filter are considered insignificant. However, the extended nature of this source is robust, because of the consistent appearance at both epochs, including the position angle, and the significantly larger FWHM than found for the calibrator stars. The circumstellar environment of RMC 71 appears in agreement with a broad circular shell and an asymmetric core, which is offset from the centre by ∼0.2″ toward the north-west (Fig. 3).

thumbnail Fig. 3.

VISIR images of RMC 71 (top panels), and the standard calibrators observed on the same nights (bottom panels, with identifications given in the lower right corners). North is up and east is left. The colour-scale is in arbitrary units.

2.2. Archival data

We consulted the infrared catalogues with the CDS VizieR (Ochsenbein et al. 2000) and IRSA Gator2 tools and retrieved photometry from the 2MASS Point Source Catalogue (Cutri et al. 2003), the Spitzer SAGE legacy survey (Data Release 3/final Meixner et al. 2006), the “Optically bright post-AGB population of LMC” (van Aarle et al. 2011) and “Spitzer Atlas of Stellar Spectra” (Ardila et al. 2010) catalogues, the AKARI IRC all-sky survey and LMC point source catalogues (Ishihara et al. 2010; Kato et al. 2012), WISE surveys (Cutri et al. 2012), the MSX Point Source Catalogue (Egan et al. 2003), and Herschel HERITAGE survey (Meixner et al. 2013). We evaluated all stars for possible contamination by neighbouring sources by visually inspecting the infrared images when available (for each source, these are shown in Figs. B.1B.6).

We checked for ALMA observations of the sample. ALMA Band 7 (343 GHz) observations of RMC 127, RMC 143 and LHA 120-S 61 were previously reported in Agliozzo et al. (2017a,b, 2019). Unpublished ALMA/ACA data of a mosaic region containing RMC 14 are available in the archive, in Bands 3 (∼93 GHz; project ID: 2018.A.00049.S) and 6 (∼220 GHz; project ID: 2017.A.00054.S). RMC 14 is not detected, with 3σ upper limits of 1.41 mJy at 93 GHz and 6.0 mJy at 220 GHz. We also made use of the SPT+Planck maps of the MCs at 1.4 mm (Crawford et al. 2016) for the stacking analysis in Sect. 5.2. None of the sources is detected in these maps obtained with the large beam of 1.5′.

Finally, we used the available Infrared Spectrograph (IRS; Houck et al. 2004) spectra from the Combined Atlas of Sources with Spitzer IRS Spectra (CASSIS; Lebouteiller et al. 2011, 2015) or from the IRSA data archive (as SSC Enhanced Products). RMC 127, LHA 120-S 61, and RMC 84 were observed under Program Name/Id: STELLARATLAS/485 and published by Ardila et al. (2010). The RMC 66, HD 38489, HD 34664, and LHA 115-S 18 spectra were published in Kastner et al. (2010), although we found that the spectrum of LHA 115-S 18 is not on source. The spectra of RMC 14 and RMC 40 were discussed in Ruffle et al. (2015), and those of RMC 110 and S Dor in Jones et al. (2017). Data of RMC 99 and Sk –69 279 are available in the archive: we are not aware of any previous publications. The IRS spectra are shown with grey continuous lines in the flux density distribution plots in Appendix B. RMC 71 and RMC 66 also have Infrared Space Observatory (ISO) data3 and low resolution Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) SED-mode spectra (presented in Kemper et al. 2010; van Loon et al. 2010, which also include RMC 126) but they are not considered in this work. The IRS spectra are mainly used to guide the interpretation of the infrared photometry when no measurement at 24 μm is possible due to confusion.

3. Data analysis

3.1. Photometry from Spitzer MIPS24/MIPS70

A small number of sources are not included in the Spitzer SAGE catalogues, or they have poor catalogue photometry, most likely because of the strict requirements adopted to identify point-sources (Meixner et al. 2006). Indeed we find that these sources are usually extended or embedded in a confused region, requiring a manual extraction of the flux measurement. The most remarkable case is Sk –69 279, which appears as a well detected extended bubble in the MIPS24 map, but has no catalogue photometry. We performed aperture photometry for all sources that were either lacking measurements or are extended in the images. As a cross-check, we also measured all targets with catalogue values. Following this verification, we chose to use the existing catalogue photometry for the non-extended sources, unless indicated otherwise.

Our aperture photometry was performed with the astropy package “photoutils”. The aperture radii, sky annuli, and aperture corrections, as recommended in the MIPS instrument handbook4, are listed in Tables 3 and 4, together with the measured flux densities, comparison with the catalogue value when available, and relevant comments on the sources and their environments. The reported flux density errors only account for the uncertainty in determining the background (defined as the standard deviation in the background annulus integrated over the aperture), which generally dominates over the absolute flux calibration error and statistical noise in the source aperture. For non-detected sources, we provide 3σ upper limits.

Table 3.

MIPS-24 aperture photometry.

Table 4.

MIPS-70 aperture photometry and 3σ upper limits.

We checked the quality of our photometry against the catalogues (see Tables 3 and 4): the results are consistent with the published values. In bold we highlight our measurements that are used in the later analysis instead of the catalogue values.

3.2. Photometry and upper limits from Herschel

Only a few sources were well detected with the Photodetector Array Camera and Spectrometer (PACS, Poglitsch et al. 2010) at 100 and 160 μm, and their flux densities were taken from the HERITAGE catalogue (Meixner et al. 2013). Additionally, we report in Table 5 new PACS100 photometry for LHA 120-S 119 and HD 38489 (although the latter is only a marginal detection). To maximise the signal-to-noise ratio, we adopted a procedure equivalent to Gaussian fitting. Local peaks were searched for in boxes with size equivalent to the beam. The regions used to evaluate the sky contribution were separated by one or two beams. The 1σ total uncertainty is determined by summing in quadrature the standard deviation of the background and the flux calibration uncertainty of 10%. Because of contamination of S Dor from a strong nearby nebular region, we reduced the background spacing to one beam. Nevertheless, the source was not unambiguously detected. For the sources with no detection (either in this work or in the catalogue), we report the 3σ upper limit at 100 μm determined as described above.

Table 5.

PACS-100 point source extraction and 3σ upper limits.

Finally, in PACS160 we only find a single new detection, LHA 120-S 61, which has a flux density of 60 ± 16 mJy. We do not report any new detection obtained at longer wavelengths with the Spectral and Photometric Imaging Receiver (SPIRE) instrument (Griffin et al. 2010). The upper limits at 160 μm and longer wavelengths are not provided as they do not add useful information owing to the larger beam sizes and high interstellar dust brightness at such wavelengths.

4. The LBV SED classification scheme and colour-magnitude diagram

We first need to ensure to limit our sample to LBV stars and find a classification criterion that does not rely on the spectroscopic and photometric variations only, but also on properties of the stellar ejecta. Table 1 summarises the relevant information on our sample, including the stellar classification from the literature. Stellar parameters, extinction, and nebular size when known are summarised in Table A.1. Unfortunately, we cannot provide information on the binary status, as to our knowledge this is known only for RMC 14 in the SMC and RMC 81 in the LMC. The complete set of photometry from this work and from the catalogues are in Tables C.1C.3. For all sources, the flux density distribution starting from 1 μm are shown in Figs. B.1B.6, where in addition to the photometric points the stellar continuum (on the basis of stellar parameters as in Table A.1) and the IRS spectra are also included when available.

In their sample of 1750 massive stars in the LMC, Bonanos et al. (2009) included six LBVs (S Dor, RMC 127, RMC 71, RMC 110, RMC 85, HD 269582) and concluded that these stars do not represent a homogeneous class of objects on the basis of their mid-infrared colours. We identify four classes based on the individual infrared spectral energy distributions (SED):

  1. Sources with a dusty nebula peaking in the mid- to far-infrared (Class 1). These can be further divided into two groups: with or without a significant free-free excess above the stellar photosphere (Classes 1a and 1b, respectively). An example can be found in the top-right panel in Fig. 4, where the object shows a moderate free-free continuum up to 10 μm, and at longer wavelengths a dominant infrared excess due to dust in the extended nebula, previously imaged by Weis et al. (2003) and Agliozzo et al. (2012). By contrast, the top left panel shows an example of strong free-free emission.

  2. Sources with only a free-free emission excess (Class 2). An example is given in the middle-left panel in Fig. 4. For these sources the MIPS24 flux density is consistent with a power law describing the free-free continuum.

  3. Sources dominated at all wavelengths ≳2 μm by a hot dusty component, suggestive of a disc close to the star (Class 3). In the literature, these are all classified as B[e] supergiants (see Table A.1). An example of a Class 3 SED is HD 34664, shown in the second-right panel of Fig. 4.

  4. Featureless sources with a single component consistent with the stellar photosphere at all wavelengths (Class 4). An example is presented in the bottom panel of Fig. 4.

thumbnail Fig. 4.

Examples of each of the four LBV SED classes defined in Sect. 4.

Table 1 summarises the classification of the sources in the scheme proposed here. Among them, only one object cannot be classified, RMC 85. This is an active LBV (van Genderen et al. 1998). The star is embedded in a very crowded region. Because of the absence of reliable photometry longward of ∼20 μm (confirming Bonanos et al. 2009), we are not able to distinguish whether it is SED Class 1 or Class 2 (see Fig. B.4).

The dust composition of sources in Class 1 typically comprises amorphous silicates, although a mixed chemistry is sometimes reported, like in the case of RMC 71 in the LMC (Guha Niyogi et al. 2014) or HD 168625 in the Milky Way (Umana et al. 2010). Crystalline silicates are also found in LBVs, albeit rarely (e.g. RMC 71 and HR Car, Voors et al. 1999; Umana et al. 2009, respectively), while they seem to be a common feature of sgB[e] stars (e.g. Kastner et al. 2006, 2010; de Wit et al. 2014).

In Fig. 5 we present the [8.0] versus [8.0]–[24] colour-magnitude diagram (CMD) of stars in the LMC, depicting the position of the LBV sample analysed in this work (filled circles of various colours representing the different LBV SED classes as described in the legend); early and late supergiants as in Bonanos et al. (2009) (downward pointing triangle symbols); grey diamonds are low- and intermediate-mass, but also some high mass stars, planetary nebulae and SN remnants from Jones et al. (2015); AGBs and RSGs from Riebel et al. (2012) and Young Stellar Objects (YSOs) from Whitney et al. (2008). Stars with SED Classes 2, 3 and 4 occupy the same region as RSGs and other types of massive stars in the diagram, which Bonanos et al. found indistinguishable from AGBs and young stellar objects, highlighting the importance of using wide-baseline photometry. However, we notice some differences: Class 2 objects tend to be fainter than RSGs at 8 μm, due to lacking a hot dusty envelope; Class 3 are usually brighter than RSGs, due to their hot dusty disc; the only Class 4 object with MIPS24 photometry is the faintest in this region of the diagram, consistent with no excess above the photosphere. For stars in the SED Class 1, the [8.0]−[24] colours coincide with the sample of stars in Jones et al. However, the nature of 7 out of their 35 objects was not spectroscopically confirmed. These could all be sources with a dusty gaseous nebula (e.g. SN remnants, planetary nebulae, LBV nebulae). For stars in the SMC we refer to the colour-diagram by Bonanos et al. (2010), who included all the LBVs in the census (Richardson & Mehner 2018). Bonanos et al. found that LBVs are among the most luminous stars at 24 μm, but they are not clearly separated from other types of star.

thumbnail Fig. 5.

Mid-infrared colour-magnitude diagram for stars in the LMC. The LBVs are shown with filled circles in different colours, depending on the SED class defined in this work. Several types of supergiant stars (RSGs, sgB[e], supergiants A and B, six LBVs and other massive stars) as in Bonanos et al. (2009) are displayed with triangles of different colours as identified in the legend. Dusty evolved stars from Jones et al. (2015) are shown with diamond symbols (these include post-AGBs, planetary nebulae, sgB[e] stars, one LBV, one WR and two SN remnants). Pixel symbols indicate the AGBs and RSGs from Riebel et al. (2012) and point symbols are confirmed and candidate YSOs from Whitney et al. (2008).

4.1. Class 1: sources with dusty nebulae peaking in the mid- to far-infrared

This class contains well-known extended MCs LBV nebulae (LBVNe), extensively studied in the literature. Some of them are associated with very active LBVs (Walborn et al. 2017). Some stars are classified as candidate LBVs: these are stars that do not satisfy all LBV classification criteria (both observed spectroscopic and photometric variability; Humphreys & Davidson 1994), however the presence of a circumstellar nebula around some of them suggests that they are ex/dormant LBVs.

Depending on whether or not they have an excess of near-infrared emission due to the ionised stellar wind, stars in this class can be further divided into: Class 1a, which are objects that present free-free emission from the wind; Class 1b, which are objects with no detected free-free emission from the wind. For two objects, Sk –69 279 and RMC 78, we report the first identification of a dusty circumstellar nebula. The dust component in these sources peaks between 24 and 100 μm, similarly to several Galactic objects with extended nebulae (typically with sizes of the order of 0.1−1 pc, e.g. G79.29+0.46, Pistol Star, AG Car, AFGL2298, HR Car; Agliozzo et al. 2014; Lau et al. 2014; Vamvatira-Nakou et al. 2015; Buemi et al. 2010, 2017, and references therein). Among the sources in this class, five (RMC 127, RMC 71, S Dor, LHA 120-S 61, and RMC 40) have an IRS spectrum and all of them show, together with the thermal continuum, the silicate bump around 10−12 μm.

The objects classified as belonging to Class 1 are briefly introduced (see Figs. B.1 and B.2):

S Dor is the prototype of the most common LBV variability (a summary of the pioneering works can be found in Wolf 1989). The infrared excess from a dusty nebula was reported in Bonanos et al. (2009) from mid-infrared photometry and in Jones et al. (2017) from the IRS spectrum. The power-law behaviour in the IRS data below 10 μm is more likely consistent with strong free-free emission in the ionised stellar wind rather than the stellar photosphere. The bright surrounding environment does not allow for the disentanglement of the nebula at 70 μm and beyond. The SED Class of this object is 1a.

RMC 127: this star belongs to the group of active LBVs (Walborn et al. 2008, 2017). This is another object with SED Class 1a. See Agliozzo et al. (2017b) for a multiwavelength study of its nebula.

LHA 120-S 61: candidate LBV (Humphreys & Davidson 1994) with SED Class 1a, its dusty nebula was modelled by Agliozzo et al. (2017a).

Sk –69 279: candidate LBV with SED Class 1a, the nebula was previously imaged by HST (Weis & Duschl 2002). We find that at 24 μm it appears as a bubble, resembling the morphology revealed in the optical. To our knowledge, this is the first report of dust in this object.

LHA 120-S 119: candidate LBV (van Genderen 2001) with a known nebula (SED Class 1a) that was previously imaged at different wavelengths (Agliozzo et al. 2012, and references therein). In the infrared it is possible to notice a thermal component of dust, detected up to 100 μm. The nebula is marginally resolved at 24 μm, revealing a marked asymmetry similar to the optical and radio images. This asymmetry was explained as a bow-shock, given the run-away nature of its associated star (Danforth & Chu 2001).

RMC 81: infrared emission from cool dust around this candidate LBV (van Genderen 2001) was first reportedly detected by IRAS at 60 μm (Beichman et al. 1988; Tubbesing et al. 2002). However, inspection of the MIPS70 image, comparison with the MIPS photometry, the lack of 25 or 12 μm IRAS detections, and considering the 1σ position offset, it appears that this flux was greatly elevated by ISM confusion. The object is detected in both MIPS bands, even at 70 μm where the surrounding environment is bright. The star is an eclipsing binary (Stahl et al. 1987). The SED Class of this object is 1a.

LHA 120-S 18: at 24 μm and beyond, the environment around this candidate LBV (van Genderen 2001) is bright and the object is not distinguishable from the background. The SED Class of this object is 1a.

RMC 143: the dusty and massive nebula (SED Class 1b) was extensively studied by Agliozzo et al. (2019), who also confirmed the LBV status of this star.

RMC 71: this star also belongs to the group of very active LBVs (Mehner et al. 2013; Walborn et al. 2017), with SED Class 1b. It is possible to observe variability in the infrared SED (compare the mid-infrared photometry with the IRS spectrum), as the data were acquired at different epochs during the most recent stellar outburst (Mehner et al. 2017). Here we show a similar dataset, but extended to longer wavelengths. The entire SED was previously modelled by Guha Niyogi et al. (2014).

RMC 78: at 24 μm the object (a candidate LBV, van Genderen & Sterken 1999) appears as a double source (SED Class 1b).

RMC 40: this star belongs to the group of confirmed LBVs (Szeifert et al. 1993) in the SMC. The first detection at 24 μm was reported by Bonanos et al. (2010). Here we add a new clear detection at 70 μm. The SED Class of this object is 1b. An IRS spectrum is available and suggests that the dust emission peaks between the two MIPS bands.

HD 269216: a confirmed LBV (Prinja & Schild 1991), this object is well detected and marginally resolved in MIPS 24. We add a new measurement as the catalogue PSF fitting may have missed a large fraction of its flux density. The SED Class of this object is 1b.

RMC 110: this is a confirmed LBV manifesting the eruptive phenomenon a few times in the past decades (Stahl et al. 1990; Campagnolo et al. 2018). Because of confusion at wavelengths beyond ∼20 μm it is not possible to obtain reliable photometry of dust in its ejecta (SED Class 1b).

4.2. Class 2: sources with free-free excess only

This class consists of the following objects (see Fig. B.3):

RMC 116 is a strongly active S Dor variable (van Genderen 2001). Its SED shows an ionised stellar wind in the form of a free-free excess above the stellar photosphere detected up to 24 μm.

RMC 123 and RMC 99 belong to a sub-group proposed by Stahl et al. (1984) and Stahl & Wolf (1987), characterised by an excess of free-free emission from a gaseous disc, rich in forbidden lines, but with no dust features. The free-free emission of RMC 123 is detected up to 24 μm. At longer wavelengths, the SED of RMC 99 is dominated by contamination from neighbouring sources so that the presence of an infrared excess is questionable. Both stars are catalogued as weakly active S Dor variables in van Genderen (2001) and are considered candidate LBVs. RMC 99 was observed with IRS, but the spectrum is likely dominated by interstellar gas and dust.

HD 269582 is a very active LBV (Walborn et al. 2017) with a strong ionised stellar wind detected up to 24 μm. Differences between 2MASS and IRAC data suggest that the star was active (highly variable) in the years of the infrared observations.

RMC 74 is a candidate LBV (van Genderen 2001) with a SED similar to that of RMC 116.

RMC 14, a well-known multiple-system containing an active LBV in the SMC (e.g. Barba et al. 1996), shows an SED similar to HD 269582, with a strong free-free excess above the continuum (a similar SED can also be found in Bonanos et al. 2010). This is also seen in two other active LBVs (RMC 127, Agliozzo et al. 2017a, and S Dor). The system was observed with IRS, but, because of the large slit width, the spectrum is likely dominated by interstellar gas and dust.

4.3. Class 3: sources dominated by a hot dusty disc

RMC 66, RMC 84, HD 38489, HD 34664, RMC 126, RMC 4, and LHA 115-S 18: These stars are all known in the literature as supergiant B[e] stars (e.g. Zickgraf et al. 1996a; van Genderen & Sterken 1999; Aret et al. 2012; Bonanos et al. 2009, 2010; Kastner et al. 2006, 2010), but sometimes they are also listed in the sample of candidate LBVs because of other observational commonalities (e.g. variability of emission lines and P Cyg profiles, van Genderen 2001). Their SEDs are suggestive of a different dust formation mechanism. Dust in these stars forms in the equatorial plane (perhaps as a consequence of fast rotation or binarity), and resides in a circumstellar disc-like configuration, resembling the mid-infrared spectra of T Tauri and disc-enshrouded Herbig Ae/Be stars (Kastner et al. 2006, 2010). Furthermore, observations suggest the silicates to be crystalline, which requires long-lived dust particles permanently close to the star, as opposed to outflowing material like in the case of LBV eruptions (for a review, see de Wit et al. 2014). The IRS spectrum, and IRAC and MIPS images of HD 38489, RMC 66, LHA 115-S 18, RMC 126 and HD 34664 first appeared in Kastner et al. (2006, 2010).

4.4. Class 4: sources with no evident free-free or dust excess

This class comprises the following objects:

HD 269604 was classified as A1 Ia-0 by Osmer (1973). It appears in several works as a candidate LBV, but van Genderen (2001) flagged it as non-candidate. We do not find evidence of free-free and/or dust infrared excess.

RMC 128 (B2 Ia with N weak, Fitzpatrick 1991) was identified as non-candidate by van Genderen, but it is mentioned as a candidate LBV in several other works. RMC 128 is only 20″ away from RMC 127, which is much more luminous and the star is hard to discern in the mid-infrared. In the near-infrared, there is no evidence of free-free and dust excess. We do not find any observational argument to keep the star in the list of candidate LBVs, although its spatial proximity to the very active RMC 127 certainly makes for an interesting target (e.g. if the two stars are co-eval).

RMC 149 was characterised by van Genderen (2001) as a weakly active S Dor variable (with low-amplitude variations over decades). However, Fariña et al. (2009) re-classified it as O8.5 II((f)). We do not find typical signatures of LBV nebular features (ionised gas and/or dust excess). The infrared images reveal that the source is embedded in a crowded environment.

Sk –69 271 is a yellow supergiant (YSG) (Neugent et al. 2012). It appeared in the list of candidate LBVs because of nearby nebular material resembling a half shell, most likely of interstellar origin (Weis et al. 1997). This source is particularly faint already at 8 μm and does not show any features at longer wavelengths. For lack of other evidence to keep these stars in the LBV sample, we discard them from further analysis and also suggest their removal from the LBV census.

5. Modelling the infrared spectral energy distributions with grey-body fitting

5.1. Individual stars

Several dusty LBV nebulae are extended. Optical instruments detected Hα emission, almost unextinguished by the dust grains (Agliozzo et al. 2012). Indeed the extinction AHα was mapped in three nebulae, revealing relatively low values and thus suggesting an optically thin configuration for the dust (Agliozzo et al. 2017a,b, 2019). For such sources we can assume that radiative transfer effects are negligible. Among the stars classified as Class 1, five have a detection in at least one PACS filter, namely RMC 71, RMC 127, RMC 143, LHA 120-S 61, and LHA 120-S 119. With the exception of LHA 120-S 119, the dust mass for these sources has been previously determined in the literature (see Table 6). For LHA 120-S 119, we thus fit the infrared SED with a greybody to derive the characteristic temperature, dust mass and β parameter (Fig. 6), using an absorption coefficient value κ850 μm of 1.7 cm2 g−1 (Sodroski et al. 1987), the same as Agliozzo et al. (2017a,b, 2019).

Table 6.

Summary of estimated dust properties.

thumbnail Fig. 6.

Example of single grey-body fit of the 2MASS, WISE, Spitzer, and Herschel data, similarly to Agliozzo et al. (2017a,b, 2019).

We assume a distance of 49.97 kpc for the LMC (Pietrzyński et al. 2013) for all sources. The distances to individual sources will vary around this nominal value due to the inclined geometry and finite thickness of the LMC, but the errors due to this single distance assumption are expected to be less than 2 kpc or 4% (van der Marel & Cioni 2001, and considering a maximum separation of any source from the LMC centre of 2°). This uncertainty (σD ∼ 1 kpc ∼ 2%) is an almost negligible contribution to the dust mass uncertainty compared to the systematic uncertainty in κ, and the typical measurement uncertainties of the far-infrared photometry.

Radiative transfer modelling was adopted to fit the data of RMC 71 by Guha Niyogi et al. (2014). For comparison, we also fit the data of RMC 71 with a grey-body and find a dust mass which agrees with the literature value. Three other sources (RMC 40, RMC 81, HD 269216) are detected in both the 24 and 70 μm MIPS bands but not in the PACS bands. Their measurements suggest that the thermal emission peaks between 24 and 70 μm. We fit these data by fixing the parameter β = 1.3 (the median of the β values determined for the five LMC LBVs above, see also Table 6). The distance assumed for the SMC LBV RMC 40 is 62.44 kpc (Graczyk et al. 2020). We note that the uncertainty of HD 269216 is nearly 100% due to the marginal detection in MIPS70. For the other Class 1 sources the lack of detection longward of 24 μm does not allow to constrain the dust parameters. The results of this analysis are summarised in Table 6.

5.2. Stack analysis of LBVs in the LMC

Given the low detection rates of individual sources at 70 μ m and longer, largely due to confusion, we have also carried out a stacking analysis to provide a statistical measure of total dust content. The stacking analysis adopted here consists in co-adding the images, at a given wavelength, of individual sources with known coordinates, in order to reduce the noise by averaging, resulting in an aggregate flux density measurement. This technique is particularly useful when the target sources are not detected. The stacking procedure is described in Kurczynski & Gawiser (2010), and references therein. Here, we do not adopt any deblending method described by the authors, as the confusion in our case is dominated by diffuse ISM emission, of which a-priori positional information is difficult to obtain. For the stacking analysis we considered all confirmed and candidate LBVs (SED Classes 1 and 2), thus excluding sgB[e] stars and the sources discarded in Sect. 4. The LMC stack sample thus comprises the first 18 entries in Table 1. In the SMC there are only two sources to consider (ID number 28 and 29 in the same table), of which only one clearly shows a dust component (number 29). For this reason, we perform the stacking analysis on the LBVs (Class 1 and 2 in Table 1) of the LMC only.

To obtain stacked images, for each source we created cutouts of the SAGE and HERITAGE mosaic of the LMC centred, to the nearest pixel, at the sky coordinates as in Table 1, after cross-checking the coordinates with the actual detections in the map. As a Wide-Field Infrared Survey Explorer (WISE) image of the whole LMC galaxy is not available, in this case we queried the WISE image service in the IRSA website and downloaded images centred at the sky coordinates of each source. For none of the instruments it was necessary to re-sample the images as they originate from the same grid. Images of the following bands were stacked: WISE W3, MIPS24, MIPS70, PACS100, PACS160, and SPIRE250. For each band we summed the cutouts of all sources. The resulting stacked images are shown in Fig. 7, while the associated photometry is summarised in Table 7. In the same figure, we also show the same stack images without RMC 143, which although the most massive LBV nebula recorded so far in the LMC (Agliozzo et al. 2019), is embedded in a very bright and crowded region, just south of the 30 Doradus star forming complex. The comparison of the resulting photometry allows us to assess the relative contribution of the other stars to the infrared emission.

thumbnail Fig. 7.

Left column: stacked images of LMC LBVs (SED Classes 1 and 2 in Table 1) and residuals after Gaussian fitting. Right column: same as before but without RMC143. The field of view corresponds approximately to 70″ × 70″. The stack image at 160 μm on the bottom right was derived with a smaller number of sources (IDs: 2, 5, 6, 7, 8, 10, 12, 13, 14, 15, 16, 17, 18) and used as a lower limit.

Table 7.

Flux densities from the stacked data.

Following the same method, we also stacked cutout images extracted from the LMC 1.4 mm SPT+Planck map of Crawford et al. (2016). However, we did not obtain any detection above the noise (∼250 mJy beam−1) and so do not show the resulting image here.

We then determined the aggregate flux density of the combined source by performing a 2D Gaussian fitting in each image. Peak height, x and y centre positions and background level are fitted parameters; the width parameters σa and σb were fixed to the instrumental beam σ for MIPS70, PACS100, and PACS160 images, while they were left as free parameters (along with a position angle) for the WISE W3 and MIPS24 fitting because at these wavelengths sources were slightly resolved. We used the non-linear least-squares method of the Scipy Optimize package (Virtanen et al. 2020) for minimising the fitting function residuals. To account for systematic errors in this procedure (e.g. missing flux in sidelobes), we corrected the values by a factor determined in each filter by comparing our own estimates of individual test sources with catalogue measurements. Figure 7 shows the stacked images and their residuals after subtraction of the 2D Gaussian. We do not include SPIRE 250 because there is no source detected, owing to source confusion over the large beam. For the stack analysis at the IRAC bands, we only summed the photometry from the catalogue, as all the individual sources were well detected. We expect the amount of flux density missed by the small catalogue aperture to be negligible since the IRAC images do not show hints of resolved nebular emission at these wavelengths.

For each passband, the resulting flux densities from the stack analysis are plotted in the SEDs in Fig. 8 as blue circles. The top figure visualises the sum of all 18 candidate and confirmed LBVs in the LMC, while the bottom shows the same without RMC 143. In the top figure the pink symbols represent lower limits. The 70 μm limit is given by the measurement in the stacked image without RMC 143, where the central source is better detected without the contamination by the 30 Dor region (as is evident from the comparison of the 70 μm error bars in the top and bottom panels). The 250 μm limit is due to a source detected at the central position of a stack image built from the only four well detected sources in this filter (namely RMC 127, LHA 120-S 119, LHA 120-S 61, and RMC 71). In the lower plot, the 160 μm lower limit is set based on the detection in the stack image obtained from fewer 13 sources (IDs: 2, 5, 6–8, 10, 12–18) that are less affected by crowding.

thumbnail Fig. 8.

Flux density distribution of stacks of the LBVs in the LMC (SED Classes 1 and 2 in Table 1), with (top) and without (bottom) RMC 143. 2MASS and IRAC values are the sum of the individual source fluxes, while WISE-W3, MIPS and PACS values are from photometry performed on the stacked images (Fig. 7). Parameters from the fitting procedure are shown (α and c are the spectral index and offset of the best-fit power-law describing the near-infrared data; Td, Md and β are the dust temperature, mass and opacity power law index of the best-fit greybody describing the thermal emission from dust). The lower limits (red points and arrow symbols) for MIPS70 and SPIRE250 in the upper plot and PACS160 in the lower plot are from stacks of subsets of sources (see text for the subsets in each case), as there is not a reliable detection in the stacked images of the whole sample.

The total SED from about 1 to 250 μm resembles that of LBV stars with dense ionised stellar winds, surrounded by extended dusty nebulae (Class 1a). We initially fit the 2MASS J, H, Ks and IRAC 3.6, 4.5 and 5.8 μm data with a power law representing the sum of stellar winds and photospheres (power-law free-free emission from ionised stellar winds, Panagia & Felli 1975, and Fν ∝ λ−2 Rayleigh-Jeans photosphere emission), and then we added the resulting best fit and uncertainties to the grey-body function built to fit the WISE W3, MIPS24, MIPS70, PACS100, and PACS160 measurements. In the figure, the best-fit total model is given by a black continuous line. The power-law fit of the wind component is the blue dashed line. The grey shaded region represents all the grey-body components that fit the data producing a χ2 which is twice the minimum value of the best fit. Because of the large uncertainty of the 70 μm produced by the crowded region around RMC 143, we set as lower limit the photometric measurement at 70 μm without RMC 143 to constrain the range of possible models. We also introduce another constraint at longer wavelengths by using the lower limit at 250 μm. The resulting dust temperature, mass, and parameter β are listed in the figures. The dust mass of the aggregate source is computed by assuming the flux density extracted from the best fit at 850 μm. The adopted absorption coefficient κ850 μm is 1.7 cm2 g−1 (Sodroski et al. 1987).

5.3. Stack analysis of other evolved stars in the LMC

The same stacking analysis was performed for the sample of RSGs and AGBs in the catalogue of Riebel et al. (2012), which contains 33 501 sources (6709 carbon-rich AGB, 19 566 oxygen-rich AGB, 5874 RSG and 1352 extreme “x-AGB” stars). Previously Dharmawardena (2019) stacked the PACS100 images of the sources in Riebel et al., split up by mass-loss rate bin and chemical class, and found some detections from carbon-rich and oxygen-rich AGB stars with the highest mass-loss rates. The wavelengths covered here are those of the four IRAC filters at 3.6, 4.5, 5.8 and 8 μm, the two MIPS filters at 24 and 70 μm, and the PACS filters at 100 and 160 μm. For completeness, we also considered the SPIRE image at 250 μm, although with low expectation of obtaining a detection because of source confusion at those wavelengths. We cross-matched the list of LBVs and the catalogues of Riebel et al. and found two common sources (RMC 66 and RMC 84), which we removed from their sample (they were catalogued as RSGs). Various attempts were made to obtain a detection at wavelengths longer than 24 μm. Initially, the images of individual categories of stars were stacked separately, with and without applying a confusion threshold on the basis of the standard deviation evaluated from all the pixels in each field (measuring about 50 × 50 pixels). Finally, all the 33 501 sources were summed together. In the 70 μm and the stack images at longer wavelength, we did not obtain any statistically meaningful detection for the RSGs, the C- and O-rich AGBs, or the whole sample.

In the case of the “extreme” AGBs, which are considered the main dust producers among all AGBs and RSGs in the MCs (Boyer et al. 2012; Riebel et al. 2012), we evaluated possible detections at 24, 100, 160, 250, 350 and 500 μm. We noticed that at far-infrared wavelengths the stacked images are dominated by confusion when the following fields are included: J053929.34–694614.8, J053911.41–690824.7, J054515.83–694648.1, J051913.89–693818.3, J045506.54–691708.6. We also found that five other sources have been identified as YSOs or other types of stars (J053840.77–690603.3, J053856.58–690417.3, J053839.68–690537.9, J045400.16–691155.4, J053238.60–682522.1, Evans et al. 2011; Walborn et al. 2013; van Gelder et al. 2020; Oliveira et al. 2019). We thus excluded these ten sources from further analysis. The final sample of extreme AGB stars from Riebel et al. (2012) that we stacked comprises 1342 objects. We obtained significant detections up to 100 μm, and a marginal detection at 160 μm, as shown in Fig. 9. To produce the composite continuum spectrum in the figure, no weighting based on the stellar luminosity or normalisation was applied. The spectrum simply comprises the aggregated flux densities measured from the stacked images at each wavelength.

thumbnail Fig. 9.

Stacked images of the extreme-AGBs from Riebel et al. (2012) and residuals after Gaussian source extraction. The field of view corresponds approximately to 70″ × 70″. Ten out of 1352 sources were excluded (see text), because the emission in their field is dominated by bright large-scale structures or because they are another type of stellar objects (e.g. YSOs).

We extracted the photometry as described in the preceding subsection. At 24 μm we compared our photometry extracted from the stack image with the stacked photometry (obtained by summing the photometry in Riebel et al. 2012) and found a consistent value. We fitted the data with two components: a black body of fixed temperature (3000 K) describing the stellar photosphere, and a grey body with dust spectral index β = 0. Because the near-infrared data of AGBs are extremely reddened by the optically thick dust, we also fitted the extinction by using the extinction curve of Fitzpatrick (1999), with selective parameter RV = 3.1. Figure 10 shows the photometry and resulting fit with and without extinction correction. For comparison, the integrated flux density at 100 μm from 18 LBVs is of the same order of magnitude of that from 1342 extreme-AGB stars (see Table 7). The dust mass derived from this fitting assumes that the dust is optically thin (at the wavelength at which κ is defined), and thus must be considered a lower limit.

thumbnail Fig. 10.

Flux density distribution of the 1342 extreme x-AGBs in the LMC from the catalogue by Riebel et al. (2012), constructed from the Gaussian fitting of the central detection in the stacked images (at wavelengths shorter than ∼10 μm the value is just the sum of the catalogue measurements). Blue points: observed phtometry; green points: de-reddened photometry; red arrow: 3σ PACS160 upper limit.

6. Discussion

6.1. LBV dust masses versus stellar luminosity and nebular size

The uncertainty on the distance is only ∼2% for the Magellanic sources (see Sect. 5.1), and thus their stellar luminosity is well constrained (∼4% uncertainty due to the distance). This is a great opportunity to investigate any possible dependence of the dust mass on the stellar parameters and on the nebular size. Figure 11 shows the dust masses estimated in Sect. 5.1 for 7 Class 1 LBVs in the LMC and RMC 40 in the SMC versus the stellar bolometric luminosity as in the literature (Table A.1). In the same plot, dust masses for some well-known nebulae around candidate and confirmed LBVs in the Milky Way are displayed. Dust production in these stars does not seem correlated to the stellar properties.

thumbnail Fig. 11.

Estimated dust masses for individual stars as a function of the stellar bolometric luminosity. Red circles: LMC LBVs (this work, see Table 6). Blue circles: Galactic LBVs, specifically AG Car (Vamvatira-Nakou et al. 2015); G24.73+0.69, G25.5+0.2, HD 168625, Hen 3–519, AFGL 2298 (Clark et al. 2003, and references therein); G26.47+0.02 (Umana et al. 2012); G79.49+0.26 (Agliozzo et al. 2014); HR Car (Machado et al. 2002; Buemi et al. 2017); Pistol star (Lau et al. 2014); Wra 751 (Vamvatira-Nakou et al. 2013). Error bars are shown when available from the literature.

For five objects the nebula size is known in the literature and one has been determined in this work. Sizes span from 0.3 to 2.2 pc (Table A.1). We do not observe a systematic trend in the dust mass versus the nebular size.

6.2. Total dust mass in the LMC

Since the ultimate goal of this work is to compare the LBV dust masses with the entire LMC dust yield we also fit the stack data with the simple grey model that Gordon et al. (2014) used to derive the global dust properties of the ISM in the LMC and SMC. In particular, we adopt their absorption coefficient κλ (calibrated with ultraviolet and optical gas-phase absorption measurements in the diffuse ISM of the Galaxy, which provide the expected amount of dust and are independent from the far-infrared and submillimetre dust emission). At 160 μm, Gordon et al. derived κ160 = 9.6 ± 0.4 ± 2.5 where the first error is statistical and the second is the systematic uncertainty due to the assumption that the dust in the MCs is like that in the Galaxy. With this κ160, our stack PACS160 flux and our stack SED fit temperature (top panel in Fig. 8), we obtain for the 18 LMC (Class 1 and 2) LBVs M dust = 0 . 11 0.03 + 0.06 M $ M_{\mathrm{dust}} = 0.11^{+0.06}_{-0.03}\,M_{\odot} $, and an average of M dust = 0 . 006 0.002 + 0.003 M $ \langle M_{\mathrm{dust}}\rangle = 0.006^{+0.003}_{-0.002}\,M _{\odot} $ per source. The uncertainty of about 25% in κλ due to the systematic error is not taken into account.

Table 8 summarises the dust masses derived from grey body fitting of the stack data. We notice that the integrated dust mass from 18 sources is of the same order of magnitude of the dust mass estimated in the Homunculus nebula of η Carinae (Morris et al. 2017). This difference may be the result of different metallicities. It may also suggest that several stars have still not experienced the required mass loss to form dust (giant eruptions) or that they do not have colliding winds within a close binary. It should be noted that only a handful of sources contribute significantly to the stacked dust mass (compare Tables 6 and 8). Additionally, the dust masses estimated here could be a lower limit, as cooler dust may have escaped the observations. Sensitive sub-millimetre observations with ALMA of all Class 1 and 2 objects would be an asset for further studies.

Table 8.

Summary of estimated LBV dust masses in the LMC.

For the LMC, Gordon et al. (2014) derived an integrated dust mass of (7.3 ± 1.7) × 105M (averaging the results from three different models, or (8.1 ± 2.8) × 105M with the simple grey body model only). A more recent analysis based on the same dataset reports estimates about a factor of two lower (Chastenet et al. 2017). In the classical scenario, the LBV phase is estimated to last around 104 − 105 yr, which is a small fraction compared to the age of the LMC (about 15 Gyr, Harris & Zaritsky 2009). To evaluate the potential contribution of LBVs to the ISM dust yield, we must account for the entire population of LBVs over the lifetime of the LMC.

This is not an easy task, because stellar evolution models are not able to predict how many BSGs become LBVs (or manifest the LBV phenomenon). Ideally, we would like to compare the ratio of LBVs and BSGs of similar luminosity, but finding a list near to completeness is challenging.

Urbaneja et al. (2017) analysed 90 LMC BSGs to study the gravity-luminosity relationship, but this number is probably a lower limit. In fact, in Bonanos et al. (2009) we count 133 stars earlier than B4 and of luminosity classes Ia and Iab in the LMC, excluding four objects that overlap with the LBV sample. However, Bonanos et al. (2009) could not claim completeness, as they targeted O, B supergiant stars with both accurate coordinates and spectral classification. We queried in Simbad for stars with spectral types earlier than B4 and of luminosity class I, within 2.5 deg from the LMC centre, and we find around 230 BSGs (several of these objects appear in Massey 2002). Massey (2010) estimated ∼6100 unevolved massive stars with M >  20 M, and ∼700 with M >  50 M.

6.2.1. Simple case: Constant number of LBVs over time

The simplest estimate we can make of the total dust mass produced by LBVs in the LMC over all time is by assuming that the number of LBVs has on average been the same as today. This simple assumption avoids having to consider the OB to LBV rate, but it does depend directly on the assumed LBV phase duration. For the range of accepted LBV phase durations of 104 − 105 yr, the LMC age of 15 Gyr, and the mass produced by the current LBV sample of 0 . 11 0.03 + 0.06 M $ 0.11^{+0.06}_{-0.03}\,M_\odot $, the range of total dust mass is from 1.2 × 104M to 1.7 × 105M, with a nominal value of 6 . 6 5.4 + 10.4 × 10 4 M $ 6.6^{+10.4}_{-5.4}\times 10^4\,M_\odot $ for an assumed LBV phase duration of 25 kyr.

6.2.2. Accounting for star formation history and initial mass function

We assume that the LBV population closely followed the fluctuations in the star formation history (SFH) of the LMC. The latter was computed by Harris & Zaritsky (2009), who derived star formation rates (SFR) in about 1400 LMC sub-regions in time bins of Δlog(t)=0.2. We integrate the SFR over these sub-regions, in order to obtain the total SFH for the LMC as in Fig. 11 of Harris & Zaritsky (2009). For the Initial Mass Function (IMF) we adopt the power law slope derived by Schneider et al. (2018), who found an excess of massive stars by 30% in the 30 Dor star forming region. We also use the Salpeter IMF as a possible lower limit in the LMC. Finally, we derive the total dust yield by LBVs as,

( i SFR ( t i ) Δ t i ) × f LBV , MS × 25 M 100 M ξ ( m ) d m 0.1 M 100 M m ξ ( m ) d m × M dust , $$ \begin{aligned} \frac{\left(\displaystyle \sum _{i} \mathrm{SFR}(t_{i})\Delta t_i\right) \times f_{\rm LBV, MS}\times \displaystyle \int _{25\,M_\odot }^{100\,M_\odot } \xi (m) \,\mathrm{d}m}{\displaystyle \int _{0.1\,M_\odot }^{100\,M_\odot } m \xi (m) \,\mathrm{d}m} \times \langle M_{\rm dust}\rangle , \end{aligned} $$

where ξ(m)=ξ0 × mγ is the IMF, SFR(ti) is the SFR in the time bin i, ⟨Mdust⟩ is the mean dust mass produced per LBV, and fLBV, MS is the fraction of stars in mass range 25−100 M experiencing the LBV phenomenon. The IMF is normalised to the total mass of stars with masses between 0.1 and 100 M. Assuming the slope γ = −1.9 from Schneider et al. (2018) (Case 1) and γ = −2.35 from Salpeter (Case 2), we find the following ranges of dust mass: 17 . 1 3.7 + 14.0 × 10 4 × f LBV , MS M $ 17.1^{+14.0}_{-3.7}\times 10^{4}\times f_{\mathrm{LBV, MS}}\,M_{\odot} $ (Case 1) and 4 . 4 1.0 + 3.6 × 10 4 × f LBV , MS M $ 4.4^{+3.6}_{-1.0}\times 10^{4}\times f_{\mathrm{LBV, MS}}\,M_{\odot} $ (Case 2). The error bars include the uncertainty on the dust mass and on the SFR. A major uncertainty in this relation lies in the fraction of massive main-sequence stars that will become LBV stars (fLBV, MS). Obviously, the number cannot exceed unity. Because the conditions leading to the LBV phenomenon are unknown, it is unconstrained by evolutionary models.

The observed fraction with respect to BSGs, fLBV, obs = nLBV, obs/nBSG, obs, needs to be corrected for incompleteness and, more importantly, for the ratio of the duration of the BSG and LBV phases, tBSG/tLBV. Evolutionary models suggest tBSG ≤ 106 yr for stars more massive than 25 M in the He-burning phase, and tMS ≤ 107 yr for stars in the H-burning phase (Maeder & Meynet 1989). While the number of stars in the LMC currently exhibiting the LBV phenomenon should be reasonably correct (nLBV, obs), it should be subject to fluctuations, and some current LBVs might be inactive. The number of BSGs in the LMC that may become LBVs (nBSG, obs ∼ 230, see above) is also uncertain. The number of MS O,B stars (6700) may be more reliable (Massey 2002, 2010). From the currently available information, it seems plausible that fLBV, MS ranges between about 30%5 and 100%6.

We repeat a similar calculation for the AGB/RSG stars. We use the average dust mass per star determined from the stack of all extreme-AGBs ( 9 3.0 + 2.2 × 10 8 M $ 9^{+2.2}_{-3.0}\times 10^{-8}\,M_{\odot} $) and we integrate over stellar masses between 1 and 25 M to also include RSGs. We find a total dust mass of 136 35 + 56 M $ 136^{+56}_{-35}\,M_{\odot} $ (Case 1) or 144 37 + 59 M $ 144^{+59}_{-37}\,M_{\odot} $ (Case 2), which are much lower values than LBVs, but maybe consistent with the long dust injection timescales of AGBs (∼2 × 109 yr−1Micelotta et al. 2018). Moreover, the dust lifetime against processing and destruction by supernova shocks in the ISM is shorter, around 2 × 107 yr (Temim et al. 2015). The inefficiency of existing AGB stars to produce the dust mass observed in the ISM was already pointed out by Matsuura et al. (2009) for the LMC and by Boyer et al. (2012) and Matsuura et al. (2013) for the SMC. Dust injection timescales from LBVs (and subsequent SNe) are much shorter (≤10−100 Myr). In their study of the total dust input from stars in the SMC Boyer et al. (2012) estimated a dust production rate of ∼10−6M yr−1 from LBVs, matching the input from all cool stars, but concluded that the grains would not survive the subsequent supernova explosion.

6.3. Potential dust survival from SN explosions

Now we evaluate how likely this LBV dust is to survive the SN forward shock (blast wave) which will follow the LBV or WR evolutionary stage. The case of ISM grains processing by the forward shock was analysed by Barlow (1978), Nozawa et al. (2006), and Bocchio et al. (2014). The processes of interest are generally sputtering and shattering. There are a few arguments to mention in favour of dust survival in the SN blast wave in LBVNe. (1) LBV eruptions are a source of large grains (> 0.1 μm) due to the favourable conditions for particle growth (Kochanek 2011a). The relatively low β values found in this work are also consistent with large grains. However, they could also be low due to a temperature gradient over the nebula. (2) The material in LBVNe seems often distributed in a clumpy configuration, spread at large distances from the star. (3) Such nebulae can be asymmetric (e.g. RMC 143, LHA 120-S 119), subtending solid angles of less than 4π so that some of the SN energy will leak out without interacting with the nebula.

Such properties reduce the importance of the destructive action by the shock, as for example evidenced in studies of the newly formed dust in SN ejecta impacted by the reverse shock. In particular, survival rates vary between 40 and 98% (Nath et al. 2008; Silvia et al. 2010, 2012; Biscaro & Cherchneff 2016; Gall & Hjorth 2018) for dust reached by the reverse shock, and around 85% for large grains in the ISM (e.g. Nozawa et al. 2006). More results are summarised in Table 2 of Micelotta et al. (2018). Some studies also show that shocks can induce dust (re)formation and growth (e.g. Kochanek 2011b; Matsuura et al. 2019). Recent 3D hydrodynamical simulations report that the wind-driven bubbles created during the entire life of a massive star create the conditions for the dust to largely survive (Martínez-González et al. 2019), or that under favourable ambient conditions implantation and trapping of ions in dust grains can counteract grain destruction by sputtering (Kirchschlager et al. 2020). Dust destruction is particularly negligible if the pre-existing shell is massive (Martínez-González et al. 2019), which is the case of LBV nebulae.

From the observational point of view, a great example is SN1987A, the progenitor of which may have been an LBV (as speculated by Smith 2007). Along with the discovery of dust formation in the SN ejecta a few years after the explosion (Danziger et al. 1989), followed up by far-infrared and sub-millimetre observations about three decades later (which revealed large amounts ∼0.5 M of cold dust, Matsuura et al. 2011; Indebetouw et al. 2014; Wesson et al. 2015), another dust component distributed in an equatorial ring lost by the progenitor 20−40 kyr ago was found (Arnett et al. 1989; McCray 1993; Bouchet et al. 2006). The fastest part of the forward shock has passed the ring (Fransson et al. 2015), offering unique insights into the resulting dust destructive processing. The silicate grains in the ring are collisionally heated by the SN blast wave (Bouchet et al. 2006). Arendt et al. (2016) proposed that some dust grains were destroyed by sputtering, supporting the theory (Dwek et al. 2008, 2010). Subsequently, Matsuura et al. (2019) reported a 30−70 μm excess at day ∼10 000 on top of the earlier modelled hot and warm component (in the ring) and cold component (in the SN ejecta), and a dust mass about a factor of 10 larger than previous estimates at day ∼8000. They identified two possible mechanisms behind this infrared excess: dust re-formation (dust grains formed from the gas phase) and dust growth (survived grains accreting atoms from the gas phase). Dust re-formation was also suggested in the cool dense shell created between the forward and reverse shocks as the SN ejecta plows into the existing circumstellar material. Several studies of type IIn SNe (which, in some cases, may have an LBV as progenitor) are a good example (e.g. Smith et al. 2008a; Fox et al. 2009; Gall et al. 2014; Andrews et al. 2016; Chugai 2018). Grain growth may also be a viable mechanism as explained by Kirchschlager et al. (2020).

7. Summary and conclusions

We have analysed the infrared images of Magellanic LBVs available in the archives, adding for a few sources new photometry that either was missing or not accurate because of extended objects or crowded fields. We have revised the Magellanic LBV sample on the basis of a new infrared SED-based classification scheme and literature information. The sources can be divided in four classes. All strongly active LBVs show either a cool dust excess peaking in the mid-infrared (Class 1a and 1b) or an excess of near-infrared emission due to free-free emission from the ionised stellar wind (Class 2), or both (Class 1a). These features are also often present in several candidate LBVs, which appear in the literature as ex-/dormant LBVs. Four sources can be removed from our initial sample, as we find that in past studies these stars did not satisfy the LBV requirements and, in addition, they do not present any of the features mentioned above (Class 4). Six other sources are known as sgB[e] (Class 3): their flux density distribution does not resemble those of confirmed LBVs, consistent with a different mechanism of dust formation, as suggested in the literature. For this reason we excluded them from further analysis.

We have employed a simple grey-body fitting method to model the infrared SED of individual sources. In the LMC, which has the most numerous list of LBVs, large amounts of dust are observed (∼10−3 − 10−2M), similar to Galactic LBVNe. We stacked the infrared images of the LMC LBVs and extracted the photometry of the resulting source, detected up to 160−250 μm. The integrated SED from the stacks resembles that of LBVs with a strong ionised stellar wind and an extended dusty nebula. The SED can be fitted with only two components: a power-law describing the free-free spectrum of ionised stellar winds and the stellar photosphere, and a single-component grey-body for the dust. For the grey-body we adopted two different values of the κ parameter, including the value determined by Gordon et al. (2014) to fit the integrated ISM SED of the LMC. A significant contribution to the stack SED comes from a few sources, the most important one is RMC143. This was already identified as a massive nebula (Agliozzo et al. 2019). We obtain an integrated present dust mass of 0 . 11 0.03 + 0.06 M $ 0.11^{+0.06}_{-0.03}\,M_{\odot} $. We have repeated a similar analysis on the sample of AGBs and RSGs by Riebel et al. (2012). We obtain a detection in the stacked images only when considering the “extreme”-AGBs. We find that the integrated 160 μm emission of 1342 extreme-AGBs is of the same order of magnitude as that of 18 LBVs. The integrated dust mass from these sources is 1 . 2 0.4 + 0.3 × 10 4 M $ 1.2^{+0.3}_{-0.4}\times 10^{-4}\,M_{\odot} $. We do not find any correlation between the dust masses and the stellar luminosities. This could be due to the fact that such stars have different evolutionary histories or that the dust production mechanism does not depend on the initial mass of the star. Most likely we are also unable to detect the lowest-mass nebulae. To estimate the total dust mass produced by LBVs in the LMC during its full lifetime, we consider two cases: constant number of LBVs across time and a case accounting for IMF and SFH. The uncertainty on the duration of LBV phase in the first case, or on the population of LBVs in the second case, add a significant uncertainty in the total mass produced by LBVs.

The range of total dust masses produced by LBVs over the lifetime of the LMC spans between 104 and 105M, suggesting that LBVs are potentially significant dust factories in the LMC. The LBV dust masses per source of 10−3 − 10−2M are at least one order of magnitude smaller than the amount of dust (0.1−1 M) postulated in SN ejecta by Todini & Ferrara (2001) and Nozawa et al. (2003). However, those works adopted ideal assumptions for the seed material condensing into dust grains. Mid-infrared observations of CC-SNe suggest smaller dust masses (10−4 − 10−2M, Kotak et al. 2009; Fox et al. 2011, although estimates from single epoch observations of the rapidly evolving SN ejecta may not be representative of the total dust produced or mid-infrared observations may have not captured the emission from colder dust). On the other hand ALMA sub-millimetre observations revealed larger amounts of dust in both SN1987A (0.2 M, Indebetouw et al. 2014) and LBV RMC143 (0.055 M, Agliozzo et al. 2019), however these works assume different opacities at 850 μm, making the comparison difficult. Still the LBV phenomenon has the potential to be the second most important source of dust, before low- and intermediate-mass stars. This was also the conclusion by Boyer et al. (2012) in the SMC. It should be noted that if, on the one hand, the LBV dust mass has a large uncertainty after analysing a relatively large sample, in the case of CC-SNe the dust mass is also very uncertain and based on observations of only a handful of objects. Destruction of dust in LBVNe by a possible subsequent SN blast wave remains to be determined. A significant fraction of this dust mass may survive thanks to several properties of these nebulae, such as large nebular masses, big grains, clumpiness and asymmetry. Dust grain (re)formation or growth may also occur, like in the case of SN1987A (Matsuura et al. 2019).


5

Comparing the LBVs with the MS OBs, taking 107 yr as the MS lifetime and an upper limit of 105 yr for the LBV phase duration gives, fLBV = (18 LBVs/6700 MS OBs) × (107 yr/105 yr)=27%.

6

Comparing the LBVs with the BSGs, taking 106 yr as the BSG phase duration and 105 yr for the LBV phase duration gives, (18 LBVs/230 BSGs) × (106 yr/105 yr)=78%. As the LBV phase may be an order of magnitude shorter than this, it suggests that the number of 230 BSGs may be incomplete.

Acknowledgments

D.A. acknowledges funding through the European Union’s Horizon 2020 and Innovation programme under the Marie Skłodowska-Curie grant agreement no. 793499 (DUSTDEVILS). F.K. is supported by the Ministry of Science and Technology (MoST) of Taiwan, under grant number MOST107-2119-M-001-031-MY3, and also by Academia Sinica, in the form of Investigator Award AS-IA-106-MO3. GP is supported by ANID – Millennium Science Initiative – ICN12_009. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2017.A.00054.S and ADS/JAO.ALMA#2018.A.00049.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan) and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This paper also includes data collected: at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programmes 096.D-0047(A), 097.D-0612(A,B), and 0100.D-0469(A,B). This work made use of PyAstronomy. This research has made use of the International Variable Star Index (VSX) database, operated at AAVSO, Cambridge, Massachusetts, USA, and the VizieR catalogue access tool, CDS, Strasbourg, France. The original description of the VizieR service was published in A&AS 143, 23. This research has made use of the NASA/IPAC Infrared Science Archive, which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology. This research made use of Photutils, an Astropy package for detection and photometry of astronomical sources (Bradley et al. 2021). The Cornell Atlas of Spitzer/IRS Sources (CASSIS) is a product of the Infrared Science Center at Cornell University, supported by NASA and JPL.

References

  1. Agliozzo, C., Umana, G., Trigilio, C., et al. 2012, MNRAS, 426, 181 [NASA ADS] [CrossRef] [Google Scholar]
  2. Agliozzo, C., Noriega-Crespo, A., Umana, G., et al. 2014, MNRAS, 440, 1391 [NASA ADS] [CrossRef] [Google Scholar]
  3. Agliozzo, C., Nikutta, R., Pignata, G., et al. 2017a, MNRAS, 466, 213 [NASA ADS] [CrossRef] [Google Scholar]
  4. Agliozzo, C., Trigilio, C., Pignata, G., et al. 2017b, ApJ, 841, 130 [NASA ADS] [CrossRef] [Google Scholar]
  5. Agliozzo, C., Mehner, A., Phillips, N. M., et al. 2019, A&A, 626, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Andrews, J. E., Krafton, K. M., Clayton, G. C., et al. 2016, MNRAS, 457, 3241 [NASA ADS] [CrossRef] [Google Scholar]
  7. Ardila, D. R., Van Dyk, S. D., Makowiecki, W., et al. 2010, ApJS, 191, 301 [NASA ADS] [CrossRef] [Google Scholar]
  8. Arendt, R. G., Dwek, E., Bouchet, P., et al. 2016, AJ, 151, 62 [NASA ADS] [CrossRef] [Google Scholar]
  9. Aret, A., Kraus, M., Muratore, M. F., & Borges Fernandes, M. 2012, MNRAS, 423, 284 [NASA ADS] [CrossRef] [Google Scholar]
  10. Arneson, R. A., Shenoy, D., Smith, N., & Gehrz, R. D. 2018, ApJ, 864, 31 [NASA ADS] [CrossRef] [Google Scholar]
  11. Arnett, W. D., Bahcall, J. N., Kirshner, R. P., & Woosley, S. E. 1989, ARA&A, 27, 629 [NASA ADS] [CrossRef] [Google Scholar]
  12. Barba, R. H., Morrell, N. I., Niemela, V. S., et al. 1996, in Revista Mexicana de Astronomia y Astrofisica Conference Series, eds. V. Niemela, N. Morrell, P. Pismis, & S. Torres-Peimbert, 5, 85 [NASA ADS] [Google Scholar]
  13. Barlow, M. J. 1978, MNRAS, 183, 367 [NASA ADS] [CrossRef] [Google Scholar]
  14. Beichman, C. A., Neugebauer, G., Habing, H. J., Clegg, P. E., & Chester, T. J. 1988, Infrared Astronomical Satellite (IRAS) Catalogs and Atlases. Volume 1: Explanatory Supplement (The Joint IRAS Science Working Group), 1 [Google Scholar]
  15. Biscaro, C., & Cherchneff, I. 2016, A&A, 589, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Bocchio, M., Jones, A. P., & Slavin, J. D. 2014, A&A, 570, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Boffin, H. M. J., Rivinius, T., Mérand, A., et al. 2016, A&A, 593, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Boian, I., & Groh, J. H. 2018, A&A, 617, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Bonanos, A. Z., Massa, D. L., Sewilo, M., et al. 2009, AJ, 138, 1003 [Google Scholar]
  20. Bonanos, A. Z., Lennon, D. J., Köhlinger, F., et al. 2010, AJ, 140, 416 [NASA ADS] [CrossRef] [Google Scholar]
  21. Bouchet, P., Dwek, E., Danziger, J., et al. 2006, ApJ, 650, 212 [NASA ADS] [CrossRef] [Google Scholar]
  22. Boyer, M. L., Srinivasan, S., Riebel, D., et al. 2012, ApJ, 748, 40 [Google Scholar]
  23. Bradley, L., Sipőcz , B., Robitaille, T., et al. 2021, https://www.doi.org/10.5281/zenodo.4624996 [Google Scholar]
  24. Buemi, C. S., Umana, G., Trigilio, C., Leto, P., & Hora, J. L. 2010, ApJ, 721, 1404 [NASA ADS] [CrossRef] [Google Scholar]
  25. Buemi, C. S., Trigilio, C., Leto, P., et al. 2017, MNRAS, 465, 4147 [NASA ADS] [CrossRef] [Google Scholar]
  26. Campagnolo, J. C. N., Borges Fernandes, M., Drake, N. A., et al. 2018, A&A, 613, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157 [Google Scholar]
  28. Chastenet, J., Bot, C., Gordon, K. D., et al. 2017, A&A, 601, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Chugai, N. N. 2018, MNRAS, 481, 3643 [CrossRef] [Google Scholar]
  30. Clark, J. S., Egan, M. P., Crowther, P. A., et al. 2003, A&A, 412, 185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Clark, J. S., Bartlett, E. S., Coe, M. J., et al. 2013, A&A, 560, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Cohen, M., Walker, R. G., Carter, B., et al. 1999, AJ, 117, 1864 [Google Scholar]
  33. Conti, P. S. 1984, in Observational Tests of the Stellar Evolution Theory, eds. A. Maeder, & A. Renzini, 105, 233 [NASA ADS] [CrossRef] [Google Scholar]
  34. Crowther, P. A. 2007, ARA&A, 45, 177 [Google Scholar]
  35. Crowther, P. A., & Smith, L. J. 1997, A&A, 320, 500 [NASA ADS] [Google Scholar]
  36. Crawford, T. M., Chown, R., Holder, G. P., et al. 2016, ApJS, 227, 23 [NASA ADS] [CrossRef] [Google Scholar]
  37. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  38. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2012, VizieR Online Data Catalog: II/281 [Google Scholar]
  39. Cutri, R. M., Wright, E. L., Conrow, T., et al. 2021, VizieR Online Data Catalog: II/328 [Google Scholar]
  40. Danforth, C. W., & Chu, Y.-H. 2001, ApJ, 552, L155 [CrossRef] [Google Scholar]
  41. Danziger, I. J., Gouiffes, C., Bouchet, P., & Lucy, L. B. 1989, IAU Circ., 4746, 1 [NASA ADS] [Google Scholar]
  42. de Wit, W. J., Oudmaijer, R. D., & Vink, J. S. 2014, Adv. Astron., 2014, 270848 [CrossRef] [Google Scholar]
  43. Dharmawardena, T. E. 2019, PhD Thesis, National Central University, Taiwan [Google Scholar]
  44. Dole, H., Lagache, G., Puget, J. L., et al. 2006, A&A, 451, 417 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Draine, B. T. 2009, in Cosmic Dust – Near and Far, eds. T. Henning, E. Grün, & J. Steinacker, ASP Conf. Ser., 414, 453 [Google Scholar]
  46. Dwek, E., Arendt, R. G., Bouchet, P., et al. 2008, ApJ, 676, 1029 [NASA ADS] [CrossRef] [Google Scholar]
  47. Dwek, E., Arendt, R. G., Bouchet, P., et al. 2010, ApJ, 722, 425 [NASA ADS] [CrossRef] [Google Scholar]
  48. Egan, M. P., Clark, J. S., Mizuno, D. R., et al. 2002, ApJ, 572, 288 [Google Scholar]
  49. Egan, M. P., Price, S. D., Kraemer, K. E., et al. 2003, VizieR Online Data Catalog: V/114 [Google Scholar]
  50. Elias-Rosa, N., Benetti, S., Cappellaro, E., et al. 2018, MNRAS, 475, 2614 [NASA ADS] [CrossRef] [Google Scholar]
  51. Enard, D. 1991, J. Opt., 22, 33 [NASA ADS] [CrossRef] [Google Scholar]
  52. Evans, C. J., Crowther, P. A., Fullerton, A. W., & Hillier, D. J. 2004, ApJ, 610, 1021 [NASA ADS] [CrossRef] [Google Scholar]
  53. Evans, C. J., Taylor, W. D., Hénault-Brunet, V., et al. 2011, A&A, 530, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Fariña, C., Bosch, G. L., Morrell, N. I., Barbá, R. H., & Walborn, N. R. 2009, AJ, 138, 510 [CrossRef] [Google Scholar]
  55. Feast, M. W., Thackeray, A. D., & Wesselink, A. J. 1960, MNRAS, 121, 337 [NASA ADS] [CrossRef] [Google Scholar]
  56. Fitzpatrick, E. L. 1991, PASP, 103, 1123 [NASA ADS] [CrossRef] [Google Scholar]
  57. Fitzpatrick, E. L. 1999, PASP, 111, 63 [Google Scholar]
  58. Fox, O., Skrutskie, M. F., Chevalier, R. A., et al. 2009, ApJ, 691, 650 [NASA ADS] [CrossRef] [Google Scholar]
  59. Fox, O. D., Chevalier, R. A., Skrutskie, M. F., et al. 2011, ApJ, 741, 7 [NASA ADS] [CrossRef] [Google Scholar]
  60. Fransson, C., Larsson, J., Migotto, K., et al. 2015, ApJ, 806, L19 [NASA ADS] [CrossRef] [Google Scholar]
  61. Gaia Collaboration 2018, VizieR Online Data Catalog: I/345 [Google Scholar]
  62. Gal-Yam, A., & Leonard, D. C. 2009, Nature, 458, 865 [NASA ADS] [CrossRef] [Google Scholar]
  63. Gall, C., & Hjorth, J. 2018, ApJ, 868, 62 [CrossRef] [Google Scholar]
  64. Gall, C., Hjorth, J., Watson, D., et al. 2014, Nature, 511, 326 [CrossRef] [Google Scholar]
  65. Gallagher, J. S. 1989, in IAU Colloq. 113: Physics of Luminous Blue Variables, eds. K. Davidson, A. F. J. Moffat, & H. J. G. L. M. Lamers, Astrophys. Space Sci. Lib., 157, 185 [CrossRef] [Google Scholar]
  66. Gordon, K. D., Roman-Duval, J., Bot, C., et al. 2014, ApJ, 797, 85 [NASA ADS] [CrossRef] [Google Scholar]
  67. Graczyk, D., Pietrzyński, G., Thompson, I. B., et al. 2020, ApJ, 904, 13 [Google Scholar]
  68. Gräfener, G., Owocki, S. P., & Vink, J. S. 2012, A&A, 538, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Grassitelli, L., Langer, N., Mackey, J., et al. 2021, A&A, 647, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3 [EDP Sciences] [Google Scholar]
  71. Groh, J. H. 2014, A&A, 572, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Groh, J. H., Damineli, A., Hillier, D. J., et al. 2009, ApJ, 705, L25 [NASA ADS] [CrossRef] [Google Scholar]
  73. Groh, J. H., Meynet, G., & Ekström, S. 2013, A&A, 550, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Guha Niyogi, S., Min, M., Meixner, M., et al. 2014, A&A, 569, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Gvaramadze, V. V., Kniazev, A. Y., & Fabrika, S. 2010, MNRAS, 405, 1047 [NASA ADS] [Google Scholar]
  76. Gvaramadze, V. V., Kniazev, A. Y., Maryeva, O. V., & Berdnikov, L. N. 2018, MNRAS, 474, 1412 [NASA ADS] [CrossRef] [Google Scholar]
  77. Hainich, R., Rühling, U., Todt, H., et al. 2014, A&A, 565, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Harris, J., & Zaritsky, D. 2009, AJ, 138, 1243 [NASA ADS] [CrossRef] [Google Scholar]
  79. Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 [NASA ADS] [CrossRef] [Google Scholar]
  80. Humphreys, R. M., & Davidson, K. 1994, PASP, 106, 1025 [NASA ADS] [CrossRef] [Google Scholar]
  81. Hutsemekers, D. 1997, in Luminous Blue Variables: Massive Stars in Transition, eds. A. Nota, & H. Lamers, ASP Conf. Ser., 120, 316 [NASA ADS] [Google Scholar]
  82. Indebetouw, R., Matsuura, M., Dwek, E., et al. 2014, ApJ, 782, L2 [NASA ADS] [CrossRef] [Google Scholar]
  83. Ishihara, D., Onaka, T., Kataza, H., et al. 2010, A&A, 514, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Jones, O. C., Meixner, M., Sargent, B. A., et al. 2015, ApJ, 811, 145 [NASA ADS] [CrossRef] [Google Scholar]
  85. Jones, O. C., Woods, P. M., Kemper, F., et al. 2017, MNRAS, 470, 3250 [NASA ADS] [CrossRef] [Google Scholar]
  86. Kastner, J. H., Buchanan, C. L., Sargent, B., & Forrest, W. J. 2006, ApJ, 638, L29 [NASA ADS] [CrossRef] [Google Scholar]
  87. Kastner, J. H., Buchanan, C., Sahai, R., Forrest, W. J., & Sargent, B. A. 2010, AJ, 139, 1993 [NASA ADS] [CrossRef] [Google Scholar]
  88. Kato, D., Ita, Y., Onaka, T., et al. 2012, AJ, 144, 179 [CrossRef] [Google Scholar]
  89. Käufl, H. U., Kerber, F., Asmus, D., et al. 2015, The Messenger, 159, 15 [NASA ADS] [Google Scholar]
  90. Kemper, F., Woods, P. M., Antoniou, V., et al. 2010, PASP, 122, 683 [NASA ADS] [CrossRef] [Google Scholar]
  91. Kerber, F., Käufl, H. U., Tristram, K., et al. 2016, SPIE Conf. Ser., 9908, 99080D [NASA ADS] [Google Scholar]
  92. Kirchschlager, F., Barlow, M. J., & Schmidt, F. D. 2020, ApJ, 893, 70 [CrossRef] [Google Scholar]
  93. Kochanek, C. S. 2011a, ApJ, 743, 73 [NASA ADS] [CrossRef] [Google Scholar]
  94. Kochanek, C. S. 2011b, ApJ, 741, 37 [NASA ADS] [CrossRef] [Google Scholar]
  95. Koenigsberger, G., Auer, L. H., Georgiev, L., & Guinan, E. 1998, ApJ, 496, 934 [NASA ADS] [CrossRef] [Google Scholar]
  96. Koenigsberger, G., Morrell, N., Hillier, D. J., et al. 2014, AJ, 148, 62 [Google Scholar]
  97. Kotak, R., & Vink, J. S. 2006, A&A, 460, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Kotak, R., Meikle, W. P. S., Farrah, D., et al. 2009, ApJ, 704, 306 [NASA ADS] [CrossRef] [Google Scholar]
  99. Kurczynski, P., & Gawiser, E. 2010, AJ, 139, 1592 [Google Scholar]
  100. Lagage, P. O., Pel, J. W., Authier, M., et al. 2004, The Messenger, 117, 12 [NASA ADS] [Google Scholar]
  101. Laporte, N., Ellis, R. S., Boone, F., et al. 2017, ApJ, 837, L21 [CrossRef] [Google Scholar]
  102. Lau, R. M., Herter, T. L., Morris, M. R., & Adams, J. D. 2014, ApJ, 785, 120 [NASA ADS] [CrossRef] [Google Scholar]
  103. Lau, R. M., Eldridge, J. J., Hankins, M. J., et al. 2020, ApJ, 898, 74 [CrossRef] [Google Scholar]
  104. Lebouteiller, V., Barry, D. J., Spoon, H. W. W., et al. 2011, ApJS, 196, 8 [NASA ADS] [CrossRef] [Google Scholar]
  105. Lebouteiller, V., Barry, D. J., Goes, C., et al. 2015, ApJS, 218, 21 [Google Scholar]
  106. Machado, M. A. D., de Araújo, F. X., Pereira, C. B., & Fernandes, M. B. 2002, A&A, 387, 151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  107. Maeder, A. 1989, in IAU Colloq. 113: Physics of Luminous Blue Variables, eds. K. Davidson, A. F. J. Moffat, & H. J. G. L. M. Lamers, Astrophys. Space Sci. Lib., 157, 15 [NASA ADS] [CrossRef] [Google Scholar]
  108. Maeder, A., & Meynet, G. 1989, A&A, 210, 155 [Google Scholar]
  109. Margutti, R., Milisavljevic, D., Soderberg, A. M., et al. 2014, ApJ, 780, 21 [Google Scholar]
  110. Martínez-González, S., Wünsch, R., Silich, S., et al. 2019, ApJ, 887, 198 [CrossRef] [Google Scholar]
  111. Massey, P. 2002, ApJS, 141, 81 [NASA ADS] [CrossRef] [Google Scholar]
  112. Massey, P. 2010, in Hot and Cool: Bridging Gaps in Massive Star Evolution, eds. C. Leitherer, P. D. Bennett, P. W. Morris, & J. T. Van Loon, ASP Conf. Ser., 425, 3 [NASA ADS] [Google Scholar]
  113. Massey, P., McNeill, R. T., Olsen, K. A. G., et al. 2007, AJ, 134, 2474 [NASA ADS] [CrossRef] [Google Scholar]
  114. Matsuura, M., Barlow, M. J., Zijlstra, A. A., et al. 2009, MNRAS, 396, 918 [NASA ADS] [CrossRef] [Google Scholar]
  115. Matsuura, M., Dwek, E., Meixner, M., et al. 2011, Science, 333, 1258 [NASA ADS] [CrossRef] [Google Scholar]
  116. Matsuura, M., Woods, P. M., & Owen, P. J. 2013, MNRAS, 429, 2527 [NASA ADS] [CrossRef] [Google Scholar]
  117. Matsuura, M., Dwek, E., Barlow, M. J., et al. 2015, ApJ, 800, 50 [NASA ADS] [CrossRef] [Google Scholar]
  118. Matsuura, M., De Buizer, J. M., Arendt, R. G., et al. 2019, MNRAS, 482, 1715 [NASA ADS] [CrossRef] [Google Scholar]
  119. Mauerhan, J. C., Smith, N., Filippenko, A. V., et al. 2013, MNRAS, 430, 1801 [NASA ADS] [CrossRef] [Google Scholar]
  120. McCray, R. 1993, ARA&A, 31, 175 [NASA ADS] [CrossRef] [Google Scholar]
  121. McKinnon, R., Torrey, P., & Vogelsberger, M. 2016, MNRAS, 457, 3775 [CrossRef] [Google Scholar]
  122. Mehner, A., Baade, D., Rivinius, T., et al. 2013, A&A, 555, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  123. Mehner, A., Baade, D., Groh, J. H., et al. 2017, A&A, 608, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Meixner, M., Gordon, K. D., Indebetouw, R., et al. 2006, AJ, 132, 2268 [NASA ADS] [CrossRef] [Google Scholar]
  125. Meixner, M., Panuzzo, P., Roman-Duval, J., et al. 2013, AJ, 146, 62 [NASA ADS] [CrossRef] [Google Scholar]
  126. Micelotta, E. R., Matsuura, M., & Sarangi, A. 2018, Space Sci. Rev., 214, 53 [NASA ADS] [CrossRef] [Google Scholar]
  127. Morris, P. W., Gull, T. R., Hillier, D. J., et al. 2017, ApJ, 842, 79 [NASA ADS] [CrossRef] [Google Scholar]
  128. Muratorio, G., & Friedjung, M. 1988, A&A, 190, 103 [NASA ADS] [Google Scholar]
  129. Nath, B. B., Laskar, T., & Shull, J. M. 2008, ApJ, 682, 1055 [NASA ADS] [CrossRef] [Google Scholar]
  130. Neugent, K. F., Massey, P., Skiff, B., & Meynet, G. 2012, ApJ, 749, 177 [NASA ADS] [CrossRef] [Google Scholar]
  131. Nozawa, T., Kozasa, T., Umeda, H., Maeda, K., & Nomoto, K. 2003, ApJ, 598, 785 [NASA ADS] [CrossRef] [Google Scholar]
  132. Nozawa, T., Kozasa, T., & Habe, A. 2006, ApJ, 648, 435 [NASA ADS] [CrossRef] [Google Scholar]
  133. Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  134. Oliveira, J. M., van Loon, J. T., Sewiło, M., et al. 2019, MNRAS, 490, 3909 [NASA ADS] [CrossRef] [Google Scholar]
  135. Osmer, P. S. 1973, ApJ, 186, 459 [NASA ADS] [CrossRef] [Google Scholar]
  136. Owocki, S. P., Townsend, R. H. D., & Quataert, E. 2017, MNRAS, 472, 3749 [Google Scholar]
  137. Panagia, N., & Felli, M. 1975, A&A, 39, 1 [Google Scholar]
  138. Pasquali, A., Schmutz, W., Nota, A., & Origlia, L. 1997, A&A, 327, 265 [NASA ADS] [Google Scholar]
  139. Pastorello, A., Kochanek, C. S., Fraser, M., et al. 2017, MNRAS, 474, 197 [Google Scholar]
  140. Pietrzyński, G., Graczyk, D., Gieren, W., et al. 2013, Nature, 495, 76 [Google Scholar]
  141. Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  142. Portegies Zwart, S. F., & van den Heuvel, E. P. J. 2016, MNRAS, 456, 3401 [NASA ADS] [CrossRef] [Google Scholar]
  143. Prinja, R. K., & Schild, H. 1991, A&A, 249, 428 [NASA ADS] [Google Scholar]
  144. Reguitti, A., Pastorello, A., Pignata, G., et al. 2019, MNRAS, 482, 2750 [NASA ADS] [Google Scholar]
  145. Richardson, N. D., & Mehner, A. 2018, Res. Notes Am. Astron. Soc., 2, 121 [Google Scholar]
  146. Riebel, D., Srinivasan, S., Sargent, B., & Meixner, M. 2012, ApJ, 753, 71 [NASA ADS] [CrossRef] [Google Scholar]
  147. Riechers, D. A., Leung, T. K. D., Ivison, R. J., et al. 2017, ApJ, 850, 1 [Google Scholar]
  148. Rieke, G. H., Young, E. T., Engelbracht, C. W., et al. 2004, ApJS, 154, 25 [Google Scholar]
  149. Ruffle, P. M. E., Kemper, F., Jones, O. C., et al. 2015, MNRAS, 451, 3504 [NASA ADS] [CrossRef] [Google Scholar]
  150. Schmutz, W., Leitherer, C., Hubeny, I., et al. 1991, ApJ, 372, 664 [NASA ADS] [CrossRef] [Google Scholar]
  151. Schneider, F. R. N., Sana, H., Evans, C. J., et al. 2018, Science, 359, 69 [NASA ADS] [CrossRef] [Google Scholar]
  152. Silvia, D. W., Smith, B. D., & Shull, J. M. 2010, ApJ, 715, 1575 [NASA ADS] [CrossRef] [Google Scholar]
  153. Silvia, D. W., Smith, B. D., & Shull, J. M. 2012, ApJ, 748, 12 [NASA ADS] [CrossRef] [Google Scholar]
  154. Smith, N. 2007, AJ, 133, 1034 [NASA ADS] [CrossRef] [Google Scholar]
  155. Smith, N. 2010, MNRAS, 402, 145 [NASA ADS] [CrossRef] [Google Scholar]
  156. Smith, N. 2014, ARA&A, 52, 487 [NASA ADS] [CrossRef] [Google Scholar]
  157. Smith, N., Li, W., Foley, R. J., et al. 2007, ApJ, 666, 1116 [NASA ADS] [CrossRef] [Google Scholar]
  158. Smith, N., Foley, R. J., & Filippenko, A. V. 2008a, ApJ, 680, 568 [Google Scholar]
  159. Smith, N., Chornock, R., Li, W., et al. 2008b, ApJ, 686, 467 [NASA ADS] [CrossRef] [Google Scholar]
  160. Sodroski, T. J., Dwek, E., Hauser, M. G., & Kerr, F. J. 1987, ApJ, 322, 101 [NASA ADS] [CrossRef] [Google Scholar]
  161. Srinivasan, S., Meixner, M., Leitherer, C., et al. 2009, AJ, 137, 4810 [NASA ADS] [CrossRef] [Google Scholar]
  162. Srinivasan, S., Boyer, M. L., Kemper, F., et al. 2016, MNRAS, 457, 2814 [NASA ADS] [CrossRef] [Google Scholar]
  163. Stahl, O., & Wolf, B. 1987, A&A, 181, 293 [NASA ADS] [Google Scholar]
  164. Stahl, O., Wolf, B., Zickgraf, F. J., et al. 1983, A&A, 120, 287 [NASA ADS] [Google Scholar]
  165. Stahl, O., Wolf, B., Leitherer, C., et al. 1984, A&A, 140, 459 [NASA ADS] [Google Scholar]
  166. Stahl, O., Wolf, B., & Zickgraf, F. J. 1987, A&A, 184, 193 [NASA ADS] [Google Scholar]
  167. Stahl, O., Wolf, B., Klare, G., Juettner, A., & Cassatella, A. 1990, A&A, 228, 379 [NASA ADS] [Google Scholar]
  168. Sugerman, B. E. K., Ercolano, B., Barlow, M. J., et al. 2006, Science, 313, 196 [Google Scholar]
  169. Sundqvist, J. O., Björklund, R., Puls, J., & Najarro, F. 2019, A&A, 632, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  170. Szeifert, T., Stahl, O., Wolf, B., et al. 1993, A&A, 280, 508 [NASA ADS] [Google Scholar]
  171. Tamura, Y., Mawatari, K., Hashimoto, T., et al. 2019, ApJ, 874, 27 [NASA ADS] [CrossRef] [Google Scholar]
  172. Tartaglia, L., Elias-Rosa, N., Pastorello, A., et al. 2016, ApJ, 823, L23 [NASA ADS] [CrossRef] [Google Scholar]
  173. Temim, T., Dwek, E., Tchernyshyov, K., et al. 2015, ApJ, 799, 158 [NASA ADS] [CrossRef] [Google Scholar]
  174. Todini, P., & Ferrara, A. 2001, MNRAS, 325, 726 [NASA ADS] [CrossRef] [Google Scholar]
  175. Trundle, C., Kotak, R., Vink, J. S., & Meikle, W. P. S. 2008, A&A, 483, L47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  176. Tubbesing, S., Kaufer, A., Stahl, O., et al. 2002, A&A, 389, 931 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  177. Umana, G., Buemi, C. S., Trigilio, C., et al. 2009, ApJ, 694, 697 [NASA ADS] [CrossRef] [Google Scholar]
  178. Umana, G., Buemi, C. S., Trigilio, C., Leto, P., & Hora, J. L. 2010, ApJ, 718, 1036 [NASA ADS] [CrossRef] [Google Scholar]
  179. Umana, G., Ingallinera, A., Trigilio, C., et al. 2012, MNRAS, 427, 2975 [NASA ADS] [CrossRef] [Google Scholar]
  180. Urbaneja, M. A., Kudritzki, R. P., Gieren, W., et al. 2017, AJ, 154, 102 [NASA ADS] [CrossRef] [Google Scholar]
  181. Vamvatira-Nakou, C., Hutsemékers, D., Royer, P., et al. 2013, A&A, 557, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  182. Vamvatira-Nakou, C., Hutsemékers, D., Royer, P., et al. 2015, A&A, 578, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  183. van Aarle, E., van Winckel, H., Lloyd Evans, T., et al. 2011, A&A, 530, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  184. van der Marel, R. P., & Cioni, M.-R. L. 2001, AJ, 122, 1807 [Google Scholar]
  185. Van Dyk, S. D., & Matheson, T. 2012, ApJ, 746, 179 [NASA ADS] [CrossRef] [Google Scholar]
  186. van Gelder, M. L., Kaper, L., Japelj, J., et al. 2020, A&A, 636, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  187. van Genderen, A. M. 2001, A&A, 366, 508 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  188. van Genderen, A. M., & Sterken, C. 1996, A&A, 308, 763 [NASA ADS] [Google Scholar]
  189. van Genderen, A. M., & Sterken, C. 1999, A&A, 349, 537 [NASA ADS] [Google Scholar]
  190. van Genderen, A. M., Sterken, C., & de Groot, M. 1998, A&A, 337, 393 [NASA ADS] [Google Scholar]
  191. van Loon, J. T., Oliveira, J. M., Gordon, K. D., et al. 2010, AJ, 139, 68 [NASA ADS] [CrossRef] [Google Scholar]
  192. Vink, J. S., & de Koter, A. 2002, A&A, 393, 543 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  193. Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 1999, A&A, 350, 181 [Google Scholar]
  194. Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nat. Meth., 17, 261 [Google Scholar]
  195. Voors, R. H. M., Waters, L. B. F. M., Morris, P. W., et al. 1999, A&A, 341, L67 [NASA ADS] [Google Scholar]
  196. Walborn, N. R., Stahl, O., Gamen, R. C., et al. 2008, ApJ, 683, L33 [NASA ADS] [CrossRef] [Google Scholar]
  197. Walborn, N. R., Barbá, R. H., & Sewiło, M. M. 2013, AJ, 145, 98 [Google Scholar]
  198. Walborn, N. R., Gamen, R. C., Morrell, N. I., et al. 2017, AJ, 154, 15 [NASA ADS] [CrossRef] [Google Scholar]
  199. Weis, K. 2003, A&A, 408, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  200. Weis, K., & Duschl, W. J. 2002, A&A, 393, 503 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  201. Weis, K., Chu, Y. H., Duschl, W. J., & Bomans, D. J. 1997, A&A, 325, 1157 [NASA ADS] [Google Scholar]
  202. Weis, K., Duschl, W. J., & Bomans, D. J. 2003, A&A, 398, 1041 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  203. Wesson, R., Barlow, M. J., Matsuura, M., & Ercolano, B. 2015, MNRAS, 446, 2089 [NASA ADS] [CrossRef] [Google Scholar]
  204. Whitney, B. A., Sewilo, M., Indebetouw, R., et al. 2008, AJ, 136, 18 [NASA ADS] [CrossRef] [Google Scholar]
  205. Wolf, B. 1989, A&A, 217, 87 [NASA ADS] [Google Scholar]
  206. Wolf, B., Stahl, O., de Groot, M. J. H., & Sterken, C. 1981, A&A, 99, 351 [NASA ADS] [Google Scholar]
  207. Wolf, B., Stahl, O., & Seifert, W. 1987, A&A, 186, 182 [NASA ADS] [Google Scholar]
  208. Zickgraf, F. J., Wolf, B., Stahl, O., Leitherer, C., & Appenzeller, I. 1986, A&A, 163, 119 [NASA ADS] [Google Scholar]
  209. Zickgraf, F. J., Kovacs, J., Wolf, B., et al. 1996a, A&A, 309, 505 [Google Scholar]
  210. Zickgraf, F. J., Humphreys, R. M., Lamers, H. J. G. L. M., et al. 1996b, A&A, 315, 510 [Google Scholar]

Appendix A: Table of stellar and nebular parameters

We use stellar parameters collected from the literature. Only values corresponding to the quiescent state of the stars were taken into account. To compensate for the lack of error estimates, we take the mean value and the corresponding standard deviation when more than one estimate is available. In this case all the relevant works are mentioned. We also include information on the nebula size when known from the literature and from this work.

Table A.1.

Stellar and nebular parameters.

Appendix B: IR images and SEDs

Flux density distribution plots of individual sources and selected space telescope images are provided here. The photometric data are corrected for interstellar dust extinction, using the extinction curve from Fitzpatrick (1999), E(B-V) as in Table A.1 and RV = 3.1. For RMC 143 an RV = 4.0 was adopted as previously determined by Agliozzo et al. (2019). The images are useful for assessing the confusion and thus the reliability of the photometry in each band. The field-of-view in each image is 240″ × 240″. The position of the source is highlighted with a cyan cross-hair, which is rotated accordingly with the sky projection. The colour-scale is squared-root in order to highlight the local environment around each source.

thumbnail Fig. B.1.

Confirmed and candidate LBVs with a SED characterised by a free-free continuum from their stellar wind and by a thermal continuum due to dust in their circumstellar ejecta (see Class 1a sources in Sect. 4). The legend for the photometric points can be found in Fig. 5 of Sect. 4.

thumbnail Fig. B.1.

continued.

thumbnail Fig. B.2.

Confirmed and candidate LBVs with a SED characterised by a thermal continuum due to dust in their circumstellar ejecta and relatively weak or absent free-free continuum from their stellar wind (see Class 1b sources in Sect. 4). The legend for the photometric points can be found in Fig. 5 of Sect. 4.

thumbnail Fig. B.2.

continued.

thumbnail Fig. B.3.

Confirmed and candidate LBVs with a SED dominated by free-free continuum up to 24 μm and by confusion at longer wavelengths (Class 2 sources in Sect. 4). For HD 269582 the jump between 2MASS and other photometry is due to real source variability. The legend for the photometric points can be found in Fig. 5 of Sect. 4.

thumbnail Fig. B.3.

continued. We note that photometry for RMC 14 has not been extinction corrected due to no AV value in the literature.

thumbnail Fig. B.4.

Unclassified LBV RMC 85. No reliable photometry could be obtained longward of 20 μm due to confusion with bright ISM dust. The legend for the photometric points can be found in Fig. 5 of Sect. 4.

thumbnail Fig. B.5.

sgB[e] stars as in literature (see Class 3 sources in Sect. 4). The legend for the photometric points can be found in Fig. 5 of Sect. 4.

thumbnail Fig. B.5.

continued. We note that photometry for LHA 115-S 18 has not been extinction corrected due to no AV value in the literature.

thumbnail Fig. B.6.

Featureless sources with no recognised LBV membership (see Class 4 sources in Sect. 4). We note that photometry for Sk −69 271 has not been extinction corrected due to no AV value in the literature. The legend for the photometric points can be found in Fig. 5 of Sect. 4.

Appendix C: Tables of infrared photometry

In this appendix the adopted photometry for all sources from the near-infrared to far-infrared are provided. In the following tables all photometry is provided as frequency flux density in units of mJy. Where catalogue values were provided as magnitudes, these have been converted to flux densities using zero point (zero magnitude) fluxes which are indicated in the table header.

All 2MASS photometry comes from the 2MASS Point Source Catalogue (Cutri et al. 2003), and were converted to flux densities using the zero point fluxes from Sect. 6.4a of the corresponding Explanatory Supplement7. AKARI N3, S7, S11, L15 and L24 photometry come from the AKARI-LMC Point-source catalog (Kato et al. 2012), and were converted to flux densities using the zero point fluxes from Kato et al. (2012) (which, although of low precision, are the logical choice to use, as they explain that the photometry was first calibrated into units of Jy and then converted to the magnitudes in the catalogue using these zero point fluxes). AKARI S9W and L18W photometry come from the AKARI/IRC mid-infrared all-sky Survey (Ishihara et al. 2010), originally in Jy units. WISE photometry almost exclusively comes from The AllWISE data Release (updated version, 16-Feb-2021) (Cutri et al. 2021). WISE W4 values for RMC 116 (ID 6) and Sk −69 271 (ID 24) come from The WISE All-Sky data Release (Cutri et al. 2012) as they are not detected in ALLWISE (both are faint and hard to distinguish from the background emission). WISE magnitudes were converted to flux densities using zero point fluxes for a Fν ∝ ν−2 spectrum from Sect. IV.4.i.1 of the Explanatory Supplement to the WISE All-Sky Data Release Products8. All MSX photometry comes from the MSX Point Source Catalog Version 2.3 (Egan et al. 2003), originally in Jy units.

The Spitzer and Herschel photometry are discussed in the main text, with some values newly extracted in this work.

For completeness, catalogue values that we have identified as incorrect are included, with their value suffixed by a flag character. A ‘c’ denotes that the value is most likely affected by confusion with nearby bright emission, and a ‘r’ denotes that the value is likely affected by the source being resolved. These were assessed with the aid of the flux density distributions and images in Appendix B.

Table C.1.

Photometry at infrared wavelengths shorter than 5 μm.

Table C.2.

Photometry at wavelengths 5−20 μm.

Table C.3.

Photometry at wavelengths longer than 20 μm.

All Tables

Table 1.

Sources, celestial coordinates (Gaia Collaboration 2018), LBV status (previous: Richardson & Mehner 2018; new: this work), and SED-based classification proposed in this work.

Table 2.

Journal of VISIR observations.

Table 3.

MIPS-24 aperture photometry.

Table 4.

MIPS-70 aperture photometry and 3σ upper limits.

Table 5.

PACS-100 point source extraction and 3σ upper limits.

Table 6.

Summary of estimated dust properties.

Table 7.

Flux densities from the stacked data.

Table 8.

Summary of estimated LBV dust masses in the LMC.

Table A.1.

Stellar and nebular parameters.

Table C.1.

Photometry at infrared wavelengths shorter than 5 μm.

Table C.2.

Photometry at wavelengths 5−20 μm.

Table C.3.

Photometry at wavelengths longer than 20 μm.

All Figures

thumbnail Fig. 1.

Cutout of the SAGE MIPS24 image (Meixner et al. 2006) of the LMC and positions of the candidate and confirmed LBVs analysed in this work. The numbers are the IDs in Table 1, and the red circles have a radius of 30″. North is towards the top right.

In the text
thumbnail Fig. 2.

Cutout of the SAGE MIPS24 image (Meixner et al. 2006) of the SMC and positions of the candidate and confirmed LBVs analysed in this work. The numbers are the IDs in Table 1, and the red circles have a radius of 30″. North is up and east is left.

In the text
thumbnail Fig. 3.

VISIR images of RMC 71 (top panels), and the standard calibrators observed on the same nights (bottom panels, with identifications given in the lower right corners). North is up and east is left. The colour-scale is in arbitrary units.

In the text
thumbnail Fig. 4.

Examples of each of the four LBV SED classes defined in Sect. 4.

In the text
thumbnail Fig. 5.

Mid-infrared colour-magnitude diagram for stars in the LMC. The LBVs are shown with filled circles in different colours, depending on the SED class defined in this work. Several types of supergiant stars (RSGs, sgB[e], supergiants A and B, six LBVs and other massive stars) as in Bonanos et al. (2009) are displayed with triangles of different colours as identified in the legend. Dusty evolved stars from Jones et al. (2015) are shown with diamond symbols (these include post-AGBs, planetary nebulae, sgB[e] stars, one LBV, one WR and two SN remnants). Pixel symbols indicate the AGBs and RSGs from Riebel et al. (2012) and point symbols are confirmed and candidate YSOs from Whitney et al. (2008).

In the text
thumbnail Fig. 6.

Example of single grey-body fit of the 2MASS, WISE, Spitzer, and Herschel data, similarly to Agliozzo et al. (2017a,b, 2019).

In the text
thumbnail Fig. 7.

Left column: stacked images of LMC LBVs (SED Classes 1 and 2 in Table 1) and residuals after Gaussian fitting. Right column: same as before but without RMC143. The field of view corresponds approximately to 70″ × 70″. The stack image at 160 μm on the bottom right was derived with a smaller number of sources (IDs: 2, 5, 6, 7, 8, 10, 12, 13, 14, 15, 16, 17, 18) and used as a lower limit.

In the text
thumbnail Fig. 8.

Flux density distribution of stacks of the LBVs in the LMC (SED Classes 1 and 2 in Table 1), with (top) and without (bottom) RMC 143. 2MASS and IRAC values are the sum of the individual source fluxes, while WISE-W3, MIPS and PACS values are from photometry performed on the stacked images (Fig. 7). Parameters from the fitting procedure are shown (α and c are the spectral index and offset of the best-fit power-law describing the near-infrared data; Td, Md and β are the dust temperature, mass and opacity power law index of the best-fit greybody describing the thermal emission from dust). The lower limits (red points and arrow symbols) for MIPS70 and SPIRE250 in the upper plot and PACS160 in the lower plot are from stacks of subsets of sources (see text for the subsets in each case), as there is not a reliable detection in the stacked images of the whole sample.

In the text
thumbnail Fig. 9.

Stacked images of the extreme-AGBs from Riebel et al. (2012) and residuals after Gaussian source extraction. The field of view corresponds approximately to 70″ × 70″. Ten out of 1352 sources were excluded (see text), because the emission in their field is dominated by bright large-scale structures or because they are another type of stellar objects (e.g. YSOs).

In the text
thumbnail Fig. 10.

Flux density distribution of the 1342 extreme x-AGBs in the LMC from the catalogue by Riebel et al. (2012), constructed from the Gaussian fitting of the central detection in the stacked images (at wavelengths shorter than ∼10 μm the value is just the sum of the catalogue measurements). Blue points: observed phtometry; green points: de-reddened photometry; red arrow: 3σ PACS160 upper limit.

In the text
thumbnail Fig. 11.

Estimated dust masses for individual stars as a function of the stellar bolometric luminosity. Red circles: LMC LBVs (this work, see Table 6). Blue circles: Galactic LBVs, specifically AG Car (Vamvatira-Nakou et al. 2015); G24.73+0.69, G25.5+0.2, HD 168625, Hen 3–519, AFGL 2298 (Clark et al. 2003, and references therein); G26.47+0.02 (Umana et al. 2012); G79.49+0.26 (Agliozzo et al. 2014); HR Car (Machado et al. 2002; Buemi et al. 2017); Pistol star (Lau et al. 2014); Wra 751 (Vamvatira-Nakou et al. 2013). Error bars are shown when available from the literature.

In the text
thumbnail Fig. B.1.

Confirmed and candidate LBVs with a SED characterised by a free-free continuum from their stellar wind and by a thermal continuum due to dust in their circumstellar ejecta (see Class 1a sources in Sect. 4). The legend for the photometric points can be found in Fig. 5 of Sect. 4.

In the text
thumbnail Fig. B.2.

Confirmed and candidate LBVs with a SED characterised by a thermal continuum due to dust in their circumstellar ejecta and relatively weak or absent free-free continuum from their stellar wind (see Class 1b sources in Sect. 4). The legend for the photometric points can be found in Fig. 5 of Sect. 4.

In the text
thumbnail Fig. B.3.

Confirmed and candidate LBVs with a SED dominated by free-free continuum up to 24 μm and by confusion at longer wavelengths (Class 2 sources in Sect. 4). For HD 269582 the jump between 2MASS and other photometry is due to real source variability. The legend for the photometric points can be found in Fig. 5 of Sect. 4.

In the text
thumbnail Fig. B.3.

continued. We note that photometry for RMC 14 has not been extinction corrected due to no AV value in the literature.

In the text
thumbnail Fig. B.4.

Unclassified LBV RMC 85. No reliable photometry could be obtained longward of 20 μm due to confusion with bright ISM dust. The legend for the photometric points can be found in Fig. 5 of Sect. 4.

In the text
thumbnail Fig. B.5.

sgB[e] stars as in literature (see Class 3 sources in Sect. 4). The legend for the photometric points can be found in Fig. 5 of Sect. 4.

In the text
thumbnail Fig. B.5.

continued. We note that photometry for LHA 115-S 18 has not been extinction corrected due to no AV value in the literature.

In the text
thumbnail Fig. B.6.

Featureless sources with no recognised LBV membership (see Class 4 sources in Sect. 4). We note that photometry for Sk −69 271 has not been extinction corrected due to no AV value in the literature. The legend for the photometric points can be found in Fig. 5 of Sect. 4.

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.