Discovery of a complex linearly polarized spectrum of Betelgeuse dominated by depolarization of the continuum Based on observations obtained at the Télescope Bernard Lyot (TBL) at Observatoire du Pic du Midi, CNRS/INSU and Université de Toulouse, France.

# Discovery of a complex linearly polarized spectrum of Betelgeuse dominated by depolarization of the continuum ††thanks: Based on observations obtained at the Télescope Bernard Lyot (TBL) at Observatoire du Pic du Midi, CNRS/INSU and Université de Toulouse, France.

M. Aurière Université de Toulouse, UPS-OMP, Institut de Recherche en Astrophysique et Planétologie, Toulouse, France CNRS, UMR5277, Institut de Recherche en Astrophysique et Planétologie, 57 Avenue d’Azereix,65000 Tarbes, France    A. López Ariste Université de Toulouse, UPS-OMP, Institut de Recherche en Astrophysique et Planétologie, Toulouse, France CNRS, UMR5277, Institut de Recherche en Astrophysique et Planétologie, 14 Avenue Edouard Belin, 31400 Toulouse, France    P. Mathias Université de Toulouse, UPS-OMP, Institut de Recherche en Astrophysique et Planétologie, Toulouse, France CNRS, UMR5277, Institut de Recherche en Astrophysique et Planétologie, 57 Avenue d’Azereix,65000 Tarbes, France    A. Lèbre LUPM - UMR5299 - CNRS and Université Montpellier - Place E. Bataillon, 34090 Montpellier, France    E. Josselin LUPM - UMR5299 - CNRS and Université Montpellier - Place E. Bataillon, 34090 Montpellier, France    M. Montargès Institut de Radioastronomie Millimétrique, 300 rue de la Piscine, 38406, Saint Martin d’Hères, France LESIA, Observatoire de Paris, PSL Research University, CNRS UMR 8109, Sorbonne Universités, UPMC, Université Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, F-92195 Meudon, France    P. Petit Université de Toulouse, UPS-OMP, Institut de Recherche en Astrophysique et Planétologie, Toulouse, France CNRS, UMR5277, Institut de Recherche en Astrophysique et Planétologie, 14 Avenue Edouard Belin, 31400 Toulouse, France
A. Chiavassa
Laboratoire Lagrange, Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Blvd de L’Observatoire CS 34229,06004 Nice Cedex 4, France
F. Paletou Université de Toulouse, UPS-OMP, Institut de Recherche en Astrophysique et Planétologie, Toulouse, France CNRS, UMR5277, Institut de Recherche en Astrophysique et Planétologie, 14 Avenue Edouard Belin, 31400 Toulouse, France    N. Fabas Kiepenheuer-Institut for Sonnenphysik, 79 104 Freiburg, Germany    R. Konstantinova-Antova Institute of Astronomy and NAO, Bulgarian Academy of Sciences, 72 Tsarigradsko shose, 1784 Sofia, Bulgaria    J.-F. Donati Université de Toulouse, UPS-OMP, Institut de Recherche en Astrophysique et Planétologie, Toulouse, France CNRS, UMR5277, Institut de Recherche en Astrophysique et Planétologie, 14 Avenue Edouard Belin, 31400 Toulouse, France    J.H. Grunhut European Southern Observatories, Karl-Schwarzschild-Str. 2, D-85748, Garching, Germany    G.A. Wade Department of Physics, Royal Military College of Canada, PO Box 17000, Station ‘Forces’, Kingston, Ontario, Canada K7K 4B4    F. Herpin Université de Bordeaux, LAB, UMR 5804, F-33270, Floirac, France CNRS, LAB, UMR 5804, F-33270, Floirac, France    P. Kervella LESIA, Observatoire de Paris, PSL Research University, CNRS UMR 8109, Sorbonne Universités, UPMC, Université Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, F-92195 Meudon, France Unidad Mixta Internacional Franco-Chilena de Astronomía (UMI 3386), CNRS/INSU, France & Departamento de Astronomía, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile.    G. Perrin LESIA, Observatoire de Paris, PSL Research University, CNRS UMR 8109, Sorbonne Universités, UPMC, Université Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, F-92195 Meudon, France    B. Tessore LUPM - UMR5299 - CNRS and Université Montpellier - Place E. Bataillon, 34090 Montpellier, France
###### Key Words.:
stars: individual: Betelgeuse – stars: polarization – stars: late-type – stars: supergiants
###### Abstract

Context:Betelgeuse is an M supergiant that harbors spots and giant granules at its surface and presents linear polarization of its continuum.

Aims:We have previously discovered linear polarization signatures associated with individual lines in the spectra of cool and evolved stars. Here, we investigate whether a similar linearly polarized spectrum exists for Betelgeuse.

Methods:We used the spectropolarimeter Narval, combining multiple polarimetric sequences to obtain high signal-to-noise ratio spectra of individual lines, as well as the least-squares deconvolution (LSD) approach, to investigate the presence of an averaged linearly polarized profile for the photospheric lines.

Results:We have discovered the existence of a linearly polarized spectrum for Betelgeuse, detecting a rather strong signal (at a few times 10 of the continuum intensity level), both in individual lines and in the LSD profiles. Studying its properties and the signal observed for the resonant Na i D lines, we conclude that we are mainly observing depolarization of the continuum by the absorption lines. The linear polarization of the Betelgeuse continuum is due to the anisotropy of the radiation field induced by brightness spots at the surface and Rayleigh scattering in the atmosphere. We have developed a geometrical model to interpret the observed polarization, from which we infer the presence of two brightness spots and their positions on the surface of Betelgeuse. We show that applying the model to each velocity bin along the Stokes and profiles allows the derivation of a map of the bright spots. We use the Narval linear polarization observations of Betelgeuse obtained over a period of 1.4 years to study the evolution of the spots and of the atmosphere.

Conclusions:Our study of the linearly polarized spectrum of Betelgeuse provides a novel method for studying the evolution of brightness spots at its surface and complements quasi-simultaneous observations obtained with PIONIER at the VLTI.

## 1 Introduction

The study of the linearly polarized spectrum of the Sun enabled the discovery of a very rich variety of structures which opened up ‘a new window for diagnostics of the Sun’ (Stenflo, 1997). This spectrum is composed of structures resulting partly from intrinsic polarization of the lines and also from the depolarization of the continuum by the absorption lines. The existence of this ‘second solar spectrum’ is due to the existence of anisotropy of the radiation field, induced by limb darkening. In the stellar case, it is considered to be very difficult to observe the second solar spectrum, since in general distant stars present a symmetric aspect (e.g. López Ariste et al., 2011). However, in the case of cool and evolved low- and intermediate-mass stars, linear polarization has already been reported in individual lines of variable Mira stars: in the Balmer hydrogen lines (McLean & Coyne, 1978; Fabas et al., 2011; Fabas, 2011) and in the resonance line of Ca i at 422.6 nm (Boyle et al., 1986). In addition Lèbre et al. (2014) discovered as yet unexplained nonzero Stokes and signals when averaging a mean profile about 10 000 atomic lines of the Mira star $χ$ Cyg observed at a luminosity maximum, using the least squares deconvolution (LSD) approach (Donati et al., 1997). Similar LSD signals as well as intrinsic polarization in individual lines such as Na i D2 and Sr i at 460.7 nm were discovered in other Miras and RV Tauri stars (Lèbre et al., 2015). It therefore appears that anisotropy of the radiation field may also exist in very evolved cool stars, and that a rich linearly polarized spectrum does exist in these pulsating variable stars.

We report here the subsequent study of the linearly polarized spectrum of Betelgeuse ($α$ Ori, HD 39801) which is a nearby M2 Iab supergiant, and which may be considered as a high-mass counterpart to the evolved AGB stars (as are the Miras). Being one of the stars with the largest apparent diameters, it has naturally deserved direct imaging and interferometric studies of its surface, which was found to deviate from circular symmetry and to present variable behavior (Wilson et al., 1997; Haubois et al., 2009; Kervella et al., 2015). A surface magnetic field has also been detected which appears to vary on a timescale of weeks/months (Aurière et al., 2010; Petit et al., 2013; Bedecarrax et al., 2013, Mathias et al. in preparation). In this paper, we demonstrate that Stokes observations reveal a rich linearly polarized spectrum suitable for analysis using both individual and averaged spectral line profiles.

In Section 2 we present the Narval spectropolarimetric observations, in Section 3 we describe the detection of the linear polarization signal and we show that the detected signal is of stellar origin. In Section 4 we infer the origin of the polarization together with some of its properties. Finally, we present a geometric model in Section 5, study the variations of the linear polarization with time in Section 6, and give our conclusions in Section 7.

## 2 Observations with Narval

Observations of Betelgeuse were obtained at the 2m Télescope Bernard Lyot (TBL) using the Narval spectropolarimeter, which is a twin of the ESPaDOnS instrument at the Canada-France-Hawaii Telescope (CFHT) (Donati, 2003; Donati et al., 2006). The observations that include linear polarization span from November 2013 to April 2015. A standard polarization observation consists of a series of four sub-exposures between which the half-wave retarders (Fresnel rhombs) are rotated to exchange the paths of the orthogonally-polarized beams within the whole instrument (and therefore the positions of the two spectra on the Charge-Couple Device (CCD), thereby reducing spurious polarization signatures. To avoid saturation of the CCD, we used 3 s exposure times for each sub-exposure. The extraction of the spectra, including wavelength calibration, correction to the heliocentric frame and continuum normalisation, was performed using the Libre-ESpRIT package (Donati et al., 1997), installed both at CFHT and TBL. The extracted spectra are in ASCII format, and consist of the normalized Stokes () spectrum and Stokes (, ) parameters as a function of wavelength, together with their associated uncertainties (where represents the continuum intensity). In the Libre-ESpRIT reduction the continuum polarization level is automatically removed. In the case of Betelgeuse for which a polarized continuum does exist (see Section 4.1), the level of polarization of the continuum cannot be recovered accurately from our spectra since Narval is not designed for this purpose. Also included in the output are ‘diagnostic null’ spectra , which are in principle featureless, and therefore serve to diagnose the presence of spurious contributions to the polarized spectra. Each single spectrum used in this work has a peak signal-to-noise ratio (S/N) in Stokes per 1.8 km s spectral bin between 1700 and 2100. Details of the observation and reduction procedures are provided by Donati et al. (1997) and Aurière et al. (2010, 2009). Table 1 reports the log of observations.

To obtain a high precision diagnosis of the spectral line polarization, the least-squares deconvolution approach (LSD, Donati et al., 1997) was applied to each reduced Stokes , and spectrum. We used a solar abundance line mask calculated from data provided by the Vienna Atomic Line Database (VALD, Kupka et al., 1999) for an effective temperature of 3750 K, surface gravity , and a microturbulence of 4.0 km s, consistent with the physical parameters of Betelgeuse (Josselin & Plez, 2007; Lambert et al., 1984). The complete mask contains about 15 000 atomic lines with a central depth greater than 40% of the continuum (which are below 60% of the continuum). Specific elements with transitions expected to trace a shocked region, a chromosphere or a circumstellar medium (such as H , Ca i, Na i) have not been considered in the mask. As supported by the works of Semel et al. (2009) and Paletou (2012) for linear polarization due to diffusion processes, we adopted an equal weight for each atomic line considered in the mask. This allows very strong detections of linear polarization in each LSD profile. On the other hand, to study the polarization of individual lines required the averaging of 8 and 8 spectra.

In the final eight columns of Table 1 we give the polarimetric parameters related to two bright spots inferred from our measurements (see Section 5.): observed maximum of linear polarization , polarization angle , position angle , and projection angle to disk center , as described in Section 5.

## 3 Stellar origin of the linearly polarized spectrum of Betelgeuse

Figure 1 shows the averaged LSD profiles obtained in November-December 2013 and January 2014 (hereafter named Set 1), and September 2014 (hereafter named Set 2). Note that both sets consist of the average of eight Stokes and profiles (Table 1). Figure 1 shows that the Stokes and signals full amplitude are up to the 7-9 level and centered on the mean profile. For clarity, we also show the heliocentric radial velocity of Betelgeuse (HRV hereafter), about 21 km s(Mathias et al., in preparation). For comparison, the circular polarization Stokes signal obtained during the same nights is about ten times weaker than the linear polarization (Mathias et al. in preparation, see also Aurière et al., 2010). Finally, Fig. 1 presents, at the same scale, the null polarization profiles : they do not show any signal above the noise level of a few . Since these profiles (including Stokes ) are obtained from different combinations of spectra acquired through the same rhombs but with different orientations, it demonstrates that the linear polarization observed in Betelgeuse is not a spurious signal due to Narval. This statement is also supported by the non detection of linear polarization in AGB stars which is described later in this section.

Figure 2 illustrates, for the same sets of data shown in Fig. 1, the linear polarization spectrum in the region of the D1 and D2 Na i lines (we recall that for all the presented polarized spectra the continuum polarization level is set to zero). In addition to the Na i D lines, this figure also shows a blend of lines of Ni I and Fe I (at 589.3 nm, between D2 and D1) and a line of Ti i at 589.9 nm. Other examples of photospheric lines are presented in Fig. 3. The structure of the individual polarized lines is similar to that of the LSD profiles (the case of the resonance lines of Na i will be studied in Section 4.2). Slow time variations from November 2013 to April 2015 of individual and LSD profiles are observed (see Section 6.1 and Fig. 8). These observations suggest that the polarized signal is related to the lines of Betelgeuse.

The linear polarization spectrum of cool evolved stars has already been detected using the LSD method in a number of Miras (Lèbre et al., 2014, 2015) and RV Tauri stars (Sabin et al., 2015; Lèbre et al., 2015). In these cases, the polarization seems to be linked to the atmospheric shock wave (Fabas et al., 2011) since its shape changes with pulsation phase, being maximum near maximum brightness, which corresponds to the shock passage through the atmosphere. The and profiles change from one Mira star to another and, for a given star, from one cycle to another. On the other hand, no linear polarization was detected at the 5 x level on four magnetic, non-pulsating, AGB stars of about the same effective temperature (namely $β$ Peg, 15 Tri, 8 And & $β$ And) studied with Narval (Lèbre et al. in preparation). For these stars, a significant Zeeman Stokes signal has been observed previously (Konstantinova-Antova et al., 2010, 2013).

The observations reported above show that the linearly polarized spectrum of Betelgeuse, as well as that observed in Miras and RV Tauri stars, is not due to a spurious signal from the polarimeters but has a stellar origin. In addition, all stars for which linear polarization in the photospheric lines has been detected belong to a category of cool stars for which linear polarization of the continuum has also been observed.

At this point we can argue that this linear polarization is not of Zeeman origin: it is much stronger than the circular polarization observed in Betelgeuse, whereas Zeeman linear polarization is typically ten times smaller than the associated circular polarization (e.g. Wade et al., 2000; Rosén et al., 2013, 2015). We will also show in Section 4 that the polarization decreases with wavelength, a phenomenon that is also in contradiction with a Zeeman origin.

## 4 Interpretation of the linear polarization in the atomic lines of Betelgeuse

### 4.1 Linear polarization of the continuum of Betelgeuse and its depolarization by atomic lines

Betelgeuse is known to present linear polarization of its continuum beyond the 0.5 % level in the blue (e.g. Hayes, 1984; Clarke & Schwarz, 1984; Nordsieck et al., 1994; Magalhaes & Nordsieck, 2000, for the most recent studies). This is a clear sign of the existence of anisotropy in the radiation field: it is generally believed that the linear polarization of the continuum is mainly due to photospheric Rayleigh scattering and Mie scattering by circumstellar dust with an anisotropy induced by one or several brightness spots (Schwarz & Clarke, 1984; Clarke & Schwarz, 1984; Doherty, 1986; Magalhaes & Nordsieck, 2000).

From our high S/N, high resolution linearly polarized spectra of Betelgeuse, we observe significant linear polarization signal associated with most of the relatively strong lines. This behavior is reminiscent of that observed on the Sun and called the ‘second solar spectrum’, except that for the latter the anisotropy of the radiation field is due to limb darkening (Stenflo et al., 2000b). The second solar spectrum has two main contributions (Stenflo et al., 1983a): i) depolarization of the continuum by the lines; ii) intrinsic polarization of individual lines. While the majority of the lines are expected to depolarize the continuum (Fluri & Stenflo, 2003), only a small number of them are expected to provide significant intrinsic polarization contributions (Stenflo et al., 1983a, b). While explaining the origin of the intrinsic polarization of individual lines is the subject of many publications, the case of the depolarization of the continuum is less well documented in the literature. This is because each intrinsic line polarization is unique and has to be explained by atomic physics; furthermore the involved coherent scattering is sensitive to magnetic field through the Hanlé effect (e.g. Stenflo, 2009) which has been applied to the investigation of weak fields on the Sun. On the other hand the physics involved in the process of line depolarization is the same for all the lines. We propose two levels of interpretation: Stenflo et al. (1983a) present empirical relationships which show that ‘the polarization line profiles of depolarizing lines appear to have approximately the same shape as their corresponding intensity profiles’; then Fluri & Stenflo (2003) present a theoretical analysis of the depolarizing lines; their properties may be also inferred from the work of Landi Degl’Innocenti & Landolfi (2004) from which it can be derived that the depolarization is mainly due to absorption and spontaneous re-emission of any photon by the atom at the origin of the spectral line: this is an incoherent process which emits an unpolarized photon, whatever the incident photon is polarized or unpolarized. In the present work the process that we call ‘depolarization of the continuum’ of Betelgeuse is the same as that described above for the Sun.

In all our observations of Betelgeuse the Na i D lines appear strongly polarized and D1 is as polarized, or even more polarized, than D2 (see Set 2 of Fig. 2). However, the intrinsic polarizability of the D1 line is much smaller than that of D2 (Landi Degl’Innocenti & Landolfi, 2004). The strong polarization of D1 in the spectrum of Betelgeuse, comparable to that of D2, cannot be intrinsic. We propose that it originates from the depolarization of the continuum. Since in our Narval data we set the (unmeasured) continuum polarization level to zero, any depolarization appears as a net signal, as observed and will be illustrated later in Section 5.3 . Next to D1, D2 also depolarizes the continuum and, in addition to that signal, the intrinsic polarization of this line is contributed. This supports the notion that the signals in the two lines are not identical. Despite these differences, the similar signal amplitudes in D1 and D2 points to the fact that depolarization dominates the signal in D2. From this observation, we can generalize and assume that the linear polarization signal present in the photospheric lines of the spectrum of Betelgeuse is dominated by the depolarization of the continuum.

Figure 2 shows that metallic lines in the region of D1 and D2 are polarized (see Section 3), having the same shape than the LSD profiles presented in Fig 1. Figure 3 shows for the same two dates an extraction of the linearly polarized spectra around three lines (Fe i 558.7 nm, Ni i 558.8 nm, Ca i 558.9 nm) which are also known to depolarize the continuum of the solar spectrum (Gandorfer, 2000). We see that all the prominent lines visible in this region are polarized and present shapes similar to that of the LSD profiles in Fig. 1. Stenflo et al. (1983a) show that in the Sun "most spectral lines have a very small intrinsic polarization, their main effect being to depolarize the continuum polarization". For the depolarizing lines, the shape of the polarized profiles appears to be statistically identical for all absorption lines (relationship 3.11 of Stenflo et al., 1983a). Then this common shape would lead to an averaged profile in the LSD procedure of the and profiles of Betelgeuse, and would explain why we obtain a very significant signal. However in the case of Betelgeuse, the linearly polarized profiles of the lines are much more complex than in the case of the Sun and in general they do not mimic the intensity profiles. Actually, they appear to present shapes varying with time, and this complexity is due to the formation process of the polarization of the continuum (see Section 5.). Hence, depolarization of the continuum seems to be an important contribution to the linearly polarized spectrum of Betelgeuse.

### 4.2 Properties of the linearly polarized spectrum

In order to investigate the hypothesis of depolarization of the continuum as the origin of the linearly polarized spectrum of Betelgeuse, we have studied some properties of its LSD averaged profiles (as introduced in Fig. 1), making sub-masks selected for their mean depth or their mean wavelength, from a mask similar to the one described in Section 2 but including many more atomic lines (all lines with a central depth greater than 0.01 of the continuum). Depending on the required conditions, the number of lines in a given sub-mask is between 2 000 and 15 000. Changing the number of lines showed that this parameter does not significantly affect the profiles obtained using the sub-masks. We have checked the variations of the strength of the Stokes and signals and of the linear polarization with respect to depth and wavelength of the lines and found strong correlations. In the following, we concentrate only on data from Set 1, but fully compatible results have also been obtained for Set 2.

In order to investigate the dependence of the linear polarization with the depth of the lines, we have made ten sub-masks, restricting the individual spectral intensity line depths from 0.01-0.1 of the continuum (the weakest lines, i.e., those mainly formed in the lower part of the atmosphere) to 0.9-1.0 of the continuum (the strongest lines, formed higher in the atmosphere). Figure 4 shows the mean polarization profiles and the maxima of (i.e., main peak maxima, represented by dots for the blue lobe and squares for the red lobe) for the 10 mean depths: the polarization strongly increases with line depth, that is, with altitude. We checked the mean wavelength of the ten sub-masks and found that they were included in a 100 nm range, meaning that they did not differ significantly. Also, in Fig. 4 the velocities corresponding to the maxima of deviate more and more from the photospheric velocity (quoted in the Figures as Betelgeuse’s heliocentric radial velocity) with increasing line depth (i.e. height in the atmosphere). We interpret this trend as due to the expansion of the atmosphere which will be discussed again in Section 5. We obtain similar results for both lobes.

In addition to this variation of the linear polarization with line depth, we also searched for a wavelength effect. Indeed, the linear polarization of the continuum in Betelgeuse is known to increase steeply towards the blue (e.g. Clarke & Schwarz, 1984), almost as steeply as the law. Figure 5 illustrates the variations of corresponding to seven sub-masks with wavelength-bandpasses of 350-420 nm, 420-490 nm, 490-560 nm, 560-630 nm, 630-700 nm, 700-770 nm, and 770-1050 nm. The corresponding average wavelengths are given in the figure. The measurements were normalized to a common mean depth of the Stokes I profiles (that corresponding to the reddest sub-mask).

Figure 5 presents the variation of the linear polarization with wavelength accross LSD profiles, and the decrease from blue to red is confirmed. However we find that it is consistent with a law. In order to interpret this phenomenon, we explored both the diffusion and limb darkening mechanisms.

Figure 6 (left panel) presents the computation of the expected wavelength dependence of Rayleigh scattering by hydrogen atoms as (Kramers, 1924), where is the Lyman  wavelength ( nm), and its usual approximation which closely matches it (blue and green curves, respectively). However, this basic wavelength dependence of Rayleigh scattering is modified in stellar atmospheres by the chromatic variation of limb darkening. The red curve in the left plot shows the case of the Sun as studied and described by Fluri & Stenflo (1999). Therefore, due to limb darkening, the polarization of the continuum in the solar case can fall faster than the 4 power of wavelength.

In the right plot of Fig. 6, we followed Landi Degl’Innocenti & Landolfi (2004) to simulate the wavelength dependence of this depolarization signal in the lines. It is proportional to the Rayleigh coefficient , as is the continuum, but inversely proportional to the line opacity and the emission coefficients of the lines. These three combined contributions result in the green curve labelled achromatic illumination since we have not considered any wavelength dependence in the anisotropy of the light due to limb darkening. The blue curve of this right plot (Fig. 6) introduces limb darkening, again for the solar case which is used here as an example and a first approximation. We note that in both approximations, with and without chromatic limb darkening, the decrease of the depolarization amplitude in the lines is moderate. For each wavelength we have also represented on the same plot the values of the different maxima of the main peaks shown in Fig. 5 and an estimate of their errors. Superposed to these points we plot as a reference the law (red curve), that qualitatively reproduces the wavelength dependence of our measurements of Betelgeuse, hence the depolarization signal. Given that limb darkening in Betelgeuse and in the Sun are expected to be similar since the main opacity source is the same (), one thus expects a dependence of the depolarization with wavelength of the order of . A more detailed analysis of the line formation conditions in the atmosphere of Betelgeuse would likely give a better fit to this observation, but it is beyond the scope of this paper. A very interesting consequence is that such observations would provide new constraints on the limb darkening of Betelgeuse. Our present investigation shows that the difference of the wavelength dependence observed for the continuum (about law, e.g., Clarke & Schwarz, 1984) and the lines (about law) can be explained if we deal with depolarization of the continuum in the lines, taking into account the diffusion and limb darkening mechanisms.

In summary, the linear polarization observed in the atomic lines of Betelgeuse appears to be due to the depolarization of the continuum. Betelgeuse is the first star in which the depolarization of the continuum is observed in the Na i D1 line, since it is neither observed in the Sun (Stenflo et al., 2000a) nor in the Mira star $χ$ Cyg, nor in the RV Tauri star R Sct (Lèbre et al., 2015). Also, no prominent intrinsic polarization is observed for the Na iD lines. However the polarized profiles are different for D1 and for D2 both in Stokes and . As stated above, it is believed that strong polarization in D1 can only originate from depolarization of the continuum, so we can infer that in D2 we also deal with intrinsic polarization which cancels, in part, the effect of depolarization of the continuum. The polarization difference can reach about 10 whereas the intrinsic polarization of D2 observed in the polarized solar spectrum is 3 x10 (Belluzzi & Landi Degl’Innocenti, 2009).

## 5 An analytic model for mapping bright spots at the surface of Betelgeuse using its linearly polarized spectrum

The analysis of our spectra of Betelgeuse described in Section 4 shows that we are mainly observing the depolarization of the continuum. In this section we use the averaged and LSD profiles. Although they do not provide the local polarization of the continuum, the information contained in these spectra allows the inference of the surface location and relative brightness of the bright spots inducing anisotropies in the radiation field.

### 5.1 Polarization of the continuum of Betelgeuse

The observed polarization of the continuum may be explained as a combination of Rayleigh and Mie scattering (e.g. Doherty, 1986). However, in order to depolarize the continuum, the observed photospheric lines must be formed at or above the height at which the continuum itself becomes polarized. Mie polarization by molecules or dust above the photosphere is thus formed too high and can be excluded from our picture. This does not mean that there is no Mie polarization in the observed continuum, but just that Mie scattering takes place farther from the star, above where the photospheric lines have formed and have depolarized the continuum, which at this point can only carry polarization due to Rayleigh scattering. Since the average line depolarization is of a few times of the continuum flux, this is the minimum contribution of Rayleigh scattering to the total polarization of the continuum.

Scattering polarization is linear, and orthogonal to the scattering plane. At any position angle of the disk of Betelgeuse, the polarization due to scattering will be at a tangent to the local limb, orthogonally to photons coming from the deep photosphere and scattered towards us. As we change azimuth, the plane of polarization will rotate. If the star at the height of its photospheric layers were a perfect sphere and if the amount of polarization were azimuthally symmetric, the integration over the stellar disk would amount to zero and no depolarization signal would be detected. As a consequence, two possibilities appear: Betelgeuse is not a sphere, and/or its brightness is not azimuthally symmetric. Most observations and analysis conclude that Betelgeuse conserves a spherical symmetry at the height associated with our observations (e.g., observations by SPHERE, Kervella et al., 2015). We are left with the second possibility: the locally-induced polarization is not azimuthally symmetric. There are several manners in which the azimuthal symmetry may be broken, among them are: clouds over the star that only let us see some regions of the lower atmosphere; plumes of material rising high above the photospheric level with higher radiation anisotropy and higher polarization level; dark regions of the photosphere emitting less light and hence inducing less polarization; brighter, hotter regions of the photosphere emitting more light and more polarization. In the following section, we adopt the latter explanation for its simplicity and since bright spots have been observed by interferometry on Betelgeuse (e.g. Haubois et al., 2009) and on Antares (Tuthill et al., 1997). In addition, numerical simulations (Chiavassa et al., 2009, 2010, 2011) also indicate that bright convective spots are expected in the atmospheres of red supergiants.

### 5.2 A two spots hypothesis for interpreting the polarized spectrum

We assume the presence of bright spots in the deep photosphere of Betelgeuse, located at particular positions over the disk. We suppose in the following that the (Rayleigh) scattering matter is located vertically above these spots. The observed LSD profiles shown in Fig. 1 are composed of two lobes, a blue one (of positive polarity in both Stokes and in both Sets), and a red one (of negative polarity in both Stokes and , also for both Sets), with the blue lobe being in general stronger than the red one. With respect to the center of the intensity mean profile, the two lobes are located respectively at about +5 and +25 km s  for both sets. To explain this behavior we consider the simple scenario of two photospheric bright spots referred to above and in Section 4: one spot would be related to the structures that move toward the observer (blue features, spot1 hereafter), the other would be linked to structures moving away (red features, spot2 hereafter). Actually, the observed radial velocities are the combinations of those of the layer where the lines are formed in the deep photosphere and that of the layer higher in the atmosphere where the light is scattered (because of a very long rotational period (e.g., Uitenbroek et al., 1998), contributions from the rotational velocity of the star are negligible). Both layers are able to contribute to velocities that can be blueshifted or redshifted with respect to the central part of the line: ascending versus sinking elements exist in the photosphere, and we assume the outer layers of Betelgeuse are in expansion. In front of the disk the expanding shells are blueshifted, whereas beyond the limb the expanding matter is redshifted. The observed intensity profiles are consistent with a model in which radiation from the stellar photosphere is scattered and reemitted by an expanding envelope(s) of radius larger than the photosphere radius (Goldberg et al., 1975). Therefore, the blue component of the LSD polarimetric profile could be the consequence of scattering by the approaching half of the expanding envelope, while the corresponding red component could be due to the scattering of the same expanding envelope, radiated from its receding half. The expansion of the atmosphere is illustrated in Fig. 4 where the strongest lines, formed higher in the atmosphere, present larger expansion velocities.

### 5.3 An analytic model for mapping bright spots at the surface of Betelgeuse

Following this scenario, we can propose an analytical model that links the observed polarizations to the positions and relative brightness of the bright photospheric spots on the visible stellar disk. The two polar coordinates that we derive from the polarimetric data are the position angle and the projection angle . Figure 7 presents what is observed in the case of one bright spot. In the left of the figure, a cartoon of the stellar disk shows the polar coordinates for one bright spot at three locations, with the same position angle but at different projected distances to the disk center. On top of each spot the arrow shows both the orientation and the relative polarization expected from the spot alone. Around the disk the dashed arrows show the definition of the Stokes parameters as measured in the spectral lines and described in Section 5.3.1. In the right-hand part of the figure, the leftmost two plots show Stokes Q (top) and U (bottom) for the three spots as they would be seen on a resolved stellar disk. The continuum is polarized and the lines depolarize at the wavelength corresponding to the projected radial velocity. The two rightmost plots show what would be observed with Narval: the continuum polarization is set to zero and what was a depolarization signal becomes a net signal over zero. The radial velocity scale is also shifted to the of Betelgeuse.

#### 5.3.1 Getting the position angle χ from our Narval data

First, each spot will lead to an excess of polarization perpendicular to the radius joining the center of the disk of Betelgeuse with the spot, this radius forming an angle with the North-South meridian of the disk, its position angle. With respect to this North-South meridian we expect an excess of polarization with orientation , the polarization angle. Conversely, observing a polarization at position, we can infer as 90 that is, an ambiguity of . Using the classical linear polarization relations, as described by Bagnulo et al. (2009), and taking into account that with Narval (as for ESPaDOnS), positive is defined in the North-South direction and positive is defined at 45 (counterclockwise) relative to positive , the polarization angle is given by

 θ=0.5arctan(U/Q),

(if and ;+180 if and ; +90 if ).

The position angle is therefore . In the next sections we choose the option . This choice will be discussed in Section 6.2. Figure 8 shows the LSD profiles of Stokes and and for our 12 observations. From the relationships above we can determine the maxima of linear polarization and its associated , from which we can infer .

In order to map the spots at the surface of Betelgeuse we now compute another coordinate in addition to , the projection angle from the disk center, and brightness ratios among the different spots.

#### 5.3.2 Projection angle μ and brightness ratios

To infer the projection angle we make the approximation that all photospheric motions and the expansion velocity of the scattering layer are radial and combine at the polarization level to an identical velocity . We represent the velocities as in the case of an expanding atmosphere (as referred to in Section 5.2). Then the blueshifted motions are related to the visible (approaching) hemisphere of the star and redshifted signals are related to the opposite (receding) hemisphere of the star. As we explain below we assign an ad hoc value  km s. The polarization of a spot at a projection angle from the disk center will acquire a line-of-sight velocity of .

The sign of is negative or positive respectively for blueshifted or redshifted motions with respect to the of Betelgeuse (about 21 km s).

By fixing we can determine on the visible/opposite hemisphere from the position of each observed polarization lobe. Having thus determined and , we have the coordinates of the spots on the disk.

The total polarization can be approximated to be due exclusively to the intrinsic brightness of the photospheric spot, corrected by the angular dependence of Rayleigh scattering. The exact dependence of the depolarization signal on the brightness is complex. To simplify the model we assume that all other contributing factors apart from brightness and position will be identical from one spot to the other. Hence we can deduce brightness ratios among the different spots by assigning the brightest one an ad hoc brightness .

Hence

### 5.4 Mapping the bright spots

The above computations can be generalized at every single 1.8 km s velocity bin along the Stokes and profiles, enabling us to derive a map of the involved bright spots on Betelgeuse.

In the case when the peaks in and are at the same velocity, as during 2013 and in January 2014, the peak of coincides with them (e.g., as shown in Fig. 8). However in general the peaks are not aligned. The model will interpret this non-coincidence as an extended spot spanning from the position of the peak in to the position of the peak in . This is due to a combination of stellar phenomena broadening spectral lines, along with the intrinsic spectral resolution of our instrument. In order to take into account the wings of these broadened signals with our model, we must interpret not only the peaks but also the wings of these broadened profiles, implying that we must use a value of which encompasses both the physical velocities in the photosphere and the broadening of the profiles. To determine the latter, we fit a Gaussian to each of the observed lobes in and . The standard deviation of these Gaussians was found to be between 5 and 10 km s in roughly 75 % of the cases. To take into account this width we have to set a value of = 50 km s. This is not a physical value for the velocity of the photosphere, but rather only a numerical value that allows us to implement our model including the broadening of the signal by our instrument and by stellar phenomena neglected in our model. We conducted several trials and saw that the result of the mapping does not depend much on the value of . We can apply the model to all our data and map the spots at the surface of Betelgeuse. This can be seen in Fig. 9 where one map per observation date is shown. For each observation date, Fig. 9 shows two planispheres (with the same orientation on the sky): the left one shows spot1 (corresponding to blueshifted motions), the right one spot2 (corresponding to the redshifted motions). In our present interpretation, the redshifted signals corresponding to spot2 arise from the invisible hemisphere of Betelgeuse. In the case of an interpretation with the redshifted signals being of photospheric origin, the spot would have the same and coordinates but would be on the visible hemisphere. The model automatically maps the two spots inferred from the profiles and for which some derived characteristic numbers are listed at the end of Section 2 and presented in Table 1. In general two maxima are determined (as seen in Fig. 8 and Fig. 9). In this case the numbers in Table 1 are the outputs of the model. In the case of 16 and 23 October 2014, when the two maxima appeared to merge , the numbers were calculated ‘by hand’.

## 6 Variations of the linear polarization with time and evolution of the spots at the surface of Betelgeuse

### 6.1 Variations of the linear polarization and of the modeled spots

As described in Section 2 and Table 1, we observed Betelgeuse on 12 nights between November 2013 and April 2015. Table 1 gives the variations of the parameters defined above. This shows that the blue-shifted lobe (subscript 1 in Table 1) is more polarized than the red-shifted lobe (subscript 2 in Table 1). For both lobes the polarization is relatively stable during the time-span, but with some decrease. Between November 2013 and January 2014 the polarization angles of the lobes were stable, which we interpret as a stable position angle of the spots. Then and increased smoothly, corresponding to a motion of spot1 towards the South. After 2015, the motion of spot1 accelerates and spot1 reaches the south pole in April. During this same period, spot2 moves the same angle towards the north. These motions of the spots are illustrated in Fig. 9.

Changing shapes of the spots are seen on the maps. The two main sources of errors that we consider as intrinsic to the model are the already-mentioned PSF broadening of the profiles (with standard deviation of 5-10 km s) and the photon noise over the profiles (measured to be of the order of of the total intensity over the LSD profiles). It appears, in the end, that our spatial resolution at the surface of the star is roughly square degrees. In Fig. 9 we have made the size of our markers equal to this error.

With respect to linear polarization measurements in the continuum, our geometric model allows us to disentangle several spots with different and to infer their positions on the disk of Betelgeuse. However, it does not provide the local polarization of the continuum since we do not know the depolarization coefficient of the lines. An order of magnitude can be derived for averaged lines. Published polarization measurements of the continuum of Betelgeuse are dominated by Mie polarization (e.g. Hayes, 1984; Clarke & Schwarz, 1984; Kervella et al., 2015). Predictions for the Rayleigh contribution (Doherty, 1986; Josselin et al., 2015) give about  % at maximum in the blue. The order of magnitude of the depolarization factor by the lines is therefore about 0.1 in the blue. We show in Section 3.2 that the depolarization varies with wavelength.

### 6.2 Comparison with quasi-simultaneous observations of VLTI/PIONIER and with older polarimetric observations

Betelgeuse has been surveyed once per year with VLTI/PIONIER since 2012 (Montargès, 2014; Montargès et al., 2015) and a large hot spot has been resolved near the limb. This hot spot appears very large compared to those which have been detected so far on the disk of Betelgeuse (e.g., Haubois et al., 2009). Two PIONIER observations (Montargès et al., 2015) are quasi-simultaneous with our Stokes survey of Betelgeuse, namely on 11 January 2014 and 21 November 2014 (see Table 1 for very near Narval’s dates). The data concerning the three first PIONIER observations (January 2012, February 2013, January 2014) have been satisfactorily fitted with a model of disk with limb darkening (hereafter LDD model) and a Gaussian hot spot (Montargès, 2014). The spot is centered near the limb of Betelgeuse and its FWHM is about one radius of the star. The center of the spot has ranging from about 80 to about 110 which corresponds to the E-SE part of Betelgeuse’s limb. For PIONIER data of 21 November 2014, the situation appears different and the fit of the model LDD Gaussian spot is worse than for the previous dates. In fact, a significantly better fit to this observation is obtained for a two-spot model. In this case the spot in the NW quadrant is the strongest one and has a Gaussian shape similar to that of the previous observations; the spot in the SE quadrant has a negligible diameter (Montargès et al., 2015). Comparing the results from Table 1 (with ) and from the PIONIER observations, we see that the position of our spot1 in Set1 is similar to that of the large hot spot in January 2012, February 2013 and January 2014. For the November 2014 positions to be consistent, we have to choose for the polarimetric observation the location with 180 away the value reported in Table 1. Figure 10 shows the PIONIER intensity maps presented by Montargès et al. (2015) and the maps obtained from polarimetry at the dates of quasi-simultaneous observations. This comparison suggests that the large hot spot observed by PIONIER in 2012, 2013 and January 2014 is responsible for the stable polarization of our blue lobe from November 2013 to January 2014 and coincides with spot1. The evolution of the PIONIER hot spot observed in November 2014 could have the same origin as the evolution of spot1 observed from October 2014 to April 2015.

Our spot2 has not been detected with PIONIER. Our hypothesis that it is a spot located on the hidden hemisphere is one possible explanation. Also, Dupree (2010) proposed that the hot spots on Betelgeuse are composed of photospheric and chromospheric features. Dupree also proposed that these chromospheric features, are not detected by interferometry working in the infrared, as was the hot spot imaged by the HST/FOC (Gilliland & Dupree, 1996; Uitenbroek et al., 1998). Remarkably, the polarization angle of our redshifted lobe (linked to spot2) in 2015 is consistent with that (about 150) corresponding to the maximum of polarization observed by Hayes (1984) and Nordsieck et al. (1994), as reported by Uitenbroek et al. (1998).

## 7 Conclusions

We have discovered that most of the photospheric lines of Betelgeuse present a linearly polarized spectrum. Together with the strong polarization of the Na i D1 line, and including the polarization laws with both wavelength and depth, we show that this linear polarization originates from the depolarization of the continuum. Taking into account the discovery of linearly polarized spectra in the Mira star  Cyg and in the RV Tauri star R Sct (Lèbre et al., 2015), in which intrinsic polarization of resonance lines (similar to the Sun) is detected, all these observations of cool evolved stars point to a kind of second stellar spectrum.

We interpret our results for Betelgeuse as being due to anisotropies of the radiation field induced by hot spots at the surface and Rayleigh scattering in the atmospheric envelope. We then propose an analytical model to interpret the observed polarization and we infer the presence of two hot spots and their positions on Betelgeuse. We show that applying the model to each velocity bin along the Stokes and profiles allows us to derive a map of the involved bright spots on Betelgeuse. We compare our observations to quasi-simultaneous VLTI/PIONIER observations (Montargès et al., 2015) and older measurements of the polarization of the continuum. Our linear polarimetric observations suggest that we detected a hot spot seen by quasi-simultaneous PIONIER observations, and another one whose associated polarization angle implies a position near where the largest percentage of polarization of the continuum occurs (as reported by Uitenbroek et al., 1998).

###### Acknowledgements.
We thank the TBL team for providing service observing with Narval. We acknowledge the use of the database VALD (Vienna, Austria) and financial support from the "Programme National de Physique Stellaire" (PNPS) of CNRS/INSU, France. R.K.-A. acknowledges support by the Bulgarian-French mobility contract DRILA 01/3. GAW acknowledges Discovery Grant support from the Natural Sciences and Engineering Research Council (NSERC) of Canada.

## References

• Aurière et al. (2010) Aurière, M., Donati, J.-F., Konstantinova-Antova, R. et al. 2010, A&A, 516, L2
• Aurière et al. (2009) Aurière, M., Wade, G.A., Konstantinova-Antova, R. et al. 2009, A&A, 504, 231
• Bagnulo et al. (2009) Bagnulo, S., Landolfi M., Landstreet J.D. et al. 2009, PASP, 121, 993
• Bedecarrax et al. (2013) Bedecarrax, I., Petit, P., Aurière, M. et al. 2013, EAS, 60, 161
• Belluzzi & Landi Degl’Innocenti (2009) Belluzzi, L. & Landi Degl’Innocenti 2009, A&A, 495, 577
• Boyle et al. (1986) Boyle, R.P., Aspin, C., Coyne, G.V. & McLean, I.S. 1986, A&A, 164, 310
• Chiavassa et al. (2009) Chiavassa, A., Plez B., Josselin E., Freytag B. 2009, A&A, 506, 1351
• Chiavassa et al. (2010) Chiavassa, A., Haubois X., Young J.S. et al. 2010, A&A, 515, 12
• Chiavassa et al. (2011) Chiavassa, A., Pasquato E., Jorissen A. et al. 2011, A&A, 528, A120
• Clarke & Schwarz (1984) Clarke, D. & Schwarz, H.E. 1984, A&A, 132, 375
• Doherty (1986) Doherty, L.R. 1986, APJ, 307, 261
• Donati (2003) Donati, J.-F. 2003, ASPC, 307, 41D
• Donati et al. (1997) Donati, J.-F., Semel, M., Carter, B.D. et al. 1997, MNRAS, 291, 658
• Donati et al. (2006) Donati, J.-F., Catala, C., Landstreet, J., Petit, P. 2006, in Casini R., Lites B., eds, Solar Polarization Workshop n4, Vol.358 of ASPC series, 362
• Dupree (2010) Dupree, A. K. 2010, Proceedings IAU Symposium No. 273, D.P. Choudhary, K.G. Strassmeier, eds., 188
• Fabas et al. (2011) Fabas, N., Lèbre, A. & Gillet, D. 2011, A&A, 535, A12
• Fabas (2011) Fabas, N. 2011, Ph D Thesis, Université de Montpellier II
• Fluri & Stenflo (1999) Fluri, D. M. & Stenflo, J.O. 1999, A&A, 341, 902
• Fluri & Stenflo (2003) Fluri, D. M. & Stenflo, J.O. 2003, A&A, 398, 763
• Gandorfer (2000) Gandorfer, A. 2000, The Second Solar Spectrum: A High Resolution Polarimetric Survey of Scattering Polarization at the Solar Limb in Graphical Representation, Vol I: 4625 Åto 6995 Å, ISBN: 3-7281-2764-7
• Gilliland & Dupree (1996) Gilliland, R.L., Dupree A.K. 1996, ApJ, 463, L29
• Goldberg et al. (1975) Goldberg, L., Ramsey, L., Testerman, L. & Carbon, D. 1975, ApJ, 199, 427
• Hayes (1984) Hayes, D.P. 1984, APJSS, 55, 179
• Haubois et al. (2009) Haubois, X., Perrin, G., Lacour, S. et al. 2009, A&A, 508, 923
• Josselin et al. (2015) Josselin E., Lambert J., Aurière M., Petit P. & Ryde N. 2015, in proceedings of IAUS 305
• Josselin & Plez (2007) Josselin, E., Plez, B. 2007, A&A, 469, 671
• Kervella et al. (2015) Kervella, P., Lagadec E., Montargès M. et al. 2015, A&A, 585, A28
• Konstantinova-Antova et al. (2010) Konstantinova-Antova, R., Aurière, M., Charbonnel, C. et al. 2010, A&A, 524, 57
• Konstantinova-Antova et al. (2013) Konstantinova-Antova, R., Aurière, M., Charbonnel, C. et al. 2013, BlgAJ, 19, 15
• Kramers (1924) Kramers, H. A. 1924, Natur., 113, 673
• Kupka et al. (1999) Kupka, F.G., Piskunov, N.E., Ryabchikova, T.A. et al. 1999, A&AS, 138, 119
• Lambert et al. (1984) Lambert, D.L., Brown, J.A., Hinkle, K.H., Johnson, H.R. 1984, ApJ, 284, 223
• Landi Degl’Innocenti & Landolfi (2004) Landi Degl’Innocenti, E. & Landolfi, M. 2004, Polarization in Spectral Lines, ASSL, 307, (Dordrecht:Kluwer)
• Lèbre et al. (2014) Lèbre, A., Aurière, M., Fabas, N. et al. 2014, A&A, 561, A85
• Lèbre et al. (2015) Lèbre, A., Aurière, M., Fabas, N. et al. 2015, in proceedings of IAUS 305, in press
• López Ariste et al. (2011) López Ariste, A., Asensio Ramos, A. & González Fernández, C. 2011, A&A, 527, A120
• Magalhaes & Nordsieck (2000) Magalhaes, A. M. & Nordsieck, K.H. 2000, proceedings IAUS 177, 433
• McLean & Coyne (1978) McLean, I.S. & Coyne, G.V. 1978, ApJ, 226, L145
• Montargès (2014) Montargès, M. 2014, PhD thesis, Pierre et Marie Curie University
• Montargès et al. (2015) Montargès, M. , Kervella P., Perrin G. et al. 2015, A&A, 588, A130
• Nordsieck et al. (1994) Nordsieck, K.H., Rao, L., Meade, M.R. et al. 1994, BAAS, 26, 864
• Paletou (2012) Paletou, F. 2012, A&A, 544, 4
• Petit et al. (2013) Petit, P., Aurière, M., Konstantinova-Antova, R. et al. 2013, LNP, 857, 231
• Rosén et al. (2015) Rosén, L., Kochukhov, O. & Wade, G.A. 2015, ApJ, 805, 169R
• Rosén et al. (2013) Rosén, L., Kochukhov, O. & Wade, G.A. 2013, MNRAS, 4361, 10R
• Sabin et al. (2015) Sabin, L., Wade, G. & Lèbre, A. 2015, MNRAS 446, 1988
• Schwarz & Clarke (1984) Schwarz, H.E. & Clarke, D. 1984, A&A, 132, 370
• Semel et al. (2009) Semel, M., Ramírez Vélez, J.C., Martínez Gonzáles, M.J. et al. 2009, A&A, 504, 1003
• Stenflo (1997) Stenflo, J.O. 1997, A&A, 324 344
• Stenflo (2009) Stenflo, J.O. 2009, in ”Cosmic Magnetic Fields: From Planets, to Stars and Galaxies”, IAUS 259, 211, K.G. Strassmeier, A.G. Kosovichev, J.E. Beckman, eds.
• Stenflo et al. (2000a) Stenflo, J.O., Gandorfer, A. & Keller, C.U. 2000a, A&A, 355, 781
• Stenflo et al. (2000b) Stenflo, J.O., Keller, C.U. & Gandorfer, A. 2000b, A&A, 355, 789
• Stenflo et al. (1983a) Stenflo, J.O., Twerenbold, D. & Harvey, J.W. 1983a, A&ASS, 52, 161
• Stenflo et al. (1983b) Stenflo, J.O., Twerenbold, D., Harvey, J.W. & Brault J.W. 1983b, A&ASS, 54, 505
• Tuthill et al. (1997) Tuthill P.G., Haniff C.A., Baldwin J.E. 1997, MNRAS, 285, 529
• Uitenbroek et al. (1998) Uitenbroek, H., Dupree, A.K., Gilliland, R.L. 1998, AJ, 116, 2501
• Wade et al. (2000) Wade, G.A., Donati, J.-F., Landstreet, J.D. & Shorlin, S.L.S. 2000, MNRAS, 313, 851
• Wilson et al. (1997) Wilson, R.W., Dhillon, V.S. & Haniff, C.A. 1997, MNRAS, 291, 819
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters