Search

Red Supergiants (RSGs) are cool, evolved massive stars in their final evolutionary stage before exploding as a supernova. However, the evolution and fate of the most luminous RSGs is uncertain. Observational evidence for luminous warm, post-RSG objects and the lack of luminous RSGs as supernova progenitors suggest a blueward evolution. In this work, we analyze WOH G64, considered since the 1980s as the most extreme RSG in the Large Magellanic Cloud in terms of its size, luminosity, and mass-loss rate. Time-series photometry over the last 30 years reveals a sudden, yet smooth change from semi-regular to irregular variability in 2014. Multi-epoch optical spectroscopy confirms the transition, as WOH~G64 now exhibits properties of a B[e] star in the optical, and warm-star features in the near-infrared. We discovered that WOH G64 is a rare, massive symbiotic system formed by a RSG, which transitioned to a Yellow Hypergiant, and a B-star companion. The dramatic transition can be explained by: a) binary interactions partially stripping the envelope, b) the return of WOH G64 to a quiescent state after an outstanding eruption exceeding 30 years, and c) the expulsion of its outer layers due to a pre-SN superwind phase, indicating its imminent explosion. WOH~G64 offers a unique opportunity to witness stellar evolution in real-time, providing crucial clues for the final phases of massive stars and their resulting supernovae., Comment: 19 pages, 15 figures. Submitted to a high-impact journal. This paper is under press embargo, contact the corresponding author for details

The rate and mechanism of mass loss of red supergiants (RSGs) remain poorly understood, especially at low metallicities. Motivated by the new empirical prescription by Yang et al. 2023, based on the largest and most complete sample in the Small Magellanic Cloud, we investigate the impact of different popular and recent RSG mass-loss prescriptions that span a range of RSG mass-loss rates on the evolution and observable properties of single massive stars. Our results show that higher mass-loss rates result in earlier envelope stripping and shorter RSG lifetimes, particularly for the more luminous stars, leading to a steeper luminosity function and predicting hotter final positions for the SN progenitors. None of the considered mass-loss prescriptions is fully consistent with all observational constraints, highlighting ongoing uncertainties in deriving and modeling RSGs mass loss. The mass-loss rates suggested by Kee et al. predict rapid envelope stripping, inconsistent with the observed population of luminous RSGs and SN progenitor detections, while the models implementing the commonly used de Jager et al. and the recent Beasor et al. prescriptions overestimate the number of luminous RSGs. While the increased mass-loss rates for luminous RSGs predicted by Yang et al. lead to better agreement with the observed RSG luminosity function, naturally reproducing the updated Humphreys-Davidson limit, they also produce luminous yellow supergiant progenitors not detected in nearby supernovae. We also estimate that binary interactions tend to slightly increase the formation of luminous RSGs due to mass accretion or merging. Our study examines the impact of RSG mass loss during the late stages of massive stars, highlighting the significance of using comprehensive observational data, exploring the uncertainties involved, and considering the effects of binary-induced or episodic mass loss., Comment: submitted in A&A

This study delves into [W60] B90, one of the most luminous and extreme Red Supergiants (RSGs) in the Large Magellanic Cloud (LMC), aiming to search for evidence of episodic mass loss. Our discovery of a bar-like nebular structure at 1 pc, reminiscent of the bar around Betelgeuse, raised the question of whether [W60] B90 also has a bow shock. We collected and analyzed proper motion data from Gaia, as well as new multi-epoch spectroscopic and imaging data, and archival time-series photometry in the optical and mid-infrared. We found [W60] B90 to be a walkaway star, with a supersonic peculiar velocity in the direction of the bar. We detected shocked emission between the bar and the star, based on the [S II]/H$\alpha$ > 0.4 criterion, providing strong evidence for a bow shock. The 30-year optical light curve revealed semi-regular variability, showing three similar dimming events with $\Delta V \sim 1$ mag, a recurrence of $\sim$12 yr, and a rise time of 400 d. We found the mid-IR light curve to vary by 0.51 mag and 0.37 mag in the WISE1 and WISE2 bands, respectively, and by 0.42 mag and 0.25 mag during the last dimming event. During this event, optical spectroscopy revealed spectral variability (M3I to M4I), a correlation between the $T_{\rm eff}$ and the brightness, increased extinction, and, after the minimum, spectral features incompatible with the models. We also found a difference of >300 K between the $T_{\rm eff}$ measured from the TiO bands in the optical and the atomic lines from our $J$-band spectroscopy. We inferred that [W60] B90 is a more massive analog of Betelgeuse in the LMC and the first extragalactic single RSG with a suspected bow shock. Its high luminosity $\log(L/L_{\odot})=5.32$ dex, mass-loss rate, and mid-IR variability compared to other RSGs in the LMC, indicate that it is in an unstable evolutionary state undergoing episodes of mass loss., Comment: 28 pages, 14 figures, submitted to A&A

We present results of the NELIOTA campaign for lunar impact flashes observed with the 1.2 m Kryoneri telescope. From August 2019 to August 2023, we report 113 validated and 70 suspected flashes. For the validated flashes, we calculate the physical parameters of the corresponding projectiles, the temperatures developed during the impacts, and the expected crater sizes. For the multiframe flashes we present light curves and thermal evolution plots. Using the whole sample of NELIOTA that encompasses 192 validated flashes in total from 2017, the statistics of the physical parameters of the meteoroids, the peak temperatures of the impacts and the expected crater sizes has been updated. Using this large sample, empirical relations correlating the luminous energies per photometric band were derived and used to roughly estimate the parameters of 92 suspected flashes of the NELIOTA archive. For a typical value of the luminous efficiency, we found that more than the 75% of the impacting meteoroids have masses between 1-200 g, radii between 0.5-3 cm and produced craters up to 3.5 m. 85% of the peak temperatures of the impacts range between 2000 and 4500 K. Statistics regarding the magnitude decline and the cooling rates of the multiframe flashes are also presented. The recalculation of the appearance frequency of meteoroids, lying within the aforementioned ranges of physical parameters, on the Moon yields that the total lunar surface is bombarded with 7.4 sporadic meteoroids/ hour and up to 12.6 meteoroids/hour when the Earth-Moon system passes through a strong meteoroid stream. By extrapolating these rates on Earth, the respective rates for various distances from its surface are calculated and used to estimate the probability of an impact of a meteoroid with a hypothetical infrastructure on the Moon, or with a satellite orbiting Earth for various impact surfaces and duration times of the missions., Comment: 32 pages, accepted for publication in A&A

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&A

The high mass-loss rates of red supergiants (RSGs) drastically affect their evolution and final fate, but their mass-loss mechanism remains poorly understood. Various empirical prescriptions scaled with luminosity have been derived in the literature, yielding results with a dispersion of 2-3 orders of magnitude. We determine an accurate mass-loss rate relation with luminosity and other parameters using a large, clean sample of RSGs and explain the discrepancy between previous works. We assembled a sample of 2,219 RSG candidates in the Large Magellanic Cloud, with ultraviolet to mid-infrared photometry in up to 49 filters. We determined the luminosity of each RSG by integrating the spectral energy distribution and the mass-loss rate using the radiative transfer code DUSTY. Our derived RSG mass-loss rates range from $10^{-9} M_\odot$ yr$^{-1}$ to $10^{-5} M_\odot$ yr$^{-1}$, mainly depending on the luminosity. The average mass-loss rate is $9.3\times 10^{-7} M_\odot$ yr$^{-1}$ for $\log{(L/L_\odot)}>4$. We established a mass-loss rate relation as a function of luminosity and effective temperature. Furthermore, we found a turning point in the relation of mass-loss rate versus luminosity relation at approximately $\log{(L/L_\odot)} = 4.4$, indicating enhanced rates beyond this limit. We show that this enhancement correlates with photometric variability. Moreover, we compared our results with prescriptions from the literature, finding an agreement with works assuming steady-state winds. Additionally, we examined the effect of different assumptions on our models and found that radiatively driven winds result in mass-loss rates higher by 2-3 orders of magnitude, which are unrealistically high for RSGs. Finally, we found that 21% of our sample constitute current binary candidates. This has a minor effect on our mass-loss relation., Comment: 17 pages, 15 figures, 4 tables; Accepted for publication in A&A

The ASSESS project aims to determine the role of episodic mass-loss in the evolution of massive stars. As a first step, we construct a catalog of spectroscopically identified dusty, evolved massive stars in ten southern galaxies for which Spitzer point-source catalogs are available. We conducted multi-object spectroscopy of dusty massive star candidates in these galaxies (spanning Z = 0.06-1.6 Zo) using the VLT. We obtained 763 spectra in WLM, NGC 55, NGC 247, NGC 253, NGC 300, NGC 1313, NGC 3109, Sextans A, M83 and NGC 7793. The targets were selected using their Spitzer photometry, by prioritizing targets with a strong infrared excess. We determined a spectral classification for each target. Additionally, we used archival images from the HST to provide a visual classification for 80 targets, as a star, cluster, or galaxy. We provide a catalog of 541 spectroscopically classified sources including 185 massive stars, of which 154 are newly classified massive stars. The catalog contains 129 red supergiants, 27 blue supergiants, 10 yellow supergiants, four luminous blue variable candidates, seven supergiant B[e] stars and eight emission line objects. Evidence for circumstellar dust is found in 24% of these massive stars, based on their infrared colors. We report a success rate of 28% for identifying massive stars among our observed spectra, while the average success rate of our priority system in selecting evolved massive stars was 36%. Additionally, the catalog contains 21 background galaxies (including AGN and quasars), 10 carbon stars and 99 HII regions. We measured the line ratios [NII]/Ha and [SII]/Ha for 76 HII regions and 36 other spectra with nebular emission-lines, thereby identifying eight sources with shocked emission. We present the largest catalog of evolved massive stars and in particular of red supergiants in nearby galaxies at low Z beyond the Local Group., Comment: 47 pages, 7 figures, 9 tables, 5 appendices. Accepted for publication in A&A

Episodic mass loss is not understood theoretically, neither accounted for in state-of-the-art models of stellar evolution, which has far-reaching consequences for many areas of astronomy. We introduce the ERC-funded ASSESS project (2018-2024), which aims to determine whether episodic mass loss is a dominant process in the evolution of the most massive stars, by conducting the first extensive, multi-wavelength survey of evolved massive stars in the nearby Universe. It hinges on the fact that mass-losing stars form dust and are bright in the mid-infrared. We aim to derive physical parameters of $\sim$1000 dusty, evolved massive stars in $\sim$25 nearby galaxies and estimate the amount of ejected mass, which will constrain evolutionary models, and quantify the duration and frequency of episodic mass loss as a function of metallicity. The approach involves applying machine-learning algorithms to select dusty, luminous targets from existing multi-band photometry of nearby galaxies. We present the first results of the project, including the machine-learning methodology for target selection and results from our spectroscopic observations so far. The emerging trend for the ubiquity of episodic mass loss, if confirmed, will be key to understanding the explosive early Universe and will have profound consequences for low-metallicity stars, reionization, and the chemical evolution of galaxies., Comment: To appear in the proceedings of IAU Symposium 361: Massive stars near and far; 6 pages, 1 figure

Mass loss is one of the key parameters that determine stellar evolution. Despite the progress we have achieved over the last decades we still cannot match the observational derived values with theoretical predictions. Even worse, there are certain phases, such as the B[e] supergiants (B[e]SGs) and the Luminous Blue Variables (LBVs), where significant mass is lost through episodic or outburst activity. This leads to various structures around them that permit dust formation, making these objects bright IR sources. The ASSESS project aims to determine the role of episodic mass in the evolution of massive stars, by examining large numbers of cool and hot objects (such as B[e]SGs/LBVs). For this, we initiated a large observing campaign to obtain spectroscopic data for $\sim$1000 IR selected sources in 27 nearby galaxies. Within this project we successfully identified 7 B[e] supergiants (one candidate) and 4 Luminous Blue Variables of which 6 and 2, respectively, are new discoveries. We used spectroscopic, photometric, and light curve information to better constrain the nature of the reported objects. We particularly note the presence of B[e]SGs at metallicity environments as low as 0.14 Z$_{\odot}$., Comment: 23 pages, 10 figures, 4 tables; published in the Special Issue of Galaxies "Theory and Observation of Active B-type Stars"

We present the optical spectroscopic evolution of SN~2023ixf seen in sub-night cadence spectra from 1.18 to 14 days after explosion. We identify high-ionization emission features, signatures of interaction with material surrounding the progenitor star, that fade over the first 7 days, with rapid evolution between spectra observed within the same night. We compare the emission lines present and their relative strength to those of other supernovae with early interaction, finding a close match to SN~2020pni and SN~2017ahn in the first spectrum and SN~2014G at later epochs. To physically interpret our observations we compare them to CMFGEN models with confined, dense circumstellar material around a red supergiant progenitor from the literature. We find that very few models reproduce the blended \NC{} emission lines observed in the first few spectra and their rapid disappearance thereafter, making this a unique diagnostic. From the best models, we find a mass-loss rate of $10^{-3}-10^{-2}$ \mlunit{}, which far exceeds the mass-loss rate for any steady wind, especially for a red supergiant in the initial mass range of the detected progenitor. These mass-loss rates are, however, similar to rates inferred for other supernovae with early circumstellar interaction. Using the phase when the narrow emission features disappear, we calculate an outer dense radius of circumstellar material $R_\mathrm{CSM, out}\sim5\times10^{14}~\mathrm{cm}$ and a mean circumstellar material density of $\rho=5.6\times10^{-14}~\mathrm{g\,cm^{-3}}$. This is consistent with the lower limit on the outer radius of the circumstellar material we calculate from the peak \Halpha{} emission flux, $R_\text{CSM, out}\gtrsim9\times10^{13}~\mathrm{cm}$., Comment: Published in ApJL

We assemble the most complete and clean red supergiant (RSG) sample (2,121 targets) so far in the Small Magellanic Cloud (SMC) with 53 different bands of data to study the MLR of RSGs. In order to match the observed spectral energy distributions (SEDs), a theoretical grid of 17,820 Oxygen-rich models (``normal'' and ``dusty'' grids are half-and-half) is created by the radiatively-driven wind model of the DUSTY code, covering a wide range of dust parameters. We select the best model for each target by calculating the minimal modified chi-square and visual inspection. The resulting MLRs from DUSTY are converted to real MLRs based on the scaling relation, for which a total MLR of $6.16\times10^{-3}$ $M_\odot$ yr$^{-1}$ is measured (corresponding to a dust-production rate of $\sim6\times10^{-6}$ $M_\odot$ yr$^{-1}$), with a typical MLR of $\sim10^{-6}$ $M_\odot$ yr$^{-1}$ for the general population of the RSGs. The complexity of mass-loss estimation based on the SED is fully discussed for the first time, indicating large uncertainties based on the photometric data (potentially up to one order of magnitude or more). The Hertzsprung-Russell and luminosity versus median absolute deviation diagrams of the sample indicate the positive relation between luminosity and MLR. Meanwhile, the luminosity versus MLR diagrams show a ``knee-like'' shape with enhanced mass-loss occurring above $\log_{10}(L/L_\odot)\approx4.6$, which may be due to the degeneracy of luminosity, pulsation, low surface gravity, convection, and other factors. We derive our MLR relation by using a third-order polynomial to fit the sample and compare our result with previous empirical MLR prescriptions. Given that our MLR prescription is based on a much larger sample than previous determinations, it provides a more accurate relation at the cool and luminous region of the H-R diagram at low-metallicity compared to previous studies., Comment: 16 pages, 19 figures, accepted by A&A

We present current and future activities regarding lunar impact flash and NEO observations and satellite tracking from Kryoneri Observatory. In particular, we present results from the ESA-funded NELIOTA program, which has been monitoring the Moon for impact flashes since early 2017. Using the 1.2 m Kryoneri telescope, which is equipped with two high frame-rate cameras recording simultaneously in two optical bands, NELIOTA has recorded over 170 validated lunar impact flashes, while another ~90 have been characterized as suspected. We present statistical results concerning the sizes, the masses and the appearance frequency of the meteoroids in the vicinity of the Earth, as well as the temperatures developed during the impacts. Moreover, we present the capabilities of the Kryoneri telescope as a sensor for satellite tracking and the future plans regarding the provision of high-quality services for both the Planetary Defense activities of ESA (S2P/PDO) and the European Union's Space Surveillance and Tracking programme (EU/SST)., Comment: 14 pages, 19 figures. Proceedings of the 2nd NEO and Debris Detection Conference, Darmstadt, Germany, 24-26 January 2023

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 tables

Mass loss is a key property to understand stellar evolution and in particular for low-metallicity environments. Our knowledge has improved dramatically over the last decades both for single and binary evolutionary models. However, episodic mass loss although definitely present observationally, is not included in the models, while its role is currently undetermined. A major hindrance is the lack of large enough samples of classified stars. We attempted to address this by applying an ensemble machine-learning approach using color indices (from IR/Spitzer and optical/Pan-STARRS photometry) as features and combining the probabilities from three different algorithms. We trained on M31 and M33 sources with known spectral classification, which we grouped into Blue/Yellow/Red/B[e] Supergiants, Luminous Blue Variables, classical Wolf-Rayet and background galaxies/AGNs. We then applied the classifier to about one million Spitzer point sources from 25 nearby galaxies, spanning a range of metallicites ($1/15$ to $\sim3~Z_{\odot}$). Equipped with spectral classifications we investigated the occurrence of these populations with metallicity., Comment: 6 pages, 2 figures, submited to the proceedings of IAUS 361 Massive Stars Near and Far

(abridged) Mass loss is a key parameter in the evolution of massive stars, with discrepancies between theory and observations and with unknown importance of the episodic mass loss. To address this we need increased numbers of classified sources stars spanning a range of metallicity environments. We aim to remedy the situation by applying machine learning techniques to recently available extensive photometric catalogs. We used IR/Spitzer and optical/Pan-STARRS, with Gaia astrometric information, to compile a large catalog of known massive stars in M31 and M33, which were grouped in Blue, Red, Yellow, B[e] supergiants, Luminous Blue Variables, Wolf-Rayet, and background galaxies. Due to the high imbalance, we implemented synthetic data generation to populate the underrepresented classes and improve separation by undersampling the majority class. We built an ensemble classifier using color indices. The probabilities from Support Vector Classification, Random Forests, and Multi-layer Perceptron were combined for the final classification. The overall weighted balanced accuracy is ~83%, recovering Red supergiants at ~94%, Blue/Yellow/B[e] supergiants and background galaxies at ~50-80%, Wolf-Rayets at ~45%, and Luminous Blue Variables at ~30%, mainly due to their small sample sizes. The mixing of spectral types (no strict boundaries in their color indices) complicates the classification. Independent application to IC 1613, WLM, and Sextans A galaxies resulted in an overall lower accuracy of ~70%, attributed to metallicity and extinction effects. The missing data imputation was explored using simple replacement with mean values and an iterative imputor, which proved more capable. We also found that r-i and y-[3.6] were the most important features. Our method, although limited by the sampling of the feature space, is efficient in classifying sources with missing data and at lower metallicitites., Comment: 27 pages, 18 figures, 6 tables; accepted for publication in A&A; typos and minor editing corrections

We present a case study of using a novel method to identify red supergiant (RSG) candidates in NGC 6822, based on their 1.6 $\mu$m "H-bump". We collected 32 bands of photometric data for NGC 6822 ranging from optical to MIR. By using the theoretical spectra from MARCS, we demonstrate that there is a prominent difference around 1.6 $\mu$m ("H-bump") between low-surface-gravity (LSG) and high-surface-gravity (HSG) targets. Taking advantage of this feature, we identify efficient color-color diagrams (CCDs) of rzH and rzK to separate HSG and LSG targets from crossmatching of optical and NIR data. Moreover, synthetic photometry from ATLAS9 also give similar results. Further separating RSG candidates from the rest of the LSG candidates is done by using semi-empirical criteria on NIR CMDs and resulted in 323 RSG candidates. Meanwhile, the simulation of foreground stars from Besan\c{c}on models also indicates that our selection criteria is largely free from the contamination of Galactic giants. In addition to the "H-bump" method, we also use the traditional BVR method as a comparison and/or supplement, by applying a slightly aggressive cut to select as much as possible RSG candidates (358 targets). Furthermore, the Gaia astrometric solution is used to constrain the sample, where 181 and 193 targets were selected from the "H-bump" and BVR method, respectively. The percentages of selected targets in both methods are similar as $\sim$60\%, indicating the comparable accuracy of the two methods. In total, there are 234 RSG candidates after combining targets from both methods with 140 ($\sim$60\%) of them in common. The final RSG candidates are in the expected locations on the MIR CMDs, while the spatial distribution is also coincident with the FUV-selected star formation regions, suggesting the selection is reasonable and reliable., Comment: A&A accepted

We present a clean, magnitude-limited (IRAC1 or WISE1$\leq$15.0 mag) multiwavelength source catalog for the LMC. The catalog was built upon crossmatching ($1''$) and deblending ($3''$) between the SEIP source list and Gaia DR2, with strict constraints on the Gaia astrometric solution to remove the foreground contamination. The catalog contains 197,004 targets in 52 different bands including 2 ultraviolet, 21 optical, and 29 infrared bands. Additional information about radial velocities and spectral/photometric classifications were collected from the literature. The bright end of our sample is mostly comprised of blue helium-burning stars (BHeBs) and red HeBs with inevitable contamination of main sequence stars at the blue end. After applying modified magnitude and color cuts based on previous studies, we identify and rank 2,974 RSG, 508 YSG, and 4,786 BSG candidates in the LMC in six CMDs. The comparison between the CMDs of the LMC and SMC indicates that the most distinct difference appears at the bright red end of the optical and near-infrared CMDs, where the cool evolved stars (e.g., RSGs, AGB, and RGs) are located, which is likely due to the effect of metallicity and SFH. Further quantitative comparison of colors of massive star candidates in equal absolute magnitude bins suggests that, there is basically no difference for the BSG candidates, but large discrepancy for the RSG candidates as LMC targets are redder than the SMC ones, which may be due to the combined effect of metallicity on both spectral type and mass-loss rate, and also the age effect. The $T_{\rm eff}$ of massive star populations are also derived from reddening-free color of $(J-K_{\rm S})_0$. The $T_{\rm eff}$ ranges are $35008000$ K for BSG population, with larger uncertainties towards the hotter stars., Comment: 18 pages, 14 figures, A&A accepted

Rapid variability before and near the maximum brightness of supernovae has the potential to provide a better understanding of nearly every aspect of supernovae, from the physics of the explosion up to their progenitors and the circumstellar environment. Thanks to modern time-domain optical surveys, which are discovering supernovae in the early stage of their evolution, we have the unique opportunity to capture their intraday behavior before maximum. We present high-cadence photometric monitoring (on the order of seconds-minutes) of the optical light curves of three Type Ia and two Type II SNe over several nights before and near maximum light, using the fast imagers available on the 2.3~m Aristarchos telescope at Helmos Observatory and the 1.2~m telescope at Kryoneri Observatory in Greece. We applied differential aperture photometry techniques using optimal apertures and we present reconstructed light curves after implementing a seeing correction and the Trend Filtering Algorithm(TFA). TFA yielded the best results, achieving a typical precision between 0.01-0.04~mag. We did not detect significant bumps with amplitudes greater than 0.05~mag in any of the SNe targets in the VR-, R-, and I- bands light curves obtained. We measured the intraday slope for each light curve, which ranges between -0.37-0.36 mag/day in broadband VR, -0.19-0.31 mag/day in R band, and -0.13-0.10 mag/day in I band. We used SNe light curve fitting templates for SN 2018gv, SN 2018hgc and SN 2018hhn to photometrically classify the light curves and to calculate the time of maximum. We provide values for the maximum of SN 2018zd after applying a low-order polynomial fit and SN 2018hhn for the first time. We suggest monitoring early supernovae light curves in hotter (bluer) bands with a cadence of hours as a promising way of investigating the post-explosion photometric behavior of the progenitor stars., Comment: 22 pages, 17 figures

We present the most comprehensive RSG sample for the SMC up to now, including 1,239 RSG candidates. The initial sample is derived based on a source catalog for the SMC with conservative ranking. Additional spectroscopic RSGs are retrieved from the literature, as well as RSG candidates selected from the inspection of CMDs. We estimate that there are in total $\sim$ 1,800 or more RSGs in the SMC. We purify the sample by studying the infrared CMDs and the variability of the objects, though there is still an ambiguity between AGBs and RSGs. There are much less RSGs candidates ($\sim$4\%) showing PAH emission features compared to the Milky Way and LMC ($\sim$15\%). The MIR variability of RSG sample increases with luminosity. We separate the RSG sample into two subsamples ("risky" and "safe") and identify one M5e AGB star in the "risky" subsample. Most of the targets with large variability are also the bright ones with large MLR. Some targets show excessive dust emission, which may be related to previous episodic mass loss events. We also roughly estimate the total gas and dust budget produced by entire RSG population as $\rm \sim1.9^{+2.4}_{-1.1}\times10^{-6}~M_{\odot}/yr$ in the most conservative case. Based on the MIST models, we derive a linear relation between $T_{\rm eff}$ and observed $\rm J-K_S$ color with reddening correction for the RSG sample. By using a constant bolometric correction and this relation, the Geneva evolutionary model is compared with our RSG sample, showing a good agreement and a lower initial mass limit of $\sim$7 $\rm M_\odot $ for the RSG population. Finally, we compare the RSG sample in the SMC and the LMC. Despite the incompleteness of LMC sample in the faint end, the result indicates that the LMC sample always shows redder color (except for the $\rm IRAC1-IRAC2$ and $\rm WISE1-WISE2$ colors due to CO absorption) and larger variability than the SMC sample., Comment: 18 pages, 12 figures, A&A accepted

The detection of X-ray emission constitutes a reliable and efficient tool for the selection of Active Galactic Nuclei (AGNs), although it may be biased against the most heavily absorbed AGNs. Simple mid-IR broad-band selection criteria identify a large number of luminous and absorbed AGNs, yet again host contamination could lead to non-uniform and incomplete samples. Spectral Energy Distribution (SED) decomposition is able to decouple the emission from the AGN versus that from star-forming regions, revealing weaker AGN components. We aim to identify the obscured AGN population in the VIPERS survey in the CFHTLS W1 field through SED modelling. We construct SEDs for 6,860 sources and identify 160 AGNs at a high confidence level using a Bayesian approach. Using optical spectroscopy, we confirm the nature of ~85% of the AGNs. Our AGN sample is highly complete (~92%) compared to mid-IR colour selected AGNs, including a significant number of galaxy-dominated systems with lower luminosities. In addition to the lack of X-ray emission (80%), the SED fitting results suggest that the majority of the sources are obscured. We use a number of diagnostic criteria in the optical, infrared and X-ray regime to verify these results. Interestingly, only 35% of the most luminous mid-IR selected AGNs have X-ray counterparts suggesting strong absorption. Our work emphasizes the importance of using SED decomposition techniques to select a population of type II AGNs, which may remain undetected by either X-ray or IR colour surveys., Comment: Accepted for publication in MNRAS in May 4, 2020. 18 figures, 3 tables

Increasing the statistics of evolved massive stars in the Local Group enables investigating their evolution at different metallicities. During the late stages of stellar evolution, the physics of some phenomena, such as episodic and systematic mass loss, are not well constrained. For example, the physical properties of red supergiants (RSGs) in different metallicity regimes remain poorly understood. Thus, we initiated a systematic study of RSGs in dwarf irregular galaxies (dIrrs) in the Local Group. The target selection is based on 3.6 $\mu$m and 4.5 $\mu$m photometry from archival Spitzer Space Telescope images of nearby galaxies. We selected 46 targets in the dIrrs IC 10, IC 1613, Sextans B, and the Wolf-Lundmark-Melotte (WLM) galaxy that we observed with the GTC-OSIRIS and VLT-FORS2 instruments. We used several photometric techniques together with a spectral energy distribution analysis to derive the luminosities and effective temperatures of known and newly discovered RSGs. We identified and spectroscopically confirmed 4 new RSGs, 5 previously known RSGs, and 5 massive asymptotic giant branch (AGB) stars. We added known objects from previous observations. In total, we present spectral classification and fundamental physical parameters of 25 late-type massive stars in the following dIrrs: Sextans A, Sextans B, IC 10, IC 1613, Pegasus, Phoenix, and WLM. This includes 17 RSGs and 8 AGB stars that have been identified here and previously. Based on our observational results and PARSEC evolutionary models, we draw the following conclusions: (i) a trend to higher minimum effective temperatures at lower metallicities and (ii) the maximum luminosity of RSGs appears to be constant at $log$($L/L$$_{\odot}$) $\approx$ $5.5$, independent of the metallicity of the host environment (up to $\mathrm{[Fe/H]}$ $\approx$ $-1$ dex)., Comment: 22 pages, 17 figures, in press (A&A)

The Hubble Space Telescope (HST) has obtained multi-epoch observations providing the opportunity for a comprehensive variability search aiming to uncover new variables. We have therefore undertaken the task of creating a catalog of variable sources based on version 3 of the Hubble Source Catalog (HSC), which relies on publicly available images obtained with the WFPC2, ACS, and WFC3 instruments onboard the HST. We adopted magnitude-dependent thresholding in median absolute deviation (a robust measure of light curve scatter) combined with sophisticated preprocessing techniques and visual quality control to identify and validate variable sources observed by Hubble with the same instrument and filter combination five or more times. The Hubble Catalog of Variables (HCV) includes 84,428 candidate variable sources (out of 3.7 million HSC sources that were searched for variability) with $V \leq 27$ mag; for 11,115 of them the variability is detected in more than one filter. The data points in the light curves of the variables in the HCV catalog range from five to 120 points (typically having less than ten points); the time baseline ranges from under a day to over 15 years; while $\sim$8% of all variables have amplitudes in excess of 1 mag. Visual inspection performed on a subset of the candidate variables suggests that at least 80% of the candidate variables that passed our automated quality control are true variable sources rather than spurious detections resulting from blending, residual cosmic rays, and calibration errors. The HCV is the first, homogeneous catalog of variable sources created from archival HST data and currently is the deepest catalog of variables available. The catalog includes variable stars in our Galaxy and nearby galaxies, as well as transients and variable active galactic nuclei. (abbreviated), Comment: 33 pages, 15 figures, 11 tables; accepted in A&A. The HCV is available from the ESA Hubble Science Archive (eHST) (http://archives.esac.esa.int/ehst) at ESAC, and can be easily explored using the HCV Explorer (http://archives.esac.esa.int/hcv-explorer). The HCV is also hosted at STScI in the form of a High Level Science Product (HLSP) in MAST (https://archive.stsci.edu/hlsp/hcv)

Studies on the massive star population in galaxies beyond the Local Group are the key to understand the link between their numbers and modes of star formation in different environments. We present the analysis of the massive star population of the galaxies NGC 1326A, NGC 1425 and NGC 4548 using archival Hubble Space Telescope Wide Field Planetary Camera 2 images in the F555W and F814W filters. Through high precision point spread function fitting photometry for all sources in the three fields we identified 7640 candidate blue supergiants, 2314 candidate yellow supergiants, and 4270 candidate red supergiants. We provide an estimation the ratio of blue to red supergiants for each field as a function of galactocentric radius. Using Modules for Experiments in Stellar Astrophysics (MESA) at solar metallicity, we defined the luminosity function and estimated the star formation history of each galaxy. We carried out a variability search in the V and I filters using three variability indexes: the median absolute deviation, the interquartile range, and the inverse von-Neumann ratio. This analysis yielded 243 new variable candidates with absolute magnitudes ranging from Mv= -4 to -10 mag. We classified the variable stars based on their absolute magnitude and their position on the color-magnitude diagram using the MESA evolutionary tracks at solar metallicity. Our analysis yielded 8 candidate variable blue supergiants, 12 candidate variable yellow supergiants, 21 candidate variable red supergiants, and 4 candidate periodic variables., Comment: 28 pages, 11 figures, Astronomy & Astrophysics, accepted. arXiv admin note: text overlap with arXiv:1807.11493

We present a clean, magnitude-limited (IRAC1 or WISE1 $\leq$ 15.0 mag) multiwavelength source catalog for the SMC with 45,466 targets in total, with the purpose of building an anchor for future studies, especially for the massive star populations at low-metallicity. The catalog contains data in 50 different bands including 21 optical and 29 infrared bands, ranging from the ultraviolet to the far-infrared. Additionally, radial velocities and spectral classifications were collected from the literature, as well as infrared and optical variability statistics were retrieved from different projects. The catalog was essentially built upon a $1''$ crossmatching and a $3''$ deblending between the SEIP source list and Gaia DR2 photometric data. Further constraints on the proper motions and parallaxes from Gaia DR2 allowed us to remove the foreground contamination. We estimated that about 99.5\% of the targets in our catalog were most likely genuine members of the SMC. By using the evolutionary tracks and synthetic photometry from MIST and the theoretical $\rm J-K_S$ color cuts, we identified 1,405 RSG, 217 YSG and 1,369 BSG candidates in the SMC in five different CMDs, where attention should also be paid to the incompleteness of our sample. We ranked the candidates based on the intersection of different CMDs. A comparison between the models and observational data shows that the lower limit of initial mass for the RSGs population may be as low as 7 or even 6 $M_{\odot}$ and the RSG is well separated from the AGB population even at faint magnitude, making RSGs a unique population connecting the evolved massive and intermediate stars, since stars with initial mass around 6 to 8 $M_{\odot}$ are thought to go through a second dredge-up to become AGBs. We encourage the interested reader to further exploit the potential of our catalog., Comment: 19 pages, 17 figures, A&A accepted, catalogs are included in the ancillary files

Identifying Active Galactic Nuclei (AGNs) through their X-ray emission is efficient, but necessarily biased against X-ray-faint objects. We aim to characterize this bias by comparing X-ray-selected AGNs to the ones identified through optical variability and mid-IR colours. We present a catalogue of AGNs selected through optical variability using all publicly available z-band Hubble Space Telescope images in the GOODS-South field. For all objects in the catalogue, we compute X-ray upper limits or discuss detections in the deepest available 7 Ms Chandra Deep Field South images and present the Spitzer/IRAC mid-IR colours. For the variability study, we consider only sources observed over at least five epochs and over a time baseline of up to ten years. We adopt the elevated median absolute deviation as a variability indicator robust against individual outlier measurements and identify 113 variability-selected AGN candidates. Among these, 26 have an X-ray counterpart and lie within the conventional AGN area in the FX/Fopt diagram. The candidates with X-ray upper limits are on average optically fainter, have higher redshifts compared to the X-ray detected ones and are consistent with low luminosity AGNs. Out of 41 variable optical sources with IR detections, 13 fulfill the IR AGN colour selection criteria. Our work emphasizes the importance of optical variability surveys for constructing complete samples of AGNs including the ones that remain undetected even by the deepest X-ray and IR surveys., Comment: 22 pages, 23 figures, accepted for publication in MNRAS

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes., Comment: accepted for publication in A&A, 24 pages, 10 figures, tables with the data will be available electronically

We present the technical specifications and first results of the ESA-funded, lunar monitoring project "NELIOTA" (NEO Lunar Impacts and Optical TrAnsients) at the National Observatory of Athens, which aims to determine the size-frequency distribution of small Near-Earth Objects (NEOs) via detection of impact flashes on the surface of the Moon. For the purposes of this project a twin camera instrument was specially designed and installed at the 1.2 m Kryoneri telescope utilizing the fast-frame capabilities of scientific Complementary Metal-Oxide Semiconductor detectors (sCMOS). The system provides a wide field-of-view (17.0' $\times$ 14.4') and simultaneous observations in two photometric bands (R and I), reaching limiting magnitudes of 18.7 mag in 10 sec in both bands at a 2.5 signal-to-noise level. This makes it a unique instrument that can be used for the detection of NEO impacts on the Moon, as well as for any astronomy projects that demand high-cadence multicolor observations. The wide field-of-view ensures that a large portion of the Moon is observed, while the simultaneous, high-cadence, monitoring in two photometric bands makes possible, for the first time, the determination of the temperatures of the impacts on the Moon's surface and the validation of the impact flashes from a single site. Considering the varying background level on the Moon's surface we demonstrate that the NELIOTA system can detect NEO impact flashes at a 2.5 signal-to-noise level of ~12.4 mag in the I-band and R-band for observations made at low lunar phases ~0.1. We report 31 NEO impact flashes detected during the first year of the NELIOTA campaign. The faintest flash was at 11.24 mag in the R-band (about two magnitudes fainter than ever observed before) at lunar phase 0.32. Our observations suggest a detection rate of $1.96 \times 10^{-7}$ events $km^{-2} h^{-1}$., Comment: Accepted for publication in A&A

We analyzed the massive star population of the Virgo Cluster galaxy NGC 4535 using archival Hubble Space Telescope Wide Field Planetary Camera 2 images in filters F555W and F814W, equivalent to Johnson V and Kron-Cousins I. We performed high precision point spread function fitting photometry of 24353 sources including 3762 candidate blue supergiants, 841 candidate yellow supergiants and 370 candidate red supergiants. We estimated the ratio of blue to red supergiants as a decreasing function of galactocentric radius. Using Modules for Experiments in Stellar Astrophysics isochrones at solar metallicity, we defined the luminosity function and estimated the star formation history of the galaxy over the last 60 Myrs. We conducted a variability search in the V and I filters using three variability indexes: the median absolute deviation, the interquartile range and the inverse von-Neumann ratio. This analysis yielded 120 new variable candidates with absolute magnitudes ranging from M$_{V}$ = $-$4 to $-$11 mag. We used the MESA evolutionary tracks at solar metallicity, to classify the variables based on their absolute magnitude and their position on the color-magnitude diagram. Among the new candidate variable sources are eight candidate variable red supergiants, three candidate variable yellow supergiants and one candidate luminous blue variable, which we suggest for follow-up observations., Comment: Accepted by A&A, 7 pages, 7 Tables, 53 figures

The characteristics of IR properties and MIR variability of RSGs in the LMC are analyzed based on 12 bands of NIR to MIR co-added data from 2MASS, Spitzer and WISE, and $\sim$6.6 years of MIR time-series data collected by the ALLWISE and NEOWISE-R projects. 773 RSGs candidates were compiled from the literature and verified by using the CMD, SED and MIR variability. About 15\% of valid targets in the $IRAC1-IRAC2$/$IRAC2-IRAC3$ diagram may show PAH emission. We show that arbitrary dereddening Q parameters related to the IRAC4, S9W, WISE3, WISE4, and MIPS24 bands could be constructed based on a precise measurement of MIR interstellar extinction law. Several peculiar outliers in our sample are discussed, in which one outlier might be a RSG right before the explosion or an x-AGB star in the very late evolutionary stage based on the MIR spectrum and photometry. There are 744 identified RSGs in the final sample having both the WISE1- and WISE2-band time-series data. The results show that the MIR variability is increasing along with the increasing of brightness. There is a relatively tight correlation between the MIR variability, MLR, and the warm dust or continuum, where the MIR variability is evident for the targets with $K_S-WISE3>1.0~mag$ and $WISE4<6.5~mag$, while the rest of the targets show much smaller MIR variability. The MIR variability is also correlated with the MLR for which targets with larger variability also show larger MLR with an approximate upper limit of $-6.1~M_\odot/yr^{-1}$. Both the variability and the luminosity may be important for the MLR since the WISE4-band flux is increasing exponentially along with the degeneracy of luminosity and variability. The identified RSG sample has been compared with the theoretical evolutionary models and shown that the discrepancy between observation and evolutionary models can be mitigated by considering both variability and extinction., Comment: 24 pages, 22 figures, A&A accepted

Principal Component Analysis (PCA) is being extensively used in Astronomy but not yet exhaustively exploited for variability search. The aim of this work is to investigate the effectiveness of using the PCA as a method to search for variable stars in large photometric data sets. We apply PCA to variability indices computed for light curves of 18152 stars in three fields in M 31 extracted from the Hubble Source Catalogue. The projection of the data into the principal components is used as a stellar variability detection and classification tool, capable of distinguishing between RR Lyrae stars, long period variables (LPVs) and non-variables. This projection recovered more than 90% of the known variables and revealed 38 previously unknown variable stars (about 30% more), all LPVs except for one object of uncertain variability type. We conclude that this methodology can indeed successfully identify candidate variable stars., Comment: 23 pages, 18 figures, MNRAS accepted

The Hubble Source Catalog (HSC) combines lists of sources detected on images obtained with the WFPC2, ACS and WFC3 instruments aboard the Hubble Space Telescope (HST) available in the Hubble Legacy Archive. The catalog contains time-domain information with about two million of its sources detected with the same instrument and filter in at least five HST visits. The Hubble Catalog of Variables (HCV) project aims to identify HSC sources showing significant brightness variations. A magnitude-dependent threshold in the median absolute deviation of photometric measurements (an outlier-resistant measure of lightcurve scatter) is adopted as the variability-detection statistic. It is supplemented with a cut in $\chi_{\rm red}^2$ that removes sources with large photometric errors. A pre-processing procedure involving bad image identification, outlier rejection and computation of local magnitude zero-point corrections is applied to HSC lightcurves before computing the variability detection statistic. About 52000 HSC sources are identified as candidate variables, among which 7800 show variability in more than one filter. Visual inspection suggests that $\sim 70\%$ of the candidates detected in multiple filters are true variables while the remaining $\sim 30\%$ are sources with aperture photometry corrupted by blending, imaging artifacts or image processing anomalies. The candidate variables have AB magnitudes in the range 15-27$^{m}$ with the median 22$^{m}$. Among them are the stars in our own and nearby galaxies as well as active galactic nuclei., Comment: 4 pages, 1 figure, proceedings of the IAU Symposium 339 Southern Horizons in Time-Domain Astronomy, 13-17 November 2017, Stellenbosch, South Africa

This paper presents a detailed analysis of the light and radial velocity curves of the semi-detached eclipsing binary system OGLE-LMC-ECL-09937. The system is composed of a hot, massive and luminous primary star of a late-O spectral type, and a more evolved, but less massive and luminous secondary, implying an Algol-type system that underwent a mass transfer episode. We derive masses of 21.04 +/- 0.34 M_Sun and 7.61 +/- 0.09 M_Sun and radii of 9.93 +/- 0.06 R_Sun and 9.18 +/- 0.04 R_Sun, for the primary and the secondary component, respectively, which make it the most massive known Algol-type system with masses and radii of the components measured with <2% accuracy. Consequently, the parameters of OGLE-LMC-ECL-09937 provide an important contribution to the sparsely populated high-mass end of the stellar mass distribution, and an interesting object for stellar evolution studies, being a possible progenitor of a binary system composed of two neutron stars., Comment: 16 pages, 7 figures, 4 Tables, published in Acta Astronomica

The Hubble Catalog of Variables (HCV) project aims to identify the variable sources in the Hubble Source Catalog (HSC), which includes about 92 million objects with over 300 million measurements detected by the WFPC2, ACS and WFC3 cameras on board of the Hubble Space Telescope (HST), by using an automated pipeline containing a set of detection and validation algorithms. All the HSC sources with more than a predefined number of measurements in a single filter/instrument combination are pre-processed to correct systematic effect and to remove the bad measurements. The corrected data are used to compute a number of variability indexes to determine the variability status of each source. The final variable source catalog will contain variables stars, active galactic nuclei (AGNs), supernovae (SNs) or even new types of variables, reaching an unprecedented depth (V$\leq$27 mag). At the end of the project, the first release of the HCV will be available at the Mikulski Archive for Space Telescopes (MAST) and the ESA Hubble Science Archives. The HCV pipeline will be deployed at the Space Telescope Science Institute (STScI) so that an updated HCV may be generated following future releases of HSC., Comment: 4 pages, 1 figure, to appear in the proceeding of "Stellar Populations and the Distance Scale", Joseph Jensen (eds), ASPCS

The survey for DUST in Nearby Galaxies with Spitzer (DUSTiNGS) identified several candidate Asymptotic Giant Branch (AGB) stars in nearby dwarf galaxies and showed that dust can form even in very metal-poor systems (Z ~ 0.008 $Z_\odot$). Here, we present a follow-up survey with WFC3/IR on the Hubble Space Telescope (HST), using filters that are capable of distinguishing carbon-rich (C-type) stars from oxygen-rich (M-type) stars: F127M, F139M, and F153M. We include six star-forming DUSTiNGS galaxies (NGC 147, IC 10, Pegasus dIrr, Sextans B, Sextans A, and Sag DIG), all more metal-poor than the Magellanic Clouds and spanning 1 dex in metallicity. We double the number of dusty AGB stars known in these galaxies and find that most are carbon rich. We also find 26 dusty M-type stars, mostly in IC 10. Given the large dust excess and tight spatial distribution of these M-type stars, they are most likely on the upper end of the AGB mass range (stars undergoing Hot Bottom Burning). Theoretical models do not predict significant dust production in metal-poor M-type stars, but we see evidence for dust excess around M-type stars even in the most metal-poor galaxies in our sample (12+log(O/H) = 7.26-7.50). The low metallicities and inferred high stellar masses (up to ~10 $M_\odot$) suggest that AGB stars can produce dust very early in the evolution of galaxies (~30 Myr after they form), and may contribute significantly to the dust reservoirs seen in high-redshift galaxies., Comment: 18 pages, 9 figures, accepted for publication in ApJ

We report the first scientific results from the NELIOTA (NEO Lunar Impacts and Optical TrAnsients) project, which has recently begun lunar monitoring observations with the 1.2-m Kryoneri telescope. NELIOTA aims to detect faint impact flashes produced by near-Earth meteoroids and asteroids and thereby help constrain the size-frequency distribution of near-Earth objects in the decimeter to meter range. The NELIOTA setup, consisting of two fast-frame cameras observing simultaneously in the $R$ and $I-$bands, enables - for the first time - direct analytical calculation of the flash temperatures. We present the first 10 flashes detected, for which we find temperatures in the range ~1,600-3,100 K, in agreement with theoretical values. Two of these flashes were detected on multiple frames in both filters and therefore yield the first measurements of the temperature drop for lunar flashes. In addition, we compute the impactor masses, which range between ~100 g and ~50 kg., Comment: 6 pages, 1 figure, 1 table, accepted for publication in A&A

We have characterized the spectroscopic orbit of the TWA 3A binary and provide preliminary families of probable solutions for the TWA 3A visual orbit as well as for the wide TWA 3A--B orbit. TWA 3 is a hierarchical triple located at 34 pc in the $\sim$10 Myr old TW Hya association. The wide component separation is 1."55; the close pair was first identified as a possible binary almost 20 years ago. We initially identified the 35-day period orbital solution using high-resolution infrared spectroscopy which angularly resolved the A and B components. We then refined the preliminary orbit by combining the infrared data with a re-analysis of our high-resolution optical spectroscopy. The orbital period from the combined spectroscopic solution is $\sim$35 days, the eccentricity is $\sim$0.63, and the mass ratio is $\sim$0.84; although this high mass ratio would suggest that optical spectroscopy alone should be sufficient to identify the orbital solution, the presence of the tertiary B component likely introduced confusion in the blended optical spectra. Using millimeter imaging from the literature, we also estimate the inclinations of the stellar orbital planes with respect to the TWA 3A circumbinary disk inclination and find that all three planes are likely misaligned by at least $\sim$30 degrees. The TWA 3A spectroscopic binary components have spectral types of M4.0 and M4.5; TWA 3B is an M3. We speculate that the system formed as a triple, is bound, and that its properties were shaped by dynamical interactions between the inclined orbits and disk., Comment: Accepted to ApJ

The Hubble Catalog of Variables (HCV) is a 3 year ESA funded project that aims to develop a set of algorithms to identify variables among the sources included in the Hubble Source Catalog (HSC) and produce the HCV. We will process all HSC sources with more than a predefined number of measurements in a single filter/instrument combination and compute a range of lightcurve features to determine the variability status of each source. At the end of the project, the first release of the Hubble Catalog of Variables will be made available at the Mikulski Archive for Space Telescopes (MAST) and the ESA Science Archives. The variability detection pipeline will be implemented at the Space Telescope Science Institute (STScI) so that updated versions of the HCV may be created following the future releases of the HSC., Comment: 4 pages, 5 figures, To appear in the conference proceedings of the IAU Symposium 325 AstroInformatics (2016 October 20-24, Sorrento, Italy)

The evolution of massive stars surviving the red supergiant (RSG) stage remains unexplored due to the rarity of such objects. The yellow hypergiants (YHGs) appear to be the warm counterparts of post-RSG classes located near the Humphreys-Davidson upper luminosity limit, which are characterized by atmospheric instability and high mass-loss rates. We aim to increase the number of YHGs in M33 and thus to contribute to a better understanding of the pre-supernova evolution of massive stars. Optical spectroscopy of five dust-enshrouded YSGs selected from mid-IR criteria was obtained with the goal of detecting evidence of extensive atmospheres. We also analyzed BVI photometry for 21 of the most luminous YSGs in M33 to identify changes in the spectral type. To explore the properties of circumstellar dust, we performed SED-fitting of multi-band photometry of the 21 YSGs. We find three luminous YSGs in our sample to be YHG candidates, as they are surrounded by hot dust and are enshrouded within extended, cold dusty envelopes. Our spectroscopy of star 2 shows emission of more than one H$\alpha$ component, as well as emission of CaII, implying an extended atmospheric structure. In addition, the long-term monitoring of the star reveals a dimming in the visual light curve of amplitude larger than 0.5 mag that caused an apparent drop in the temperature that exceeded 500 K. We suggest the observed variability to be analogous to that of the Galactic YHG $\rho$ Cas. Five less luminous YSGs are suggested as post-RSG candidates showing evidence of hot or/and cool dust emission. We demonstrate that mid-IR photometry, combined with optical spectroscopy and time-series photometry, provide a robust method for identifying candidate YHGs. Future discovery of YHGs in Local Group galaxies is critical for the study of the late evolution of intermediate-mass massive stars., Comment: 24 pages, 12 figures, 7 Tables. A&A in press

We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr., Comment: Accepted for Publication in A&A, 16 pages, 17 figures and 4 tables

Photometric measurements are prone to systematic errors presenting a challenge to low-amplitude variability detection. In search for a general-purpose variability detection technique able to recover a broad range of variability types including currently unknown ones, we test 18 statistical characteristics quantifying scatter and/or correlation between brightness measurements. We compare their performance in identifying variable objects in seven time series data sets obtained with telescopes ranging in size from a telephoto lens to 1m-class and probing variability on time-scales from minutes to decades. The test data sets together include lightcurves of 127539 objects, among them 1251 variable stars of various types and represent a range of observing conditions often found in ground-based variability surveys. The real data are complemented by simulations. We propose a combination of two indices that together recover a broad range of variability types from photometric data characterized by a wide variety of sampling patterns, photometric accuracies, and percentages of outlier measurements. The first index is the interquartile range (IQR) of magnitude measurements, sensitive to variability irrespective of a time-scale and resistant to outliers. It can be complemented by the ratio of the lightcurve variance to the mean square successive difference, 1/h, which is efficient in detecting variability on time-scales longer than the typical time interval between observations. Variable objects have larger 1/h and/or IQR values than non-variable objects of similar brightness. Another approach to variability detection is to combine many variability indices using principal component analysis. We present 124 previously unknown variable stars found in the test data., Comment: 29 pages, 8 figures, 7 tables; accepted to MNRAS; for additional plots, see http://scan.sai.msu.ru/~kirx/var_idx_paper/

We present a detailed analysis of the pre-main-sequence (PMS) population of the young star cluster Westerlund~2 (Wd2), the central ionizing cluster of the HII region RCW 49, using data from a high resolution multi-band survey with the Hubble Space Telescope. The data were acquired with the Advanced Camera for Surveys in the F555W, F814W, and F658N filters and with the Wide Field Camera 3 in the F125W, F160W, and F128N filters. We find a mean age of the region of 1.04+-0.72 Myr. The combination of dereddened F555W and F814W photometry in combination with F658N photometry allows us to study and identify stars with H_alpha excess emission. With a careful selection of 240 bona-fide PMS H_alpha excess emitters we were able to determine their H_alpha luminosity, which has a mean value L(H_alpha)=1.67 x 10^{-31} erg s^{-1}. Using the PARSEC 1.2S isochrones to obtain the stellar parameters of the PMS stars we determined a mean mass accretion rate \dot M_acc=4.43 x 10^{-8} M_sun yr^{-1} per star. A careful analysis of the spatial dependence of the mass-accretion rate suggests that this rate is ~25% lower in center of the two density peaks of Wd2 in close proximity to the luminous OB stars, compared to the Wd2 average. This rate is higher with increasing distance from the OB stars, indicating that the PMS accretion disks are being rapidly destroyed by the far-ultra-violet radiation emitted by the OB population., Comment: 21 pages, 17 figures, 4 tables

We present $Spitzer$ IRAC mid-infrared point source catalogs for mosaics covering the fields of the nearby ($\lesssim$4 Mpc) galaxies NGC 55, NGC 253, NGC 2366, NGC 4214, and NGC 5253. We detect a total of 20159 sources in these five fields. Point spread function photometry was performed on sources detected in both $Spitzer$ IRAC 3.6 $\mu$m and 4.5 $\mu$m bands at greater than 3$\sigma$ above background. These data were then supplemented by aperture photometry in the IRAC 5.8 $\mu$m and 8.0 $\mu$m bands conducted at the positions of the shorter wavelength sources. For sources with no detected object in the longer wavelengths, we estimated magnitude limits based on the local sky background. The individual galaxy point source breakdown is the following: NGC 55, 8746 sources; NGC 253, 9001 sources; NGC 2366, 505 sources; NGC 4214, 1185 sources; NGC 5253, 722 sources. The completeness limits of the full catalog vary with bandpass and were found to be $m_{3.6}=18.0$, $m_{4.5}=17.5$, $m_{5.8}=17.0$, and $m_{8.0}=16.5$ mag. For all galaxies, this corresponds to detection of point sources brighter than $M_{3.6}=-10$. These catalogs can be used as a reference for stellar population investigations, individual stellar object studies, and in planning future mid-infrared observations with the James Webb Space Telescope., Comment: 20 pages, 14 figures, accepted to Astronomy and Astrophysics

NELIOTA is a new ESA activity launched at the National Observatory of Athens in February 2015 aiming to determine the distribution and frequency of small near-earth objects (NEOs) via lunar monitoring. The project involves upgrading the 1.2m Kryoneri telescope at the National Observatory of Athens, procuring two fast-frame cameras, and developing a software system, which will control the telescope and the cameras, process the images and automatically detect NEO impacts. NELIOTA will provide a web-based user interface, where the impact events will be reported and made available to the scientific community and the general public. The objective of this 3.5 year activity is to design, develop and implement a highly automated lunar monitoring system, which will conduct an observing campaign for 2 years in search of NEO impact flashes on the Moon. The impact events will be verified, characterised and reported. The 1.2m telescope will be capable of detecting flashes much fainter than current, small-aperture, lunar monitoring telescopes. NELIOTA is therefore expected to characterise the frequency and distribution of NEOs weighing as little as a few grams., Comment: 3 pages, 1 figure, IAU GA Symposium 318 proceedings

Increasing the statistics of spectroscopically confirmed evolved massive stars in the Local Group enables the investigation of the mass loss phenomena that occur in these stars in the late stages of their evolution. We aim to complete the census of luminous mid-IR sources in star-forming dwarf irregular (dIrr) galaxies of the Local Group. To achieve this we employed mid-IR photometric selection criteria to identify evolved massive stars, such as red supergiants (RSGs) and luminous blue variables (LBVs), by using the fact that these types of stars have infrared excess due to dust. The method is based on 3.6 $\mu$m and 4.5 $\mu$m photometry from archival ${\it Spitzer}$ Space Telescope images of nearby galaxies. We applied our criteria to 4 dIrr galaxies: Pegasus, Phoenix, Sextans A, and WLM, selecting 79 point sources, which we observed with the VLT/FORS2 spectrograph in multi-object spectroscopy mode. We identified 13 RSGs, of which 6 are new discoveries, also 2 new emission line stars, and 1 candidate yellow supergiant. Among the other observed objects we identified carbon stars, foreground giants, and background objects, such as a quasar and an early-type galaxy that contaminate our survey. We use the results of our spectroscopic survey to revise the mid-IR and optical selection criteria for identifying RSGs from photometric measurements. The optical selection criteria are more efficient in separating extragalactic RSGs from foreground giants than mid-IR selection criteria, however the mid-IR selection criteria are useful for identifying dusty stars in the Local Group. This work serves as a basis for further investigation of the newly discovered dusty massive stars and their host galaxies., Comment: 19 pages, 7 tables, 11 figures. Accepted for publication in A&A

We present a study of the properties of the O-type, massive eclipsing binary 2MASS J13130841-6239275 located in the outskirts of the Danks 2 cluster in the G305 star-forming complex, using near-infrared spectroscopy from VLT/ISAAC. We derive the masses and radii to be 24.5$\pm$0.9 M$_{\odot}$ and 9.2$\pm$0.1 R$_{\odot}$ for the primary and 21.7$\pm$0.8 M$_{\odot}$ and 8.7$\pm$0.1 R$_{\odot}$ for the secondary component. In addition, we evaluate the sensitivity of our parameters to the choice of the spectral features used to determine the radial velocities. Both components appear to be main-sequence O6.5$-$O7 type stars at an age of $\sim$5 Myr, which is in agreement with the age of the cluster. A high visual extinction of A$_{5495}$=11.9$\pm$0.1 mag is reported, which is likely attributed to the cold molecular gas contaminating the north-east region of the cluster. By fitting the spectral energy distribution of the system to the available $BVI_{c}JHK_{s}$ photometry, we determine a distance to the system of 3.52$\pm$0.08 kpc with a precision of 2$\%$, which is the most well-determined distance to the Danks 2 cluster and the host complex reported in the literature., Comment: 13 pages, 9 figures, 6 tables. Accepted for publication in A&A

We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. We estimate $T_{\rm eff}=27$ kK, $\log{(g)}=2.9$ (cgs) for the VFTS 450 primary spectrum (34kK, 3.6: for the secondary spectrum); and $T_{\rm eff} = 22$kK, $\log{(g)}=2.8$ for the VFTS 652 primary spectrum (35kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89d and 8.59d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range $\sim$20--50M$_\odot$. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.

We present results of high-cadence monitoring of the optical light curve of the nearby, Type Ia SN 2014J in M82 using the 2.3m Aristarchos telescope. $B$ and $V$-band photometry on days 15-18 after $t_{max}(B)$, obtained with a cadence of 2 min per band, reveals evidence for rapid variability at the 0.02-0.05 mag level on timescales of 15-60 min on all four nights, taking the red noise estimation at face value. The decline slope was measured to be steeper in the $B$-band than in $V$-band, and to steadily decrease in both bands from 0.15 mag/day (night 1) to 0.04 mag/day (night 4) in V and from 0.19 mag/day (night 1) to 0.06 mag/day (night 4) in B, corresponding to the onset of the secondary maximum. We propose that rapid variability could be due to one or a combination of the following scenarios: the clumpiness of the ejecta, their interaction with circumstellar material, the asymmetry of the explosion, or the mechanism causing the secondary maximum in the near-infrared light curve. We encourage the community to undertake high-cadence monitoring of future, nearby and bright supernovae to investigate the intraday behavior of their light curves., Comment: 9 pages, 6 figures, 2 tables. Accepted for publication in A&A

We present first results from a high resolution multi-band survey of the Westerlund 2 region with the Hubble Space Telescope. Specifically, we imaged Westerlund 2 with the Advanced Camera for Surveys through the $F555W$, $F814W$, and $F658N$ filters and with the Wide Field Camera 3 in the $F125W$, $F160W$, and $F128N$ filters. We derive the first high resolution pixel-to-pixel map of the color excess $E(B-V)_g$ of the gas associated with the cluster, combining the H$\alpha$ ($F658N$) and Pa$\beta$ ($F128N$) line observations. We demonstrate that, as expected, the region is affected by significant differential reddening with a median of $E(B-V)_g=1.87$~mag. After separating the populations of cluster members and foreground contaminants using a $(F814W-F160W)$ vs. $F814W$ color-magnitude diagram, we identify a pronounced pre-main-sequence population in Westerlund 2 showing a distinct turn-on. After dereddening each star of Westerlund 2 individually in the color-magnitude diagram we find via over-plotting PARSEC isochrones that the distance is in good agreement with the literature value of $\sim4.16 \pm 0.33$~kpc. With zero-age-main-sequence fitting to two-color-diagrams, we derive a value of total to selective extinction of $R_V=3.95 \pm 0.135$. A spatial density map of the stellar content reveals that the cluster might be composed of two clumps. We estimate the same age of 0.5-2.0 Myr for both clumps. While the two clumps appear to be coeval, the northern clump shows a $\sim 20 \%$ lower stellar surface density., Comment: 24 pages, 27 figures, 7 tables; Accepted for publication to The Astronomical Journal

We present a quantitative spectroscopic study of twenty-seven red supergiants in the Sculptor Galaxy NGC 300. J-band spectra were obtained using KMOS on the VLT and studied with state of the art synthetic spectra including NLTE corrections for the strongest diagnostic lines. We report a central metallicity of [Z]= -0.03 +/- 0.05 with a gradient of -0.083 +/- 0.014 [dex/kpc], in agreement with previous studies of blue supergiants and H II-region auroral line measurements. This result marks the first application of the J-band spectroscopic method to a population of individual red supergiant stars beyond the Local Group of galaxies and reveals the great potential of this technique., Comment: 11 pages

We present an analysis of archival Spitzer images and new ground-based and Hubble Space Telescope (HST) near-infrared (IR) and optical images of the field of M83 with the goal of identifying rare, dusty, evolved massive stars. We present point source catalogs consisting of 3778 objects from $Spitzer$ Infrared Array Camera (IRAC) Band 1 (3.6 $\mu$m) and Band 2 (4.5 $\mu$m), and 975 objects identified in Magellan 6.5m FourStar near-IR $J$ and $K_{\rm s}$ images. A combined catalog of coordinate matched near- and mid-IR point sources yields 221 objects in the field of M83. Using this photometry we identify 185 massive evolved stellar candidates based on their location in color-magnitude and color-color diagrams. We estimate the background contamination to our stellar candidate lists and further classify candidates based on their appearance in $HST$ Wide Field Camera 3 (WFC3) observations of M83. We find 49 strong candidates for massive stars which are very promising objects for spectroscopic follow-up. Based on their location in a $B-V$ versus $V-I$ diagram, we expect at least 24, or roughly 50%, to be confirmed as red supergiants., Comment: 32 pages, 23 figures, accepted for publication in A&A