Broad Balmer Wings in BA Hyper/Supergiants Distorted by Diffuse Interstellar Bands: Five Examples in the 30 Doradus Region from the VLT-FLAMES Tarantula Survey
Extremely broad emission wings at H and H have been found in VFTS data for five very luminous BA supergiants in or near 30 Doradus in the Large Magellanic Cloud. The profiles of both lines are extremely asymmetrical, which we have found to be caused by very broad diffuse interstellar bands (DIBs) in the longward wing of H and the shortward wing of H. These DIBs are well known to interstellar but not to many stellar specialists, so that the asymmetries may be mistaken for intrinsic features. The broad emission wings are generally ascribed to electron scattering, although we note difficulties for that interpretation in some objects. Such profiles are known in some Galactic hyper/supergiants and are also seen in both active and quiescent Luminous Blue Variables. No prior or current LBV activity is known in these 30 Dor stars, although a generic relationship to LBVs is not excluded; subject to further observational and theoretical investigation, it is possible that these very luminous supergiants are approaching the LBV stage for the first time. Their locations in the HRD and presumed evolutionary tracks are consistent with that possibility. The available evidence for spectroscopic variations of these objects is reviewed, while recent photometric monitoring does not reveal variability. A search for circumstellar nebulae has been conducted, with an indeterminate result for one of them.
A property of Luminous Blue Variables both in outburst and quiescense is the occurrence of extremely broad Balmer emission wings, usually ascribed to electron scattering in dense extended atmospheres or envelopes (Bernat & Lambert 1978; Hubeny & Leitherer 1989; Santolaya-Rey et al. 1997; Najarro & Hillier 2012). In some cases, a narrow P Cygni profile is superimposed, producing the “Prussian Helmet” morphology. Examples may be found in Stahl et al. (1985), Hutsemékers & Van Drom (1991), and Walborn & Fitzpatrick (2000). Similar profiles are observed in some very luminous Galactic B supergiants with no evidence of previous LBV activity (e.g., Lennon et al. 1992; Marco & Negueruela 2009; Clark et al. 2012), although it would not be surprising if these objects were approaching the LBV stage. A puzzling property of the H wings is their apparent asymmetry, more extended and shallow shortward but curtailed and steeper longward. As emphasized by Hutsemékers & Van Drom, that morphology cannot be produced by electron scattering; moreover, it is not shared by H. Here we propose a likely resolution of this discrepancy.
We report the discovery of these broad Balmer wings in five BA hyper/supergiants in or near the 30 Doradus Nebula of the Large Magellanic Cloud (LMC), but again, with no record or known evidence of outbursts. The only confirmed LBV in the 30 Dor region is Radcliffe (R) 143 (Feast et al. 1960), located in an evolved association immediately SE of the main Tarantula Nebula (Walborn & Blades 1997). Feast et al. reported an F-supergiant spectrum, but Parker et al. (1993) discovered it two magnitudes fainter with a late-B spectrum (see also Walborn 1997), establishing its LBV nature. Smith et al. (1998; see also Weis 2003) subsequently found ejected nitrogen-rich circumstellar nebulosity, another frequent characteristic of LBVs and indicative of a larger, earlier eruption. The spectra, locations, and known variability of these five 30 Dor objects, one of which may be associated with R143, are discussed below, as well as a search for circumstellar nebulae.
The primary spectroscopic data used here are from the VLT-FLAMES Tarantula Survey (VFTS; ESO Large Programme 182.D-0222, PI CJE). Full details of the instrumental parameters and data reductions were provided by Evans et al. (2011, 2015 and references therein). In brief, the observations were obtained in the Medusa–Giraffe configuration of the Fibre Large Array Multi-Element Spectrograph (FLAMES; Pasquini et al. 2002) at the Very Large Telescope (VLT) on Cerro Paranal, Chile. Each target was observed with the standard LR02 and LR03 settings of the Giraffe spectrograph, providing coverage of 3960–5071 Å at a spectral resolving power (R) of 7000–8500. At least six observations were obtained with the LR02 setting to support the investigation of spectroscopic binaries (Sana et al. 2013; Evans et al. 2015). In addition, the H region was observed with the HR15N setting at R of 16,000. Most of the observations were done during 2008 October through 2009 February, with a final epoch in 2009 October to extend the binary period sensitivity.
The locations of the five VFTS objects in the field of 30 Dor are shown in Figure 1, along with that of R143. Coordinates and photometry are provided by Evans et al. (2011). Their full blue-violet spectra are displayed in Figure 2, while enlargements of their H and H profiles are shown in Figure 3. Several interesting common properties of these Balmer profiles are discussed here, while further details of the individual objects and spectra are considered in following subsections.
The similar, marked asymmetry of the broad H wings in all the objects is remarkable and characteristic: the shortward wings are more extended and shallower, while the longward ones are shorter and steeper (Figure 3). The broad emission has a full width at zero intensity of 50 Å or 3000 km s, far too large to be ascribed to a Doppler effect, which is why such profiles have been interpreted in terms of electron scattering. (Typical wind terminal velocities of Galactic B1–2 hypergiants are a few hundred km s; Crowther et al. 2006, Clark et al. 2012.) In luminous late-B and A supergiants with lower mass-loss rates, peculiar Balmer emission profiles can be produced by NLTE effects (Hubeny & Leitherer 1989; Santolaya-Rey et al. 1997; Puls et al. 1998). The observed peak intensities of the broad emission wings are 1–3% relative to the continuum.
The H wings are entirely different from those of H, with a broad absorption depression in the shortward one. Our first thought on the latter was a feature related to the stellar winds, although its similarity across the entire B-type spectral range (Figure 2) was disconcerting, as was its extent relative to the stellar-wind velocities of such objects. Subsequently, we realized that the depths of these features are correlated with the reddenings of the stars, immediately suggesting an origin related to extinction and consultation of Herbig (1995), who indeed lists a broad, shallow Diffuse Interstellar Band (DIB) centered at 6533 Å, with a FWHM of 21 Å, in good agreement with the data. An analogous origin of the H asymmetry was then suggested by the extreme case of Cygnus OB2-12, as shown in Figure 4. Here the DIB responsible is centered at 4882 Å with FWHM 25 Å according to Herbig, but Maíz Apellániz et al. (2014) have decomposed it into two features, a narrower one centered at 4880 and a broader one at 4887 Å.
3.1 Vfts 003 = Hd 38029b
This luminous B supergiant is located only 28 from and is most likely associated in a small cluster with the WC4 + O6-6.5 III system HD 38029A = VFTS 002 (Evans et al. 2011), which is currently much fainter in the optical but likely also descends from a very massive progenitor. To our knowledge, VFTS 003 lacks accurate, spatially resolved optical photometry; in fact, there is confusion in the literature about the magnitudes of the two stars. Remarkably, Sanduleak (1970) resolved them and correctly gave their relative positions and magnitudes; in particular, he lists m of 11.8 for the B supergiant and 14.2 for the WC system. However, Breysacher et al. (1999) quote a magnitude of 11.75 for the “WC star,” which clearly corresponds primarily to the B supergiant instead. The same remark applies to the magnitudes listed by SIMBAD, except for . In contrast, the issue was correctly resolved by Prevot-Burnichon et al. (1981). In any event, this association should eventually yield valuable mutual clues to the basic parameters of both evolved objects. The VFTS 002/003 system is located at the edge of the survey field, far SW of the Tarantula Nebula, and is thus not directly associated with it (Fig. 1).
The H profile (Fig. 3) shows the classic “Prussian Helmet” morphology, with a narrow P Cygni profile superimposed on the broad emission wings. Since this star is not near strong H ii emission, a significant nebular contribution is not expected.
The spectral type of B1 Iap is the earliest among the present sample; it is based upon the Si iv/Si iii and Si iii/He i absorption-line ratios (Fig. 2) for the temperature and luminosity classes, respectively, allowing for LMC metal deficiency in the latter (Fitzpatrick 1991). The “p” is due to the H profile. The N ii spectrum is very strong, morphologically normal for the spectral type, which is expected to correspond to a significant physical enhancement. Indeed, McEvoy et al. (2015) have derived log(N/H)+12 (hereafter N/H) = 8.0 for VFTS 003; for comparison, Kurt & Dufour (1998) give an LMC baseline N/H of 6.9 determined from H ii regions.
3.2 Vfts 533 = R142
The H profile again appears to show the classic “Prussian Helmet” morphology, although in this case there is a substantial, blended nebular contribution to the narrow emission. Unlike in VFTS 003, the narrow H structure is very similar to that of H. For reference in the context of the DIBs, .
The spectral type is B1.5 Iap Nwk. Remarkably for a very high-luminosity object possibly nearing the LBV phase, the nitrogen lines are weaker than those of oxygen and carbon; while the N ii lines are well marked, they are too weak for the spectral type, and the O ii spectrum is much stronger. Compared with VFTS 003 discussed above, the slightly later type of VFTS 533 places it nearer to the normal N ii maximum at B2, exacerbating the discrepancy. Fitzpatrick (1991) found the same anomalous nitrogen deficiency in the LMC B1–2 hypergiants Sanduleak (1970, Sk) , , and , as well as in several slightly less luminous supergiants. Another striking comparison of N-weak and N-strong objects at this type from the VFTS sample is shown by Evans et al. (2015, their Fig. 1). McEvoy et al. (2015) have derived N/H of 7.4 for VFTS 533, consistent with the morphological appearance of the spectrum. With respect to the LMC baseline abundance cited above, the VFTS 533 value is certainly enhanced, but far less than that of other hyper/supergiants, e.g., VFTS 003.
Parker (1993) derived a spectral type of B0 Ia, and Walborn & Blades (1997; data from 1982) B0.5-B0.7 I; while the quality of the present data is far superior, nevertheless both of those classifications imply Si iv line strengths much greater than seen here. Thus a real spectral variation during an interval of 26 years appears likely. Moreover, the depth of the absorption component in the narrow P Cyg profile at H increased by a factor of two to 20% of the continuum between VFTS epochs of 2008 October and December. During the same interval the metallic lines underwent a radial-velocity decrease of 12 km s and a small intensity decrease, perhaps consistent with atmospheric motions as in the Galactic supergiant Sher 25 (Taylor et al. 2014).
As shown in Fig. 1, this star is located just N(E) of R136 at the center of the nebula. However, it lies among the older stars identified as a separate cluster by Sabbi et al. (2012). Indeed, this spectral type is unlikely to be associated with the 1–2 Myr old, extremely early O and supermassive WNh stars in R136 (Crowther et al. 2010; Crowther et al. in prep.) Nevertheless, it is an extremely luminous star: with , it is one magnitude brighter than the typical Ia (Walborn 1972). With the bolometric corrections listed by Walborn et al. (2008, Table 2), the corresponding is 10.
3.3 Vfts 424 = R138
In this case, the narrow H emission appears to be nebular, centered within a stellar absorption feature that in turn reverses the broad emission. It is probably not a coincidence that the spectral type is much later than those of VFTS 003 and 533, i.e., B9.5 Iap with very strong Si ii, Mg ii, and Fe ii absorption lines. The presence of N ii 3995 at such a late type, albeit weak, may indicate an overabundance. The H emission profile is also a stellar/nebular composite; in this case .
VFTS 424 is also very near R136, although on the opposite (N)W side from VFTS 533. Still, it is most likely associated with the older population related to the Sabbi et al. (2012) cluster; see also Walborn & Blades (1997) and Selman et al. (1999) for more on the distribution of this population. VFTS 424 is also highly luminous, again with , albeit with only 8.4 because of the lower temperature.
3.4 Vfts 458 = P662
This star was suggested as an LBV candidate by Walborn & Blades (1997), albeit for different reasons; their data did not cover H. Although the S/N was low, N ii 3995 appeared to be present with a strength greater than that of the adjacent He i 4009, an anomaly for the spectral type of B6 Iap. There also appeared to be a broad absorption consistent with an Fe ii complex at 4550–4650 Å. Ironically, neither of those features has the same appearance in the superior VFTS data (Fig. 2), although 3995 and other N ii lines are still too strong for the spectral type. The apparent broad feature is resolved into N ii and Fe ii lines in the VFTS data. McEvoy et al. derive N/H = 8.0. The triple H structure is very similar to that of VFTS 424 and likely has similar origins. The H emission is again a stellar/nebular composite, and is the largest in this sample consistent with the strongest 6533 DIB.
VFTS 458 = Parker (1993) 662 is located about 2 (30 pc in projection) south of R136, near the previously known LBV R143 and within the “R143 Association” first distinguished and discussed by Walborn & Blades (1997). It also has (as does R143 itself in its current state), and of 8.6.
3.5 VFTS 739 = Sk
The H profile of this star consists of only two components: the very broad emission wings and an approximately central stellar absorption reversal (Fig. 3). Note that the lack of a central emission feature definitely eliminates electron scattering as the origin of the broad wings in this profile (cf. Hutsemékers & Van Drom 1991). Remarkably, its H profile displays a second set of stronger, Be-like emission wings with a peak-to-peak separation of only 3.3 Å or 150 km s. Such a profile in a supergiant does not necessarily imply the presence of a disk, as it can be produced by NLTE radiative transfer effects, as can the broad wings (Hubeny & Leitherer 1989; Puls et al. 1998). As further discussed in the next paragraph, this star is far from any nebulosity, thus supporting the interpretation of the profiles in the previous two spectra. The spectral type of A0 Iap is the latest of the five objects discussed here, while .
VFTS 739 does not belong to 30 Dor per se but is located even further to the south, at the northern edge of a large region of other discrete OB associations and scattered “field” stars. The comparably bright star to its immediate SE in Fig. 1 is VFTS 764 = Sk of type O9.7 Ia Nstr (Walborn et al. 2014), also indicative of a fairly advanced evolutionary stage. McEvoy et al. derive N/H = 7.8 for VFTS 764. The of VFTS 739 is 7.2 ( 7.3), while that of VFTS 764 is 7.4 ( 10.2).
In addition to spectroscopic characteristics, photometric histories and circumstellar nebulae may provide information about the status of candidate or precursor LBVs. Here we briefly describe the results of some inquiries in those areas. Some Galactic counterparts to these VFTS objects and LBVs with similar Balmer profiles are also listed.
4.1 Photometric Variability
We investigated our targets in the ASAS and the OGLE shallow databases. The ASAS data show apparent variability of VFTS 458 at the 1 mag level and of VFTS 739 at 0.3 mag. However, we are concerned about the bright nebulosity associated with the former (next subsection) and the latter range is near the accuracy limit. Indeed, the more precise OGLE data exclude variations at those levels for both stars; and data for VFTS 003, 424, 458, and 739 between HJD 2,453,000 and 2,455,000 show only a few points deviating by the order of 0.1 mag. (VFTS 533 was not reported.) Thus, we do not find LBV-like variability in light among the available data for these stars.
4.2 Circumstellar Nebulae
It is well known that the relatively rare giant eruptions of LBVs (as opposed to their more modest “outbursts” or photospheric expansions) eject nitrogen- and sometimes dust-rich circumstellar nebulae that may persist for thousands of years (Weis 2003; Walborn et al. 2008 and references therein). Hence, any such nebulae associated with our objects would be relevant to their evolutionary status. Accordingly, we have investigated available images from the Hubble Space Telescope cameras.
Sabbi et al. (2013) have described a massive Treasury Program in 30 Doradus using the WFC3 and ACS cameras, including NIR observations with the former and H with the latter. The H bandpass is 72 Å, which includes the adjacent [N ii] lines. The optical image clearly shows the structured, asymmetrical ejected nebulosity associated with R143, discovered by Smith et al. (1998) and further discussed by Weis (2003).
The only one of our objects with a possibly associated nebulosity is
VFTS 458, which appears centered in a relatively isolated, complex, asymmetrical (bipolar?) structure in both the optical and -band images (Figure 5). We have further investigated this structure with two new optical spectroscopic datasets: three longslit position angles from the LCO 2.5m Boller & Chivens cassegrain spectrograph (R ) passing through the star and bright nebular knots; and three pointings with the ESO La Silla 2.2m FEROS echelle spectrograph (R = 48,000) both on the star and a few arcsec N or S of it. The principal result is even more extreme complexity of the region, with multiple nebular line components varying in both velocity and relative intensity on arcsec scales, which these data are inadequate to fully characterize; further observations (preferably with IFUs) will be required to do so. Most importantly, we have not encountered any point with an [N ii] 6584/H intensity ratio greater than 10%, while ratios in the range 0.5 to 1 are typical of LBV and related circumstellar nebulae (references in the previous paragraph). Of course, the R143 nebula provides a cautionary note: most of the nebulosity projected near the star is ambient Hii, but there are also compact knots of ejected processed material present. We note that the Spitzer/SAGE project (Meixner et al. 2006) found very bright  magnitudes of 2.4–2.8 for VFTS 458;
4.3 Galactic Counterparts
Several Galactic hypergiants or supergiants with the same peculiar H profiles found in these LMC objects are known. To promote future comparative analyses, some details are provided here for several published cases. We are also aware of a number of unpublished ones that will appear in due course.
HD 190603, B1.5 Ia, has a Prussian Helmet profile, while HD 4841, B5 Ia, has the broad, asymmetrical wings with no central emission (Lennon et al. 1992, Figs. 8 and 16, respectively). This atlas also displays many spectra throughout the B Ia sequence without such H profiles.
Scorpii (HD 152236), B1.5 Ia, has a Prussian Helmet as well, which may be variable (Clark et al. 2012, Fig A.1). This spectrum is very similar to that of VFTS 533, albeit with much stronger nitrogen.
HD 80077, B2 Ia, also has a Prussian Helmet (Marco & Negueruela 2009, Fig. 14). See also van Genderen (2001).
Cygnus OB2-12, B4 Ia, as classified here from the Si iii/Si ii ratio in the high-S/N GOSSS spectrogram shown in Fig. 4, which shows its marked Prussian Helmet H profile (see also Clark et al. 2012 again showing possible variability).
4.4 Anomalous H Profiles in LBVs
Here we list several Luminous Blue Variables with references demonstrating the presence of H profiles similar to those in the VFTS objects discussed. It should be recalled that the intensities of these broad wings are only a few percent above the continuum, so high S/N is required to detect them, as well as adequate intensity and wavelength scales to display them. They may be lost in some published LBV spectrograms, especially the extensive early photographic material.
HR Carinae (HD 90177) is the poster child for this phenomenon in LBVs, because of both its prominence and the fact that the discussion by Hutsemékers & Van Drom (1991) is the only previous one addressing it to our knowledge. They show the asymmetrical H profiles at both minimum and maximum phases of the LBV variation (their Fig. 2). It can also be seen at minimum in Fig. 8 of Walborn & Fitzpatrick (2000). Although at the limit of the data quality, that figure likely also shows the same asymmetrical H wings in the neighboring LBV candidates HD 168607 and 168625, sans any central emission or P Cyg profile in the latter.
AG Carinae (HD 94910) likely also has this kind of H profile at maximum in Fig. 1 of Hutsemékers & Kohoutek (1988), although near the limit of the data and display.
R127 (HDE 269858) in the LMC shows the same asymmetrical H wings upon entering its current extended outburst in Fig. 3 of Stahl et al. (1983).
We have discussed peculiar H profiles in the spectra of five luminous hypergiants or supergiants in the 30 Doradus region, found in data from the VLT-FLAMES Tarantula Survey. These profiles contain extremely broad wings (total width 3000 km s), which are markedly asymmetrical, with a more extended, shallow slope shortward and a less extended, steeper slope longward. At the earlier B types a narrow, central P Cygni feature is also present (Prussian Helmet profile), whereas at mid-B to A types only the broad wings are seen. Identical profiles have been known for some time in a number of Galactic hypergiants. These H profiles are quite different from those of H in the same spectra, which have more symmetrical emission wings except for a broad absorption depression in the shortward one.
Most importantly, we have proposed explanations of these peculiar, discrepant Balmer profiles in terms of broad Diffuse Interstellar Bands curtailing the longward wing of H and depressing the shortward one of H.
Further, we have drawn attention to similar H profiles in a number of both active and quiescent Luminous Blue Variables or candidates thereto. In Figure 6, the locations in a theoretical Hertzsprung-Russell Diagram of several hypergiants/supergiants with such profiles, including those discussed here, are shown in comparison with some LBVs. Note that the objects below the Humphreys-Davidson (HD) Limit may be returning from the red supergiant region. It is reasonable to hypothesize that the hyper/supergiants may be approaching the LBV state for the first time, which would be of considerable interest to the understanding of late massive stellar evolution if substantiated by further investigation. In fact, Cyg OB2-12 lies beyond the HD Limit on the assumption of association membership, which may be consistent with the presence of very massive O stars (Wright et al. 2015 and references therein). We note that our VFTS and related Galactic objects are unlikely to already be quiescent or post-LBVs, because established members of those categories have additional spectral peculiarities such as Fe ii and [Fe ii] emission lines, whereas the only peculiarities in these spectra are the broad Balmer wings.
- affiliation: Operated by AURA, Inc., under NASA contract NAS5-26555.
- affiliation: Operated by AURA, Inc., under NASA contract NAS5-26555.
- Bernat, A. P., & Lambert, D. L. 1978, PASP, 90, 520
- Breysacher, J., Azzopardi, M., & Testor, G. 1999, A&AS, 137, 117
- Brott, I., de Mink, S. E., Cantiello, M., et al. 2011, A&A, 530, A115
- Clark, J. S., Najarro, F., Negueruela, I., et al. 2012, A&A, 541, A145
- Crowther, P. A., Lennon, D. J., & Walborn, N. R. 2006, A&A, 446, 279
- Crowther, P. A., Schnurr, O., Hirschi, R., et al. 2010, MNRAS, 408, 731
- Evans, C. J., Taylor, W. D., Hénault-Brunet, V., et al. 2011, A&A, 530, A108
- Evans, C. J., Kennedy, M. B., Dufton, P. L., et al. 2015, A&A, 574, A13
- Feast, M. W., Thackeray, A. D., & Wesselink, A. J. 1960, MNRAS, 121, 337
- Fitzpatrick, E. L. 1991, PASP, 103, 1123
- Herbig, G. H. 1995, ARA&A, 33, 19
- Hubeny, I., & Leitherer, C. 1989, PASP, 101, 114
- Humphreys, R. M., & Davidson, K. 1994, PASP, 106, 1025
- Hutsemékers, D., & Kohoutek, L. 1988, A&AS, 73, 217
- Hutsemékers, D., & Van Drom, E. 1991, A&A, 248, 141
- Kurt, C. M., & Dufour, R. J. 1998, RMxAC, 7, 202
- Lennon, D. J., Dufton, P. L., & Fitzsimmons, A. 1992, A&AS, 94, 569
- Maíz Apellániz, J. 2013, http://adsabs.harvard.edu/abs/2013hsa7.conf..583M
- Maíz Apellániz, J., Sota, A., Barbá, R. H., et al. 2014, IAU Symp, 297, 117
Maíz Apellániz, J., Sota, A., Walborn, N. R., et al. 2011,
- Marco, A., & Negueruela, I. 2009, A&A, 493, 79
- McEvoy, C. M., Dufton, P. L., Evans, C. J., et al. 2015, A&A, 575, A70
- Meixner, M., Gordon, K. D., Indebetouw, R., et al. 2006, AJ, 132, 2268
- Mel’nik, A. M., & Dambis, A. K. 2009, MNRAS, 400, 518
- Najarro, F.(P.), & Hillier, D. J. 2012, ASSL, 384, 67
- Parker, J. Wm. 1993, AJ, 106, 560
- Parker, J. Wm., Clayton, G. C.; Winge, C., et al. 1993, ApJ, 409, 770
- Pasquini, L., Avila, G., Blecha, A., et al. 2002, ESO Messenger, 110, 1
- Prevot-Burnichon, M. L., Prevot, L., Rebeirot, E., et al. 1981, A&A, 103, 83
- Puls, J., Kudritzki, R.-P., Santolaya-Rey, A. E., et al. 1998, ASPC, 131, 245
- Sabbi, E., Anderson, J., Lennon, D. J., et al. 2013, AJ, 146, 53
- Sabbi, E., Lennon, D. J., Gieles, M., et al. 2012, ApJ, 754, L37
- Sana, H., de Koter, A., de Mink, S. E., et al. 2013, A&A, 550, A107
- Sanduleak, N. 1970, CTIO Contribution, No. 89
- Santolaya-Rey, A. E., Puls, J., & Herrero, A. 1997, A&A, 323, 488
- Selman, F., Melnick, J., Bosch, G., et al. 1999, A&A, 347, 532
- Smith, L. J., Nota, A., Pasquali, A., et al. 1998, ApJ, 503, 278
- Stahl, O., Wolf, B., Klare, G., et al. 1983, A&A, 127, 49
- Stahl, O., Wolf, B., Leitherer, C., et al. 1985, A&AS, 61, 237
- Taylor, W. D., Evans, C. J., Simón-Díaz, S., et al. 2014, MNRAS, 442, 1483
- van Genderen, A. M. 2001, A&A, 366, 508
- Walborn, N. R. 1972, AJ, 77, 312
- Walborn, N. R. 1997, ASPC, 120, 8
- Walborn, N. R., & Blades, J. C. 1997, ApJS, 112, 457
- Walborn, N. R., & Fitzpatrick, E. L. 2000, PASP, 112, 50
- Walborn, N. R., & Howarth, I. D. 2000, PASP, 112, 1446
- Walborn, N. R., Stahl, O., Gamen, R. C., et al. 2008, ApJ, 683, L33
- Walborn, N. R., Sana, H., Simón-Díaz, S., et al. 2014, A&A, 564, A40
- Weis, K. 2003, A&A, 408, 205
- Wright, N. J., Drew, J. E., & Mohr-Smith, M. 2015, MNRAS, 449, 741