We study the global dust and (molecular) gas content in the Lyman Alpha Reference Sample (LARS), i.e. 14 local star-forming galaxies. We characterize their interstellar medium and relate newly derived properties to quantities relevant for Ly$\alpha$ escape. We observed LARS galaxies with Herschel/PACS, SOFIA/FIFI-LS, the IRAM 30m telescope and APEX, targeting far-infrared (FIR) continuum and emission lines of [C II]158$\mu$m, [O I]63$\mu$m, [O III]88$\mu$m and low-J CO lines. Using Bayesian methods we derive dust model parameters and estimate total gas masses for all LARS galaxies, taking into account a metallicity-dependent gas-to-dust ratio. Star formation rates were estimated from FIR, [C II]158$\mu$m and [O I]63$\mu$m luminosities. LARS covers a wide dynamic range in the derived properties, with FIR-based star formation rates from $\sim$0.5-100 $M_{\odot}\ yr^{-1}$, gas fractions between $\sim$15-80% and gas depletion times ranging from a few hundred Myr up to more than 10 Gyr. The distribution of LARS galaxies in the $\Sigma_{gas}$ vs. $\Sigma_{SFR}$ (Kennicutt-Schmidt plane) is thus quite heterogeneous. However, we find that LARS galaxies with the longest gas depletion times, i.e. relatively high gas surface densities ($\Sigma_{gas}$) and low star formation rate densities ($\Sigma_{SFR}$), have by far the highest Ly$\alpha$ escape fraction. A strong $\sim$linear relation is found between Ly$\alpha$ escape fraction and the total gas (HI+H$_2$) depletion time. We argue that the Ly$\alpha$ escape in those galaxies is driven by turbulence in the star-forming gas that shifts the Ly$\alpha$ photons out of resonance close to the places where they originate. We further report on an extreme [C II]158$\mu$m excess in LARS 5, corresponding to $\sim$14$\pm$3% of the FIR luminosity, i.e. the most extreme [C II]-to-FIR ratio observed in a non-AGN galaxy to date.

Understanding the production and escape of Lyman $\alpha$ (Ly$\alpha$) radiation from star-forming galaxies is a long standing problem in astrophysics. The ability to predict the Ly$\alpha$ luminosity of galaxies would open up new ways of exploring the Epoch of Reionization (EoR), and to estimate Ly$\alpha$ emission from galaxies in cosmological simulations where radiative transfer calculations cannot be done. We apply multivariate regression methods to the Lyman Alpha Reference Sample dataset to obtain a relation between the galaxy properties and the emitted Ly$\alpha$. The derived relation predicts the Ly$\alpha$ luminosity of our galaxy sample to good accuracy, regardless of whether we consider only direct observables (root-mean-square (RMS) dispersion around the relation of $\sim 0.19$ dex) or derived physical quantities (RMS $\sim 0.27$ dex). We confirm the predictive ability on a separate sample of compact star-forming galaxies and find that the prediction works well, but that aperture effects on measured Ly$\alpha$ luminosity may be important, depending on the redshift of the galaxy. We apply statistical feature selection techniques to determine an order of importance of the variables in our dataset, enabling future observations to be optimized for predictive ability. When using physical variables, we are able to determine that the most important predictive parameters are, in order, star formation rate, dust extinction, compactness and the gas covering fraction. We discuss the application of our results in terms of studying the EoR and intensity mapping experiments.

Using MUSE AO data, we probe the ISM in the local spiral galaxy NGC 7793 at a spatial resolution of $\sim$ 10 pc. We identify HII regions and compile a catalogue of supernova remnants (SNRs), planetary nebulae (PNe) and Wolf Rayet stars (WR). We compute electron densities and temperatures from the [SII]6716/6731 and [SIII]6312/9069 line ratios. We study the properties of the ionised gas through BPT diagrams combined with gas velocity dispersion. We spectroscopically confirm 2 previously detected WR stars and a SNR and report the discovery of 7 WR, 1 SNR, and 2 PNe. The diffuse ionized gas (DIG) fraction is between $\sim$ 27 and 42% depending on the method used to define the HII region boundaries. In agreement with previous studies, we find that the DIG exhibits enhanced [SII]/H$\alpha$ and [NII]/H$\alpha$ ratios and a median temperature that is $\sim$ 3000 K higher than in HII regions. We also observe an apparent inverse correlation between temperature and H$\alpha$ surface brightness. Overall, the observed [SII]6716/6731 ratio is consistent within 1$\sigma$ with $n_e$ < 30 cm$^{-3}$, with an almost identical distribution for the DIG and HII regions. The velocity dispersion of the ionised gas indicates that the DIG has a higher degree of turbulence than the HII regions. Comparison with photoionisation and shock models reveals that the diffuse component can only partially be explained via shocks and that it is most likely consistent with photons leaking from density bounded HII regions or with radiation from evolved field stars. [abridged]

The James Webb Space Telescope will provide deep imaging and spectroscopy for sources at redshifts above 6, covering the Epoch of Reionization (EoR, 6 < z < 10). The Mid-IR instrument (MIRI) integral field spectrograph (MRS) will be the only instrument on board JWST able to observe the brightest optical emission lines H$\alpha$ and [OIII]0.5007$\mu$m at redshifts above 7 and 9, respectively. This paper presents a study of the H$\alpha$ fluxes predicted by FIRSTLIGHT cosmological simulations for galaxies at redshifts of 6.5 to 10.5, and its detectability with MIRI. Deep (40 ks) spectroscopic integrations with MRS will be able to detect (S/N > 5) EoR sources at redshifts above 7 with intrinsic star formation rates of more than 2 M$_{\odot}$ yr$^{-1}$, and stellar masses above 4-9 $\times$ 10$^7$ M$_{\odot}$. In addition, the paper presents realistic MRS simulated observations of the expected (rest-frame) optical and near-infrared spectra for some spectroscopically confirmed EoR sources detected by ALMA as [OIII]88$\mu$m emitters. The MRS simulated spectra cover a wide range of low metallicities from about 0.2 to 0.02Z$_{\odot}$, and different [OIII]88$\mu$m/[OIII]0.5007$\mu$m line ratios. The simulated 10ks MRS spectra show S/N in the range of 5 to 90 for H$\beta$, [OIII]0.4959,0.5007$\mu$m, H$\alpha$ and HeI1.083$\mu$m emission lines of MACS1149-JD1 at z = 9.11, independent of metallicity. In addition, deep 40 ks simulated spectra of the luminous merger candidate B14-65666 at z=7.15 shows the MRS capabilities of detecting, or putting strong upper limits, on the [NII]0.6584$\mu$m, [SII]0.6717,0.6731$\mu$m, and [SIII]0.9069,0.9532$\mu$m emission lines. In summary, MRS will enable the detailed study of key physical properties like internal extinction, instantaneous star formation, hardness of the ionising continuum, and metallicity, in bright (intrinsic or lensed) EoR sources.

We have used the capability of the MUSE instrument to explore the impact of stellar feedback at large scales in Haro 11, a galaxy under extreme starburst condition and one of the first galaxies where Lyman continuum (LyC) has been detected. Using Ha, [OIII] and [OI] emission lines from deep MUSE observations, we have constructed a sequence of velocity-dependent maps of the Ha emission, the state of the ionised gas and a tracer of fast shocks. These allowed us to investigate the ionisation structure of the galaxy in 50 kms^2 bins over a velocity range of -400 to 350 kms. The ionised gas in Haro 11 is assembled by a rich arrangement of structures, such as superbubbles, filaments, arcs and galactic ionised channels, whose appearances change drastically with velocity. The central star forming knots and the star forming dusty arm are the main engines that power the strong mechanical feedback in this galaxy, although with different impact on the ionisation structure. Haro 11 appears to leak LyC radiation in many directions. We found evidence of a kpc-scale fragmented superbubble, that may have cleared galactic-scale channels in the ISM. Additionally, the southwestern hemisphere is highly ionised in all velocities, hinting at a density bound scenario. A compact kpc-scale structure of lowly ionised gas coincides with the diffuse Lya emission and the presence of fast shocks. Finally, we find evidence that a significant fraction of the ionised gas mass may escape the gravitational potential of the galaxy.

We study young star-forming clumps on physical scales of 10-500 pc in the Lyman-Alpha Reference Sample (LARS), a collection of low-redshift (z = 0.03-0.2) UV-selected star-forming galaxies. In each of the 14 galaxies of the sample, we detect clumps for which we derive sizes and magnitudes in 5 UV-optical filters. The final sample includes $\sim$1400 clumps, of which $\sim$600 have magnitude uncertainties below 0.3 in all filters. The UV luminosity function for the total sample of clumps is described by a power-law with slope $\alpha = -2.03^{+0.11}_{-0.13}$. Clumps in the LARS galaxies have on average $\Sigma_{SFR}$ values higher than what observed in HII regions of local galaxies and comparable to typical SFR densities of clumps in z = 1-3 galaxies. We derive the clumpiness as the relative contribution from clumps to the UV emission of each galaxy, and study it as a function of galactic-scale properties, i.e. $\Sigma_{SFR}$ and the ratio between rotational and dispersion velocities of the gas ($v_s/\sigma_0$). We find that in galaxies with higher $\Sigma_{SFR}$ or lower $v_s/\sigma_0$, clumps dominate the UV emission of their host systems. All LARS galaxies with Ly$\alpha$ escape fractions larger than 10% have more than 50% of the UV luminosity from clumps. We tested the robustness of these results against the effect of different physical resolutions. At low resolution, the measured clumpiness appears more elevated than if we could resolve clumps down to single clusters. This effect is small in the redshift range covered by LARS, thus our results are not driven by the physical resolution.

We constrain the stellar population properties of a sample of 52 massive galaxies, with stellar mass log Ms>10.5, over the redshift range 0.5<z<2 by use of observer-frame optical and near-infrared slitless spectra from HST's ACS and WFC3 grisms. The deep exposures (~100 ks) allow us to target individual spectra of massive galaxies to F160W=22.5AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via an MCMC method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Dt) that is found to dominate the variance of the retrieved parameters according to Principal Component Analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Dt appears to weakly anti-correlate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

Mrk 71/NGC 2366 is the closest Green Pea (GP) analog and candidate Lyman Continuum (LyC) emitter. Recently, 11 LyC-leaking GPs have been detected through direct observations of the ionizing continuum, making this the most abundant class of confirmed LyC-emitters at any redshift. High resolution, multi-wavelength studies of GPs can lead to an understanding of the method(s), through which LyC escapes from these galaxies. The proximity of Mrk 71/NCG 2366 offers unprecedented detail on the inner workings of a GP analog, and enables us to identify the mechanisms of LyC escape. We use 5825-7650{\AA} integral field unit PMAS observations to study the kinematics and physical conditions in Mrk 71. An electron density map is obtained from the [S II] ratio. A fortuitous second order contamination by the [O II]3727 doublet enables the construction of an electron temperature map. Resolved maps of sound speed, thermal broadening, "true" velocity dispersion, and Mach number are obtained and compared to the high resolution magneto-hydrodynamic SILCC simulations. Two regions of increased velocity dispersion indicative of outflows are detected to the north and south of the super star cluster, knot B, with redshifted and blueshifted velocities, respectively. We confirm the presence of a faint broad kinematical component, which is seemingly decoupled from the outflow regions, and is fainter and narrower than previously reported in the literature. Within uncertainties, the low- and high-ionization gas move together. Outside of the core of Mrk 71, an increase in Mach numbers is detected, implying a decrease in gas density. Simulations suggest this drop in density can be as high as ~4 dex, down to almost optically thin levels, which would imply a non-zero LyC escape fraction along the outflows... [abridged]

Stellar feedback strongly affects the interstellar medium (ISM) of galaxies. Stellar feedback in the first galaxies likely plays a major role in enabling the escape of LyC photons, which contribute to the re-ionization of the Universe. Nearby starburst galaxies serve as local analogues allowing for a spatially resolved assessment of the feedback processes in these galaxies. We characterize the feedback effects from the star clusters in the local high-redshift analogue ESO 338-IG04 on the ISM and compare the results with the properties of the most massive clusters. We use high quality VLT/MUSE optical integral field data to derive the physical properties of the ISM such as ionization, density, shocks, and perform new fitting of the spectral energy distributions of the brightest clusters in ESO 338-IG04 from HST imaging. ESO 338-IG04 has a large ionized halo which we detect to a distance of 9 kpc. We identify 4 Wolf-Rayet (WR) clusters based on the blue and red WR bump. We follow previously identified ionization cones and find that the ionization of the halo increases with distance. Analysis of the galaxy kinematics shows two complex outflows driven by the numerous young clusters in the galaxy. We find a ring of shocked emission traced by an enhanced [OI]/H$\alpha$ ratio surrounding the starburst and at the end of the outflow. Finally we detect nitrogen enriched gas associated with the outflow, likely caused by the WR stars in the massive star clusters. Photo-ionization dominates the central starburst and sets the ionization structure of the entire halo, resulting in a density bounded halo, facilitating the escape of LyC photons. Outside the central starburst, shocks triggered by an expanding super bubble become important. The shocks at the end of the outflow suggest interaction between the hot outflowing material and the more quiescent halo gas.

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH.

It has recently been established that the properties of young star clusters (YSCs) can vary as a function of the galactic environment in which they are found. We use the cluster catalogue produced by the Legacy Extragalactic UV Survey (LEGUS) collaboration to investigate cluster properties in the spiral galaxy M51. We analyse the cluster population as a function of galactocentric distance and in arm and inter-arm regions. The cluster mass function exhibits a similar shape at all radial bins, described by a power law with a slope close to $-2$ and an exponential truncation around $10^5\ \rm{M}_{\odot}$ . While the mass functions of the YSCs in the spiral arm and inter-arm regions have similar truncation masses, the inter-arm region mass function has a significantly steeper slope than the one in the arm region; a trend that is also observed in the giant molecular cloud mass function and predicted by simulations. The age distribution of clusters is dependent on the region considered, and is consistent with rapid disruption only in dense regions, while little disruption is observed at large galactocentric distances and in the inter-arm region. The fraction of stars forming in clusters does not show radial variations, despite the drop in the $H_2$ surface density measured as function of galactocentric distance. We suggest that the higher disruption rate observed in the inner part of the galaxy is likely at the origin of the observed flat cluster formation efficiency radial profile.

The Lyman Alpha Reference Sample (LARS) of 14 star-forming galaxies offers a wealth of insight into the workings of these local analogs to high-redshift star-forming galaxies. The sample has been well-studied in terms of LyA and other emission line properties, such as HI mass, gas kinematics, and morphology. We analyze deep surface photometry of the LARS sample in UBIK broadband imaging obtained at the Nordic Optical Telescope and the Canada-France-Hawaii Telescope, and juxtaposition their derived properties with a sample of local high-redshift galaxy analogs, namely, with blue compact galaxies (BCGs). We construct radial surface brightness and color profiles with both elliptical and isophotal integration, as well as RGB images, deep contours, color maps, a burst fraction estimate, and a radial mass-to-light ratio profile for each LARS galaxy. Standard morphological parameters like asymmetry, clumpiness, the Gini and M20 coefficients are [...] analyzed, as well as isophotal asymmetry profiles for each galaxy. [...] We compare the LARS to the properties of the BCG sample and highlight the differences. Several diagnostics indicate that the LARS galaxies have highly disturbed morphologies even at the level of the faintest isophotes [...]. The ground-based photometry [...] reveals previously unexplored isophotes [...]. The burst fraction estimate suggests a spatially more extended burst region in LARS than in the BCGs. [...] The galaxies in the LARS sample appear to be in earlier stages of a merger event compared to the BCGs. Standard morphological diagnostics like asymmetry, clumpiness, Gini and M20 coefficients cannot separate the two samples, although an isophotal asymmetry profile successfully captures the average difference in morphology. These morphological diagnostics do not show any correlation with the equivalent width or the escape fraction of Lyman Alpha. [abridged]

We present the results of an unbiased search for Ly{\alpha} emission from continuum-selected 6 < z < 8 galaxies. Our dataset consists of 160 orbits of G102 slitless grism spectroscopy obtained with the Hubble Space Telescope (HST) Wide Field Camera 3 as part of the Faint Infrared Grism Survey (FIGS; PI: Malhotra), which obtains deep slitless spectra of all sources in four fields, and was designed to minimize contamination in observations of previously-identified high-redshift galaxy candidates. The FIGS data can potentially spectroscopically confirm the redshifts of galaxies, and as Ly{\alpha} emission is resonantly scattered by neutral gas, FIGS can also constrain the ionization state of the intergalactic medium (IGM) during the epoch of reionization. These data have sufficient depth to detect Ly{\alpha} emission in this epoch, as Tilvi et al. (2016) have published the FIGS detection of previously known (Finkelstein et al. 2013) Ly{\alpha} emission at z = 7.51. The FIGS data use five separate roll-angles of HST to mitigate the contamination by nearby galaxies. We created a method that accounts for and removes the contamination from surrounding galaxies, and also removes any dispersed continuum light from each individual spectrum (Pirzkal et al. 2017). We searched for significant (> 4{\sigma}) emission lines using two different automated detection methods, free of any visual inspection biases. Applying these methods on photometrically-selected high-redshift candidates between 6 < z < 8 we find two emission lines, one previously published by Tilvi et al. (2016), and a new line at 1.028{\mu}m. We identify this lines as Ly{\alpha} at z = 7.452 +/- 0.003. This newly spectroscopically confirmed galaxy has the highest Ly{\alpha} rest-frame equivalent width (EW Ly{\alpha}) yet published at z > 7 (140.3 +/- 19.0{\AA}).

We examine the dust geometry and Ly{\alpha} scattering in the galaxies of the Lyman Alpha Reference Sample (LARS), a set of 14 nearby (0.02 < $z$ < 0.2) Ly{\alpha} emitting and starbursting systems with Hubble Space Telescope Ly{\alpha}, H{\alpha}, and H{\beta} imaging. We find that the global dust properties determined by line ratios are consistent with other studies, with some of the LARS galaxies exhibiting clumpy dust media while others of them show significantly lower Ly{\alpha} emission compared to their Balmer decrement. With the LARS imaging, we present Ly{\alpha}/H{\alpha} and H{\alpha}/H{\beta} maps with spatial resolutions as low as $\sim$ 40 pc, and use these data to show that in most galaxies, the dust geometry is best modeled by three distinct regions: a central core where dust acts as a screen, an annulus where dust is distributed in clumps, and an outer envelope where Ly{\alpha} photons only scatter. We show that the dust that affects the escape of Ly{\alpha} is more restricted to the galaxies' central regions, while the larger Ly{\alpha} halos are generated by scattering at large radii. We present an empirical modeling technique to quantify how much Ly{\alpha} scatters in the halo, and find that this "characteristic" scattering distance correlates with the measured size of the Ly{\alpha} halo. We note that there exists a slight anti-correlation between the scattering distance of Ly{\alpha} and global dust properties.

We report on the discovery of ionised gas filaments in the circum-galactic halo of the extremely metal-poor compact starburst SBS 0335-052E in a 1.5h integration with the MUSE integral-field spectrograph. We detect these features in H${\alpha}$ and [OIII] emission down to surface-brightness levels of $5 \times 10^{-19}$erg s$^{-1}$cm$^{-2}$arcsec$^{-2}$. The filaments have projected diameters of 2.1 kpc and extend more than 9 kpc to the north and north-west from the main stellar body. We also detect extended nebular HeII $\lambda$4686 emission that brightens towards the north-west at the rim of a star-burst driven super-shell, suggestive of a locally enhanced UV radiation field due to shocks. We also present a velocity field of the ionised gas. The filaments appear to connect seamlessly in velocity space to the kinematical disturbances caused by the shell. Similar to high-$z$ star-forming galaxies, the ionised gas in this galaxy is dispersion dominated. We argue that the filaments were created via feedback from the starburst and that these ionised structures in the halo may act as escape channels for Lyman continuum radiation in this gas-rich system.

Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used to study the young star cluster (YSC) population of this interacting system. Our newly extracted source catalogue contains 2834 cluster candidates, morphologically classified to be compact and uniform in colour, for which ages, masses and extinction are derived. In this first work we study the main properties of the YSC population of the whole galaxy, considering a mass-limited sample. Both luminosity and mass functions follow a power law shape with slope -2, but at high luminosities and masses a dearth of sources is observed. The analysis of the mass function suggests that it is best fitted by a Schechter function with slope -2 and a truncation mass at $1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm this result and link the shape of the luminosity function to the presence of a truncation in the mass function. A mass limited age function analysis, between 10 and 200 Myr, suggests that the cluster population is undergoing only moderate disruption. We observe little variation in the shape of the mass function at masses above $1\times10^4\ M_\odot$, over this age range. The fraction of star formation happening in the form of bound clusters in M51 is $\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over the whole range from 1 to 200 Myr.

Low- and intermediate mass galaxies are widely discussed as cause of reionization at redshift $z\sim10-6$. However, observational proof of galaxies that are leaking ionizing radiation (Lyman continuum; LyC) is a currently ongoing challenge and the list of LyC emitting candidates is still short. Tololo 1247-232 is among those very few galaxies with recently reported leakage. We performed intermediate resolution ultraviolet (UV) spectroscopy with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope and confirm ionizing radiation emerging from Tololo 1247-232. Adopting an improved data reduction procedure, we find that LyC escapes from the central stellar clusters, with an escape fraction of 1.5$\pm$0.5% only, i.e. the lowest value reported for the galaxy so far. We further make use of FUV absorption lines of Si II and Si IV as a probe of the neutral and ionized interstellar medium. We find that most of the ISM gas is ionized, likely facilitating LyC escape from density bounded regions. Neutral gas covering as a function of line-of-sight velocity is derived using the apparent optical depth method. The ISM is found to be sufficiently clumpy, supporting the direct escape of LyC photons. We further report on broadband UV and optical continuum imaging as well as narrowband imaging of Ly$\alpha$, H$\alpha$ and H$\beta$. Using stellar population synthesis, a Ly$\alpha$ escape fraction of 8% was derived. We also performed VLA 21cm imaging. The hydrogen hyperfine transition was not detected, but a deep upper limit atomic gas mass of $\lesssim10^9 M_{\odot}$ could be derived. The upper limit gas fraction defined as $\frac{M_{H I}}{M_*}$ is only 20%. Evidence is found that the H I gas halo is relatively small compared to other Lyman Alpha emitters.

Star forming galaxies are believed to be a major source of Lyman Continuum (LyC) radiation responsible for reionizing the early Universe. Direct observations of escaping ionizing radiation have however been few and with low escape fractions. In the local Universe, only ~10 emitters have been observed, with typical escape fractions of a few percent. The mechanisms regulating this escape need to be strongly evolving with redshift in order to account for the Epoch of Reionization. Gas content and star formation feedback are among the main suspects, known to both regulate neutral gas coverage and evolve with cosmic time. In this paper, we reanalyze HST-COS spectrocopy of the first detected local LyC leaker, Haro 11. We examine the connection between LyC leakage and Lyman-$\alpha$ line shape, and feedback-influenced neutral ISM properties like kinematics and gas distribution. We discuss the two extremes of an optically thin, density bounded ISM and a riddled, optically thick, ionization bounded ISM, and how Haro 11 fits into their theoretical predictions. We find that the most likely ISM model is a clumpy neutral medium embedded in a highly ionized medium with a combined covering fraction of unity and a residual neutral gas column density in the ionized medium high enough to be optically thick to Lyman-$\alpha$, but low enough to be at least partly transparent to Lyman continuum and undetected in Si II. This suggests that SF feedback and galaxy-scale interaction events play a major role in opening passageways for ionizing radiation through the neutral medium.

While the Lyman-$\alpha$ ($\mathrm{Ly}\alpha$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$\alpha$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$\alpha$ equivalent widths, EW(Ly$\alpha$) $\gtrsim 150$ \AA{}, at $z>6$ are good candidates for Ly$\alpha$ follow-up observations, and can stand out in multiband imaging surveys because of their unusual colors. We have conducted a photometric search for such objects in the Cluster Lensing And Supernova survey with Hubble (CLASH), and report here the identification of three likely gravitationally-lensed images of a single LAE candidate at $z\sim6.3$, behind the galaxy cluster Abell 2261($z = 0.225$). In the process, we also measured with Keck/MOSFIRE the first spectroscopic redshift of a multiply-imaged galaxy behind Abell 2261, at $z = 3.337$. This allows us to calibrate the lensing model, which in turn is used to study the properties of the candidate LAE. Population III galaxy spectral energy distribution (SED) model fits to the CLASH broadband photometry of the possible LAE provide a slightly better fit than Population I/II models. The best fitted model suggests intrinsic EW(Ly$\alpha$) $\approx 160$ \AA{} after absorption in the interstellar and intergalactic medium. Future spectroscopic observations will examine this prediction as well as shed more light on the morphology of this object, which indicates it may be a merger of two smaller galaxies.

We present the first robust detection of HI 21 cm emission in the blue compact galaxy Haro 11 using the 100m Robert C. Byrd Green Bank Telescope (GBT). Haro 11 is a luminous blue compact galaxy with emission in both Lyman Alpha and the Lyman continuum. We detect (5.1 $\pm$ 0.7 $\times$10$^8$) M$_{\odot}$ of HI gas at an assumed distance of 88 Mpc, making this galaxy HI deficient compared to other local galaxies with similar optical properties. Given this small HI mass, Haro 11 has an elevated M$_{H2}$/M$_{HI}$ ratio and a very low gas fraction compared to most local galaxies, and contains twice as much mass in ionized hydrogen as in neutral hydrogen. The HI emission has a linewidth of 71 kms$^{-1}$ and is offset 60 kms$^{-1}$ redward of the optical line center. It is undergoing a starburst after a recent merger which has elevated the star formation rate, and will deplete the gas supply in $<$ 0.2 Gyr. Although this starburst has elevated the SFR compared to galaxies with similar HI masses and linewidths, Haro 11 matches a trend of lower gas fractions toward higher star formation rates and is below the general trend of increasing HI mass with increasing luminosity. Taken together, our results paint Haro 11 as a standard low-mass galaxy that is undergoing an unusually efficient star formation episode.

We report results from a new HST study of the OVI 1032,1038\AA\ doublet in emission around intensely star-forming galaxies. The programme aims to characterize the energy balance in starburst galaxies and gas cooling in the difficult-to-map coronal temperature regime of 2-5 x $10^5$K. We present the first resolved image of gas emission in the OVI line. Our target, SDSS J1156+5008, is very compact in the continuum but displays OVI emission to radii of 23 kpc. The surface brightness profile is well fit by an exponential with a scale of 7.5kpc. This is ten times the size of the photoionized gas, and we estimate that 1/6 the total OVI luminosity comes from resonantly scattered continuum radiation. Spectroscopy - which closely resembles a stacked sample of archival spectra - confirms the OVI emission, and determines the column density and outflow velocity from blueshifted absorption. The combination of measurements enables several new calculations with few assumptions. The OVI regions fill only ~$10^{-3}$ of the volume. By comparing the cooling time with the cloud sound-crossing time, the cooling distance with the size, and the pressure in the OVI and nebular gas, we conclude that the OVI-bearing gas cannot have been lifted to the scale height at this temperature, and must be cooling in situ through this coronal temperature regime. The coronal phase contains ~1% of the ionized mass, and its kinetic energy is currently ~1% of the budget set by supernova feedback. However a much larger amount of the gas must have cooled through this phase during the star formation episode. The outflow exceeds the escape velocity and the gas may become unbound, but it will recombine before it escapes and become visible to Lyman (and OI) spectroscopy. The mapping of this gas represents a crucial step in further constraining galaxy formation scenarios and guiding the development of future satellites.

We test the predictions of spectral synthesis models based on seven different massive-star prescriptions against Legacy ExtraGalactic UV Survey (LEGUS) observations of eight young massive clusters in two local galaxies, NGC 1566 and NGC 5253, chosen because predictions of all seven models are available at the published galactic metallicities. The high angular resolution, extensive cluster inventory and full near-ultraviolet to near-infrared photometric coverage make the LEGUS dataset excellent for this study. We account for both stellar and nebular emission in the models and try two different prescriptions for attenuation by dust. From Bayesian fits of model libraries to the observations, we find remarkably low dispersion in the median E(B-V) (~0.03 mag), stellar masses (~10^4 M_\odot) and ages (~1 Myr) derived for individual clusters using different models, although maximum discrepancies in these quantities can reach 0.09 mag and factors of 2.8 and 2.5, respectively. This is for ranges in median properties of 0.05-0.54 mag, 1.8-10x10^4 M_\odot and 1.6-40 Myr spanned by the clusters in our sample. In terms of best fit, the observations are slightly better reproduced by models with interacting binaries and least well reproduced by models with single rotating stars. Our study provides a first quantitative estimate of the accuracies and uncertainties of the most recent spectral synthesis models of young stellar populations, demonstrates the good progress of models in fitting high-quality observations, and highlights the needs for a larger cluster sample and more extensive tests of the model parameter space.

While numerical simulations suggest that the strength of the Lyman alpha (Lya) line of star-forming disk galaxies strongly depends on the inclination at which they are observed (i.e. from edge-on to face-on, we expect to see a change from an attenuated Lya line to a strong Lya emission line), recent observations with the Hubble space telescope (HST) have highlighted few low-redshift highly inclined (edge-on) disk galaxies that breaks this trend. We aim to understand how a strong Lya emission line is able to escape from one of those inclined disk galaxies, named Mrk1486 (z=0.0338). For that purpose we used a large set of HST imaging and spectroscopic data to investigate both the ISM structure and the dominant source of Lya radiation inside Mrk1486. Moreover, we used a 3D Monte Carlo Lya radiation transfer code to study the radiative transfer of Lya and UV continuum photons inside a 3D geometry of neutral hydrogen (HI) and dust that models the ISM structure at the galaxy center. The analysis of IFU Halpha spectroscopic data of Mrk1486 indicates the presence of two bipolar galactic winds of HI gas above and bellow the disk plane of Mrk1486. Furthermore, comparing different diagnostic diagrams (such as [OIII]5007/Hbeta versus [OI]6300/Halpha) to photo- and shock-ionization models, we find that the Lya production of Mrk1486 is dominated by photoionization inside the galaxy disk. From this perspective, our numerical simulations succeed in reproducing the strength and spectral shape of the observed Lya line of Mrk1486 by assuming a scenario in which the Lya photons are produced inside the disk, travel along the galactic winds and scatter on cool HI materials toward the observer. As bipolar galactic winds are ubiquitous in star-forming disk galaxies, this mechanism may explain the origin of strong Lya emission lines recently observed from highly inclined galaxies at high-redshift.

Despite strong interest in the starburst (hereafter SB) phenomenon, the concept remains ill-defined. We use a strict definition of SB to examine the statistical properties of local SB and post-starburst (hereafter PB) galaxies. We also seek relationships to active galaxies. Potential SB galaxies are selected from the SDSS DR7 and their stellar content is analysed. We apply an age dependent dust attenuation correction and derive star formation rates (SFR), ages and masses of the young and old populations. The photometric masses nicely agree with dynamical masses derived from the H-alpha emission line width. To select SB galaxies, we use the birthrate parameter b=SFR/<SFR>, requiring b>=3. The PB sample is selected from the citerion EW(Hdelta_abs)>=6 A. Only 1% of star-forming galaxies are found to be SB galaxies. They contribute 3-6% to the stellar production and are therefore unimportant for the local star formation activity. The median SB age is 70 Myr, roughly independent of mass. The b-parameter strongly depends on burst age. Values close to b=60 are found at ages ~10 Myr, while almost no SBs are found at ages >1 Gyr. The median baryonic burst mass fraction of sub-L* galaxies is 5%, decreasing slowly with mass. The median mass fraction of the recent burst in the PB sample is 5-10%. The age-mass distribution of the progenitors of the PBs is bimodal with a break at log(M)~10.6 above which the ages are doubled. The SB and PB luminosity functions (hereafter LFs) follow each other closely until M_r~-21, when AGNs begin to dominate. The PB LF continues to follow the AGN LF while SB loose significance. This suggests that the number of luminous SBs is underestimated by about one dex at high luminosities, due to large amounts of dust and/or AGN blending. It also indicates that the SB phase preceded the AGN phase. We also discuss the conditions for global gas outflow caused by stellar feedback.

We present integral field spectroscopic observations with the Potsdam Multi Aperture Spectrophotometer of all 14 galaxies in the $z\sim 0.1$ Lyman Alpha Reference Sample (LARS). We produce 2D line of sight velocity maps and velocity dispersion maps from the Balmer $\alpha$ (H$\alpha$) emission in our data cubes. These maps trace the spectral and spatial properties of the LARS galaxies' intrinsic Ly$\alpha$ radiation field. We show our kinematic maps spatially registered onto the Hubble Space Telescope H$\alpha$ and Lyman $\alpha$ (Ly$\alpha$) images. Only for individual galaxies a causal connection between spatially resolved H$\alpha$ kinematics and Ly$\alpha$ photometry can be conjectured. However, no general trend can be established for the whole sample. Furthermore, we compute non-parametric global kinematical statistics -- intrinsic velocity dispersion $\sigma_0$, shearing velocity $v_\mathrm{shear}$, and the $v_\mathrm{shear}/\sigma_0$ ratio -- from our kinematic maps. In general LARS galaxies are characterised by high intrinsic velocity dispersions (54\,km\,s$^{-1}$ median) and low shearing velocities (65\,km\,s$^{-1}$ median). $v_\mathrm{shear}/\sigma_0$ values range from 0.5 to 3.2 with an average of 1.5. Noteworthy, five galaxies of the sample are dispersion dominated systems with $v_\mathrm{shear}/\sigma_0 <1$ and are thus kinematically similar to turbulent star forming galaxies seen at high redshift. When linking our kinematical statistics to the global LARS Ly$\alpha$ properties, we find that dispersion dominated systems show higher Ly$\alpha$ equivalent widths and higher Ly$\alpha$ escape fractions than systems with $v_\mathrm{shear}/\sigma_0 > 1$. Our result indicates that turbulence in actively star-forming systems is causally connected to interstellar medium conditions that favour an escape of Ly$\alpha$ radiation.

The leakage of Lyman continuum photons from star forming galaxies is an elusive parameter. When observed, it provides a wealth of information on star formation in galaxies and the geometry of the interstellar medium, and puts constraints on the role of star forming galaxies in the reionization of the universe. H-alpha-selected galaxies at z~2 trace the highest star formation population at the peak of cosmic star formation history, providing a base for directly measuring Lyman continuum escape. Here we present this method, and highlight its benefits as well as caveats. We also use the method on 10 H-alpha emitters in the Chandra Deep Field South at z=2.2, also imaged with the Hubble Space Telescope in the ultraviolet. We find no individual Lyman continuum detections, and our stack puts a 5 sigma upper limit on the average absolute escape fraction of <24%, consistent with similar studies. With future planned observations, the sample sizes would rapidly increase and the method presented here would provide very robust constraints on the escape fraction.

The stellar content of young massive star clusters emit large amounts of Lyman continuum photons and inject momentum into the inter stellar medium (ISM) by the strong stellar winds of the most massive stars in the cluster. When the most massive stars explode as supernovae, large amounts of mechanical energy are injected in the ISM. A detailed study of the ISM around these massive cluster provides insights on the effect of cluster feedback. We present high quality integral field spectroscopy taken with VLT/MUSE of two starburst galaxies: ESO 338-IG04 and Haro 11. Both galaxies contain a significant number of super star clusters. The MUSE data provide us with an unprecedented view of the state and kinematics of the ionized gas in the galaxy allowing us to study the effect of stellar feedback on small and large spatial scales. We present our recent results on studying the ISM state of these two galaxies. The data of both galaxies show that the mechanical and ionization feedback of the super star clusters in the galaxy modify the state and kinematics of the ISM substancially by creating highly ionized bubbles around the cluster, making the central part of the galaxy highly ionized. This shows that the HII regions around the individual clusters are density bounded, allowing the ionizing photons to escape and ionize the ISM further out.

(abridged) Luminous blue compact galaxies are among the most active galaxies in the local universe in terms of their star formation rate per unit mass. They may be seen as the local analogs of higher redshift Lyman Break Galaxies. Studies of their kinematics is key to understanding what triggers their unusually active star formation In this work we investigate the kinematics of stars and ionised gas in Haro11, one of the most luminous blue compact galaxies in the local universe. Previous works have indicated that many such galaxies may be triggered by galaxy mergers. We have employed Fabry-Perot interferometry, long-slit spectroscopy and Integral Field Unit (IFU) spectroscopy to explore the kinematics of Haro11. We target the near infrared Calcium triplet to derive the stellar velocity field and velocity dispersion. Ionised gas is analysed through emission lines from hydrogen, [OIII] , and [SIII]. When spectral resolution and signal to noise allows we investigate the the line profile in detail and identify multiple velocity components when present. We find that to first order, the velocity field and velocity dispersions derived from stars and ionised gas agree. Hence the complexities reveal real dynamical disturbances providing further evidence for a merger in Haro11. Through decomposition of emission lines we find evidence for kinematically distinct components, for instance a tidal arm behind the galaxy. The ionised gas velocity field can be traced to large galactocentric radii, and shows significant velocity dispersion even far out in the halo. We discuss the origin of the line width, and interpreted as virial motions it indicates a mass of ~1E11 M_sun. Morphologically and kinematically Haro11 shows many resemblances with the famous Antennae galaxies, but is much denser which is the likely explanation for the higher star formation efficiency in Haro11.

We present ultraviolet (UV) observations of six nearby Type~Ia supernovae (SNe~Ia) obtained with the Hubble Space Telescope, three of which were also observed in the near-IR (NIR) with Wide-Field Camera~3. UV observations with the Swift satellite, as well as ground-based optical and near-infrared data provide complementary information. The combined data-set covers the wavelength range $0.2$--$2~\mu$m. By also including archival data of SN 2014J, we analyse a sample spanning observed colour excesses up to $E(B-V)=1.4~$mag. We study the wavelength dependent extinction of each individual SN and find a diversity of reddening laws when characterised by the total-to-selective extinction $R_V$. In particular, we note that for the two SNe with $E(B-V)\gtrsim1~$mag, for which the colour excess is dominated by dust extinction, we find $R_V=1.4\pm0.1$ and $R_V=2.8\pm0.1$. Adding UV photometry reduces the uncertainty of fitted $R_V$ by $\sim50\,$% allowing us to also measure $R_V$ of individual low-extinction objects which point to a similar diversity, currently not accounted for in the analyses when SNe~Ia are used for studying the expansion history of the universe.

Context. Lyman alpha emitting galaxies (LAEs) are used to probe the distant universe and are therefore important for galaxy evolution studies and for providing clues to the nature of the epoch of reionization, but the exact circumstances under which Lyman alpha escapes a galaxy are still not fully understood. Aims. The Trident project is designed to simultaneously examine Lyman alpha, H-alpha and Lyman Continuum emission from galaxies at redshift z~2, thus linking together these three aspects of ionising radiation in galaxies. In this paper, we outline the strategy of this project and examine the properties of LAEs in the GOODS North field. Methods. We performed a narrowband LAE survey in GOODS North using existing and two custom made filters at the Nordic Optical Telescope with MOSCA. We use complementary broad band archival data in the field to make a careful candidate selection and perform optical to near-IR SED fitting. We also estimate far-infrared luminosities by matching our candidates to detections in Spitzer/MIPS 24{\mu}m and Herschel/PACS catalogs. Results. We find a total of 25 LAE candidates, probing mainly the bright end of the LAE luminosity function with L_Ly {\alpha} ~ 1-15e42 erg/s. They display a range of masses of ~0.5-50e9 M_solar, and average ages from a few tens of Myr to 1 Gyr when assuming a constant star formation history. The majority of our candidates also show signs of recent elevated star formation. Three candidates have counterparts in the GOODS-Herschel far-IR catalogue, with luminosities consistent with ultra-luminous infrared galaxies (ULIRGs). Conclusions. The wide range of parameters derived from our SED fitting, as well as part of our sample being detected as ULIRGs, seems to indicate that at these Lyman alpha luminosities, LAEs do not necessarily have to be young dwarfs, and that a lack of dust is not required for Lyman alpha to escape.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

The Ly$\alpha$ line is an important diagnostic for star formation at high redshift, but interpreting its flux and line profile is difficult because of the resonance nature of Ly$\alpha$. Trends between the escape of Ly$\alpha$ photons and dust and properties of the interstellar medium (ISM) have been found, but detailed comparisons between Ly$\alpha$ emission and the properties of the gas in local high-redshift analogs are vital for understanding the relation between Ly$\alpha$ emission and galaxy properties. For the first time, we can directly infer the properties of the ionized gas at the same location and similar spatial scales of the extended Ly$\alpha$ halo around ESO 338-IG04. We obtained VLT/MUSE integral field spectra. We used ionization parameter mapping of the [SII]/[OIII] line ratio and the kinematics of H$\alpha$ to study the ionization state and kinematics of the ISM of ESO338-IG04. The velocity map reveals two outflows. The entire central area of the galaxy is highly ionized by photons leaking from the HII regions around the youngest star clusters. Three highly ionized cones have been identified, of which one is associated with an outflow. We propose a scenario where the outflows are created by mechanical feedback of the older clusters, while the highly ionized gas is caused by the hard ionizing photons emitted by the youngest clusters. A comparison with the Ly$\alpha$ map shows that the (approximately bipolar) asymmetries observed in the Ly$\alpha$ emission are consistent with the base of the outflows detected in H$\alpha$. No clear correlation with the ionization cones is found. The mechanical and ionization feedback of star clusters significantly changes the state of the ISM by creating ionized cones and outflows. The comparison with Ly$\alpha$ suggests that especially the outflows could facilitate the escape of Ly$\alpha$ photons [Abridged].

We present high-resolution far-UV spectroscopy of the 14 galaxies of the Lyman Alpha Reference Sample; a sample of strongly star-forming galaxies at low redshifts ($0.028 < z < 0.18$). We compare the derived properties to global properties derived from multi band imaging and 21 cm HI interferometry and single dish observations, as well as archival optical SDSS spectra. Besides the Lyman $\alpha$ line, the spectra contain a number of metal absorption features allowing us to probe the kinematics of the neutral ISM and evaluate the optical depth and and covering fraction of the neutral medium as a function of line-of-sight velocity. Furthermore, we show how this, in combination with precise determination of systemic velocity and good Ly$\alpha$ spectra, can be used to distinguish a model in which separate clumps together fully cover the background source, from the "picket fence" model named by Heckman et al. (2011). We find that no one single effect dominates in governing Ly$\alpha$ radiative transfer and escape. Ly$\alpha$ escape in our sample coincides with a maximum velocity-binned covering fraction of $\lesssim 0.9$ and bulk outflow velocities of $\gtrsim 50$ km s$^{-1}$, although a number of galaxies show these characteristics and yet little or no Ly$\alpha$ escape. We find that Ly$\alpha$ peak velocities, where available, are not consistent with a strong backscattered component, but rather with a simpler model of an intrinsic emission line overlaid by a blueshifted absorption profile from the outflowing wind. Finally, we find a strong anticorrelation between H$\alpha$ equivalent width and maximum velocity-binned covering factor, and propose a heuristic explanatory model.

Population III galaxies are predicted to exist at high redshifts and may be rendered sufficiently bright for detection with current telescopes when gravitationally lensed by a foreground galaxy cluster. Population III galaxies that exhibit strong Lya emission should furthermore be identifiable from broadband photometry because of their unusual colors. Here, we report on a search for such objects at z > 6 in the imaging data from the Cluster Lensing And Supernova survey with Hubble (CLASH), covering 25 galaxy clusters in 16 filters. Our selection algorithm returns five singly-imaged candidates with Lya-like color signatures, for which ground-based spectroscopy with current 8-10 m class telescopes should be able to test the predicted strength of the Lya line. None of these five objects have been included in previous CLASH compilations of high-redshift galaxy candidates. However, when large grids of spectral synthesis models are applied to the study of these objects, we find that only two of these candidates are significantly better fitted by Population III models than by more mundane, low-metallicity stellar populations.

The Lyman-alpha reference sample (LARS) is a program with the Hubble Space Telescope (HST) that provides a sample of local universe laboratory galaxies in which to study the astrophysics of the visibility and strength of the Lyman-alpha (Lya) line of hydrogen. This article presents an overview of the survey, its selection function and HST imaging observations. The sample was selected from the GALEX+SDSS catalogue at z=0.028-0.19, in order to allow Lya to be captured with combinations of long pass filters in the Solar Blind Channel (SBC) of HST/ACS. In addition, LARS utilises Halpha and Hbeta narrow, and U, B, i broad-band imaging with ACS and WFC3. In order to study galaxies in which large numbers of Lya photons are produced we demanded an Halpha equivalent width > 100{\AA}. The sample of 14 galaxies covers far UV (FUV) luminosities that overlaps with those of high-z Lya emitters and Lyman Break Galaxies (LBGs). We present the reduction steps used to obtain the Lya images, including our LARS eXtraction software (LaXs) which utilises pixel-by-pixel spectral synthesis fitting to subtract the continuum at Lya. We demonstrate that the use of SBC long pass filters dramatically increase the signal to noise compared to the nominal Lya filter in SBC. To exemplify, we present results for Mrk 259. This irregular galaxy shows extended but strongly asymmetric Lya emission. Spectroscopy from HST/COS show a moderate outflow in the neutral ISM and an asymmetric Lya profile. Radiative transfer modeling is able to reproduce the Lya line profile, confirming the presence of an outflow. The integrated photometry provides a Lya luminosity of 1.3E42 erg/s an equivalent width W(Lya) =45{\AA} and a FUV absolute magnitude -19.2. Mrk 259 would hence be detectable in high-z Lya and LBG surveys. The total Lya escape fraction is 12%, similar to that observed at z>4. (shortened)

We present new H I imaging and spectroscopy of the 14 UV-selected star-forming galaxies in the Lyman Alpha Reference Sample (LARS), aimed for a detailed study of the processes governing the production, propagation, and escape of Ly$\alpha$ photons. New H I spectroscopy, obtained with the 100m Green Bank Telescope (GBT), robustly detects the H I spectral line in 11 of the 14 observed LARS galaxies (although the profiles of two of the galaxies are likely confused by other sources within the GBT beam); the three highest redshift galaxies are not detected at our current sensitivity limits. The GBT profiles are used to derive fundamental H I line properties of the LARS galaxies. We also present new pilot H I spectral line imaging of 5 of the LARS galaxies obtained with the Karl G. Jansky Very Large Array (VLA). This imaging localizes the H I gas and provides a measurement of the total H I mass in each galaxy. In one system, LARS 03 (UGC 8335 or Arp 238), VLA observations reveal an enormous tidal structure that extends over 160 kpc from the main interacting systems and that contains $>$10$^9$ M$_{\odot}$ of H I. We compare various H I properties with global Ly$\alpha$ quantities derived from HST measurements. The measurements of the Ly$\alpha$ escape fraction are coupled with the new direct measurements of H I mass and significantly disturbed H I velocities. Our robustly detected sample reveals that both total H I mass and linewidth are tentatively correlated with key Ly$\alpha$ tracers. Further, on global scales, these data support a complex coupling between Ly$\alpha$ propagation and the H I properties of the surrounding medium.

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

We report upon new results regarding the Lya output of galaxies, derived from the Lyman alpha Reference Sample (LARS), focusing on Hubble Space Telescope imaging. For 14 galaxies we present intensity images in Lya, Halpha, and UV, and maps of Halpha/Hbeta, Lya equivalent width (EW), and Lya/Halpha. We present Lya and UV light profiles and show they are well-fitted by S\'ersic profiles, but Lya profiles show indices systematically lower than those of the UV (n approx 1-2 instead of >~4). This reveals a general lack of the central concentration in Lya that is ubiquitous in the UV. Photometric growth curves increase more slowly for Lya than the FUV, showing that small apertures may underestimate the EW. For most galaxies, however, flux and EW curves flatten by radii ~10 kpc, suggesting that if placed at high-z, only a few of our galaxies would suffer from large flux losses. We compute global properties of the sample in large apertures, and show total luminosities to be independent of all other quantities. Normalized Lya throughput, however, shows significant correlations: escape is found to be higher in galaxies of lower star formation rate, dust content, mass, and several quantities that suggest harder ionizing continuum and lower metallicity. Eight galaxies could be selected as high-z Lya emitters, based upon their luminosity and EW. We discuss the results in the context of high-z Lya and UV samples. A few galaxies have EWs above 50 AA, and one shows f_escLya of 80%; such objects have not previously been reported at low-z.

Aims: We aim to study the effects of an inhomogeneous interstellar medium (ISM) on the strength and the shape of the Lyman alpha (Lya) line in starburst galaxies. Methods: Using our 3D Monte Carlo Lya radiation transfer code, we study the radiative transfer of Lya, UV and optical continuum photons in homogeneous and clumpy shells of neutral hydrogen and dust surrounding a central source. Our simulations predict the Lya and continuum escape fraction, the Lya equivalent width EW(Lya), the Lya line profile and their dependence on the gas geometry and the main input physical parameters. Results: The ISM clumpiness is found to have a strong impact on the Lya line radiative transfer, entailing a strong dependence of the emergent features of the Lya line (escape fraction, EW(Lya)) on the ISM morphology. Although a clumpy and dusty ISM appears more transparent to radiation (both line and continuum) compared to an equivalent homogeneous ISM of equal dust optical depth, we find that the Lya photons are, in general, still more attenuated than UV continuum radiation. As a consequence, the observed Lya equivalent width (EWobs(Lya)) is lower than the intrinsic one (EWint(Lya)) for nearly all clumpy ISM configurations considered. There are, however, special conditions under which Lya photons escape more easily than the continuum, resulting in an enhanced EWobs(Lya). The requirement for this to happen is that the ISM is almost static (galactic outflows < 200 km/s), extremely clumpy (with density contrasts >10^7 in HI between clumps and the interclump medium), and very dusty (E(B-V) > 0.30). When these conditions are fulfilled the emergent Lya line profile shows no velocity shift and little asymmetry. Given the asymmetry and velocity shifts generally observed in starburst galaxies with Lya emission, we conclude that clumping is unlikely to significantly enhance their relative Lya/UV transmission.

Context. The Lyman alpha emission line of galaxies is an important tool for finding galaxies at high redshift, and thus probe the structure of the early universe. However, the resonance nature of the line and its sensitivity to dust and neutral gas is still not fully understood. Aims. We present measurements of the velocity, covering fraction and optical depth of neutral gas in front of two well known local blue compact galaxies that show Lyman alpha in emission: ESO 338-IG 04 and Haro 11. We thus test observationally the hypothesis that Lyman alpha can escape through neutral gas by being Doppler shifted out of resonance. Methods. We present integral field spectroscopy from the GIRAFFE/Argus spectrograph at VLT/FLAMES in Paranal, Chile. The excellent wavelength resolution allows us to accurately measure the velocity of the ionized and neutral gas through the H-alpha emission and Na D absorption, which traces the ionized medium and cold interstellar gas, respectively. We also present independent measurements with the VLT/X-shooter spectrograph which confirm our results. Results. For ESO 338-IG04, we measure no significant shift of neutral gas. The best fit velocity is -15 (16) km/s. For Haro 11, we see an outflow from knot B at 44 (13) km/s and infalling gas towards knot C with 32 (12) km/s. Based on the relative strength of the Na D absorption lines, we estimate low covering fractions of neutral gas (down to 10%) in all three cases. Conclusions. The Na D absorption likely occurs in dense clumps with higher column densities than where the bulk of the Ly-alpha scattering takes place. Still, we find no strong correlation between outflowing neutral gas and a high Lyman alpha escape fraction. The Lyman alpha photons from these two galaxies are therefore likely escaping due to a low column density and/or covering fraction.

Using high-resolution Hubble Space Telescope data, we reexamine the fundamental properties (ages, masses and extinction values) of the rich star cluster population in the dwarf starburst galaxy NGC 5253. The gain in resolution compared to previous studies is of order a factor of two in both spatial dimensions, while our accessible wavelength range transcends previous studies by incorporation of both near-ultraviolet and near-infrared (IR) passbands. We apply spectral synthesis treatments based on two different simple stellar population model suites to our set of medium-, broad-band and H-alpha images to gain an improved physical understanding of the IR-excess flux found for a subset of young clusters (30 of 149). With the caveat that our models are based on fully sampled stellar mass functions, the NGC 5253 cluster population is dominated by a significant number of relatively low-mass (M_cl <= a few 10^4 Msun) objects with ages ranging from a few x 10^6 to a few x 10^7 yr, which is in excellent agreement with the starburst age of the host galaxy. The IR-excess clusters are almost all found in this young age range and have masses of up to a few x 10^4 Msun. The IR excess in the relatively low-mass NGC 5253 clusters is most likely caused by a combination of stochastic sampling effects and colour variations due to the presence of either luminous red or pre-main-sequence stars. We also find a small number of intermediate-age (~ 1 Gyr-old), ~10^5 Msun clusters, as well as up to a dozen massive, ~10 Gyr-old globular clusters. Their presence supports the notion that NGC 5253 is a very active galaxy that has undergone multiple episodes of star cluster formation.

We report on new imaging observations of the Lyman alpha emission line (Lya), performed with the Hubble Space Telescope, that comprise the backbone of the Lyman alpha Reference Sample (LARS). We present images of 14 starburst galaxies at redshifts 0.028 < z < 0.18 in continuum-subtracted Lya, Halpha, and the far ultraviolet continuum. We show that Lya is emitted on scales that systematically exceed those of the massive stellar population and recombination nebulae: as measured by the Petrosian 20 percent radius, RP20, Lya radii are larger than those of Halpha by factors ranging from 1 to 3.6, with an average of 2.4. The average ratio of Lya-to-FUV radii is 2.9. This suggests that much of the Lya light is pushed to large radii by resonance scattering. Defining the "Relative Petrosian Extension" of Lya compared to Halpha, \xi_ext = RP20_Lya / RP20_Ha, we find \xi_ext to be uncorrelated with total Lya luminosity. However \xi_ext is strongly correlated with quantities that scale with dust content, in the sense that a low dust abundance is a necessary requirement (although not the only one) in order to spread Lya photons throughout the interstellar medium and drive a large extended Lya halo.

It has been suggested that radiative transfer effects may explain the unusually high equivalent widths (EWs) of the Lya line, observed occasionally from starburst galaxies, especially at high redshifts. If the dust is locked up inside high-density clouds dispersed in an empty intercloud medium, the Lya photons could scatter off of the surfaces of the clouds, effectively having their journey confined to the dustless medium. The continuum radiation, on the other hand, does not scatter, and would thus be subject to absorption inside the clouds. This scenario is routinely invoked when Lya EWs higher than what is expected theoretically are observed, although the ideal conditions under which the results are derived usually are not considered. Here we systematically examine the relevant physical parameters in this idealized framework, testing whether any astrophysically realistic scenarios may lead to such an effect. It is found that although clumpiness indeed facilitates the escape of Lya, it is highly unlikely that any real interstellar media should result in a preferential escape of Lya over continuum radiation. Other possible causes are discussed, and it is concluded that the observed high EWs are more likely to be caused by cooling radiation from cold accretion and/or anisotropic escape of the Lya radiation.

[Abridged] We present deep optical and near-infrared UBVRIHKs imaging data for 24 blue compact galaxies (BCGs). The sample contains luminous dwarf and intermediate-mass BCGs which are predominantly metal-poor, although a few have near-solar metallicities. We have analyzed isophotal and elliptical integration surface brightness and color profiles, extremely deep (mu_B<29 mag arcsec^{-2}) contour maps and RGB images for each galaxy in the sample. The colors are compared to different spectral evolutionary models. We detect extremely extended low surface brightness (LSB) components dominant beyond the Holmberg radius as well as optical bridges between companion galaxies at the mu_V~28th mag arcsec^{-2} isophotal level. The central surface brightness mu_0 and scale length h_r are derived from two radial ranges typically assumed to be dominated by the underlying host galaxy. We find that mu_0 and h_r of the BCGs host deviate from those of dwarf ellipticals (dE) and dwarf irregulars (dI) solely due to a strong burst contribution to the surface brightness profile almost down to the Holmberg radius. Structural parameters obtained from a fainter region, mu_B=26-28 mag arcsec^{-2}, are consistent with those of true LSB galaxies for the starbursting BCGs in our sample, and with dEs and dIs for the BCGs with less vigorous star formation.

We present deep surface photometry of a volume--limited sample of 21 UM emission line galaxies in broadband optical UBVRI and near infra-red (NIR) HKs filters. The sample comprises 19 blue compact galaxies (BCGs) and two spirals. For some targets the exposure times are the deepest to date. For the BCG UM462 we observe a previously undetected second disk component beyond a surface brightness level of mu_B=26 mag arcsec^{-2}. This is a true low surface brightness component with central surface brightness mu_0=24.1 mag arcsec^{-2} and scale length h_r=1.5 kpc. All BCGs are dwarfs, with M_B>=-18, and very compact, with an average scale length of h_r~1 kpc. We separate the burst and host populations for each galaxy and compare them to stellar evolutionary models with and without nebular emission contribution. We also measure the A_{180} asymmetry in all filters and detect a shift from optical to NIR in the average asymmetry of the sample. This shift seems to be correlated with the morphological class of the BCGs. Using the color-asymmetry relation, we identify five BCGs in the sample as mergers, which is confirmed by their morphological class. Though clearly separated from normal galaxies in the concentration-asymmetry parameter space, we find that it is not possible to distinguish luminous starbursting BCGs from the merely star forming low luminosity BCGs.

We present supernova rate measurements at redshift 0.1-1.0 from the Stockholm VIMOS Supernova Survey (SVISS). The sample contains 16 supernovae in total. The discovered supernovae have been classified as core collapse or type Ia supernovae (9 and 7, respectively) based on their light curves, colour evolution and host galaxy photometric redshift. The rates we find for the core collapse supernovae are 3.29 (-1.78,-1.45)(+3.08,+1.98) x 10^-4 yr^-1 Mpc^-3 h70^3 (with statistical and systematic errors respectively) at average redshift 0.39 and 6.40 (-3.12,-2.11)(+5.30,+3.65) x 10^-4 yr^-1 Mpc^-3 h70^3 at average redshift 0.73. For the type Ia supernovae we find a rate of 1.29 (-0.57,-0.28)(+0.88,+0.27) x 10^-4 yr^-1 Mpc^-3 h70^3 at average redshift 0.62. All of these rate estimates have been corrected for host galaxy extinction, using a method that includes supernovae missed in infrared bright galaxies at high redshift. We use Monte Carlo simulations to make a thorough study of the systematic effects from assumptions made when calculating the rates and find that the most important errors come from misclassification, the assumed mix of faint and bright supernova types and uncertainties in the extinction correction. We compare our rates to other observations and to the predicted rates for core collapse and type Ia supernovae based on the star formation history and different models of the delay time distribution. Overall, our measurements, when taking the effects of extinction into account, agree quite well with the predictions and earlier results. Our results highlight the importance of understanding the role of systematic effects, and dust extinction in particular, when trying to estimate the rates of supernovae at moderate to high redshift.

We present optical (UBVJ), ultraviolet (FUV, NUV), and high resolution atomic hydrogen (HI) observations of the nearby blue compact dwarf (BCD), VII Zw 403. We find that VII Zw 403 has a relatively high HI mass-to-light ratio for a BCD. The rotation velocity is nominally 10-15 km/s, but rises to ~20 km/s after correction for the ~8-10 km/s random motions present in the gas. The velocity field is complex; including a variation in the position angle of the major axis going from the NE to the SW parts of the galaxy. Our high resolution HI maps reveal structure in the central gas, including a large, low-density HI depression or hole between the southern and northern halves of the galaxy, coincident with an unresolved x-ray source. Although interactions have been proposed as the triggering mechanism for the vigorous star formation occurring in BCDs, VII Zw 403 does not seem to have been tidally triggered by an external interaction, as we have found no nearby possible perturbers. It also doesn't appear to fall in the set of galaxies that exhibit a strong central mass density concentration, as its optical scale length is large in comparison to similar systems. However, there are some features that are compatible with an accretion event: optical/HI axis misalignment, a change in position angle of the kinematic axis, and a complex velocity field.

We present the late-time optical light curve of the ejecta of SN 1987A measured from HST imaging observations spanning the past 17 years. We find that the flux from the ejecta declined up to around year 2001, powered by the radioactive decay of 44Ti. Then the flux started to increase, more than doubling by the end of 2009. We show that the increase is the result of energy deposited by X-rays produced in the interaction with the circumstellar medium. We suggest that the change of the dominant energy input to the ejecta, from internal to external, marks the transition from supernova to supernova remnant. The details of the observations and the modelling are described in the accompanying supplementary information.

Supernova surveys can be used to study a variety of subjects such as: (i) cosmology through type Ia supernovae (SNe), (ii) star-formation rates through core-collapse SNe, and (iii) supernova properties and their connection to host galaxy characteristics. The Stockholm VIMOS Supernova Survey (SVISS) is a multi-band imaging survey aiming to detect supernovae at redshift ~0.5 and derive thermonuclear and core-collapse supernova rates at high redshift. In this paper we present the supernovae discovered in the survey along with light curves and a photometric classification into thermonuclear and core-collapse types. To detect the supernovae in the VLT/VIMOS multi-epoch images, we used difference imaging and a combination of automatic and manual source detection to minimise the number of spurious detections. Photometry for the found variable sources was obtained and careful simulations were made to estimate correct errors. The light curves were typed using a Bayesian probability method and Monte Carlo simulations were used to study misclassification. We detected 16 supernovae, nine of which had a core-collapse origin and seven had a thermonuclear origin. The estimated misclassification errors are quite small, in the order of 5%, but vary with both redshift and type. The mean redshift of the supernovae is 0.58. Additionally, we found a variable source with a very extended light curve that could possibly be a pair instability supernova.

We study the multi-wavelength properties of a set of 171 Ly-alpha emitting candidates at redshift z = 2.25 found in the COSMOS field, with the aim of understanding the underlying stellar populations in the galaxies. We especially seek to understand what the dust contents, ages and stellar masses of the galaxies are, and how they relate to similar properties of Ly-alpha emitters at other redshifts. The candidates here are shown to have different properties from those of Ly-alpha emitters found at higher redshift, by fitting the spectral energy distributions (SEDs) using a Monte-Carlo Markov-Chain technique and including nebular emission in the spectra. The stellar masses, and possibly the dust contents, are higher, with stellar masses in the range log M_* = 8.5 - 11.0 M_sun and A_V = 0.0 - 2.5 mag. Young population ages are well constrained, but the ages of older populations are typically unconstrained. In 15% of the galaxies only a single, young population of stars is observed. We show that the Ly-alpha fluxes of the best fit galaxies are correlated with their dust properties, with higher dust extinction in Ly-alpha faint galaxies. Testing for whether results derived from a light-weighted stack of objects correlate to those found when fitting individual objects we see that stellar masses are robust to stacking, but ages and especially dust extinctions are derived incorrectly from stacks. We conclude that the stellar properties of Ly-alpha emitters at z = 2.25 are different from those at higher redshift and that they are diverse. Ly-alpha selection appears to be tracing systematically different galaxies at different redshifts.

Some of the first stars could be cooler and more massive than standard stellar models would suggest, due to the effects of dark matter annihilation in their cores. It has recently been argued that such objects may attain masses in the 10^4--10^7 solar mass range, and that such supermassive dark stars should be within reach of the upcoming James Webb Space Telescope. Notwithstanding theoretical difficulties with this proposal, we argue here that some of these objects should also be readily detectable with both the Hubble Space Telescope and ground-based 8--10 m class telescopes. Existing survey data already place strong constraints on 10^7 solar mass dark stars at z~10. We show that such objects must be exceedingly rare or short-lived to have avoided detection.

The first stars in the history of the Universe are likely to form in the dense central regions of 10^5-10^6 Msolar cold dark matter halos at z=10-50. The annihilation of dark matter particles in these environments may lead to the formation of so-called dark stars, which are predicted to be cooler, larger, more massive and potentially more long-lived than conventional population III stars. Here, we investigate the prospects of detecting high-redshift dark stars with the upcoming James Webb Space Telescope (JWST). We find that dark stars at z>6 are intrinsically too faint to be detected by JWST. However, by exploiting foreground galaxy clusters as gravitational telescopes, certain varieties of cool (Teff < 30000 K) dark stars should be within reach at redshifts up to z=10. If the lifetimes of dark stars are sufficiently long, many such objects may also congregate inside the first galaxies. We demonstrate that this could give rise to peculiar features in the integrated spectra of galaxies at high redshifts, provided that dark stars make up at least 1 percent of the total stellar mass in such objects.

We present the integrated Halpha emission line profile for 157 HII regions in the central 3.4' x 3.4' of the galaxy M 83 (NGC 5236). Using the Fabry-Perot interferometer GHaFaS, on the 4.2 m William Herschel on La Palma, we show the importance of a good characterization of the instrumental response function for the study of line profile shapes. The luminosity-velocity dispersion relation is also studied, and in the log(L)-log(sigma) plane we do not find a linear relation, but an upper envelope with equation log(L)=0.9 *log(sigma)+38.1. For the adopted distance of 4.5 Mpc, the upper envelope appears at the luminosity L=10^38.5 ergs, in full agreement with previous studies of other galaxies, reinforcing the idea of using HII regions as standard candles.

Previous studies of the low surface brightness host of the blue compact galaxy (BCG) Haro 11 have suggested an abnormally red color of V-K=4.2+-0.8 for the host galaxy. This color is inconsistent with any normal stellar population over a wide range of stellar metallicities (Z=0.001-0.02). Similar though less extreme host colors have been measured for other BCGs and may be reconciled with population synthesis models, provided that the stellar metallicity of the host is higher than that of the ionized gas in the central starburst. We present the deepest V and K band observations to date of Haro 11 and derive a new V-K color for the host galaxy. Our new data suggest a far less extreme colour of V-K=2.3+-0.2, which is perfectly consistent with the expectations for an old host galaxy with the same metallicty as that derived from nebular emission lines in the star-forming center.

We have performed a deep multi-band photometric analysis of the star cluster population of Haro 11. This starburst galaxy (log L_FUV = 10.3 L_sun) is considered a nearby analogue of Lyman break galaxies (LBGs) at high redshift. The study of the numerous star clusters in the systems is an effective way to investigate the formation and evolution of the starburst phase. In fact, the SED fitting models have revealed a surprisingly young star cluster population, with ages between 0.5 and 40 Myr, and estimated masses between 10^3 and 10^7 solar masses. An independent age estimation has been done with the EW(Halpha) analysis of each cluster. This last analysis has confirmed the young ages of the clusters. We noticed that the clusters with ages between 1 and 10 Myr show a flux excess in H (NIC3/F160W) and/or I (WFPC2/F814W) bands with respect to the evolutionary models. Once more Haro 11 represents a challenge to our understanding.

Previous studies of the host galaxy of Haro 11 have suggested an extreme red color of V-K=4.2+-0.8 which cannot be reconciled with any normal stellar population of reasonable metallicity. We present the deepest V and K band data to date of the blue compact galaxy (BCG) Haro 11 and derive a new V-K color for the host. Our new data suggest a much more modest value of V-K=2.3+-0.2, which is in the same range as V-K colors measured for several other BCGs. The new Haro 11 color is not abnormally red and can be attributed to an old metal-poor stellar population with a Salpeter initial mass function.

We present preliminary results from a study of the compact star cluster populations in three local luminous blue compact galaxies: ESO 185-IG 013, ESO 350-IG 038 (a.k.a. Haro 11), and MRK 930. These systems show peculiar morphologies and the presence of hundreds of SCs that have been produced by the past, recent, and/or current starburst phases. We use a complete set of HST images ranging from the UV to IR for each galaxy. Deep images in V (WFPC2/f606w) and I (WFPC2/f814w) are used to capture most of the star cluster candidates up to the old ones (fainter) which have had, in the past, less possibility to be detected. The other bands are used in the SED fitting technique for constraining ages and masses. Our goals are to investigate the evolution of these three blue compact galaxies and the star cluster formation impact on their star formation history.

In a programme of observations of local luminous blue compact galaxies (BCGs), we are investigating kinematics by using tracers of both stars and ionized gas. Here we summarise our program and present new data on the local Lyman break galaxy analogue Haro 11. From spatially-resolved spectroscopy around the near-infrared Ca II triplet, we find that its stars and ionized gas have similar velocity fields. Our programme so far indicates however that emission line velocities can differ locally by a few tens of km/s from the Ca II values. Comparing our data to simple stellar population models, we assess which stellar population the Ca II triplet traces and its potential beyond the local universe.

Intermediate mass galaxies (logM(Msun)>10) at z~0.6 are the likeliest progenitors of the present-day numerous population of spirals. There is growing evidence that they have evolved rapidly since the last 6 to 8 Gyr ago, and likely have formed a significant fraction of their stellar mass, often showing perturbed morphologies and kinematics. We have gathered a representative sample of 88 such galaxies and have provided robust estimates of their gas phase metallicity. For doing so, we have used moderate spectral resolution spectroscopy at VLT/FORS2 with unprecedented high S/N allowing to remove biases coming from interstellar absorption lines and extinction to establish robust values of R23=([OII]3727 + [OIII]4959,5007)/Hbeta. We definitively confirm that the predominant population of z~0.6 starbursts and luminous IR galaxies (LIRGs) are on average, two times less metal rich than the local galaxies at a given stellar mass. We do find that the metal abundance of the gaseous phase of galaxies is evolving linearly with time, from z=1 to z=0 and after comparing with other studies, from z=3 to z=0. Combining our results with the reported evolution of the Tully Fisher relation, we do find that such an evolution requires that ~30% of the stellar mass of local galaxies have been formed through an external supply of gas, thus excluding the close box model. Distant starbursts & LIRGs have properties (metal abundance, star formation efficiency & morphologies) similar to those of local LIRGs. Their underlying physics is likely dominated by gas infall probably through merging or interactions. Our study further supports the rapid evolution of z~0.4-1 galaxies. Gas exchanges between galaxies is likely the main cause of this evolution.

The aim of the work presented in this paper is to test and optimise supernova detection methods based on the optimal image subtraction technique. The main focus is on applying the detection methods to wide field supernova imaging surveys and in particular to the Stockholm VIMOS Supernova Survey (SVISS). We have constructed a supernova detection pipeline for imaging surveys. The core of the pipeline is image subtraction using the ISIS 2.2 package. Using real data from the SVISS we simulate supernovae in the images, both inside and outside galaxies. The detection pipeline is then run on the simulated frames and the effects of image quality and subtraction parameters on the detection efficiency and photometric accuracy are studied. The pipeline allows efficient detection of faint supernovae in the deep imaging data. It also allows controlling and correcting for possible systematic effects in the SN detection and photometry. We find such a systematic effect in the form of a small systematic flux offset remaining at the positions of galaxies in the subtracted frames. This offset will not only affect the photometric accuracy of the survey, but also the detection efficiencies. Our study has shown that ISIS 2.2 works well for the SVISS data. We have found that the detection efficiency and photometric accuracy of the survey are affected by the stamp selection for the image subtraction and by host galaxy brightness. With our tools the subtraction results can be further optimised, any systematic effects can be controlled and photometric errors estimated, which is very important for the SVISS, as well as for future SN searches based on large imaging surveys such as Pan-STARRS and LSST.

Luminous blue compact galaxies, common at z~1 but now relatively rare, show disturbed kinematics in emission lines. As part of a programme to understand their formation and evolution, we have investigated the stellar dynamics of a number of nearby objects in this class. We obtained long-slit spectra with VLT/FORS2 in the spectral region covering the near-infrared calcium triplet. In this paper we focus on the well-known luminous blue compact galaxy ESO 338-IG04 (Tololo 1924-416). A previous investigation, using Fabry-Perot interferometry, showed that this galaxy has a chaotic H-alpha velocity field, indicating that either the galaxy is not in dynamical equilibrium or that H-alpha does not trace the gravitational potential due to feedback from star formation. Along the apparent major axis, the stellar and ionised gas velocities for the most part follow each other. The chaotic velocity field must therefore be a sign that the young stellar population in ESO 338-IG04 is not in dynamical equilibrium. The most likely explanation, which is also supported by its morphology, is that the galaxy has experienced a merger and that this has triggered the current starburst. Summarising the results of our programme so far, we note that emission-line velocity fields are not always reliable tracers of stellar motions, and go on to assess the implications for kinematic studies of similar galaxies at intermediate redshift.

(abbreviated) We present the first results of the ESO large program, ``IMAGES'' which aims at obtaining robust measurements of the kinematics of distant galaxies using the multi-IFU mode of GIRAFFE on the VLT. 3D spectroscopy is essential to robustly measure the often distorted kinematics of distant galaxies (e.g., Flores et al. 2006). We derive the velocity fields and $\sigma$-maps of 36 galaxies at 0.4<z<0.75 from the kinematics of the [OII] emission line doublet, and generate a robust technique to identify the nature of the velocity fields based on the pixels of the highest signal-to-noise ratios (S/N). We have gathered a unique sample of 63 velocity fields of emission line galaxies (W0([OII]) > 15 A) at z=0.4-0.75, which are a representative subsample of the population of M_stellar>1.5x10^{10} M_sun emission line galaxies in this redshift range, and are largely unaffected by cosmic variance. Taking into account all galaxies -with or without emission lines- in that redshift range, we find that at least 41+/-7% of them have anomalous kinematics, i.e., they are not dynamically relaxed. This includes 26+/-7% of distant galaxies with complex kinematics, i.e., they are not simply pressure or rotationally supported. Our result implies that galaxy kinematics are among the most rapidly evolving properties, because locally, only a few percent of the galaxies in this mass range have complex kinematics.

We study the two-dimensional distribution and kinematics of the stellar and gaseous components in the centre of the blue compact dwarf galaxy, He 2-10. The aim is to compare the kinematics of gas and stars in order to determine whether they are consistent with one another, or if stars and gas can be decoupled due to gravitational perturbations and feedback from star formation. We have used the integral field unit ARGUS, part of FLAMES on the European Southern Observatory's Very Large Telescope, to target the near-infrared Ca-triplet in the central 300 times 480 parsecs of He 2-10. The selected wavelength regime includes several prominent spectral features, including the Paschen series and an [SIII] emission-line, which we have used to derive the kinematics of the ionised interstellar medium. We find no systematic trend in the velocities of the stars over the observed field of view and conclude that the stellar kinematics is governed by random motions. This is in contrast to the motions the ionised interstellar medium, where we find spatial velocity variations up to 60 km/s. Our gas velocity field is consistent with previous studies of both the molecular gas and the feedback-driven outflow in He 2-10. We interpret the kinematic decoupling between the stars and the gas as He 2-10 being in the process of transformation to a dwarf elliptical galaxy.

We have studied the afterglow of the gamma-ray burst (GRB) of February 18, 2006. This is a nearby long GRB, with a very low peak energy, and is therefore classified as an X-ray Flash (XRF). XRF 060218 is clearly associated with a supernova -- dubbed SN 2006aj. We present early spectra for SN 2006aj as well as optical lightcurves reaching out to 50 days past explosion. Our optical lightcurves define the rise times, the lightcurve shapes and the absolute magnitudes in the U, V and R bands, and we compare these data with data for other relevant supernovae. SN 2006aj evolved quite fast, somewhat similarly to SN 2002ap, but not as fast as SN 1994I. Our spectra show the evolution of the supernova over the peak, when the U-band portion of the spectrum rapidly fades due to extensive line blanketing. We compare to similar spectra of very energetic Type Ic supernovae. Our first spectra are earlier than spectra for any other GRB-SN. The spectrum taken 12 days after burst in the rest frame is similar to somewhat later spectra of both SN 1998bw and SN 2003dh, implying a rapid early evolution. This is consistent with the fast lightcurve. From the narrow emission lines from the host galaxy we derive a redshift of z=0.0331+-0.0007. This makes XRF 060218 the second closest gamma-ray burst detected. The flux of these emission lines indicate a high-excitation state, and a modest metallicity and star formation rate of the host galaxy.

The optical/near-IR stellar continuum carries unique information about the stellar population in a galaxy, its mass function and star-formation history. Star-forming regions display rich emission-line spectra from which we can derive the dust and gas distribution, map velocity fields, metallicities and young massive stars and locate shocks and stellar winds. All this information is very useful in the dissection of the starburst phenomenon. We discuss a few of the advantages and limitations of observations in the optical/near-IR region and focus on some results. Special attention is given to the role of interactions and mergers and observations of the relatively dust-free starburst dwarfs. In the future we expect new and refined diagnostic tools to provide us with more detailed information about the IMF, strength and duration of the burst and its triggering mechanisms.

We present far infrared (FIR) spectroscopy of the luminous blue compact galaxy (BCG) Haro 11 (ESO 350-IG38) obtained with the ISO Long Wavelength Spectrometer (LWS) in low resolution mode. This metal poor dwarf merger is an extremely hot IRAS source with a high [CII]158um/CO(1-0) flux ratio. We discuss the balance between dust and line cooling in the photodissociated regions (PDR), in particular the role of the [CII] line, and derive the basic properties of the PDR gas. While samples of other starburst galaxies show a correlation between the [CII]/FIR flux ratio and the IRAS f60/f100 ratio, Haro 11 deviates significantly from this relationship being brighter in [CII] than average. We propose that the relationship is caused by an increasing optical depth with increasing IRAS temperature but that the low metallicity and the extreme starburst properties of Haro 11 allow the medium to be thin despite its high f60/f100 ratio, explaining the [CII] excess. This leaves room for a more optimistic view on the possibilities to detect massive starforming mergers at high redshifts, using the [CII] line.

Observations of six luminous blue compact galaxies (BCGs) and two star forming companion galaxies were carried out with the CIGALE scanning Fabry-Perot interferometer attached to the ESO 3.6m telescope on La Silla. The observations were made in the H-alpha emission line which is prominent in BCGs. A velocity sampling of 5 km/s and a pixel size of 0.9 arcseconds were used. In this paper we present the observations and the data together with the velocity fields and the derived rotation curves. In addition we provide rough estimates of the total dynamical mass and of the ionised gas mass for each galaxy. All galaxies display rotation, but while the companion galaxies have regular velocity fields, those of the BCGs are complex and appear perturbed. This is the most extensive study to date of the optical velocity fields of BCGs. The interpretation of these results will be presented in a forthcoming paper (Paper II).

Representing single stellar populations, globular clusters (GCs) are relatively easy to model, thus providing powerful tools for studying the evolution of galaxies. This has been demonstrated for the blue compact galaxy ESO338-IG04. GC systems in galaxies may be fossils of starbursts and mergers. Thus studies of GCs in the local universe may add to our understanding of the formation and evolution of galaxies and the distant universe.

Multicolour images of the starbursting metal poor blue compact galaxy ESO338-IG04 have been obtained with WFPC2/HST. In the images we find numerous point-like sources concentrated towards the main body of the galaxy, which we identify as globular cluster candidates. We show that these objects are physically associated with the galaxy and that they are spatially extended. Given their high intrinsic luminosities, these objects cannot be individual stars. Using photometric evolution models we show that the objects constitute a rich population of massive star clusters with ages ranging from a few Myr to 10 Gyr, and masses ranging from 10^4 to more than 10^7 solar masses. We show that these objects are most likely globular clusters of varying age. There are peaks in the age distribution of the clusters: at <30 Myr, ~100 Myr, ~600 Myr, 2.5-5 Gyr and ~10 Gyr. The youngest objects are found in the crowded starburst region. The galaxy presently appears to be involved in a merger, which is the probable cause of the present globular cluster formation. The presence of a numerous intermediate age (2.5 to 5 Gyr) population of globular clusters, suggests that a previous merger might have occurred. As the starburst fades, this galaxy will become very rich in globular clusters. Transforming all objects to an age comparable to that of Milky Way globular clusters reveals a luminosity function similar to the Galactic. We suggest that this galaxy is the result of a merger between a dwarf elliptical and a gas rich dwarf. The possibility of dating the globular clusters offers an efficient way of studying the history of violent star formation in this and similar galaxies.