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.

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.

This paper contains the results from the first 30 months of the NELIOTA project for impacts of NEOs/meteoroids on the lunar surface. Our analysis on the statistics concerning the efficiency of the campaign and the parameters of the projectiles and the impacts is presented. The parameters of the lunar impact flashes based on simultaneous observations in two wavelength bands are used to estimate the distributions of the masses, sizes and frequency of the impactors. These statistics can be used both in space engineering and science. The photometric fluxes of the flashes are measured using aperture photometry and their apparent magnitudes are calculated using standard stars. Assuming that the flashes follow a black body law of irradiation, the temperatures can be derived analytically, while the parameters of the projectiles are estimated using fair assumptions on their velocity and luminous efficiency of the impacts. 79 lunar impact flashes have been observed with the 1.2 m Kryoneri telescope in Greece. The masses of the meteoroids range between 0.7 g and 8 kg and their respective sizes between 1-20 cm depending on their assumed density, impact velocity, and luminous efficiency. We find a strong correlation between the observed magnitudes of the flashes and the masses of the meteoroids. Moreover, an empirical relation between the emitted energies of each band has been derived allowing the estimation of the physical parameters of the meteoroids that produce low energy impact flashes. The NELIOTA project has so far the highest detection rate and the faintest limiting magnitude for lunar impacts compared to other ongoing programs. Based on the impact frequency distribution on Moon, we estimate that sporadic meteoroids with typical masses less than 100 g and sizes less than 5 cm enter the mesosphere of the Earth with a rate ~108 meteoroids/hr and also impact Moon with a rate of ~8 meteoroids/hr.

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.

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).

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)

One of the most exciting scientific challenges is to detect Earth-like planets in the habitable zones of other stars in the galaxy and search for evidence of life. The ability to observe and characterise dozens of potentially Earth-like planets now lies within the realm of possibility due to rapid advances in key space and imaging technologies. The associated challenge of directly imaging very faint planets in orbit around nearby very bright stars is now well understood, with the key instrumentation also being perfected and developed. Such advances will allow us to develop large transformative telescopes, covering a broad UV-optical-IR spectral range, which can carry out the detailed research programmes designed to answer the questions we wish to answer: Carry out high contrast imaging surveys of nearby stars to search for planets within their habitable zones. Characterise the planets detected to determine masses and radii from photometric measurements. Through spectroscopic studies of their atmospheres and surfaces, search for habitability indicators and for signs of an environment that has been modified by the presence of life. Active studies of potential missions have been underway for a number of years. The latest of these is the Large UV Optical IR space telescope (LUVOIR), one of four flagship mission studies commissioned by NASA in support of the 2020 US Decadal Survey. LUVOIR, if selected, will be of interest to a wide scientific community and will be the only telescope capable of searching for and characterizing a sufficient number of exoEarths to provide a meaningful answer to the question - Are we alone?. This paper is a submission to the European Space Agency Voyage 2050 call for white papers outlining the case for an ESA contribution to a Large UVOIR telescope.

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.

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.

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.

NELIOTA (NEO Lunar Impacts and Optical TrAnsients) is an ESA-funded lunar monitoring project, 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. A prime focus, high-speed, twin-camera system providing simultaneous observations in two photometric bands at a rate of 30 frames-per-second on the 1.2 m Kryoneri telescope of the National Observatory of Athens was commissioned for this purpose. A dedicated software processes the images and automatically detects candidate lunar impact flashes, which are then validated by an expert user. The four year observing campaign began in February 2017 and has so far detected more than 40 lunar impact events. The software routinely detects satellites, which typically appear as streaks or dots crossing the lunar disk. To avoid confusing these events with real flashes, we check different available catalogs with spacecraft orbital information and exclude spacecraft identifications.

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}$.

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.

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.

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.

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.

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.

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.

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.

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.

We aim to construct an exceptionally deep (V ~< 27) catalog of variable objects in selected Galactic and extragalactic fields visited multiple times by the Hubble Space Telescope (HST). While HST observations of some of these fields were searched for specific types of variables before (most notably, the extragalactic Cepheids), we attempt a systematic study of the population of variable objects of all types at the magnitude range not easily accessible with ground-based telescopes. The variability timescales that can be probed range from hours to years depending on how often a particular field has been visited. For source extraction and cross-matching of sources between visits we rely on the Hubble Source Catalog which includes 10^7 objects detected with WFPC2, ACS, and WFC3 HST instruments. The lightcurves extracted from the HSC are corrected for systematic effects by applying local zero-point corrections and are screened for bad measurements. For each lightcurve we compute variability indices sensitive to a broad range of variability types. The indices characterize the overall lightcurve scatter and smoothness. Candidate variables are selected as having variability index values significantly higher than expected for objects of similar brightness in the given set of observations. The Hubble Catalog of Variables will be released in 2018.

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.

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$.

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.

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.

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.

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.

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.

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

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.

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.

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.

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.

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 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.

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.

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.

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.

The DUSTiNGS survey (DUST in Nearby Galaxies with Spitzer) is a 3.6 and 4.5 micron imaging survey of 50 nearby dwarf galaxies designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. Using two epochs, spaced approximately six months apart, we identify a total of 526 dusty variable AGB stars (sometimes called "extreme" or x-AGB stars; [3.6]-[4.5]>0.1 mag). Of these, 111 are in galaxies with [Fe/H] < -1.5 and 12 are in galaxies with [Fe/H] < -2.0, making them the most metal-poor dust-producing AGB stars known. We compare these identifications to those in the literature and find that most are newly discovered large-amplitude variables, with the exception of approximately 30 stars in NGC 185 and NGC 147, one star in IC 1613, and one star in Phoenix. The chemical abundances of the x-AGB variables are unknown, but the low metallicities suggest that they are more likely to be carbon-rich than oxygen-rich and comparisons with existing optical and near-IR photometry confirms that 70 of the x-AGB variables are confirmed or likely carbon stars. We see an increase in the pulsation amplitude with increased dust production, supporting previous studies suggesting that dust production and pulsation are linked. We find no strong evidence linking dust production with metallicity, indicating that dust can form in very metal-poor environments.

Nearby resolved dwarf galaxies provide excellent opportunities for studying the dust-producing late stages of stellar evolution over a wide range of metallicity (-2.7 < [Fe/H] < -1.0). Here, we describe DUSTiNGS (DUST in Nearby Galaxies with Spitzer): a 3.6 and 4.5 micron post-cryogen Spitzer Space Telescope imaging survey of 50 dwarf galaxies within 1.5 Mpc that is designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. The survey includes 37 dwarf spheroidal, 8 dwarf irregular, and 5 transition-type galaxies. This near-complete sample allows for the building of statistics on these rare phases of stellar evolution over the full metallicity range. The photometry is >75% complete at the tip of the Red Giant Branch for all targeted galaxies, with the exception of the crowded inner regions of IC 10, NGC 185, and NGC 147. This photometric depth ensures that the majority of the dust-producing stars, including the thermally-pulsing AGB stars, are detected in each galaxy. The images map each galaxy to at least twice the half-light radius to ensure that the entire evolved star population is included and to facilitate the statistical subtraction of background and foreground contamination, which is severe at these wavelengths. In this overview, we describe the survey, the data products, and preliminary results. We show evidence for the presence of dust-producing AGB stars in 8 of the targeted galaxies, with metallicities as low as [Fe/H] = -1.9, suggesting that dust production occurs even at low metallicity.

We present a variability study of 4646 massive stars in the Small Magellanic Cloud (SMC) with known spectral types from the catalog of Bonanos et al. (2010) using the light curves from the OGLE-III database. The goal is to exploit the time domain information available through OGLE-III to gain insight into the processes that govern the evolution of massive stars. This variability survey of massive stars with known spectral types is larger than any previous survey by a factor of 7. We find that 60% of our sample (2766 stars) show no significant variability and 40% (1880 stars) exhibit variability distributed as follows: 807 stars display low-amplitude stochastic variability with fluctuations in I-band of up to 0.05 mag, 443 stars present irregular variability of higher amplitude (76% of these are reported as variables for the first time), 205 are eclipsing binaries (including 101 newly discovered systems), 50 are candidate rotating variables, 126 are classical Cepheids, 188 stars exhibit short-term sinusoidal periodicity (P < 3 days) making them candidate "slowly pulsating B stars" and non-radial Be pulsators, and 61 periodic stars exhibit longer periods. We demonstrate the wealth of information provided in the time domain, by doubling the number of known massive eclipsing binary systems and identifying 189 new candidate early-type Be and 20 Oe stars in the SMC. In addition, we find that ~80% of Be stars are photometrically variable in the OGLE-III time domain and provide evidence that short-term pulsating stars with additional photometric variability are rotating close to their break-up velocity.

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.

The Magellanic Quasars Survey (MQS) has now increased the number of quasars known behind the Magellanic Clouds by almost an order of magnitude. All survey fields in the Large Magellanic Cloud (LMC) and 70% of those in the Small Magellanic Cloud (SMC) have been observed. The targets were selected from the third phase of the Optical Gravitational Lensing Experiment (OGLE-III) based on their optical variability, mid-IR and/or X-ray properties. We spectroscopically confirmed 758 (565 LMC and 193 SMC) quasars behind the Clouds, of which 94% (527 LMC and 186 SMC) are newly identified. The MQS quasars have long-term (12 years and growing for OGLE), high-cadence light curves, enabling unprecedented variability studies of quasars. The MQS quasars also provide a dense reference grid for measuring both the internal and bulk proper motions of the Clouds, and 50 quasars are bright enough (I<18 mag) for absorption studies of the interstellar/galactic (ISM/IGM) medium of the Clouds.

Only a small number of high mass stars (> 30 Mo) have fundamental parameters (i.e. masses and radii) measured with high enough accuracy from eclipsing binaries to constrain formation and evolutionary models of massive stars. This work aims to increase this limited sample, by studying the 4 massive eclipsing binary candidates discovered by Bonanos in the young massive cluster Westerlund 1. We present new follow-up echelle spectroscopy of these binaries and models of their light and radial velocity curves. We obtain fundamental parameters (i.e. masses, radii) for the 8 component stars, finding masses that span a range of 10-40 Mo, and contributing accurate fundamental parameters for 1 additional very massive star, the 33 Mo component of W13. WR77o is found to have a ~40 Mo companion, which provides a second dynamical constraint on the mass of the progenitor of the magnetar known in the cluster. We also use W13 to estimate the first, direct, eclipsing binary distance to Westerlund 1 and therefore the magnetar, and find it to be at 3.7 +/- 0.6 kpc. Our results confirm previous evidence for a high mass for the progenitor of the magnetar. In addition, the availability of eclipsing binaries with accurate parameters opens the way for direct, independent, high precision eclipsing binary distance measurements to Westerlund 1.

The nearest accreting T Tauri star, TW Hya was observed with spectroscopic and photometric measurements simultaneous with a long se gmented exposure using the CHANDRA satellite. Contemporaneous optical photometry from WASP-S indicates a 4.74 day period was present during this time. Absence of a similar periodicity in the H-alpha flux and the total X-ray flux points to a different source of photometric variations. The H-alpha emission line appears intrinsically broad and symmetric, and both the profile and its variability suggest an origin in the post-shock cooling region. An accretion event, signaled by soft X-rays, is traced spectroscopically for the first time through the optical emission line profiles. After the accretion event, downflowing turbulent material observed in the H-alpha and H-beta lines is followed by He I (5876A) broadening. Optical veiling increases with a delay of about 2 hours after the X-ray accretion event. The response of the stellar coronal emission to an increase in the veiling follows about 2.4 hours later, giving direct evidence that the stellar corona is heated in part by accretion. Subsequently, the stellar wind becomes re-established. We suggest a model that incorporates this sequential series of events: an accretion shock, a cooling downflow in a supersonically turbulent region, followed by photospheric and later, coronal heating. This model naturally explains the presence of broad optical and ultraviolet lines, and affects the mass accretion rates determined from emission line profiles.

The 2008 optical transient in NGC 300 is one of a growing class of intermediate-luminosity transients that brighten several orders of magnitude from a previously optically obscured state. The origin of their eruptions is not understood. Our multi-wavelength photometry and spectroscopy from maximum light to more than a year later provide a record of its post-eruption behavior. We describe its changing spectral-energy distribution, the evolution of its absorption- and emission line-spectrum, the development ofa bipolar outflow, and the rapid transition from a dense wind to an optically thin ionized wind. In addition to strong, narrow hydrogen lines, the F-type absorption-line spectrum of the transient is characterized by strong Ca II and [Ca II] emission. The very broad wings of the Ca II triplet and the asymmetric [Ca II] emission lines are due to strong Thomson scattering in the expanding ejecta. Post-maximum, the hydrogen and Ca II lines developed double-peaked emission profiles that we attribute to a bipolar outflow. Between approximately 60 and 100 days after maximum, the F-type absorption spectrum, formed in its dense wind, weakened and the wind became transparent to ionizing radiation. We discuss the probable evolutionary state of the transient and similar objects such as SN 2008S, and conclude thatthey were most likely post-red supergiants or post-AGB stars on a blue loop to warmer temperatures when the eruption occurred. These objects are not LBVs.

I present an overview of the techniques used for detecting and following up binaries in nearby galaxies and present the current census of extragalactic binaries, with a focus on eclipsing systems. The motivation for looking in other galaxies is the use of eclipsing binaries as distance indicators and as probes of the most massive stars.

We present preliminary results of the first near-infrared variability study of the Arches cluster, using adaptive optics data from NIRI/Gemini and NACO/VLT. The goal is to discover eclipsing binaries in this young (2.5 $\pm$ 0.5 Myr), dense, massive cluster for which we will determine accurate fundamental parameters with subsequent spectroscopy. Given that the Arches cluster contains more than 200 Wolf-Rayet and O-type stars, it provides a rare opportunity to determine parameters for some of the most massive stars in the Galaxy.

We present fundamental parameters of 4 massive eclipsing binaries in the young massive cluster Westerlund 1. The goal is to measure accurate masses and radii of their component stars, which provide much needed constraints for evolutionary models of massive stars. Accurate parameters can further be used to determine a dynamical lower limit for the magnetar progenitor and to obtain an independent distance to the cluster. Our results confirm and extend the evidence for a high mass for the progenitor of the magnetar.

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers.

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code CMFGEN, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

We discuss the nature of the brightest mid-IR point source (which we dub Object X) in the nearby galaxy M33. Although multi-wavelength data on this object have existed in the literature for some time, it has not previously been recognized as the most luminous mid-IR object in M33 because it is entirely unremarkable in both optical and near-IR light. In the Local Group Galaxies Survey, Object X is a faint red source visible in VRI and H-alpha but not U or B. It was easily seen at JHK_s in the 2MASS survey. It is the brightest point source in all four Spitzer IRAC bands and is also visible in the MIPS 24-micron band. Its bolometric luminosity is 5x10^5 L_sun. The source is optically variable on short time scales (tens of days) and is also slightly variable in the mid-IR, indicating that it is a star. Archival photographic plates (from 1949 and 1991) show no optical source, so the star has been obscured for at least half a century. Its properties are similar to those of the Galactic OH/IR star IRC+10420 which has a complex dusty circumstellar structure resulting from episodic low velocity mass ejections. We propose that Object X is a M>30 M_sun evolved star obscured in its own dust ejected during episodic mass loss events over at least half a century. It may emerge from its current ultra-short evolutionary phase as a hotter post-RSG star analogous to M33 Var A. The existence and rarity of such objects can be an important probe of a very brief yet eventful stellar evolutionary phase.

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

This paper presents the first distance measurement to the massive, semi-detached, eclipsing binary LMC-SC1-105, located in the LH 81 association of the Large Magellanic Cloud (LMC). Previously determined parameters of the system are combined with new near-infrared photometry and a new temperature analysis to constrain the reddening toward the system, and determine a distance of 50.6+/-1.6 kpc (corresponding to a distance modulus of 18.52+/-0.07 mag), in agreement with previous eclipsing binary measurements. This is the sixth distance measurement to an eclipsing binary in the LMC, although the first to an O-type system. We thus demonstrate the suitability of O-type eclipsing binaries (EBs) as distance indicators. We suggest using bright, early-type EBs to measure distances along different sight lines, as an independent way to map the depth of the LMC and resolve the controversy about its three-dimensional structure.

We present a study of the infrared properties of 4922 spectroscopically confirmed massive stars in the Large and Small Magellanic Clouds, focusing on the active OB star population. Besides OB stars, our sample includes yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables (LBVs) and supergiant B[e] stars. We detect a distinct Be star sequence, displaced to the red, and find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We also find photometric variability among the active OB population and evidence for transitions of Be stars to B stars and vice versa. We furthermore confirm the presence of dust around all the supergiant B[e] stars in our sample, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs.

We present results of our study of the infrared properties of massive stars in the Large and Small Magellanic Clouds, which are based on the Spitzer SAGE surveys of these galaxies. We have compiled catalogs of spectroscopically confirmed massive stars in each galaxy, as well as photometric catalogs for a subset of these stars that have infrared counterparts in the SAGE database, with uniform photometry from 0.3 to 24 microns in the UBVIJHKs+IRAC+MIPS24 bands. These catalogs enable a comparative study of infrared excesses of OB stars, classical Be stars, yellow and red supergiants, Wolf-Rayet stars, Luminous Blue Variables and supergiant B[e] stars, as a function of metallicity, and provide the first roadmaps for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths.

We present preliminary results of the first near-infrared variability study of the Arches cluster, using adaptive optics data from NIRI/Gemini and NACO/VLT. The goal is to discover eclipsing binaries in this young (2.5 +- 0.5 Myr), dense, massive cluster for which we will determine accurate fundamental parameters with subsequent spectroscopy. Given that the Arches cluster contains more than 200 Wolf-Rayet and O-type stars, it provides a rare opportunity to determine parameters for some of the most massive stars in the Galaxy.

Westerlund 1 is one of the most massive young clusters known in the Local Group, with an age of 3-5 Myr. It contains an assortment of rare evolved massive stars, such as blue, yellow and red supergiants, Wolf-Rayet stars, a luminous blue variable, and a magnetar, as well as 4 massive eclipsing binary systems (Wddeb, Wd13, Wd36, WR77o, see Bonanos 2007). The eclipsing binaries present a rare opportunity to constrain evolutionary models of massive stars, the distance to the cluster and furthermore, to determine a dynamical lower limit for the mass of a magnetar progenitor. Wddeb, being a detached system, is of great interest as it allows determination of the masses of 2 of the most massive unevolved stars in the cluster. We have analyzed spectra of all 4 eclipsing binaries, taken in 2007-2008 with the 6.5 meter Magellan telescope at Las Campanas Observatory, Chile, and present fundamental parameters (masses, radii) for their component stars.

Motivated by the detection of a recent outburst of the massive luminous blue variable LMC-R71, which reached an absolute magnitude M_V = -9.3 mag, we undertook a systematic study of the optical variability of 1268 massive stars in the Large Magellanic Cloud, using a recent catalog by Bonanos et al. (2009) as the input. The ASAS All Star Catalog (Pojmanski 2002) provided well-sampled light curves of these bright stars spanning 10 years. Combining the two catalogs resulted in 599 matches, on which we performed a variability search. We identified 117 variable stars, 38 of which were not known before, despite their brightness and large amplitude of variation. We found 13 periodic stars that we classify as eclipsing binary (EB) stars, eight of which are newly discovered bright, massive eclipsing binaries composed of OB type stars. The remaining 104 variables are either semi- or non-periodic, the majority (85) being red supergiants. Most (26) of the newly discovered variables in this category are also red supergiants with only three B and four O stars.

We present a catalog of 5324 massive stars in the Small Magellanic Cloud (SMC), with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 3654 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer, SAGE-SMC survey database, for which we present uniform photometry from 0.3-24 um in the UBVIJHKs+IRAC+MIPS24 bands. We compare the color magnitude diagrams and color-color diagrams to those of the Large Magellanic Cloud (LMC), finding that the brightest infrared sources in the SMC are also the red supergiants, supergiant B[e] (sgB[e]) stars, luminous blue variables, and Wolf-Rayet stars, with the latter exhibiting less infrared excess, the red supergiants being less dusty and the sgB[e] stars being on average less luminous. Among the objects detected at 24 um are a few very luminous hypergiants, 4 B-type stars with peculiar, flat spectral energy distributions, and all 3 known luminous blue variables. We detect a distinct Be star sequence, displaced to the red, and suggest a novel method of confirming Be star candidates photometrically. We find a higher fraction of Oe and Be stars among O and early-B stars in the SMC, respectively, when compared to the LMC, and that the SMC Be stars occur at higher luminosities. We estimate mass-loss rates for the red supergiants, confirming the correlation with luminosity even at the metallicity of the SMC. Finally, we confirm the new class of stars displaying composite A & F type spectra, the sgB[e] nature of 2dFS1804 and find the F0 supergiant 2dFS3528 to be a candidate luminous blue variable with cold dust.

In this paper, we present for the first time a proper-motion-selected white dwarf (WD) sample of the old Galactic open cluster M 67, down to the bottom of the WD cooling sequence (CS). The color-magnitude diagram is based on data collected with the LBC-Blue camera at the prime focus of LBT. As first epoch data, we used CFHT-archive images collected 10 years before LBC data. We measured proper motions of all the identified sources. Proper motions are then used to separate foreground and background objects from the cluster stars, including WDs. Finally, the field-object cleaned WD CS in the V vs. B-I color-magnitude diagram is compared with the models. We confirm that the age derived from the location of the bottom of the WD CS is consistent with the turn off age.

We present a catalog of 1750 massive stars in the Large Magellanic Cloud, with accurate spectral types compiled from the literature, and a photometric catalog for a subset of 1268 of these stars, with the goal of exploring their infrared properties. The photometric catalog consists of stars with infrared counterparts in the Spitzer SAGE survey database, for which we present uniform photometry from 0.3-24 microns in the UBVIJHKs+IRAC+MIPS24 bands. The resulting infrared color-magnitude diagrams illustrate that the supergiant B[e], red supergiant and luminous blue variable (LBV) stars are among the brightest infrared point sources in the Large Magellanic Cloud, due to their intrinsic brightness, and at longer wavelengths, due to dust. We detect infrared excesses due to free-free emission among ~900 OB stars, which correlate with luminosity class. We confirm the presence of dust around 10 supergiant B[e] stars, finding the shape of their spectral energy distributions (SEDs) to be very similar, in contrast to the variety of SED shapes among the spectrally variable LBVs. The similar luminosities of B[e] supergiants (log L/Lo>=4) and the rare, dusty progenitors of the new class of optical transients (e.g. SN 2008S and NGC 300 OT), plus the fact that dust is present in both types of objects, suggests a common origin for them. We find the infrared colors for Wolf-Rayet stars to be independent of spectral type and their SEDs to be flatter than what models predict. The results of this study provide the first comprehensive roadmap for interpreting luminous, massive, resolved stellar populations in nearby galaxies at infrared wavelengths.

A luminous optical transient (OT) that appeared in NGC 300 in early 2008 had a maximum brightness, M_V ~ -12 to -13, intermediate between classical novae and supernovae. We present ground-based photometric and spectroscopic monitoring and adaptive-optics imaging of the OT, as well as pre- and post-outburst space-based imaging with HST and Spitzer. The optical spectrum at maximum showed an F-type supergiant photosphere with superposed emission lines of hydrogen, Ca II, and [Ca II], similar to the spectra of low-luminosity Type IIn "supernova impostors" like SN 2008S, as well as cool hypergiants like IRC +10420. The emission lines have a complex, double structure, indicating a bipolar outflow with velocities of ~75 km/s. The luminous energy released in the eruption was ~10^47 ergs, most of it emitted in the first 2 months. By registering new HST images with deep archival frames, we have precisely located the OT site, and find no detectable optical progenitor brighter than broad-band V magnitude 28.5. However, archival Spitzer images reveal a bright, non-variable mid-IR pre-outburst source. We conclude that the NGC 300 OT was a heavily dust-enshrouded luminous star, of ~10-15 Msun, which experienced an eruption that cleared the surrounding dust and initiated a bipolar wind. The progenitor was likely an OH/IR source which had begun to evolve on a blue loop toward higher temperatures, but the precise cause of the outburst remains uncertain.

Do some Wolf-Rayet stars owe their strong winds to something else besides radiation pressure? The answer to this question is still not entirely obvious, especially in certain Wolf-Rayet subclasses, mainly WN8 and WC9. Both of these types of Wolf-Rayet stars are thought to be highly variable, as suggested by observations, possibly due to pulsations. However, only the WN8 stars have so far been vigorously and systematically investigated for variability. We present here the results of a systematic survey during 3 consecutive weeks of 19 Galactic WC9 stars and 1 WC8 star for photometric variability in two optical bands, V and I. Of particular interest are the correlated variations in brightness and colour index in the context of carbon-dust formation, which occurs frequently in WC9 and some WC8 stars. In the most variable case, WR76, we used this information to derive a typical dust grain size of ~ 0.1 um. However, most photometric variations occur at surprisingly low levels and in fact almost half of our sample shows no significant variability at all above the instrumental level (sigma ~ 0.005 - 0.01 mag).

This paper presents a photometric and spectroscopic study of the bright blue eclipsing binary LMC-SC1-105, selected from the OGLE catalog as a candidate host of very massive stars (>=30Mo). The system is found to be a double-lined spectroscopic binary, which indeed contains massive stars. The masses and radii of the components are M1= 30.9+/-1.0 Mo, M2= 13.0+/-0.7 Mo, and R1= 15.1+/-0.2 Ro, R2= 11.9+/-0.2 Ro, respectively. The less massive star is found to be filling its Roche lobe, indicating the system has undergone mass-transfer. The spectra of LMC-SC1-105 display the Struve-Sahade effect, with the HeI lines of the secondary appearing stronger when it is receding and causing the spectral types to change with phase (O8+O8 to O7+O8.5). This effect could be related to the mass-transfer in this system. To date, accurate (<=10%) fundamental parameters have only been measured for 15 stars with masses greater than 30 Mo, with the reported measurements contributing valuable data on the fundamental parameters of very massive stars at low metallicity. The results of this work demonstrate that the strategy of targeting the brightest blue stars in eclipsing binaries is an effective way of studying very massive stars.

We measure the projected rotational velocities of the late B-type hypervelocity stars HVS7 and HVS8 from high resolution spectroscopy to be 60 +/- 17 km/s and 260 +/- 70 km/s. The 'slow' rotation of HVS7 is in principle consistent with having originated in a binary system, assuming a high inclination angle of the stellar rotation axis. However, the fast rotation of HVS8 is more typical of single B-type stars. HVS8 could have therefore been ejected by a mechanism other than that proposed by Hills. We also estimate the effective temperatures and surface gravities for HVS7 and HVS8 and obtain an additional measurement of their radial velocities. We find evidence in support of a blue horizontal branch nature for HVS7, and a main sequence nature for HVS8.

We measure the metallicity of the unusual hypervelocity star HE 0437-5439 from high resolution spectroscopy to be half-solar. We determine a spectral type of B2 IV-III for the star and derive an effective temperature Teff= 21,500 +/- 1,000 K and a surface gravity log(g) = 3.7 +/- 0.2 (cgs). We also present BV time series photometry and find the star to be non-variable at the 0.02 mag level. We refine the magnitude of the hypervelocity star to V=16.36+/-0.04 mag, with a color B-V=-0.23+/-0.03 mag, confirming its early-type nature. Our metallicity result establishes the origin of HE 0437-5439 in the Large Magellanic Cloud and implies the existence of a massive black hole somewhere in this galaxy.

The initial-final mass relationship of white dwarfs, which is poorly constrained, is of paramount importance for different aspects in modern astrophysics. From an observational perspective, most of the studies up to now have been done using white dwarfs in open clusters. In order to improve the initial-final mass relationship we explore the possibility of deriving a semi-empirical relation studying white dwarfs in common proper motion pairs. We have acquired long-slit spectra of the white dwarf members of the selected common proper motion pairs, as well as high resolution spectra of their companions. From these observations, a full analysis of the two members of each common proper motion pair has lead to the initial and final masses of the white dwarfs. These observations have allowed us to provide updated information for the white dwarfs, since some of them were misclassified. This work is the first one in using common proper motion pairs to improve the initial-final mass relationship, and has also allowed to cover the poorly explored low-mass domain. As in the case of studies based on white dwarfs in open clusters, the distribution of the semi-empirical data presents a large scatter, which is higher than the expected uncertainties in the derived values. This suggests that the initial-final mass relationship may not be a single-valued function.

This paper presents the first optical variability study of the Westerlund 1 super star cluster in search of massive eclipsing binary systems. A total of 129 new variable stars have been identified, including the discovery of 4 eclipsing binaries that are cluster members, 1 additional candidate, 8 field binaries, 19 field delta Scuti stars, 3 field W UMa eclipsing binaries, 13 other periodic variables and 81 long period or non-periodic variables. These include the known luminous blue variable, the B[e] star, 11 Wolf-Rayet stars, several supergiants, and other reddened stars that are likely members of Westerlund 1. The bright X-ray source corresponding to the Wolf-Rayet star WR77o (B) is found to be a 3.51 day eclipsing binary. The discovery of a reddened detached eclipsing binary system implies the first identification of main-sequence stars in Westerlund 1.

We present a re-analysis of our H- and K-band photometry and light-curves for GCIRS 16SW, a regular periodic source near the Galactic center. These data include those presented by DePoy et al. (2004); we correct a sign error in their reduction, finding GCIRS 16SW to be an eclipsing binary with no color variations. We find the system to be an equal mass overcontact binary (both stars overfilling their Roche lobes) in a circular orbit with a period P=19.4513 days, an inclination angle i=71 degrees. This confirms and strengthens the findings of Martins et al. (2006) that GCIRS 16SW is an eclipsing binary composed of two ~50Msun stars, further supporting evidence of recent star formation very close to the Galactic center. Finally, the calculated luminosity of each component is close to the Eddington luminosity, implying that the temperature of 24400 K given by Najarro et al. (1997) might be overestimated for these evolved stars.

In the last decade, over 7000 eclipsing binaries have been discovered in the Local Group through various variability surveys. Measuring fundamental parameters of these eclipsing binaries has become feasible with 8 meter class telescopes, making it possible to use eclipsing binaries as distance indicators. Distances with eclipsing binaries provide an independent method for calibrating the extragalactic distance scale and thus determining the Hubble constant. This method has been used for determining distances to eclipsing binaries in the Magellanic Clouds and the Andromeda Galaxy and most recently to a detached eclipsing binary in the Triangulum Galaxy by the DIRECT Project. The increasing number of eclipsing binaries found by microlensing and variability surveys also provide a rich database for advancing our understanding of star formation and evolution.

We present the first direct distance determination to a detached eclipsing binary in M33, which was found by the DIRECT Project. Located in the OB 66 association at coordinates (alpha, delta)=(01:33:46.17,+30:44:39.9) for J2000.0, it was one of the most suitable detached eclipsing binaries found by DIRECT for distance determination, given its apparent magnitude and orbital period. We obtained follow-up BV time series photometry, JHKs photometry and optical spectroscopy from which we determined the parameters of the system. It contains two O7 main sequence stars with masses of 33.4+/-3.5 Mo and 30.0+/-3.3 Mo and radii of 12.3+/-0.4 Ro and 8.8+/-0.3 Ro, respectively. We derive temperatures of 37000+/-1500 K and 35600+/-1500 K. Using BVRJHKs photometry for the flux calibration, we obtain a distance modulus of 24.92+/-0.12 mag (964+/-54 kpc), which is ~0.3 mag longer than the Key Project distance to M33. We discuss the implications of our result and the importance of establishing M33 as an independent rung on the cosmological distance ladder.

We present near-infrared $JHK_s$ photometry and light curves of the candidate magnetic white dwarf+brown dwarf binary SDSS J121209.31+013627.7 and report on the detection of near-infrared excess and variability in the $K_s-$band. The observed near-infrared excess can be explained by the presence of an L7 brown dwarf and an extra emission source. The $J$ and $H$ light curves appear flat, which rules out eclipses deeper than 0.2 mag and the presence of an accretion hot spot on the white dwarf. From the variable $K_s$ lightcurve, we get a refined period for the system of 88$\pm$1 minutes. We show that the observed variability in $K_s-$band can be explained by cyclotron emission, which can be modeled by a small spot on the surface of the white dwarf. SDSS 1212 exhibits similarities to the ultra-short period polar EF Eridani, however the lack of evidence for Roche-lobe overflow accretion suggests it may be a pre-polar.

We present deep optical photometry of the afterglow of gamma-ray burst (GRB) 041006 and its associated hypernova obtained over 65 days after detection (55 R-band epochs on 10 different nights). Our early data (t<4 days) joined with published GCN data indicates a steepening decay, approaching F_nu ~t^{-0.6} at early times (<<1 day) and F_nu ~t^{-1.3} at late times. The break at t_b=0.16+-0.04 days is the earliest reported jet break among all GRB afterglows. During our first night, we obtained 39 exposures spanning 2.15 hours from 0.62 to 0.71 days after the burst that reveal a smooth afterglow, with an rms deviation of 0.024 mag from the local power-law fit, consistent with photometric errors. After t~4 days, the decay slows considerably, and the light curve remains approximately flat at R~24 mag for a month before decaying by another magnitude to reach R~25 mag two months after the burst. This ``bump'' is well-fitted by a k-corrected light curve of SN1998bw, but only if stretched by a factor of 1.38 in time. In comparison with the other GRB-related SNe bumps, GRB 041006 stakes out new parameter space for GRB/SNe, with a very bright and significantly stretched late-time SN light curve. Within a small sample of fairly well observed GRB/SN bumps, we see a hint of a possible correlation between their peak luminosity and their ``stretch factor'', broadly similar to the well-studied Phillips relation for the type Ia supernovae.

We present UBVRIJHK photometry and optical spectroscopy of the so-called 'peculiar' Type Ia supernova 1999by in NGC 2841. The observations began one week before visual maximum light which is well-defined by daily observations. The light curves and spectra are similar to those of the prototypical subluminous event SN 1991bg. We find that maximum light in B occurred on 1999 May 10.3 UT (JD 2,451,308.8 +/- 0.3) with B=13.66 +/- 0.02 mag and a color of B_max-V_max=0.51 +/- 0.03 mag. The late-time color implies minimal dust extinction from the host galaxy. Our photometry, when combined with the recent Cepheid distance to NGC 2841 (Macri et al. 2001), gives a peak absolute magnitude of M_B=-17.15 +/- 0.23 mag, making SN 1999by one of the least luminous Type Ia events ever observed. We estimate a decline rate parameter of dm15(B)=1.90 mag, versus 1.93 for SN 1991bg, where 1.10 is typical for so-called 'normal' events. We compare SN 1999by with other subluminous events and find that the B_max-V_max color correlates strongly with the decline rate and may be a more sensitive indicator of luminosity than the fading rate for these objects. We find a good correlation between luminosity and the depth of the spectral feature at 580 nm, which had been attributed solely to Si II. We show that in cooler photospheres the 580 nm feature is dominated by Ti II, which provides a simple physical explanation for the correlation. Using only subluminous Type Ia supernovae we derive a Hubble parameter of H_0=75 +12 -11 km/s Mpc, consistent with values found from brighter events.

We present a high-precision I-band light curve for the Wolf-Rayet binary WR 20a, obtained as a sub-project of the Optical Gravitational Lensing Experiment. Rauw et al. have recently presented spectroscopy for this system, strongly suggesting extremely large minimum masses of 70.7 +/- 4.0 Mo and 68.8 +/- 3.8 Mo for the component stars of the system, with the exact values depending strongly on the period of the system. We detect deep eclipses of about 0.4 mag in the light curve of WR 20a, confirming and refining the suspected period of P=3.686 days and deriving an inclination angle of i=74.5 +/- 2 deg. Using these photometric data and the radial velocity data of Rauw et al., we derive the masses for the two components of WR 20a to be 83.0 +/- 5.0 Mo and 82.0 +/- 5.0 Mo. Therefore, WR 20a is confirmed to consist of two extremely massive stars and to be the most massive binary known with an accurate mass determination.

Accurate masses and, in general, all fundamental parameters of distant stars can only be measured in eclipsing binaries. Several massive star candidates with masses near 200 Mo exist, however they have large uncertainties associated with them. The most massive binary ever measured accurately is WR 20a, for which we present the light curve. Measuring the period and inclination, we derive masses greater than 80 Mo for each component. Massive binaries are bound to exist in Local Group galaxies, such as M31 and M33. These can be selected from their light curves, obtained by variability studies, such as the DIRECT project. We present photometry and spectroscopy of the detached system M33A, for which we are obtaining a direct distance determination. The DIRECT project has detected several candidate massive binaries which are brighter but non-detached systems, perhaps similar to WR 20a. We plan to obtain spectra for them and measure their masses.

We have undertaken several projects with the purpose of determining accurate distances to nearby galaxies to calibrate the extragalactic distance scale. Specifically, I describe the DIRECT project which aims to derive the distance to M31 and M33 directly, using detached eclipsing binaries and the Baade-Wesselink method for Cepheids. I also present a ``hybrid'' method of discovering Cepheids with ground-based telescopes using image subtraction and then following them up with the HST to derive Cepheid period-luminosity distances.

We present the first CCD variability study of the Draco dwarf spheroidal galaxy. The data were obtained with the FLWO 1.2 m telescope on 22 nights, over a period of 10 months, covering a 22 x 22 arcmin field centered at RA=17:19:57.5, Dec=57:50:05, J2000.0. The analysis of the BVI images produced 163 variable stars, 146 of which were RR Lyrae: 123 RRab, 16 RRc, 6 RRd and one RR12. The other variables include a SX Phe star, four anomalous Cepheids and a field eclipsing binary. Using the short distance scale statistical parallax calibration of Gould & Popowski and 94 RRab stars from our field, we obtain a distance modulus of (m-M)_0=19.40 +/- 0.02 (stat) +/- 0.15 (syst) mag for Draco, corresponding to a distance of 75.7 kpc +/- 0.7 (stat) +/- 5.4 (syst) kpc. By comparing the spread in magnitudes of RRab stars in B,V and I, we find no evidence for internal dust in the Draco dwarf spheroidal galaxy. The catalog of all variables, as well as their photometry and finding charts, is available electronically via anonymous ftp and the World Wide Web. The complete set of the CCD frames is available upon request.

We apply the image subtraction method to re-analyze the ESO Very Large Telescope data on M83 (NGC 5236), obtained and analyzed by Thim et al. Whereas Thim et al. found 12 Cepheids with periods between 12-55 days, we find 112 Cepheids with periods ranging from 7-91 days, as well as ~60 other variables. These include 2 candidate eclipsing binaries, which, if confirmed, would be the first optically discovered outside the Local Group. We thus demonstrate that the image subtraction method is much more powerful for detecting variability, especially in crowded fields. However, HST observations are necessary to obtain a Cepheid period-luminosity distance not dominated by blending and crowding. We propose a ``hybrid'' approach, where numerous Cepheids are discovered and characterized using large ground-based telescopes and then followed-up with the HST to obtain precise distances.

The DIRECT Project aims to obtain direct distances to two Local Group galaxies, M31 and M33, which occupy a crucial position near the bottom of the cosmological distance ladder. The first step is to search for detached eclipsing binaries (DEBs) and Cepheids using 1 m class telescopes to select good candidates, which will be followed up spectroscopically on 6.5-10 m class telescopes. In this ninth paper we present a catalog of variable stars discovered with image subtraction in field M31Y [RA=10.97 deg, Dec=41.69 deg, J2000.0]. The data were obtained with the FLWO 1.2 m telescope on 25 nights, over a period of 6 months. In our search covering 22x22 arcmin we discovered 41 eclipsing binaries, 126 Cepheids, and 97 other periodic or nonperiodic variables, including a LBV candidate, a nova and a galactic CV. The catalog of variables, as well as their photometry and finding charts, is available electronically via anonymous ftp and the World Wide Web. The complete set of the CCD frames is available upon request.

We describe optical photometric and spectroscopic observations of the bright variable BG Geminorum. Optical photometry shows a pronounced ellipsoidal variation of the K0 I secondary, with amplitudes of ~0.5 mag at VRI and a period of 91.645 days. A deep primary eclipse is visible for wavelengths < 4400A; a shallower secondary eclipse is present at longer wavelengths. Eclipse timings and the radial velocity curve of the K0 secondary star indicate an interacting binary where a lobe-filling secondary, M_2 ~ 0.5 Msun, transfers material into a extended disk around a massive primary, M_1 ~ 4.5 Msun. The primary star is either an early B-type star or a black hole. If it did contain a black hole, BG Gem would be the longest period black hole binary known by a factor of 10, as well as the only eclipsing black hole binary system.