Time-domain studies of pre-main sequence stars have long been used to investigate star properties during their early evolutionary phases and to trace the evolution of circumstellar environments. Historically these studies have been confined to the nearest, low-density, star forming regions. We used the Wide Field Camera 3 on board of the Hubble Space Telescope to extend, for the first time, the study of pre-main sequence variability to one of the few young massive clusters in the Milky Way, Westerlund 2. Our analysis reveals that at least 1/3 of the intermediate and low-mass pre-main sequence stars in Westerlund 2 are variable. Based on the characteristics of their light curves, we classified ~11% of the variable stars as weak-line T-Tauri candidates, ~ 52% as classical T-Tauri candidates, ~ 5% as dippers and ~26% as bursters. In addition, we found that 2% of the stars below 6Mo (~6% of the variables) are eclipsing binaries, with orbital periods shorter than 80 days. The spatial distribution of the different populations of variable pre-main sequence stars suggests that stellar feedback and UV-radiation from massive stars play an important role on the evolution of circumstellar and planetary disks.

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)

We mapped the Galactic young massive star cluster Westerlund 2 (Wd2) with the integral field spectrograph MUSE (spatial resolution: 0.2arcsec/px, spectral resolution: $\Delta \lambda$ = 1.25A, wavelength range 4600-9350A) mounted on the VLT, as part of an on-going study to measure the stellar and gas kinematics of the cluster region. In this paper we present the fully reduced dataset and introduce our new Python package "MUSEpack", which we developed to measure stellar radial velocities with an absolute precision of 1-2km/s without the necessity of a spectral template library. This novel method uses the two-dimensional spectra and an atomic transition line library to create templates around strong absorption lines for each individual star. The code runs fully automatically on multi-core machines, which makes it possible to efficiently determine stellar radial velocities of a large number of stars with the necessary precision to measure the velocity dispersion of young star clusters. MUSEpack also provides an enhanced method for removing telluric lines in crowded fields without sky exposures and a Python wrapper for ESO's data reduction pipeline. We observed Wd2 with a total of 11 short and 5 long exposures to cover the bright nebular emission and OB stars, as well as the fainter pre-main sequence stars down to ~1M$_\odot$. The survey covers an area of ~11arcmin$^2$ (15.8pc$^2$). In total, we extracted 1,725 stellar spectra with a mean S/N>5 per pixel. A typical radial velocity (RV) uncertainty of 4.78km/s, 2.92km/s, and 1.1km/s is reached for stars with a mean S/N>10, S/N>20, S/N>50 per pixel, respectively. Depending on the number of spectral lines used to measure the RVs, it is possible to reach RV accuracies of 0.9km/s, 1.3km/s, and 2.2km/s with $\geq5$, 3-4, and 1-2 spectral lines, respectively. The combined statistical uncertainty on the radial velocity measurements is 1.10km/s.

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses.

Westerlund 2 (Wd2) is the central ionizing star cluster of the \ion{H}{2} region RCW~49 and the second most massive young star cluster (${\rm M} = (3.6 \pm 0.3)\times 10^4\,{\rm M}_\odot$) in the Milky Way. Its young age ($\sim2\,$Myr) and close proximity to the Sun ($\sim 4\,$kpc) makes it a perfect target to study stars emerging from their parental gas cloud, the large number of OB-stars and their feedback onto the gas, and the gas dynamics. We combine high-resolution multi-band photometry obtained in the optical and near-infrared with the \textit{Hubble} Space Telescope (HST), and VLT/MUSE integral field spectroscopy to study the gas, the stars, and their interactions, simultaneously. In this paper we focus on a small, $64\times64\,{\rm arcsec}^2$ region North of the main cluster center, which we call the Northern Bubble (NB), a circular cavity carved into the gas of the cluster region. Using MUSE data, we determined the spectral types of 17 stars in the NB from G9III to O7.5. With the estimation of these spectral types we add 2 O and 5 B-type stars to the previously published census of 37 OB-stars in Wd2. To measure radial velocities we extracted 72 stellar spectra throughout Wd2, including the 17 of the NB, and show that the cluster member stars follow a bimodal velocity distribution centered around $(8.10 \pm 1.53)\,{\rm km}\,{\rm s}^{-1}$ and $(25.41 \pm 1.57)\,{\rm km}\,{\rm s}^{-1}$ with a dispersion of $(4.52 \pm 1.78)\,{\rm km}\,{\rm s}^{-1}$ and $(3.46 \pm 1.29)\,{\rm km}\,{\rm s}^{-1}$, respectively. These are in agreement with CO($J=1$-2) studies of RCW~49 leaving cloud-cloud collision as a viable option for the formation scenario of Wd2. The bimodal distribution is also detected in the Gaia DR2 proper motions.

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.

We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques, and correlate these maps with neutral HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching.

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.

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in five bands from the near UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 Mo, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the time scale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems.

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.

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 report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.

We present a detailed analysis of the spatial distribution of the stellar population and the present-day mass function (PDMF) of the Westerlund 2 (Wd2) region using the data from our high resolution multi-band survey with the Hubble Space Telescope. We used state-of-the-art artificial star tests to determine spatially resolved completeness maps for each of the broad-band filters. We reach a level of completeness of 50 % down to F555W=24.8 mag (0.7 $M_\odot$) and F814W=23.3 mag (0.2 $M_\odot$) in the optical and F125W=20.2 mag and F160W=19.4 mag (both 0.12 $M_\odot$) in the infrared throughout the field of view. We had previously reported that the core of Wd2 consists of two clumps: namely the main cluster (MC) and the northern clump (NC). From the spatial distribution of the completeness corrected population, we find that their stellar surface densities are 1114 stars pc$^{-2}$ and 555 stars pc$^{-2}$, respectively, down to F814W=21.8 mag. We find that the present-day mass function (PDMF) of Wd2 has a slope of $\Gamma=-1.46 \pm 0.06$, which translates to a total stellar cluster mass of $(3.6 \pm 0.3) \cdot 10^4 M_\odot$. The spatial analysis of the PDMF reveals that the cluster population is mass-segregated, most likely primordial. In addition, we report the detection of a stellar population of spatially uniformly distributed low-mass (<0.15 $M_\odot$) stars, extending into the gas ridges of the surrounding gas and dust cloud, as well as a confined region of reddened stars, likely caused by a foreground CO cloud. We find hints that a cloud-cloud collision might be the origin of the formation of Wd2.

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.

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

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 detailed clustering analysis of the young stellar population across the star-forming ring galaxy NGC 6503, based on the deep HST photometry obtained with the Legacy ExtraGalactic UV Survey (LEGUS). We apply a contour-based map analysis technique and identify in the stellar surface density map 244 distinct star-forming structures at various levels of significance. These stellar complexes are found to be organized in a hierarchical fashion with 95% being members of three dominant super-structures located along the star-forming ring. The size distribution of the identified structures and the correlation between their radii and numbers of stellar members show power-law behaviors, as expected from scale-free processes. The self-similar distribution of young stars is further quantified from their autocorrelation function, with a fractal dimension of ~1.7 for length-scales between ~20 pc and 2.5 kpc. The young stellar radial distribution sets the extent of the star-forming ring at radial distances between 1 and 2.5 kpc. About 60% of the young stars belong to the detected stellar structures, while the remaining stars are distributed among the complexes, still inside the ring of the galaxy. The analysis of the time-dependent clustering of young populations shows a significant change from a more clustered to a more distributed behavior in a time-scale of ~60 Myr. The observed hierarchy in stellar clustering is consistent with star formation being regulated by turbulence across the ring. The rotational velocity difference between the edges of the ring suggests shear as the driving mechanism for this process. Our findings reveal the interesting case of an inner ring forming stars in a hierarchical fashion.

The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community.

STEP (the SMC in Time: Evolution of a Prototype interacting late-type dwarf galaxy) is a Guaranteed Time Observation survey being performed at the VST (the ESO VLT Survey Telescope). STEP will image an area of 74 deg$^2$ covering the main body of the Small Magellanic Cloud (32 deg$^2$), the Bridge that connects it to the Large Magellanic Cloud (30 deg$^2$) and a small part of the Magellanic Stream (2 deg$^2$). Our $g,r,i,H_{\alpha}$ photometry is able to resolve individual stars down to magnitudes well below the main-sequence turnoff of the oldest populations. In this first paper we describe the observing strategy, the photometric techniques, and the upcoming data products of the survey. We also present preliminary results for the first two fields for which data acquisition is completed, including some detailed analysis of the two stellar clusters IC\,1624 and NGC\,419.

We derive the star formation history in four regions of the Small Magellanic Cloud (SMC) using the deepest VI color-magnitude diagrams (CMDs) ever obtained for this galaxy. The images were obtained with the Advanced Camera for Surveys onboard the Hubble Space Telescope and are located at projected distances of 0.5-2 degrees from the SMC center, probing the main body and the wing of the galaxy. We derived the star-formation histories (SFH) of the four fields using two independent procedures to fit synthetic CMDs to the data. We compare the SFHs derived here with our earlier results for the SMC bar to create a deep pencil-beam survey of the global history of the central SMC. We find in all the six fields observed with HST a slow star formation pace from 13 to 5-7 Gyr ago, followed by a ~ 2-3 times higher activity. This is remarkable because dynamical models do not predict a strong influence of either the LMC or the Milky Way (MW) at that time. The level of the intermediate-age SFR enhancement systematically increases towards the center, resulting in a gradient in the mean age of the population, with the bar fields being systematically younger than the outer ones. Star formation over the most recent 500 Myr is strongly concentrated in the bar, the only exception being the area of the SMC wing. The strong current activity of the latter is likely driven by interaction with the LMC. At a given age, there is no significant difference in metallicity between the inner and outer fields, implying that metals are well mixed throughout the SMC. The age-metallicity relations we infer from our best fitting models are monotonically increasing with time, with no evidence of dips. This may argue against the major merger scenario proposed by Tsujimoto and Bekki 2009, although a minor merger cannot be ruled out.

The Bar is the most productive region of the Small Magellanic Cloud in terms of star formation but also the least studied one. In this paper we investigate the star formation history of two fields located in the SW and in the NE portion of the Bar using two independent and well tested procedures applied to the color-magnitude diagrams of their stellar populations resolved by means of deep HST photometry. We find that the Bar experienced a negligible star formation activity in the first few Gyr, followed by a dramatic enhancement from 6 to 4 Gyr ago and a nearly constant activity since then. The two examined fields differ both in the rate of star formation and in the ratio of recent over past activity, but share the very low level of initial activity and its sudden increase around 5 Gyr ago. The striking similarity between the timing of the enhancement and the timing of the major episode in the Large Magellanic Cloud is suggestive of a close encounter triggering star formation.

We use deep images acquired with the Advanced Camera for Surveys (ACS) on board of the Hubble Space Telescope (HST) in the filters F555W and F814W to characterize the properties of NGC 376, a young star cluster located in the wing of the Small Magellanic Cloud (SMC). Using isochrone fitting we derive for NGC 376 an age of 28+/-7 Myr, in good agreement with previous studies. The high spatial resolution ACS data allow us to determine the center of gravity of the cluster and to construct extended surface brightness and radial density profiles. Neither of these profiles can be fitted with a theoretical model, suggesting that the cluster is not in virial equilibrium. Considering the young age of the cluster, we speculate that the distortion of the radial profiles may be the result of the rapid gas dispersal that follows the initial phase of star formation. The cluster shows clear evidence of dynamical mass segregation. From the properties of the radial profiles and the present day mass function (PDMF) we conclude that NGC 376 appears to have already lost nearly 90% of its initial stellar mass, probably as a consequence of the sudden gas dispersal that follows the early phase of star formation (SF).

We present the photometric catalogs for the star-forming cluster NGC 602 in the wing of the Small Magellanic Cloud covering a range of wavelengths from optical HST/ACS (F555W, F814W) and SMARTS/ANDICAM (V, I) to infrared (Spitzer/IRAC 3.6, 4.5, 5.8, and 8 micron and MIPS 24 micron). Combining this with IRSF (InfraRed Survey Facility) near-infrared photometry (J, H, Ks), we compare the young main sequence (MS) and pre-main sequence (PMS) populations prominent in the optical with the current young stellar object (YSO) populations revealed by the infrared (IR). We analyze the MS and PMS population with isochrones in color-magnitude diagrams to derive ages and masses. The optical data reveal ~565 PMS candidates, low mass Stage III YSOs. We characterize ~40 YSOs by fitting their spectral energy distributions (SEDs) to a grid of models (Robitaille et al. 2007) to derive luminosities, masses and evolutionary phase (Stage I-III). The higher resolution HST images reveal that ~70% of the YSO candidates are either multiples or protoclusters. For YSOs and PMS sources found in common, we find a consistency in the masses derived. We use the YSO mass function to derive a present-day star-formation rate of ~0.2-1.0 Msun/yr/kpc^2, similar to the rate derived from the optical star formation history suggesting a constant star formation rate for this region. We demonstrate a progression of star formation from the optical star cluster center to the edge of the star forming dust cloud. We derive lifetimes of a few 10^5 years for the YSO Stages I and II.

We discuss the star formation history of the SMC region NGC 346 based on Hubble Space Telescope images. The region contains both field stars and cluster members. Using a classical synthetic CMD procedure applied to the field around NGC 346 we find that there the star formation pace has been rising from a quite low rate 13 Gyr ago to \approx 1.4 \times 10^{-8} Mo yr^{-1}pc^{-2} in the last 100 Myr. This value is significantly higher than in other star forming regions of the SMC. For NGC 346 itself, we compare theoretical and observed Color-Magnitude Diagrams (CMDs) of several stellar sub-clusters identified in the region, and we derive their basic evolution parameters. We find that NGC 346 experienced different star formation regimes, including a dominant and focused "high density mode", with the sub-clusters hosting both pre-main sequence (PMS) and upper main sequence (UMS) stars, and a diffuse "low density mode", as indicated by the presence of low-mass PMS sub-clusters. Quantitatively, the star formation in the oldest sub-clusters started about 6 Myr ago with remarkable synchronization, it continued at high rate (up to 2 \times 10^{-5} Mo yr^{-1} pc^{-2}) for about 3 Myr and is now progressing at a lower rate. Interestingly, sub-clusters mainly composed by low mass PMS stars seem to experience now the first episode of star formation, following multi-seeded spatial patterns instead of resulting from a coherent trigger. Two speculative scenarios are put forth to explain the deficiency of UMS stars: the first invokes under-threshold conditions of the parent gas; the second speculates that the initial mass function (IMF) is a function of time, with the youngest sub-clusters not having had sufficient time to form more massive stars.

We present a novel approach to derive the age of very young star clusters, by using the Turn-On (TOn). The TOn is the point in the color-magnitude diagram (CMD) where the pre-main sequence (PMS) joins the main sequence (MS). In the MS luminosity function (LF) of the cluster, the TOn is identified as a peak followed by a dip. We propose that by combining the CMD analysis with the monitoring of the spatial distribution of MS stars it is possible to reliably identify the TOn in extragalactic star forming regions. Compared to alternative methods, this technique is complementary to the turn-off dating and avoids the systematic biases affecting the PMS phase. We describe the method and its uncertainties, and apply it to the star forming region NGC346, which has been extensively imaged with the Hubble Space Telescope (HST). This study extends the LF approach in crowded extragalactic regions and opens the way for future studies with HST/WFC3, JWST and from the ground with adaptive optics.

We present structural parameters for the seven intermediate-age and old star clusters NGC121, Lindsay 1, Kron 3, NGC339, NGC416, Lindsay 38, and NGC419 in the Small Magellanic Cloud. We fit King profiles and Elson, Fall, and Freeman profiles to both surface-brightness and star count data taken with the Advanced Camera for Surveys aboard the Hubble Space Telescope. Clusters older than 1 Gyr show a spread in cluster core radii that increases with age, while the youngest clusters have relatively compact cores. No evidence for post core collapse clusters was found. We find no correlation between core radius and distance from the SMC center, although consistent with other studies of dwarf galaxies, some relatively old and massive clusters have low densities. The oldest SMC star cluster, the only globular NGC121, is the most elliptical object of the studied clusters. No correlation is seen between ellipticity and distance from the SMC center. The structures of these massive intermediate-age (1-8 Gyr) SMC star clusters thus appear to primarily result from internal evolutionary processes.

We observed six fields of the Small Magellanic Cloud (SMC) with the Advanced Camera for Survey on board the Hubble Space Telescope in the F555W and F814W filters. These fields sample regions characterized by very different star and gas densities, and, possibly, by different evolutionary histories. We find that the SMC was already forming stars ~12 Gyr ago, even if the lack of a clear horizontal branch suggests that in the first few billion years the star formation activity was low. Within the uncertainties of our two-band photometry, we find evidence of a radial variation in chemical enrichment, with the SMC outskirts characterized by lower metallicity than the central zones. From our CMDs we also infer that the SMC formed stars over a long interval of time until ~2-3 Gyr ago. After a period of modest activity, star formation increased again in the recent past, especially in the bar and the wing of the SMC, where we see an enhancement in the star-formation activity starting from ~500 Myr ago. The inhomogeneous distribution of stars younger than ~100 Myr indicates that recent star formation has mainly developed locally.

Deep HST/ACS photometry of the young cluster NGC 602, located in the remote low density "wing" of the Small Magellanic Cloud, reveals numerous pre-main sequence stars as well as young stars on the main sequence. The resolved stellar content thus provides a basis for studying the star formation history into recent times and constraining several stellar population properties, such as the present day mass function, the initial mass function and the binary fraction. To better characterize the pre-main sequence population, we present a new set of model stellar evolutionary tracks for this evolutionary phase with metallicity appropriate for the Small Magellanic Cloud (Z = 0.004). We use a stellar population synthesis code, which takes into account a full range of stellar evolution phases to derive our best estimate for the star formation history in the region by comparing observed and synthetic color-magnitude diagrams. The derived present day mass function for NGC 602 is consistent with that resulting from the synthetic diagrams. The star formation rate in the region has increased with time on a scale of tens of Myr, reaching $0.3-0.7 \times 10^{-3} M_\odot yr^{-1}$ in the last 2.5 Myr, comparable to what is found in Galactic OB associations. Star formation is most complete in the main cluster but continues at moderate levels in the gas-rich periphery of the nebula.

We introduce the SMC in space and time, a large coordinated space and ground-based program to study star formation processes and history, as well as variable stars, structure, kinematics and chemical evolution of the whole SMC. Here, we present the Colour-Magnitude Diagrams(CMDs) resulting from HST/ACS photometry, aimed at deriving the star formation history (SFH) in six fields of the SMC. The fields are located in the central regions, in the stellar halo, and in the wing toward the LMC. The CMDs are very deep, well beyond the oldest Main Sequence Turn-Off, and will allow us to derive the SFH over the entire Hubble time.

The young star cluster NGC 602 and its associated HII region, N90, formed in a relatively isolated and diffuse environment in the Wing of the Small Magellanic Cloud. Its isolation from other regions of massive star formation and the relatively simple surrounding HI shell structure allows us to constrain the processes that may have led to its formation and to study conditions leading to massive star formation. We use images from Hubble Space Telescope and high resolution echelle spectrographic data from the Anglo-Australian Telescope along with 21-cm neutral hydrogen (HI) spectrum survey data and the shell catalogue derived from it to establish a likely evolutionary scenario leading to the formation of NGC 602. We identify a distinct HI cloud component that is likely the progenitor cloud of the cluster and HII region which probably formed in blister fashion from the cloud's periphery. We also find that the past interaction of HI shells can explain the current location and radial velocity of the nebula. The surrounding Interstellar Medium is diffuse and dust-poor as demonstrated by a low visual optical depth throughout the nebula and an average HI density of the progenitor cloud estimated at 1 cm^-3. These conditions suggest that the NGC 602 star formation event was produced by compression and turbulence associated with HI shell interactions. It therefore represents a single star forming event in a low gas density region.

We present a photometric analysis of the star clusters Lindsay 1, Kron 3, NGC339, NGC416, Lindsay 38, and NGC419 in the Small Magellanic Cloud (SMC), observed with the Hubble Space Telescope Advanced Camera for Surveys (ACS) in the F555W and F814W filters. Our color magnitude diagrams (CMDs) extend ~3.5 mag deeper than the main-sequence turnoff points, deeper than any previous data. Cluster ages were derived using three different isochrone models: Padova, Teramo, and Dartmouth, which are all available in the ACS photometric system. Fitting observed ridgelines for each cluster, we provide a homogeneous and unique set of low-metallicity, single-age fiducial isochrones. The cluster CMDs are best approximated by the Dartmouth isochrones for all clusters, except for NGC419 where the Padova isochrones provided the best fit. The CMD of NGC419 shows several main-sequence turn-offs, which belong to the cluster and to the SMC field. We thus derive an age range of 1.2-1.6 Gyr for NGC419. Interestingly, our intermediate-age star clusters have a metallicity spread of ~0.6 dex, which demonstrates that the SMC does not have a smooth, monotonic age-metallicity relation. We find an indication for centrally concentrated blue straggler star candidates in NGC416, while for the other clusters these are not present. Using the red clump magnitudes, we find that the closest cluster, NGC419 (~50kpc), and the farthest cluster, Lindsay 38 (~67kpc), have a relative distance of ~17kpc, which confirms the large depth of the SMC.

As first Paper of a series devoted to study the old stellar population in clusters and fields in the Small Magellanic Cloud, we present deep observations of NGC121 in the F555W and F814W filters, obtained with the Advanced Camera for Surveys on the Hubble Space Telescope. The resulting color-magnitude diagram reaches ~3.5 mag below the main-sequence turn-off; deeper than any previous data. We derive the age of NGC121 using both absolute and relative age-dating methods. Fitting isochrones in the ACS photometric system to the observed ridge line of NGC121, gives ages of 11.8 +- 0.5 Gyr (Teramo), 11.2 +- 0.5 Gyr (Padova) and 10.5 +- 0.5 Gyr (Dartmouth). The cluster ridge line is best approximated by the alpha-enhanced Dartmouth isochrones. Placing our relative ages on an absolute age scale, we find ages of 10.9 +- 0.5 Gyr (from the magnitude difference between the main-sequence turn-off and the horizontal branch) and 11.5 +- 0.5 Gyr (from the absolute magnitude of the horizontal branch), respectively. These five different age determinations are all lower by 2 - 3 Gyr than the ages of the oldest Galactic globular clusters of comparable metallicity. Therefore we confirm the earlier finding that the oldest globular cluster in the Small Magellanic Cloud, NGC121, is a few Gyr younger than its oldest counterparts in the Milky Way and in other nearby dwarf galaxies such as the Large Magellanic Cloud, Fornax, and Sagittarius. If it were accreted into the Galactic halo, NGC121 would resemble the ''young halo globulars'', although it is not as young as the youngest globular clusters associated with the Sagittarius dwarf. The young age of NGC121 could result from delayed cluster formation in the Small Magellanic Cloud or result from the random survival of only one example of an initially small number star clusters.

Deep F555W and F814W Hubble Space Telescope ACS images are the basis for a study of the present day mass function (PDMF) of NGC346, the largest active star forming region in the Small Magellanic Cloud (SMC). We find a PDMF slope of Gamma=-1.43+/-0.18 in the mass range 0.8-60 Mo, in excellent agreement with the Salpeter Initial Mass Function (IMF) in the solar neighborhood. Caveats on the conversion of the PDMF to the IMF are discussed. The PDMF slope changes, as a function of the radial distance from the center of the NGC 346 star cluster, indicating a segregation of the most massive stars. This segregation is likely primordial considering the young age (~3 Myr) of NGC346, and its clumpy structure which suggests that the cluster has likely not had sufficient time to relax. Comparing our results for NGC346 with those derived for other star clusters in the SMC and the Milky Way (MW), we conclude that, while the star formation process might depend on the local cloud conditions, the IMF does not seem to be affected by general environmental effects such as galaxy type, metallicity, and dust content.

NGC 602 is a young stellar cluster located in a peripheral region of the Small Magellanic Cloud known as the wing. Far from the main body of the galaxy and abutting the Magellanic Bridge, the SMC's wing is characterized by low gas and stellar content. With deep optical imaging from the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope (HST), we have discovered an extensive pre-Main Sequence (PMS) population, with stellar masses in the range 0.6-3 Solar masses. These low mass PMS stars formed coevally with the central cluster about 4 Myr ago. Spitzer Space Telescope (Spitzer) images of the same region from the Infrared Array Camera (IRAC) also reveal a population of Young Stellar Objects (YSOs), some of which are still embedded in nebular material and most of which likely formed even more recently than the young stars detected with HST/ACS imaging. We infer that star formation started in this region about 4 Myr ago with the formation of the central cluster and gradually propagated towards the outskirts where star formation is presently ongoing.

In the quest of understanding how star formation occurs and propagates in the low metallicity environment of the Small Magellanic Cloud (SMC), we acquired deep F555W (~V), and F814W (~I) HST/ACS images of the young and massive star forming region NGC 346. These images and their photometric analysis provide us with a snapshot of the star formation history of the region. We find evidence for star formation extending from ~10 Gyr in the past until ~150 Myr in the field of the SMC. The youngest stellar population (~3 +/- 1 Myr) is associated with the NGC 346 cluster. It includes a rich component of low mass pre-main sequence stars mainly concentrated in a number of sub-clusters, spatially co- located with CO clumps previously detected by Rubio et al. (2000). Within our analysis uncertainties, these sub-clusters appear coeval with each other. The most massive stars appear concentrated in the central sub-clusters, indicating possible mass segregation. A number of embedded clusters are also observed. This finding, combined with the overall wealth of dust and gas, could imply that star formation is still active. An intermediate age star cluster, BS90, formed ~4.3 +/-0.1 Gyr ago, is also present in the region. Thus, this region of the SMC has supported star formation with varying levels of intensity over much of the cosmic time.

Color-magnitude diagrams (CMDs) are presented for the Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) survey of the Large Magellanic Cloud (LMC). IRAC and MIPS 24 um epoch one data are presented. These data represent the deepest, widest mid-infrared CMDs of their kind ever produced in the LMC. Combined with the 2MASS survey, the diagrams are used to delineate the evolved stellar populations in the Large Magellanic Cloud as well as Galactic foreground and extragalactic background populations. Some 32000 evolved stars brighter than the tip of the red giant branch are identified. Of these, approximately 17500 are classified as oxygen-rich, 7000 carbon-rich, and another 1200 as ``extreme'' asymptotic giant branch (AGB) stars. Brighter members of the latter group have been called ``obscured'' AGB stars in the literature owing to their dusty circumstellar envelopes. A large number (1200) of luminous oxygen--rich AGB stars/M supergiants are also identified. Finally, there is strong evidence from the 24 um MIPS channel that previously unexplored, lower luminosity oxygen-rich AGB stars contribute significantly to the mass loss budget of the LMC (1200 such sources are identified).

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

We report the discovery of a rich population of low mass stars in the young, massive star forming region N66/NGC346 in the Small Magellanic Cloud, from deep V, I and H alpha images taken with the HST/ACS. These stars have likely formed together with the NGC346 cluster, ~3-5 Myr ago. Their magnitude and colors are those of pre-main sequence stars in the mass range 0.6-3 Mo, mostly concentrated in the main cluster, but with secondary subclusters spread over a region across ~45 pc. These subclusters appear to be spatially coincident with previously known knots of molecular gas identified in ground based and ISO observations. We show that N66/NGC346 is a complex region, being shaped by its massive stars, and the observations presented here represent a key step towards the understanding of how star formation occurred and has progressed in this low metallicity environment.

We present the results of NTT/VLT UBV imaging of a 260 square arcmin region containing the Galactic Luminous Blue Variable WRA751, in search for its birth-cluster, i.e. a cluster of young and massive stars spatially and physically associated with it. On the basis of the classical reddening-free parameter Q, we have identified a sample of 24 early-type stars with colours typical of spectral types earlier than B3. Interestingly, these stars are clustered within a radius of 1 arcmin from WRA751, corresponding to about 1% of the imaged field. These stars tightly distribute around (B-V) = 1.67, which in turn defines a mean extinction A(V) = 6.1 mag. The 5 brighter (V > 16.2) and bluer (Q < -0.9) stars of the sample have been subsequently observed with FORS1 and classified as 3 late O- and 2 early B- stars. The absence of stars earlier than O8 indicates an age of the cluster older than 4 Myr, although it could be due to an incomplete sampling of the upper end of the main sequence. Nevertheless, the detection of OB stars of class I certainly indicates an age of a few million years. At an assumed distance of 6 kpc, we estimate a cluster radius of 3.4 pc and a total mass of 2200 solar masses. Our discovery is only the second known instance of a Galactic Luminous Blue Variable associated with its birth-cluster.

We present the first detection of 12 CO J=2->1 and 12 CO J=1->0 emission from the LBV AG Carinae. AG Carinae resides in a region which is very rich in molecular gas with complex motions. We find evidence of a slow outflow of molecular gas, expanding at ~ 7 km/s. This emission appears spatially unresolved. We argue that it is spatially localised, rather than extended, and possibly associated with the immediate circumstellar region of AG Carinae. Does it originate from a circumstellar envelope, similar to carbon stars, or from a circumstellar disk? The option of the circumstellar disk is preferable because it is consistent with additional independent indications for the existence of wind asymmetries in close proximity to the central star, found from spectropolarimetry and analysis of the UV and optical line profiles, and it provides the conditions of density and shielding necessary for the survival of the CO molecules in proximity to such a hot star (Teff ~ 14000 K - 20000 K). In the assumption that the CO emission originated when AG Carinae was in an evolved state, we derive a lower limit to the mass of molecular gas of 2.8 solar masses. This is smaller, but still comparable with the mass of ionized gas present in the circumstellar environment (4.2 solar masses), with the implication that the molecular gas fraction can contribute significantly to the overall mass lost from the central star in its post main sequence evolution.

We present high resolution HST imaging of the nebula associated with the galactic LBV candidate HD 168625, together with ISO imaging and AAT echelle spectroscopy. The overall nebular morphology is elliptical with the major axis at PA ~ 120 deg. The dimensions of the nebula are 12'' x 16''.7 at Halpha and 15''.5 x 23''.5 at 4 micron. In the HST Halpha image, the nebula is resolved into a complex structure of filaments and arcs of different brightness. The asymmetry is lost in the HST continuum image where the nebula appears more diffuse and richer in filaments and clumps with the shape of cometary tails. At 11.3 micron the nebular emission peaks in two diametrically opposite lobes, placed on the nebula boundaries and along its major axis. A very faint loop is also visible at optical wavelengths, north and south of the shell. We suggest that the nebula is an ellipsoid with projected sizes of 14'' and 9'' (0.19 pc x 0.12 pc) along the RA and DEC directions, respectively. This ellipsoid is expanding at 19 km/s and is dynamically as old as ~ 4800 yrs; it probably interacts with the stellar wind and the loop so that PAH emission is detected from its caps, i.e. the lobes seen in the ISO images. The chemistry of the loop suggests that it is composed of un-processed material, probably from the local interstellar medium swept by the stellar wind.

We use the analysis of the heavy element abundances (C, N, O, S) in circumstellar nebulae around Luminous Blue Variables to infer the evolutionary phase in which the material has been ejected. (1) We discuss the different effects that may have changed the gas composition of the nebula since it was ejected (2) We calculate the expected abundance changes at the stellar surface due to envelope convection in the red supergiant phase. If the observed LBV nebulae are ejected during the RSG phase, the abundances of the LBV nebulae require a significantly smaller amount of mass to be lost than assumed in evolutionary models. (3) We calculate the changes in the surface composition during the main sequence phase by rotation induced mixing. If the nebulae are ejected at the end of the MS-phase, the abundances in LBV nebulae are compatible with mixing times between 5 x 10^6 and 1 x 10^7 years. The existence of ON stars supports this scenario. (4) The predicted He/H ratio in the nebulae are significantly smaller than the current observed photospheric values of their central stars. Combining various arguments we show that the LBV nebulae are ejected during the blue SG phase and that the stars have not gone through a RSG phase. The chemical enhancements are due to rotation induced mixing, and the ejection is possibly triggered by near-critical rotation. During the ejection, the outflow was optically thick, which resulted in a large effective radius and a low effective temperature. This also explains the observed properties of LBV dust.

We present HST NICMOS photometry of the resolved stellar population in NGC 1569. The CMD in the F110W and F160W photometric bands contains ~2400 stars with a formal photometric error < 0.1 mag down to J~23.5 and H~22.5. We describe the data processing which is required to calibrate the instrumental peculiarities of NICMOS. Two different packages for PSF-fitting photometry are used to strengthen the photometric results in the crowded stellar field of NGC 1569. The resulting CMD is discussed in terms of the major evolutionary properties of the resolved stellar populations. For a distance modulus of 26.71 and a reddening E(B-V)=0.56, our CMD samples stars down to ~0.8 Mo, corresponding to look-back times > 15 Gyr. This is clear indication of SF activity spanning an entire Hubble time. The metallicity of the reddest RGB stars is in better agreement with Z=0.004 as measured in HII regions, than with Z=0.0004 as expected from the stellar ages. The presence of - yet undetected - very metal-poor stars embedded in the stellar distribution around J=22.75 and J-H=1.15 is, however, not ruled out. The youngest stars (< 50 Myr) are preferentially found around the two central super star clusters, whereas the oldest population has a more uniform spatial distribution. A SFR per unit area of 1 Mo yr*(-1) kpc*(-2) and a mass formed in stars of ~ 1.4x10*6 Mo in the last 50 Myr are derived from the CMD. The NIR CMD places strong constraints on the lower limit of the onset of SF in NGC 1569. The exceptionally high crowding in the NICMOS images of NGC 1569 is a challenge for the photometric analysis. As a result, optical and NIR images of NGC 1569 sample different populations and cannot be cross-correlated. Nevertheless, we demonstrate the consistency of the SF histories derived from the optical and NIR CMDs.

We present new, NTT longslit spectroscopy of the B[e] supergiant in the binary system R4 in the Small Magellanic Cloud. The data show extended, forbidden N and S emissions which are typical signatures of circumstellar matter. Their extension along the space axis of the slit defines an angular size of 8.6'' which translates into a linear size of 2.4 pc. The N emission lines also show the velocity structure of a bipolar outflow expanding at 100 km/s on average. This implies that, for a measured radius of 1.2 pc, the outflow originated about 12000 yrs ago. The line flux ratio [NII]6584/[SII]6717 indicates that the nebula is nitrogen enriched and therefore it has been ejected from the central star. This is the first bipolar, ejection nebula detected around a well-established B[e] supergiant. The bipolar morphology and the chemical enrichment shown by the nebula associated with R4 are consistent with the picture of a binary merger (Langer & Heger 1998), in which R4 was originally a system composed by a close pair and a third star (the observed A companion). The close pair merged into a single star and the merging process produced a circumstellar nebula that was later shaped by the ensueing B star wind.

We report the result of a study in which we have used very deep broadband V and I WFPC2 images of the R136 cluster in the Large Magellanic Cloud from the HST archive, to sample the luminosity function below the detection limit of 2.8 Mo previously reached. In these new deeper images, we detect stars down to a limiting magnitude of m_F555W = 24.7 (~ 1 magnitude deeper than previous works), and identify a population of red stars evenly distributed in the surrounding of the R136 cluster. A comparison of our color-magnitude diagram with recentely computed evolutionary tracks indicates that these red objects are pre-main sequence stars in the mass range 0.6 - 3 Mo. We construct the initial mass function (IMF) in the 1.35 - 6.5 Mo range and find that, after correcting for incompleteness, the IMF shows a definite flattening below ~ 2 Mo. We discuss the implications of this result for the R136 cluster and for our understanding of starburst galaxies formation and evolution in general.

We present WFPC2-HST photometry (in the F439W and F555W bands) of the resolved stars in NGC 1569. The derived color-magnitude diagram (CMD) contains 2800 objects down to the 26th mag, and is complete down to the 23rd. Adopting the literature distance modulus and reddening, the CMD samples stars more massive than 4 Mo, corresponding to a look-back time of 0.15 Gyr approximately. The data are compared to synthetic CMDs to derive the recent star formation history in NGC 1569. It is found that the observed field has experienced a global burst of star formation (SF) lasting at least 0.1 Gyr, and ended a few Myr ago. During the burst, the SF rate was approximately constant, and, if quiescent periods occurred, they lasted less than 10 Myr. The level of the SF rate was very high: 0.5 Mo/yr for a Salpeter IMF, in stars with mass ranging from 0.1 to 120 Mo. When scaled for the surveyed area, the SF rate is approximately 100 times larger than found in the most active dwarf irregulars in the Local Group. The data are consistent with a Salpeter IMF, or slightly steeper exponents. We discuss the implications of our results in the general context of the evolution of dwarf galaxies.