Free Access
Volume 583, November 2015
Article Number A98
Number of page(s) 7
Section Stellar structure and evolution
Published online 02 November 2015

© ESO, 2015

1. Introduction

The fluorescence of small neutral polycyclic aromatic hydrocarbon molecules (PAHs) in the Red Rectangle (RR) nebula in 2003, and in reflection nebulae, was discovered by Vijh et al. (2004; 2005a; 2006, respectively). As this emission is contained in a broad band, somewhere between 3550 and 4400 Å, with most of the peaks between 3700 Å, and 4050 Å, it was called blue luminescence (BL). Vijh et al. (2005a) found a very close correlation between the BL distribution and the 3.3 μm UIR-band emission. The latter is attributed to the C–H stretch band radiated by neutral PAH molecules. The central source of the RR nebula is HD 44179, a binary, and presumably a post-AGB star of spectral type A (Teff ~ 8000 K) on the verge of becoming a planetary nebula. To be neutral, PAH molecules should be shielded from direct radiation. In the case of the RR nebula, they probably lie in the shadow of the optically thick dust ring.

The present paper deals with a presumed BL effect in the L band (λ = 3680−4080 Å) of the Walraven VBLUW photometric system of the yellow hypergiant (YHG) HR 5171A = V766 Cen, reported by van Genderen (1992, Figs. 1 and 5, hereafter Paper I), made in 1971 and 1973 (and 1977, Sects. 3.1, Hereafter, these data sets are abbreviated to sets 71, 73, and 77, respectively. No proper explanation could be offered at the time. Surprisingly, that excess appeared to be absent in the VBLUW time series 19801991. We also investigated UBV photometry of HR 5171B made in 1971, and uvby photometry of HR 5171A, made in 1979 and between 19911994. The rationale was that the blue and ultraviolet bands have some overlap with the BL spectrum considering its FWHM ~ 450 Å1. Thus, if the excess was really due to BL, the location could be either Gum48d in the background, or the dense CS matter of HR 5171A.

HR 5171A is thought to be located in the old HII region Gum48d, known for its gas, dust, and molecular complexes. However, this would implicate that HR 5171A has a luminosity that is too high for a hypergiant, with the consequence that its distance should be smaller (Nieuwenhuijzen et al., in prep.). Moreover, Karr et al. (2009) concluded that the present energy budget of Gum 48d only needed one ionizing O-type star, i.e the present B-type star HR 5171B a few million years ago on the main sequence, just like HR 5171A roughly at the same time. Thus, this might be another argument that HR 5171A is no member of Gum48d, unless the contribution of the secondary’s brightness (unknown; Sect. 2) would cancel this inconsistency.

thumbnail Fig. 1

Fluxes log f0 for the VBLUW data sets as a function of the wavelength. The theoretical SED for a star with log g = 0, and T = 5000 K (blue curve) is fitted to the 71, 73 and 80–91 fluxes. The red dashed curve fitting the 77 set fluxes (apart from the L flux as it includes the suspected BL) does not agree with any model. At the bottom, the filter coverage for the VBLUW bands, and the BL spectrum for only the highest peak. The L excesses are indicated with triangular peaks.

2. Characteristics of HR 5171A and the Gum48d nebula

HR 5171A = V766 Cen = HD 119796, in 1971 of spectral type G8Ia+ (Humphreys et al. 1971), T ~ 5000 K (Sect. 3.3.) is an evolved massive star, showing various types of photometric instabilities. It is suspected to be a contact binary surrounded by a mid-IR nebula (Chesneau et al. 2014). It belongs to the very small group of yellow hypergiants (YHG; de Jager 1980; 1998; de Jager & van Genderen 1989). The distance according to Humphreys et al. (1971) is ~3.6 kpc, but may be less (Sect. 1). The star still exhibits a high mass-loss rate (Schuster et al. 2006; Schuster 2007; Karr et al. 2009; Chesneau et al. 2014). A historical photometric study from 1953 until 1992 was presented in Paper I. The star showed the signature of a star evolving to the red (which is likely not the case, Nieuwenhuijzen et al. in prep.), in view of the increasing reddening of the colour index (BV)J (UBV system) between 1960 and 1980 (Paper I; Chesneau et al. 2014).

Gum48d, containing the photo-dissociation region RCW 80 (Schuster et al. 2006; Schuster 2007), is located in the Centaurus Arm. At present, its ionization is sustained by the near-by blue star HR 5171B (Schuster 2007), a B0Ibp star (Humphreys et al. 1971), at a distance of ~10. Most of the UIR bands of Gum48d contain emission features of PAHs (Schuster et al. 2006; Schuster 2007; Karr et al. 2009). The 3.4 μm C–H stretch band emission from the WISE archive, which is the best signature for neutral PAHs (Bakes et al. 2004), turned out to be saturated at a flux of 95 Jy (Cutri et al. 2012; see also in Chesneau et al. 2014, the Appendix Table D.1). Thus, based on these observations, the nebula is a serious candidate to be the source for the BL. On the other hand, Chesneau et al. (2014), believe that HR 5171A is responsible for the high far- and mid-IR excess (see their Table D.1) as the instruments used have apertures not larger than a few arcsec, centred on HR 5171A. However, with respect to temperature, HR 5171A is too cool (~5000 K) to produce high energetic photons (3.55 eV; Witt, priv. comm.) to create PAH molecules in the highest energetic state. The presence of atmospheric instabilities, during which a huge amount of energy is released like in some other YHG (Lobel et al. 2003; Nieuwenhuijzen et al. 2012), is suspected. These processes could be well sufficient to produce these kinds of photons (Sect. 4).

HR 5171A shows a prominent 10 μm silicate feature (Humphreys et al. 1971) and a strong Na I 2.2 μm emission (Chesneau et al. 2014), also detected in other massive evolved stars like the YHG (Oudmaijer & de Wit 2013). This could point to the presence of a pseudo-photosphere in an optically thick wind, shielding the neutral sodium from ionization by direct starlight (Oudmaijer & de Wit 2013). These authors also claim that the NaI emission in other YHG originates in a region inside the dust condensation radius. Considering the controversy about the precise origin of the L excess: the CS material of HR 5171A and Gum48d, we compute the fluxes for these two alternatives. A possible instrumental effect responsible for a false emission in the L channel during the 1970s is also discussed (Sect. 3.1). We also discuss two emission features radiated by a hot chromosphere, often a feature exhibited by cool stars, which under abnormal hot and active conditions also radiate in the wavelength area of the L channel (Sect. 3.2).

3. The photometric observations

3.1. The VBLUW photometry. Other possible causes to explain the L–excess

A total of seven observations, made between 1971 and 1977, showing the strong brightness excess in the L band (λeff = 3838 Å), were obtained with the 90-cm Dutch telescope, which was then still located at the Leiden Southern Station in South Africa, Broederstroom (Hartebeestpoortdam). At the time, the telescope was equipped with the same simultaneous five-colour photometer of Walraven as between 19801991 (no L excess) after the telescope was moved to the ESO in 1979, La Silla, Chile (Lub & Pel 1977b; de Ruiter & Lub 1986; Pel & Lub 2007). It is customary to express the brightness and colour indici of the VBLUW system in log intensity scale. The observing and calibration procedures were similar for all data sets (Paper I). However, the aperture was 215 and/or 165 in the seventies (original telescope outputs are not available anymore), and 116 in the eighties. In both cases, the near-by blue star HR 5171B, in the visual 35 fainter than HR 5171A, could be satisfactorily kept out by decentring HR 5171A with 3–5.

Table 1

Observed photometric parameters of the four data sets and the brightness excesses in the L-band.

Two observations were obtained in 1971 (Pel 1976), two observations in 1973 (Paper I) and three in 1977 (Pel, priv. comm.). The time intervals between the individual observations in a set were a few days up to a few weeks. All photometric data are given in the 1980 system (Paper I). Those made in 1971 and 1973, on both sides of maximum No. 3 (see Fig. 1 of Paper I) are marked “L excess”. The 1977 observations are located in the rising branch to maximum No. 7, very near the top (not plotted in Fig. 1 of Paper I, as we obtained them from Pel only recently. The star was brighter and bluer than usual because of intrinsic changes (Sect. 3.3.).

The amount of L excess is based on the position in the theoretical two-colour diagram VB/BL (empirical version in Fig. 5 in Paper I) based on Kurucz-atmospheres (Castelli & Kurucz 2003; Castelli 2014). The excesses in BL represent the colour index difference between the observed and expected position (at roughly constant VB), if there was no excess in L.

Table 1 lists the observed photometric parameters for the three data sets of the 70s and the median values of the 8091 set. The parameters VJ and (BV)J were transformed from V and VB by formulae of Pel 1987). The last column lists the L excesses, labelled LEx.

The L excesses in the early 1970s, which is very obvious in the empirical two-colour diagram VB/BL, could not be due to any stray light from the blue nearby star HR 5171B; whatever the size of the diaphragm, as HR 5171A was always decentred in such a way that the blue star’s distance to the edge was enlarged. Moreover, HR 5171A would have shown progressively increasing excesses in U and W, if this was done insufficiently. The suggested presence of a mysterious blue star C in the diaphragm, at a 015 distance, nearly as bright as HR 5171A2, is unlikely. Otherwise, the U and W fluxes should have shown even larger excesses. Moreover, the 71/73 BU and UW indices (in 1977 the star brightened intrinsically, Sect. 3.3.) were very similar to those of 80/91 lacking any L excess!

Another possible explanation of the false L excesses, with the observed intensity, is if the “field of view” of the photomultiplier of the L channel was accidentally shifted by ~5, (with respect to the other channels) in such a manner that HR 5171B was shifted up to, or on the edge of the diaphragm. This scenario is considered almost impossible from a technical point of view (Pel, priv. comm.). Moreover, if this had been the case, it would mean that hundreds of program stars and tens of thousands of L measurements of comparison and standard stars during the seventies, were badly centred by about 5. The much less favourable seeing conditions at the location of the Southern Station, compared to the seeing at the ESO, would have entailed a much larger scatter for BL than for the three other colour indices. A quality check of, e.g. the 150 and 100 complete light and colour curves of Population I Cepheids (Pel 1976) and RR Lyrae stars (Lub 1977a), respectively, proves that this was not the case.

3.2. Other emission sources in the wavelength area of the L band

Two emission sources, which can give rise to radiation in the wavelength area of the L band, are related to a hot chromosphere. Cool stars like the twin sister of HR5171 A, ρ Cas, at times showed the signature of a hot chromosphere. One presumes that this is caused by the dissipation of mechanical energy delivered by the different types of pulsations and occasional shell ejections (Sargent 1961; Zsoldos & Percy 1991; Percy & Kolin 2000; Lobel et al. 2003). This object showed a variable Balmer continuum radiation starting at ~4200 Å, up to shorter wavelengths, in the absolute energy distribution for two nights in 1970 and two in 1974 (Joshi & Rautela 1978). The variable amplitude of this radiation at the location of the L band was ~03. As a second emission source, we have to mention the CaII H and K resonance lines, located in the red wing of the L band. Normally these resonance lines are in absorption, but in a very heated chromosphere they appear in emission, as in flare stars. There are convincing reasons to reject a hot chromosphere as the source of the L excesses in the 1970s (Sects. 3.3, 4).

3.3. On the location of the BL. The extinctions. The spectral energy distribution for the VBLUW data sets

An important issue is the location of the supposed BL. There are two possibilities. Scenario (1) in Gum48d, or (2) in the extended envelope of HR 5171A. Depending on the location, different corrections for extinction (interstellar IS and/or circumstellar CS) should be applied to the observed L fluxes.

Scenario (1): in this case the reddening of HR 5171B has been used: E(BV)J = 1.00 ± 0.10, a rounded-off average of the reddenings based on the photometry by Humphreys et al. (1971), Dean (1980) and Paper I. Using the extinction law for the diffuse IS extinction RJ = 3.2 (Cardelli et al. 1989), AVj = 3.20 ± 0.32. The equivalent values in the Walraven VBLUW system are: E(VB) = 0.43 ± 0.04 (log scale, transformation formula by Pel 1987), and RW = 3.2, in log scale: AV = 1.38 ± 015, and in the L band: AL = 1.97 ± 0.15 (Steenman & Thé 1989). The calibration constant is given by de Ruiter & Lub (1986).

Scenario (2): this scenario is of importance for the observed, non-stellar L excesses and for the stellar spectral energy distribution (SED). The total reddening for the four data sets was determined with the reddening lines in the two-colour diagrams defined by Pedicelli et al. (2008), and the position of the 8091 set as reference. This is because it was free of any L excess and if dereddened, its position in the VB/BL and BU/UW diagrams moves to T ~ 5000 K, adopting log  g = 0. This is in agreement with spectral analyses (Humphreys et al. 1971; Warren 1973; Chesneau et al. 2014; Lobel, priv. comm.; Oudmaijer, priv. comm.).

The reddening (IS and CS) in log scale was E(VB) = 0.63 ± 0.04, corresponding with E(BV)J = 1.41 ± 0.10 (in magnitudes, transformation formula from Pel 1987). It appeared that for luminous G-type stars the extinction law is RJ = 3.5 − 3.6 (Humphreys, priv. comm.) because of their extended envelopes with gas, dust, and molecular clouds. To find the R values for the three photometric systems Johnson UBV, Strömgren uvby, and Walraven VBLUW, we used observations (almost) simultaneously made in the three systems. We found RJ = 3.6, RS = 4.7, and RW = 3.3 for which we got a satisfactory match of the three SEDs3. Hence, the extinctions are AVj = 5.08 ± 0.30 (in mag) and AL = 2.96 ± 0.15 (in log scale). The extinction-free stellar fluxes log f0 resulted from the addition of the calibration constants (de Ruiter & Lub 1986) and the extinction-free brightnesses. They are listed in Table 2.

Table 2

Extinction-free brightnesses log I0 and fluxes log f0 for the five channels of the four data sets according to scenario (2).

Figure 1 shows the SEDs of HR 5171A for scenario (2) only; the green and blue symbols denote the 71, 73, and 80–91 sets, respectively, and the red symbols denote the 77 set. It is obvious that the brightness differences between the sets 71, 73 and 8091 were small, therefore, they could be used as one homogeneous set. It appears that a satisfactory match, even with the small positive Balmer jump, is obtained for V, B, (L only for the 80-91 set), and U with the model SED for 5000 K, adopting log g = 0, and vturb = 2 kms-1, using Kurucz models (Castelli & Kurucz 2003; Castelli 2014). The W flux is too high by about 0.2 log scale (~05) for all three sets, thus, possibly caused by chromospheric radiation. Obviously, these W fluxes represented very weak Balmer continuum excesses, otherwise U should show excesses as well, which should then be greater than those in L. Therefore, these L excesses should be due to some other source. The L excesses are indicated with the two triangles.

The SED of the 77 set is a different story: no match with any hotter theoretical SED is possible. The slope of the Paschen continuum (from V to B), is slightly steeper than for the other sets. Thus, a possible photospheric temperature rise is small (by a few hundred degrees), but the extreme high excess in the ultraviolet is stunning, with on top the L excess. The pulsational activity and/or outbursts of ionization energy, represented by an excess of optical light (see Fig. 1 in Paper I), heated the chromosphere presumably to such a high level that it resulted in an excessive Balmer continuum emission of up to more than 1m in U and W (Sect. 4.). It should be emphasized that the observed L brightnesses (71, 73, and 77) contain a stellar and a non-stellar contribution. The stellar contribution was likely subject to a much higher reddening than the non-stellar one (Sect. 3.3.1).

Table 3

Computation of the fluxes f(EL0) (in 10-10 Wm-2μm-1) of the BL, based on the brightness excesses in the L band after correction for interstellar extinction.

Table 4

Comparison of the observed 1979 uvby magnitudes of HR 5171A (corrected for the light of HR 5171B) on the third line with the 1991 observed uvby magnitudes on the fourth line.

3.3.1. The L excess fluxes

In Table 3, we list the date (set); the brightness of HR 5171A in the L band (taken from Table 1), including the excess LEx; the corrected brightness of HR 5171A in L without the excess (LLEx); the brightness of the excess (EL), which is the difference of the anti-logs of L and L-LEx., after transforming it into a log scale; and the unreddened fluxes, f(EL0), from scenario (1), in 10-10 Wm-2μm-1. The uncertainty is ±0.5, and, thus between 3020%, and mainly due to the 010 uncertainty in the reddening. As the reddening differences between the three data sets were small (Sect. 3.3.), the difference in the fluxes of 1971 and 1973/1977 should be considered real and suggests that the presumed BL was variable. The diameter of the apertures used were 21″5, and/or 16″5 (19711977, Sect. 3.1.), and 11″6 (19801991). These apertures are equivalent with an area with a diameter of 0.4, 0.3 and 0.2 pc in Gum48d, respectively, assuming a distance of 3.6 kpc. The effect of the Gum48d’s radiation travelling through HR 5171A’s envelope, with its much higher reddening (see scenario (2)), can be neglected because its diameter is about 3″4 (Chesneau et al. 2014), much smaller than the size of the apertures used.

Scenario (2): the dereddened fluxes f(EL0), are listed in the sixth column. It appears that the fluxes are now larger than for scenario (1) by a factor of ten (due to the higher reddening) as well as the uncertainty. Thus, 14.7 in 1971 and 21.4 ± 5 in 1973/1977.

Some doubt about the reliability of these high values, almost of the order of the stellar L flux, is appropriate. However, it is likely that the supposed BL source did not suffer from the same high extinction as the central star HR 5171A because its location could have been near the border of the envelope (e.g. facing the Earth), where the extinction is much lower. Then, the lower reddening used in scenario (1) would be much more appropriate, resulting in fluxes about ten times smaller. The envelope of HR 5171A may be more or less similar to the envelope surrounding the AGB star IRC +10216 (CW Leo, a Mira Ceta variable), modelled by McElroy et al. (2013), Li et al. (2014) and Li (2015). That model has an in the centre (for HR 5171A: ~5m), decreasing exponentially outward and being almost negligible halfway to the edge. Thus, by virtue of this model, it is much more realistic to suppose that the fluxes for the L excesses should have been near the low side of the total range 1.4 to 21. A dust model of HR 5171A by Lobel (priv. comm.) will be released soon.

Thus, then the remarkable size of the L excesses, as evident in Table 3, and in the empirical VB/BU diagram in Fig. 5 in Paper I, being almost of the order of the stellar flux, can be understood (see our comment on the size of the 1977 flux below, and in Sect. 4). Moreover, the star was, and still is, intrinsically very red.

It is evident that the flux measured in the 77 set contains the excess of the unknown source and those radiated by the hot chromosphere. Thus, the value listed in Table 3 represents the total excess flux of different sources.

3.4. The search for any emission in the UJ band (Johnson UBV system), and in the uv bands (uvby system)

Any uv observation of HR 5171A, made after the 1971–1977 time interval, and one UJ observation of its close neighbour HR 5171B, made in 1971 (by Humphreys, priv. comm.), has been carefully inspected and analysed. Despite their partial overlap with the BL spectrum no excess in UJ was found, and in uv only a small amount of excesses, which we do not consider conclusive. Therefore, they are only briefly discussed.

As the UJ magnitude of HR 5171B in 1971 was not affected by any measurable emission, we conclude that either the contribution to the L brightness of set 71 originated accidentally behind HR 5171A in a small area of Gum48d, or its origin was the CS material of HR 5171A.

Olsen (1983) obtained seven observations in the Strömgren uvby system in the period 1979, February 26 and October 6, with the Danish 50-cm reflector at the ESO in La Silla, Chile, equipped with the simultaneous uvby photometer (Grønberg et al. 1976a; Grønberg & Olsen 1976b). The aperture was 30″ and included HR 5171B, for which we had to correct before inspecting the v magnitude of HR 5171A for any brightness excess. The mean errors in the average photometric parameters appear to be very small (Olsen 1983).

As no uvby photometry has ever been made of HR 5171B, we had to use models (solar abundancy) for early B-type stars: those of Lester et al. (1986) and of Castelli & Kurucz (2003)/Castelli (2014). We adopted Teff = 26 500 K, log  g = 3 and 3.5, and a y magnitude almost equal to the observed VJ = 10.03 (Paper I). The VJ (~y) above was dereddened using E(BV)J = 1.00 ± 0.10 (Sect. 3.2), and RJ = 3.2. For the amount of extinction in the uvby channels, we consulted the tables of Steenman & Thé (1989).

Table 4 lists Olsen’s uvby observations of both stars together (first line), the mean of the two models for the blue star (second line), and then in the third line the “observed” models of HR 5171A only. Thus, the magnitudes on the second line subtracted from the magnitudes on the first line. The estimated uncertainties in the models are smaller than 01, with the exception of the u magnitude, which is ±02. The fourth line in Table 4 lists the mean of three uvby observations made in 1991 (see below).

Strömgren uvby observations were also made by the LTPV Group of Sterken (1983) from 1991 until 1994 with the 50-cm ESO telescope in La Silla, Chile (Manfroid et al. 1991; Sterken et al. 1993). The very first three uvby observations by the LTPV Group, coincided with the last six observations of the 80–91 VBLUW data set in maximum No. 17 (Fig. 1 in Paper I; van Leeuwen et al. 1998). We used these three 1991 uvby observations to compare them with the 1979 observations of Olsen (listed in Table 4, third line). Although the aperture used by the LTPV group was 20, star HR 5171B could generally be excluded by decentring HR 5171A with ~5.

thumbnail Fig. 2

Log f0, for observations made almost simultaneously in two different photometric systems (green and red symbols) as a function of λ with the indication of the λ coverage. UjBjVj for Johnson, uvby for Strömgren, VBLUW for Walraven. The blue curves are SEDs for 5000 K and 5250 K models. Vertical arrows: see the text in Sect. 3.5.

3.5. The spectral energy distribution for the 1979 and 1991 data sets in three photometric systems

We constructed an SED (Fig. 2) simultaneously using obtained BJVJ (UBV system) photometry by Dean (1980); see Fig. 1 in Paper I, together with the simultaneously observed uvby and VBLUW data from the 1991 sets. The band widths are indicated at the bottom with horizontal line pieces. The central mark indicates λeff. For v and u, two widths are indicated: the smallest was used by the LTPV group (Sterken 1983), and the largest was used by Olsen (1983). The BJ and VJ bands are only represented by their λeff. The band width (FWHM) of the BL for the RR nebula is indicated (Vijh et al. 2004; 2005b; 2006). The full line is the portion containing the highest emission peaks. The same reddening was used as in Sect. 3.3. (scenario 2), as well as the extinction laws for the three photometric systems, and the extinctions for each pass band. Calibration constants are from Cox (1999), for UBV; Helt et al. (1991); for ubvy, and de Ruiter & Lub (1986), for VBLUW.

Table 5

Reddening free fluxes of HR 5171A (Wm-2μm-1) measured in three different photometric systems.

The computed fluxes log  f0 are listed in Table 5. The obtained SEDs of the simultaneous sets 79 uvby and BJVJ, match very satisfactorily. They also match well the 91 sets uvby and VBLUW. Figure 2 shows brightness excess by about 0.10.2 (see arrows) in the u band in log f0 with respect to the combined SED. However, considering the errors, we are very hesitant to consider them as due to BL. Errors in the models and physical differences of the star between 1979 and 1991 are possible, but some contribution by Balmer continuum radiation is possible as well. The excess in W: ~03, represents a weak Balmer continuum emission, but weaker than during the 1970s and 1980s (Fig. 1: ~05, Sect. 3.3).

4. Discussion and conclusions

It is quite well possible that the brightness excesses in the L band (VBLUW system) detected in the photometry of the yellow hypergiant HR 5171A made in the 1970s were due to blue luminescence (BL) that was discovered in the RR nebula by Vijh et al. (2004; 2005a; 2006). Our suspicion is possible in view of the complete overlap of its spectrum, including the highest peaks, with the L band. We first considered various possibilities, which might have created false brightness excesses in the L band, but they turned out to be implausible. The same can be said about the following emission sources previously mentioned in Sect. 3.2.: the Balmer continuum emission due a hot chromosphere and the CaII H and K resonance lines (in the L band), which only appear in emission in extreme circumstances, when the chromosphere is extremely heated. However, in the first case one expects much higher brightness excesses in W and U than in L (see Fig. 1 in Joshi & Rautela 1978, showing an increasing Balmer continuum emission from λ4200 Å, to shorter wavelengths in the SED of ρ Cas). It is obvious in our Figs. 1 and 2 that apart from the 77 set observed during a large amplitude pulsation showing high Balmer continuum excesses in W and U, the 71 and 73 sets show no excesses in U and the excesses in W are smaller than in L. Note that in the 1980s there was no L excess at all other than a W excess, presumably due to the presence of a chromosphere. We conclude that the L excesses in 1971 and 1973 were not caused by Balmer continuum emission. Obviously, the L excess of the 77 set in Fig. 1 contains the contribution from the supposed BL source, but also should include Balmer continuum emission, and perhaps also some emission from the CaII doublet. Therefore, the L flux listed in Table 3 is too high for the supposed BL alone, indicating that this emission was declining. After all, in 1979, no excess in that wavelength area of the uvby system was recorded, apart from some excess in the u band, presumably from uncertainties in the atmosphere models and/or from some Balmer continuum radiation (Fig. 2).

The unreddened L fluxes, depending on the location of the BL source, in units of 10-10 Wm-2μm-1, varied between 1.4 and 2.2 (71 and 73 sets), when located in Gum48d (scenario 1). If the unreddened L fluxes were located in the dense centre of the CS envelope of HR 5171A (scenario 2) then the values would be a factor of 10 higher, i.e. 15–21. However, these values were much lower, approaching the lower values that we prefer, if they were located near the outskirts of the envelope.

No significant trace of an excess due to BL was detected in other bands apart from the L band. The distribution of the highest BL peaks in the direction of HR 5171A occurred in the L band, which does not contradict the possibility that BL was the cause.

As mentioned above, no L excess was present in the 1980–1991 VBLUW time series, while the colour indices BU and UW were almost the same as in the early seventies (1971/1973), which is in contrast to the BL colour index. The short duration of the supposed BL emission could mean that we were dealing with a small-scale source. The dynamical circumstances of the CS matter surrounding HR 5171A would then be a good candidate, as shielded clouds containing neutral PAH molecules can quickly be exposed to highly energetic photons. On the other hand, one may wonder whether this cool star (~5000 K) was able to produce high energetic photons at all to excite the neutral PAH molecules to the upper electronic state. As for this requirement, Gum48d, with its hot exciting star HR 5171B, was by far the best candidate. However, it is very plausible that these photons were produced by HR 5171A itself. The star showed a number of relative short-lasting “light outbursts” during the 1970s, being blue near maximum light (Fig. 1 in Paper I; see SED of maximum No. 7 in Fig. 1 of the present paper). Amplitudes were ~05 in V, with a timescale of ~12 yr. Similar type of light curve features were exhibited by the YHG ρ Cas (Lobel et al. 2003) and HR 8752 (Nieuwenhuijzen et al. 2012). These authors concluded that they showed at times strong atmospheric ionization driven instabilities. As a result, explosive “flash” outbursts occurred by the release of stored ionization energy due to H recombination as the temperature of the atmosphere declined. These processes might have delivered the necessary high energetic photons (3.55 eV, Witt, priv. comm.).


Geneva seven-colour photometry exists as well, made in the 1970s and 1980s. The catalogue only lists two averages based on an unknown number of observations, and precise dates are lacking (Rufener 1988).


(Dommanget & Nys 1994; Dommanget 2002; Mason et al. 2001, WDS13472-6235; Mason & Hartkopf 2006; ESA 1997, HIP 67261; van Leeuwen, priv. comm.). However, its existence is still a matter of debate: HST images do not show this component (Schuster & Humphreys, priv. comm.). On the other hand, the analysis of three years of Hipparcos photometry, revealed two 4d5 oscillations with a light amplitude of 005 (Eyer 1998). Thus, a variable star should be inside the aperture used by the Hipparcos satellite.


This matching eliminated small reddening differences, which are evident from the VB values between the 70s and 80s; see Table 1. We emphasize that these R values are only applicable to this star. The reason is that the cause of the peculiar and long-lasting stellar reddening is not well understood.


We are very grateful to a number of colleagues: Prof. Xander Tielens for his support, who drew our attention to the papers of Vijh et al. (2004; 2005a); Prof. Jan Willem Pel for offering his 1977 VBLUW photometry; and the following colleagues for many fruitful discussions, help, and invaluable advices: Dr. Floor van Leeuwen, Prof. Jan Willem Pel, Prof. Adolf Witt, Dr. Jan Lub, Prof. Roberta Humphreys, and Dr. René Oudmaijer, and especially to Dr. Erik-Heyn Olsen, who spent a lot of time studying his log books, making it easy for us to establish the date of the seven uvby observations. We took much advantage of his and Dr. Chris Sterken’s expertise on the Strömgren photometric system. A.L. acknowledges partial financial support by the Belgian Federal Science Policy Office under contract No. BR/143/A2/BRASS and in connection with the ESA PRODEX programmes “Gaia-DPAC QSOs” and “Binaries, extreme stars and solar system objects” (contract C90290). We are grateful to the referee for the invaluable comments and suggestions.


  1. Bakes, E. L. O., Bauschlicher, C., & Tielens, A. G. G. M. 2004, in Astrophysics of Dust, eds. A. N. Witt, G. C. Clayton, & B. T. Draine (San Fransisco: ASP), ASP Conf. Ser., 309, 731 [Google Scholar]
  2. Cox, A. N. 1999, Allen’s Astrophysical Quantities, 4th edn. (Springer) [Google Scholar]
  3. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
  4. Castelli, F. 2014, [Google Scholar]
  5. Castelli, F., & Kurucz, R. L. 2003, IAU Symp., 210, 47 [NASA ADS] [Google Scholar]
  6. Chesneau, O., Meilland, A., Chapellier, E., et al. 2014, A&A, 563, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Cutri, R. M., et al. 2012, VizieR Online Data Catalog: II/311 [Google Scholar]
  8. Dean, J. F. 1980, IBVS, 1796, 1 [NASA ADS] [Google Scholar]
  9. Dommanget, J. 2002, VizieR Online Data Catalog: I/274 [Google Scholar]
  10. Dommanget, J., & Nys, O. 1994, Cat. of Components of Double and Multiple Stars, Comm. Obs. R. Belgium, Ser. A, No. 115 [Google Scholar]
  11. ESA 1997, The Hipparcos and Tycho Catalogues, ESA SP–1200, Vols. 10−12, DC 255 [Google Scholar]
  12. Eyer, L. 1998, Ph.D. Thesis, Observatoire de Genève [Google Scholar]
  13. Grønbech, B., & Olsen, E. H. 1976, A&AS, 25, 213 [NASA ADS] [Google Scholar]
  14. Grønbech, B., Olsen, E. H., & Strömgren, B. 1976, A&AS, 26, 155 [NASA ADS] [Google Scholar]
  15. Helt, B. E., Bibo, E. A., Cuypers, J., et al. 1991, A&AS, 89, 399 [NASA ADS] [Google Scholar]
  16. Humphreys, R. M., Strecker, D. W., & Ney, E. P. 1971, ApJ, 167, L35 [NASA ADS] [CrossRef] [Google Scholar]
  17. de Jager, C. 1980, The Brightest Stars (Dordrecht, Holland: D. Reidel Publ. Co.) [Google Scholar]
  18. de Jager, C. 1998, A&ARv, 8, 145 [NASA ADS] [CrossRef] [Google Scholar]
  19. de Jager, C., & van Genderen, A. M. 1989, A&SS, 157, 127 [Google Scholar]
  20. Joshi, S. C., & Rautela, B. S. 1978, MNRAS, 183, 55 [NASA ADS] [CrossRef] [Google Scholar]
  21. Karr, J. L., Manoj, P., & Ohashi, N. 2009, ApJ, 697, 133 [NASA ADS] [CrossRef] [Google Scholar]
  22. Lester, J. B., Gray, R. O., & Kurucz, R. L. 1986, ApJS, 61, 509 [NASA ADS] [CrossRef] [Google Scholar]
  23. Li, X. 2015, Ph.D. Thesis, Sterrewacht Leiden, Leiden Observatory, Leiden University [Google Scholar]
  24. Li, X., Millar, T. J., Walsh, C., et al. 2014, A&A, 568, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Lub, J. 1977, A&AS, 29, 345 [NASA ADS] [Google Scholar]
  26. Lub, J., & Pel, J. W. 1977, A&A, 54, 137 [NASA ADS] [Google Scholar]
  27. Lobel, A., Dupree, A. K., Stefanik, R. P., et al. 2003, ApJ, 583, 923 [NASA ADS] [CrossRef] [Google Scholar]
  28. Manfroid, J., Sterken, C., Bruch, K., et al. 1991, A&AS, 87, 481 [NASA ADS] [Google Scholar]
  29. Mason, B. D., & Hartkopf, W. I. 2006, Washington Double Star Cat., USNO Double Star CD 2006.5 [Google Scholar]
  30. Mason, B. D., Wycoff, G. L., Hartkopf, W. I., et al. 2001, AJ, 122, 3466 [Google Scholar]
  31. McElroy, D., Walsh, C., Markwick, A. J., et al. 2013, A&A, 550, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Nieuwenhuijzen, H., de Jager, C., Kolka, I., et al. 2012, A&A, 546, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Olsen, E. H. 1983, A&AS, 54, 55 [NASA ADS] [Google Scholar]
  34. Oudmaijer, R. D., & de Wit, W.-J. 2013, A&A, 551, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Pedicelli, S., Lub, J., Pel, J. W., et al. 2008, Mem. Soc. Astron. Ital., 79, 539 [NASA ADS] [Google Scholar]
  36. Pel, J. W. 1976, A&AS, 24, 413 [NASA ADS] [Google Scholar]
  37. Pel, J. W. 1987, Internal Report, Leiden Observatory [Google Scholar]
  38. Pel, J. W., & Lub, J. 2007, in The Future of Photometric, Spectrophotometric and Polarimetric Standardization, ed. C. Sterken, ASP Conf. Ser., 364, 63 [Google Scholar]
  39. Percy, J. R., & Kolin, D. L. 2000, PASP, 112, 363 [NASA ADS] [CrossRef] [Google Scholar]
  40. Rufener, F. 1988, Catalogue of Stars measured in the Geneva Observatory Photometric system (Sauverny: Geneva Observatory) [Google Scholar]
  41. de Ruiter, H. R., & Lub, J. 1986, A&AS, 63, 59 [NASA ADS] [Google Scholar]
  42. Sargent, W. L. W. 1961, ApJ, 134, 142 [NASA ADS] [CrossRef] [Google Scholar]
  43. Schuster, M. T. 2007, Ph.D. Thesis, University Minnesota [Google Scholar]
  44. Schuster, M. T., Humphreys, R. M., & Marengo, M. 2006, AJ, 131, 603 [NASA ADS] [CrossRef] [Google Scholar]
  45. Steenman, H., & Thé, P. S. 1989, A&SS, 189, 189 [Google Scholar]
  46. Sterken, C. 1983, The Messenger, 33, 10 [Google Scholar]
  47. Sterken, C., Manfroid, J., Anton, K., et al. 1993, A&AS, 102, 79 [NASA ADS] [Google Scholar]
  48. van Genderen, A. M. 1992, A&A, 257, 177 (Paper I) [NASA ADS] [Google Scholar]
  49. van Leeuwen, F., van Genderen, A. M., & Zegelaar, I. 1998, A&AS, 128, 117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Vijh, U. P., Witt, A. N., & Gordon, K. D. 2004, ApJ, 606, L65 [NASA ADS] [CrossRef] [Google Scholar]
  51. Vijh, U. P., Witt, A. N., & Gordon, K. D. 2005a, ApJ, 619, 368 [NASA ADS] [CrossRef] [Google Scholar]
  52. Vijh, U. P., Witt, A. N., & Gordon, K. D. 2005b, ApJ, 633, 262 [NASA ADS] [CrossRef] [Google Scholar]
  53. Vijh, U. P., Witt, A. N., York, D. G., et al. 2006, ApJ, 653, 1336 [NASA ADS] [CrossRef] [Google Scholar]
  54. Warren, P. R. 1973, MNRAS, 161, 427 [NASA ADS] [Google Scholar]
  55. Zsoldos, E., & Percy, J. R. 1991, A&A, 246, 441 [NASA ADS] [Google Scholar]

All Tables

Table 1

Observed photometric parameters of the four data sets and the brightness excesses in the L-band.

Table 2

Extinction-free brightnesses log I0 and fluxes log f0 for the five channels of the four data sets according to scenario (2).

Table 3

Computation of the fluxes f(EL0) (in 10-10 Wm-2μm-1) of the BL, based on the brightness excesses in the L band after correction for interstellar extinction.

Table 4

Comparison of the observed 1979 uvby magnitudes of HR 5171A (corrected for the light of HR 5171B) on the third line with the 1991 observed uvby magnitudes on the fourth line.

Table 5

Reddening free fluxes of HR 5171A (Wm-2μm-1) measured in three different photometric systems.

All Figures

thumbnail Fig. 1

Fluxes log f0 for the VBLUW data sets as a function of the wavelength. The theoretical SED for a star with log g = 0, and T = 5000 K (blue curve) is fitted to the 71, 73 and 80–91 fluxes. The red dashed curve fitting the 77 set fluxes (apart from the L flux as it includes the suspected BL) does not agree with any model. At the bottom, the filter coverage for the VBLUW bands, and the BL spectrum for only the highest peak. The L excesses are indicated with triangular peaks.

In the text
thumbnail Fig. 2

Log f0, for observations made almost simultaneously in two different photometric systems (green and red symbols) as a function of λ with the indication of the λ coverage. UjBjVj for Johnson, uvby for Strömgren, VBLUW for Walraven. The blue curves are SEDs for 5000 K and 5250 K models. Vertical arrows: see the text in Sect. 3.5.

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.