Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in diverse galactic environmental conditions accessible in the local Universe. At high redshift, where conditions for cluster formation and evolution are more extreme, we are only able to observe the integrated light of a handful of objects that we believe will become globular clusters. We therefore discuss how numerical and analytical methods, informed by the observed properties of cluster populations in the local Universe, are used to develop sophisticated simulations potentially capable of disentangling the genetic map of galaxy formation and assembly that is carried by globular cluster populations.

We present recent results obtained with the VLT/MUSE Integral Field Spectrograph fed by the 4LGSF and its laser tomography adaptive optics module GALACSI. While this so-called narrow-field mode of MUSE was not designed to perform directly imaging of exoplanets and outflows, we show that it can be a game changer to detect and characterize young exoplanets with a prominent emission lines (i.e H{\alpha}, tracer of accretion), at moderate contrasts. These performances are achieved thanks to the combo of a near-diffraction limited PSF and a medium resolution spectrograph and a cross-correlation approach in post-processing . We discuss this in the context of ground and space, infrared and visible wavelengths, preparing for missions like JWST and WFIRST in great synergy and as pathfinder for future ELT/GSMT (Extremely Large and/or Giant Segmented Mirror Telescopes) instruments.

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.

We present MUSE Integral Field Unit (IFU) observations of five individual HII regions in two giant (>100 pc in radius) star-forming complexes in the low-metallicity ($Z$~0.33 $Z_{\odot}$) nearby (D ~ 2 Mpc) dwarf spiral galaxy NGC 300. We combine the IFU data with high spatial resolution HST photometry to demonstrate the extraction of stellar spectra and the classification of individual stars from ground-based data at the distance of 2 Mpc. For the two star-forming complexes, in which no O-type stars had previously been identified, we find a total of 13 newly identified O-type stars in the mass range 15-50 M$_{\odot}$, as well as 4 Wolf-Rayet stars. We use the derived massive stellar content to analyze the impact of stellar feedback on the HII regions. As already found for HII regions in the Magellanic Clouds, the dynamics of the analyzed NGC 300 HII regions are dominated by a combination of the pressure of the ionized gas and stellar winds. By comparing the derived ionized gas mass loading factors to the total gas mass loading factor across the NGC 300 disk, we find that the latter is an order of magnitude higher, either indicating very early evolutionary stages for these HII regions, or being a direct result of the multi-phase nature of feedback-driven bubbles. Moreover, we analyze the relation between the star formation rate and the pressure of the ionized gas as derived from small (<100 pc) scales, as both quantities are systematically overestimated when derived on galactic scales. With the wealth of ongoing and upcoming IFU instruments and programs, this study serves as a pathfinder for the systematic investigation of resolved stellar feedback in nearby galaxies, and it delivers the necessary analysis tools to enable massive stellar content and feedback studies sampling an unprecedented range of HII region properties across entire galaxies in the nearby Universe.

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.

Diagnosing the physical and chemical conditions within star-forming galaxies (SFGs) is of paramount importance to understanding key components of galaxy formation and evolution: star-formation, gas enrichment, outflows, and accretion. Well established optical emission-line diagnostics used to discern such properties (i.e., metal content, density, strength/shape of ionizing radiation) will be observationally inaccessible for the earliest galaxies, emphasizing the need for robust, reliable interstellar medium (ISM) diagnostics at ultraviolet (UV) wavelengths. Calibrating these UV diagnostics requires a comprehensive comparison of the UV and optical emission lines in nearby SFGs. Optical integral field unit (IFU) surveys have revealed the inhomogeneous nature of the ISM in SFGs, which leads to non-systematic biases in the interpretation of unresolved sources. Spatial variations are especially important to consider at UV wavelengths, where the strongest emission features originate from only the highest excitation regions of the nebula and are challenging to distinguish from competing high-ionization sources (e.g., shocks, AGN, etc.). Since surveys collecting large-scale optical integral field unit (IFU) spectroscopy are already underway, this white paper calls for an IFU or multi-object far-UV (FUV) spectroscopic instrument with high sensitivity, high spatial resolution, and large field of view (FoV). Given the impact of large-scale optical IFU surveys over the past decade, this white paper emphasizes the scientific need for a comparable foundation of spatially-resolved far-UV spectroscopy survey of nearby galaxies that will lay the foundation of diagnostics critical to the interpretation of the distant universe.

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 present a new Galaxy population synthesis Model (GalMod). GalMod is a star-count model featuring an asymmetric bar/bulge as well as spiral arms and related extinction. The model, initially introduced in Pasetto et al. (2016b), has been here completed with a central bar, a new bulge description, new disk vertical profiles and several new bolometric corrections. The model can generate synthetic mock catalogs of visible portions of the Milky Way (MW), external galaxies like M31, or N-body simulation initial conditions. At any given time, e.g., a chosen age of the Galaxy, the model contains a sum of discrete stellar populations, namely bulge/bar, disk, halo. These populations are in turn the sum of different components: the disk is the sum of spiral arms, thin disks, a thick disk, and various gas components, while the halo is the sum of a stellar component, a hot coronal gas, and a dark matter component. The Galactic potential is computed from these population density profiles and used to generate detailed kinematics by considering up to the first four moments of the collisionless Boltzmann equation. The same density profiles are then used to define the observed color-magnitude diagrams in a user-defined field of view from an arbitrary solar location. Several photometric systems have been included and made available on-line and no limits on the size of the field of view are imposed thus allowing full sky simulations, too. Finally, we model the extinction adopting a dust model with advanced ray-tracing solutions. The model's web page (and tutorial) can be accessed at www.GalMod.org.

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.

The JAXA Hayabusa-2 mission was approved in 2010 and launched on December 3, 2014. The spacecraft will arrive at the near-Earth asteroid 162173 Ryugu in 2018 where it will perform a survey, land and obtain surface material, then depart in Dec 2019 and return to Earth in Dec 2020. We observed Ryugu with the Herschel Space Observatory in Apr 2012 at far-IR thermal wavelengths, supported by several ground-based observations to obtain optical lightcurves. We reanalysed previously published Subaru-COMICS and AKARI-IRC observations and merged them with a Spitzer-IRS data set. In addition, we used a large set of Spitzer-IRAC observations obtained in the period Jan to May, 2013. The data set includes two complete rotational lightcurves and a series of ten "point-and-shoot" observations. The almost spherical shape of the target together with the insufficient lightcurve quality forced us to combine radiometric and lightcurve inversion techniques in different ways to find the object's key physical and thermal parameters. We find that the solution which best matches our data sets leads to this C class asteroid having a retrograde rotation with a spin-axis orientation of (lambda = 310-340 deg; beta = -40+/-15 deg) in ecliptic coordinates, an effective diameter (of an equal-volume sphere) of 850 to 880 m, a geometric albedo of 0.044 to 0.050 and a thermal inertia in the range 150 to 300 Jm-2s-0.5K-1. Based on estimated thermal conductivities of the top-layer surface in the range 0.1 to 0.6 WK-1m-1, we calculated that the grain sizes are approximately equal to between 1 and 10 mm. The finely constrained values for this asteroid serve as a `design reference model', which is currently used for various planning, operational and modelling purposes by the Hayabusa2 team.

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 report on the study of interstellar extinction across the Tarantula nebula (30 Doradus), in the Large Magellanic Cloud, using observations from the Hubble Tarantula Treasury Project in the 0.3 - 1.6 micron range. The considerable and patchy extinction inside the nebula causes about 3500 red clump stars to be scattered along the reddening vector in the colour-magnitude diagrams, thereby allowing an accurate determination of the reddening slope in all bands. The measured slope of the reddening vector is remarkably steeper in all bands than in the the Galactic diffuse interstellar medium. At optical wavelengths, the larger ratio of total-to-selective extinction, namely Rv = 4.5 +/- 0.2, implies the presence of a grey component in the extinction law, due to a larger fraction of large grains. The extra large grains are most likely ices from supernova ejecta and will significantly alter the extinction properties of the region until they sublimate in 50 - 100 Myr. We discuss the implications of this extinction law for the Tarantula nebula and in general for regions of massive star formation in galaxies. Our results suggest that fluxes of strongly star forming regions are likely to be underestimated by a factor of about 2 in the optical.

We present and describe the astro-photometric catalog of more than 800,000 sources found in the Hubble Tarantula Treasury Project (HTTP). HTTP is a Hubble Space Telescope (HST) Treasury program designed to image the entire 30 Doradus region down to the sub-solar (~0.5 solar masses) mass regime using the Wide Field Camera 3 (WFC3) and the Advanced Camera for Surveys (ACS). We observed 30 Doradus in the near ultraviolet (F275W, F336W), optical (F555W, F658N, F775W), and near infrared (F110W, F160W) wavelengths. The stellar photometry was measured using point-spread function (PSF) fitting across all the bands simultaneously. The relative astrometric accuracy of the catalog is 0.4 mas. The astro-photometric catalog, results from artificial star experiments and the mosaics for all the filters are available for download. Color-magnitude diagrams are presented showing the spatial distributions and ages of stars within 30 Dor as well as in the surrounding fields. HTTP provides the first rich and statistically significant sample of intermediate and low mass pre-main sequence candidates and allows us to trace how star formation has been developing through the region. The depth and high spatial resolution of our analysis highlight the dual role of stellar feedback in quenching and triggering star formation on the giant HII region scale. Our results are consistent with stellar sub-clustering in a partially filled gaseous nebula that is offset towards our side of the Large Magellanic Cloud.

We performed a deep wide-field (6.76 deg^2) near-infrared survey with the VISTA telescope that covers the entire extent of the Carina nebula complex (CNC). The point-source catalog created from these data contains around four million individual objects down to masses of 0.1 M_sun. We present a statistical study of the large-scale spatial distribution and an investigation of the clustering properties of infrared-excesses objects, which are used to trace disk-bearing young stellar objects (YSOs). We find that a (J - H) versus (Ks - [4.5]) color-color diagram is well suited to tracing the population of YSO-candidates (cYSOs) by their infrared excess. We identify 8781 sources with strong infrared excess, which we consider as cYSOs. This sample is used to investigate the spatial distribution of the cYSOs with a nearest-neighbor analysis. The surface density distribution of cYSOs agrees well with the shape of the clouds as seen in our Herschel far-infrared survey. The strong decline in the surface density of excess sources outside the area of the clouds supports the hypothesis that our excess-selected sample consists predominantly of cYSOs with a low level of background contamination. This analysis allows us to identify 14 groups of cYSOs outside the central area. Our results suggest that the total population of cYSOs in the CNC comprises about 164000 objects, with a substantial fraction (~35%) located in the northern, still not well studied parts. Our cluster analysis suggests that roughly half of the cYSOs constitute a non-clustered, dispersed population.

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 study of the recent star formation of 30 Doradus in the Large Magellanic Cloud (LMC) using the panchromatic imaging survey Hubble Tarantula Treasury Project (HTTP). In this paper we focus on the stars within 20 pc of the center of the massive ionizing cluster of 30 Doradus, NGC 2070. We recovered the star formation history by comparing deep optical and NIR color-magnitude diagrams (CMDs) with state-of-the-art synthetic CMDs generated with the latest PARSEC models, which include all stellar phases from pre-main sequence to post- main sequence. For the first time in this region we are able to measure the star formation using intermediate and low mass stars simultaneously. Our results suggest that NGC2070 experienced a prolonged activity. In particular, we find that the star formation in the region: i) exceeded the average LMC rate ~ 20 Myr ago; ii) accelerated dramatically ~ 7 Myr ago; and iii) reached a peak value 1-3 Myr ago. We did not find significant deviations from a Kroupa initial mass function down to 0.5 Msun. The average internal reddening E(B-V) is found to be between 0.3 and 0.4 mag.

In September 2011, the Herschel Space Observatory performed an observation campaign with the PACS photometer observing the asteroid (101955) 1999 RQ36 in the far infrared. The Herschel observations were analysed, together with ESO VLT-VISIR and Spitzer-IRS data, by means of a thermophysical model in order to derive the physical properties of 1999 RQ36. We find the asteroid has an effective diameter in the range 480 to 511 m, a slightly elongated shape with a semi-major axis ratio of a/b=1.04, a geometric albedo of 0.045 +0.015/-0.012, and a retrograde rotation with a spin vector between -70 and -90 deg ecliptic latitude. The thermal emission at wavelengths below 12 micron -originating in the hot sub-solar region- shows that there may be large variations in roughness on the surface along the equatorial zone of 1999 RQ36, but further measurements are required for final proof. We determine that the asteroid has a disk-averaged thermal inertia of Gamma = 650 Jm-2s-0.5K-1 with a 3-sigma confidence range of 350 to 950 Jm-2s-0.5K-1, equivalent to what is observed for 25143 Itokawa and suggestive that 1999 RQ36 has a similar surface texture and may also be a rubble-pile in nature. The low albedo indicates that 1999 RQ36 very likely contains primitive volatile-rich material, consistent with its spectral type, and that it is an ideal target for the OSIRIS-REx sample return mission.