A new method of measuring center-of-mass velocities of radially pulsating stars from high-resolution spectroscopy

We present a radial velocity analysis of 20 solar neighborhood RR Lyrae and 3 Population II Cepheids variables. We obtained high-resolution, moderate-to-high signal-to-noise ratio spectra for most stars and obtained spectra were covering different pulsation phases for each star. To estimate the gamma (center-of-mass) velocities of the program stars, we use two independent methods. The first, `classic' method is based on RR Lyrae radial velocity curve templates. The second method is based on the analysis of absorption line profile asymmetry to determine both the pulsational and the gamma velocities. This second method is based on the Least Squares Deconvolution (LSD) technique applied to analyze the line asymmetry that occurs in the spectra. We obtain measurements of the pulsation component of the radial velocity with an accuracy of $\pm$ 3.5 km s$^{-1}$. The gamma velocity was determined with an accuracy $\pm$ 10 km s$^{-1}$, even for those stars having a small number of spectra. The main advantage of this method is the possibility to get the estimation of gamma velocity even from one spectroscopic observation with uncertain pulsation phase. A detailed investigation of the LSD profile asymmetry shows that the projection factor $p$ varies as a function of the pulsation phase -- this is a key parameter which converts observed spectral line radial velocity variations into photospheric pulsation velocities. As a byproduct of our study, we present 41 densely-spaced synthetic grids of LSD profile bisectors that are based on atmospheric models of RR Lyr covering all pulsation phases.Comment: 17 pages, 16 figures, accepted for publication in MNRAS; doi:10.1093/mnras/stx294

Identification of red supergiants in nearby galaxies with mid-IR photometry

The role of episodic mass loss in massive star evolution is one of the most important open questions of current stellar evolution theory. Episodic mass loss produces dust and therefore causes evolved massive stars to be very luminous in the mid-infrared and dim at optical wavelengths. We aim to increase the number of investigated luminous mid-IR sources to shed light on the late stages of these objects. To achieve this we employed mid-IR selection criteria to identity dusty evolved massive stars in two nearby galaxies. The method is based on mid-IR colors, using 3.6 {\mu}m and 4.5 {\mu}m photometry from archival Spitzer Space Telescope images of nearby galaxies and J-band photometry from 2MASS. We applied our criteria to two nearby star-forming dwarf irregular galaxies, Sextans A and IC 1613, selecting eight targets, which we followed up with spectroscopy. Our spectral classification and analysis yielded the discovery of two M-type supergiants in IC 1613, three K-type supergiants and one candidate F-type giant in Sextans A, and two foreground M giants. We show that the proposed criteria provide an independent way for identifying dusty evolved massive stars, that can be extended to all nearby galaxies with available Spitzer/IRAC images at 3.6 {\mu}m and 4.5 {\mu}m.Comment: 8 pages, 4 figures, A&A in pres

Chemical abundances of solar neighborhood RR Lyrae stars

We have analysed a sample of 18 RR Lyrae stars (17 fundamental-mode - RRab - and one first overtone - RRc) and three Population II Cepheids (two BL Her stars and one W Vir star), for which high-resolution (R $\ge$30000), high signal-to-noise (S/N$\ge$30) spectra were obtained with either SARG at the Telescopio Nazionale Galileo (La Palma, Spain) or UVES at the ESO Very Large Telescope (Paranal, Chile). Archival data were also analyzed for a few stars, sampling $\gtrsim$3 phases for each star. We obtained atmospheric parameters (T$_{\rm{eff}}$, log$g$, v$_{\rm{t}}$, and [M/H]) and abundances of several iron-peak and $\alpha$-elements (Fe, Cr, Ni, Mg, Ca, Si, and Ti) for different pulsational phases, obtaining $$=+0.31$\pm$0.19 dex over the entire sample covering -2.2$<$[Fe/H]$<$-1.1 dex. We find that silicon is indeed extremely sensitive to the phase, as reported by previous authors, and cannot be reliably determined. Apart from this, metallicities and abundance ratios are consistently determined, regardless of the phase, within 0.10-0.15 dex, although caution should be used in the range $0\lesssim\phi\lesssim0.15$. Our results agree with literature determinations for both variable and non-variable field stars, obtained with very different methods, including low and high-resolution spectroscopy. W Vir and BL Her stars, at least in the sampled phases, appear indistinguishable from RRab from the spectroscopic analysis point of view. Our large sample, covering all pulsation phases, confirms that chemical abundances can be obtained for RR Lyrae with the classical EW-based technique and static model atmospheres, even rather close to the shock phases.Comment: 16 pages, 8 figures, 7 tables (of which 1 electronic), accepted for publication on MNRA

Properties of luminous red supergiant stars in the Magellanic Clouds

There is evidence that some red supergiants (RSGs) experience short lived phases of extreme mass loss, producing copious amounts of dust. These episodic outburst phases help to strip the hydrogen envelope of evolved massive stars, drastically affecting their evolution. However, to date, the observational data of episodic mass loss is limited. This paper aims to derive surface properties of a spectroscopic sample of fourteen dusty sources in the Magellanic Clouds using the Baade telescope. These properties may be used for future spectral energy distribution fitting studies to measure the mass loss rates from present circumstellar dust expelled from the star through outbursts. We apply MARCS models to obtain the effective temperature ($T_{\rm eff}$) and extinction ($A_V$) from the optical TiO bands. We use a $\chi^2$ routine to determine the best fit model to the obtained spectra. We compute the $T_{\rm eff}$ using empirical photometric relations and compare this to our modelled $T_{\rm eff}$. We have identified a new yellow supergiant and spectroscopically confirmed eight new RSGs and one bright giant in the Magellanic Clouds. Additionally, we observed a supergiant B[e] star and found that the spectral type has changed compared to previous classifications, confirming that the spectral type is variable over decades. For the RSGs, we obtained the surface and global properties, as well as the extinction $A_V$. Our method has picked up eight new, luminous RSGs. Despite selecting dusty RSGs, we find values for $A_V$ that are not as high as expected given the circumstellar extinction of these evolved stars. The most remarkable object from the sample, LMC3, is an extremely massive and luminous evolved massive star and may be grouped amongst the largest and most luminous RSGs known in the Large Magellanic Cloud (log(L$_*$/L$_{\odot})\sim$5.5 and $R = 1400 \,\ \textrm R_{\odot}$).Comment: Accepted for publication in A&A. 17 pages, 14 figures, 4 table

The IACOB project VIII. Searching for empirical signatures of binarity in fast-rotating O-type stars

The empirical distribution of projected rotational velocities (vsini) in massive O-type stars is characterized by a dominant slow velocity component and a tail of fast rotators. Binary interaction has been proposed to play a dominant role in the formation of this tail. We perform a complete and homogeneous search for empirical signatures of binarity in a sample of 54 fast-rotating stars with the aim of evaluating this hypothesis. This working sample has been extracted from a larger sample of 415 Galactic O-type stars which covers the full range of vsini values. We use new and archival multi-epoch spectra in order to detect spectroscopic binary systems. We complement this information with Gaia proper motions and TESS photometric data to aid in the identification of runaway stars and eclipsing binaries, respectively. The identified fraction of single-lined spectroscopic binary (SB1) systems and apparently single stars among the fast-rotating sample is $\sim$18% and $\sim$70%, respectively. When comparing these percentages with those corresponding to the slow-rotating sample we find that our sample of fast rotators is characterized by a slightly larger percentage of SB1 systems ($\sim$18% vs. $\sim$13%) and a considerably smaller fraction of clearly detected SB2 systems (8% vs. 33%). Overall, there seems to be a clear deficit of spectroscopic binaries (SB1+SB2) among fast-rotating O-type stars ($\sim$26% vs. $\sim$46%). On the contrary, the fraction of runaway stars is significantly higher in the fast-rotating domain ($\sim$33-50%) than among those stars with vsini < 200 km/s. Lastly, almost 65% of the apparently single fast-rotating stars are runaways. Our empirical results seem to be in good agreement with the idea that the tail of fast-rotating O-type stars (with vsini > 200 km/s) is mostly populated by post-interaction binary products.Comment: 33 pages, 16 figures, accepted for publication in "Astronomy and Astrophysics

Method of LSD profile asymmetry for estimating the center of mass velocities of pulsating stars

We present radial velocity analysis for 20 solar neighborhood RR Lyrae and 3 Population II Cepheids. High-resolution spectra were observed with either TNG/SARG or VLT/UVES over varying phases. To estimate the center of mass (barycentric) velocities of the program stars, we utilized two independent methods. First, the 'classic' method was employed, which is based on RR Lyrae radial velocity curve templates. Second, we provide the new method that used absorption line profile asymmetry to determine both the pulsation and the barycentric velocities even with a low number of high-resolution spectra and in cases where the phase of the observations is uncertain. This new method is based on a least squares deconvolution (LSD) of the line profiles in order to an- alyze line asymmetry that occurs in the spectra of pulsating stars. By applying this method to our sample stars we attain accurate measurements (+- 2 kms^-1) of the pulsation component of the radial velocity. This results in determination of the barycentric velocity to within 5 kms^-1 even with a low number of high- resolution spectra. A detailed investigation of LSD profile asymmetry shows the variable nature of the project factor at different pulsation phases, which should be taken into account in the detailed spectroscopic analysis of pulsating stars

Investigating episodic mass loss in evolved massive stars: II. Physical properties of red supergiants at subsolar metallicity

Mass loss during the red supergiant (RSG) phase plays a crucial role in the evolution of an intermediate massive star, however, the underlying mechanism remains unknown. We aim to increase the sample of well-characterized RSGs at subsolar metallicity, by deriving the physical properties of 127 RSGs in nine nearby, southern galaxies presented by Bonanos et al. For each RSG, we provide spectral types and used \textsc{marcs} atmospheric models to measure stellar properties from their optical spectra, such as the effective temperature, extinction, and radial velocity. By fitting the spectral energy distribution, we obtained the stellar luminosity and radius for 97 RSGs, finding $\sim 50\%$ with log$(L/ \rm L_{\odot}) \geq 5.0$ and 6 RSGs with $R \gtrsim 1400 \,\ \rm R_{\odot}$. We also find a correlation between the stellar luminosity and mid-IR excess of 33 dusty, variable sources. Three of these dusty RSGs have luminosities exceeding the revised Humphreys-Davidson limit. We then derive a metallicity-dependent $J-K_s$ color versus temperature relation from synthetic photometry and two new empirical $J-K_s$ color versus temperature relations calibrated on literature TiO and $J$-band temperatures. To scale our derived, cool TiO temperatures to values in agreement with the evolutionary tracks, we derive two linear scaling relations calibrated on $J$-band and $i$-band temperatures. We find that the TiO temperatures are more discrepant as a function of the mass-loss rate and discuss future prospects of the TiO bands as a mass-loss probe. Finally, we speculate that 3 hot, dusty RSGs may have experienced a recent mass ejection ($12\%$ of the K-type sample) and indicate them as candidate Levesque-Massey variables.Comment: 22 pages, 18 figures, submitted to A&