Using a sample of galaxies at $z\approx 3$ with detected Lyman Continuum (LyC) leakage in the SSA22 field, we attempt to verify a proposed indirect method for identifying cases with high LyC escape fraction \fesc\ based on measurements of the H$\beta$ equivalent width (EW) and the $\beta$ slope of the UV continuum. To this end, we present Keck/MOSFIRE H$\beta$ emission line flux measurements of LyC galaxies (LCGs) at spectroscopic redshifts $z_{\mathrm{spec}} \sim 3.3$, Lyman break galaxies (LBGs) at photometric redshifts $z_{\mathrm{phot}} = 2.7$--$3.7$, and Ly$\alpha$ emitters at $z_{\mathrm{phot}} = 3.1$. We also reconfirm the spectroscopic redshifts and measure the H$\beta$ emission line fluxes from 2 LCGs and 6 LBGs. For the LCG in our sample with the most extreme \fesc, as revealed by the direct detection of LyC photons, we find that the EW(H$\beta$)--$\beta$ method gives a broadly consistent estimate for \fesc, although the error bars remain very large. We also discuss how a combination of $f_\mathrm{esc}$ measurements based on direct and indirect methods can shed light on the LyC escape mechanism and the anisotropy of the leakage.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

We report results of a search for galaxies at z>3 with Lyman continuum (LyC) emission using a narrow-band filter NB359 with Subaru / Suprime-Cam in a ~800 arcmin$^2$ blank field around the GOODS-N. We use 103 star-forming galaxies (SFGs) and 8 AGNs with spectroscopic redshifts in a range between 3.06 and 3.5, and 157 photometrically selected z=3.1 Lyman $\alpha$ emitter (LAE) candidates as the targets. After removing galaxies spectroscopically confirmed to be contaminated by foreground sources, we found two SFGs and one AGN as candidate LyC emitting sources among the targets with spectroscopic redshifts. Among LAE candidates, five sources are detected in the NB359 image, and three among them may be contaminated by foreground sources. We compare the sample galaxies in the GOODS-N with those in the SSA22, where a prominent protocluster at z=3.1 is known and the LyC search using the same NB359 filter has been made. Frequency of galaxies with LyC leakage in the SSA22 field may be about two times higher than the galaxies in the GOODS-N with the sample UV magnitude range, although the numbers of LyC detections in these fields are too small to make a statistically significant conclusion. By combining the sample galaxies in these fields, we place the 3$\sigma$ upper limits of the observed LyC-to-UV flux density ratio and LyC escape fraction for galaxies at z=3.1 with absolute UV magnitude $M_{UV} < -18.8$ as ($f_{LyC}/f_{UV})_{obs} < 0.036$ and $f_{esc}^{abs}<8$\%, respectively.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

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]

Low redshift, spatially resolved Lyman continuum (LyC) emitters allow us to clarify the processes for LyC escape from these starburst galaxies. We use Hubble Space Telescope (HST) WFC3 and ACS imaging of the confirmed low-redshift LyC emitter Tol 1247-232 to study the ionization structure of the gas and its relation to the ionizing star clusters. We perform ionization parameter mapping (IPM) using [O III]4959, 5007 and [O II]3727 imaging as the high- and low-ionization tracers, revealing broad, large-scale, optically thin regions originating from the center, and reaching the outskirts of the galaxy, consistent with LyC escape. We carry out stellar population synthesis modeling of the 26 brightest clusters using our HST photometry. Combining these data with the nebular photometry, we find a global LyC escape fraction of f_esc = 0.12, with uncertainties also consistent with zero escape and with all measured f_esc values for this galaxy. Our analysis suggests that, similar to other candidate LyC emitters, a two-stage starburst has taken place in this galaxy, with a 12 Myr old, massive, central cluster likely having pre-cleared regions in and around the center, and the second generation of 2 - 4 Myr old clusters dominating the current ionization, including some escape from the galaxy.

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 report the detection of CO(J=2-1) coincident with the super star cluster (SSC) Mrk 71-A in the nearby Green Pea analog galaxy, NGC 2366. Our NOEMA observations reveal a compact, ~7 pc, molecular cloud whose mass (10^5 M_sun) is similar to that of the SSC, consistent with a high star-formation efficiency, on the order of 0.5. There are two, spatially distinct components separated by 11 km/s. If expanding, these could be due to momentum-driven, stellar wind feedback. Alternatively, we may be seeing the remnant infalling, colliding clouds responsible for triggering the SSC formation. The kinematics are also consistent with a virialized system. These extreme, high-density, star-forming conditions inhibit energy-driven feedback; the co-spatial existence of a massive, molecular cloud with the SSC supports this scenario, and we quantitatively confirm that any wind-driven feedback in Mrk 71-A is momentum-driven, rather than energy-driven. Since Mrk 71-A is a candidate Lyman continuum emitter, this implies that energy-driven superwinds may not be a necessary condition for the escape of ionizing radiation. In addition, the detection of the nebular continuum emission yields an accurate astrometric position for the Mrk 71-A. We also detect four other massive, molecular clouds in this giant star-forming complex.

We present the remarkable discovery that the dwarf irregular galaxy NGC 2366 is an excellent analog of the Green Pea (GP) galaxies, which are characterized by extremely high ionization parameters. The similarities are driven predominantly by the giant H II region Markarian 71 (Mrk 71). We compare the system with GPs in terms of morphology, excitation properties, specific star-formation rate, kinematics, absorption of low-ionization species, reddening, and chemical abundance, and find consistencies throughout. Since extreme GPs are associated with both candidate and confirmed Lyman continuum (LyC) emitters, Mrk 71/NGC 2366 is thus also a good candidate for LyC escape. The spatially resolved data for this object show a superbubble blowout generated by mechanical feedback from one of its two super star clusters (SSCs), Knot B, while the extreme ionization properties are driven by the <1 Myr-old, enshrouded SSC Knot A, which has ~ 10 times higher ionizing luminosity. Very massive stars (> 100 Msun) may be present in this remarkable object. Ionization-parameter mapping indicates the blowout region is optically thin in the LyC, and the general properties also suggest LyC escape in the line of sight. Mrk 71/NGC 2366 does differ from GPs in that it is 1 - 2 orders of magnitude less luminous. The presence of this faint GP analog and candidate LyC emitter (LCE) so close to us suggests that LCEs may be numerous and commonplace, and therefore could significantly contribute to the cosmic ionizing budget. Mrk 71/NGC 2366 offers an unprecedentedly detailed look at the viscera of a candidate LCE, and could clarify the mechanisms of LyC escape.

Using four different suites of cosmological simulations, we generate synthetic spectra for galaxies with different Lyman continuum escape fractions (fesc) at redshifts z=7-9, in the rest-frame wavelength range relevant for the James Webb Space Telescope (JWST) NIRSpec instrument. By investigating the effects of realistic star formation histories and metallicity distributions on the EW(Hb)-beta diagram (previously proposed as a tool for identifying galaxies with very high fesc), we find that neither of these effects are likely to jeopardize the identification of galaxies with extreme Lyman continuum leakage. Based on our models, we expect essentially all z=7-9 galaxies that exhibit rest-frame EW(Hb)< 30 {\AA} to have fesc>0.5. Incorrect assumptions concerning the ionizing fluxes of stellar populations or the dust properties of z>6 galaxies can in principle bias the selection, but substantial model deficiencies of this type will at the same time reveal themselves as an offset between the observed and simulated distribution of z>6 galaxies in the EW(Hb)-beta diagram. Such offsets would thereby allow JWST/NIRSpec measurements of these observables to serve as input for further model refinement.

We present the largest to date sample of hydrogen Lyman continuum (LyC) emitting galaxy candidates at any redshift, with 18 Ly$\alpha$ Emitters (LAEs) and 7 Lyman Break Galaxies (LBGs), obtained from the SSA22 field with Subaru/Suprime-Cam. The sample is based on the 159 LAEs and 136 LBGs observed in the field, all with spectroscopically confirmed redshifts, and these LyC candidates are selected as galaxies with counterpart in a narrow-band filter image which traces LyC at $z\geq 3.06$. Many LyC candidates show a spatial offset between the rest-frame non-ionizing ultraviolet (UV) detection and the LyC-emitting substructure or between the Ly$\alpha$ emission and LyC. The possibility of foreground contamination complicates the analysis of the nature of LyC emitters, although statistically it is highly unlikely that all candidates in our sample are contaminated by foreground sources. Many viable LyC LAE candidates have flux density ratios inconsistent with standard models, while also having too blue UV slopes to be foreground contaminants. Stacking reveals no significant LyC detection, suggesting that there is a dearth of objects with marginal LyC signal strength, perhaps due to a bimodality in the LyC emission. The foreground contamination-corrected $3\sigma$ upper limits of the observed average flux density ratios are $f_{LyC}/f_{UV}<0.08$ from stacking LAEs and $f_{LyC}/f_{UV}<0.02$ from stacking LBGs. There is a sign of a positive correlation between LyC and Ly$\alpha$, suggesting that both types of photons escape via a similar mechanism. The LyC detection rate among proto-cluster LBGs is seemingly lower compared to the field.

Subaru/SuprimeCam narrowband photometry of the SSA22 field reveals the presence of four Lyman continuum (LyC) candidates among a sample of 14 AGN. Two show offsets and likely have stellar LyC in nature or are foreground contaminants. The remaining two LyC candidates are Type I AGN. We argue that the average LyC escape fraction of high redshift low luminosity AGN is not likely to be unity, as often assumed in the literature. From direct measurement we obtain the average LyC-to-UV flux density ratio and ionizing emissivity for a number of AGN classes and find it at least a factor of two lower than values obtained assuming f_esc = 1. Comparing to recent Ly{\alpha} forest measurements, AGNs at redshift z~3 make up at most ~12% and as little as ~5% of the total ionizing budget. Our results suggest that AGNs are unlikely to dominate the ionization budget of the Universe at high redshifts.

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)

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

The stellar haloes of galaxies can currently be studied either through observations of resolved halo stars or through surface photometry. Curiously, the two methods appear to give conflicting results, as a number of surface photometry measurements have revealed integrated colours that are too red to be reconciled with the halo properties inferred from the study of resolved stars. Several explanations for this anomaly have been proposed - including dust photoluminescence, extinction of extragalactic background light and a bottom-heavy stellar initial mass function. A decisive test is, however, still lacking. Here, we explain how observations of the halo of a nearby galaxy, involving a combination of both surface photometry and bright star counts, can be used to distinguish between the proposed explanations. We derive the observational requirements for this endeavour and find that star counts in filters VI and surface photometry in filters VIJ appears to be the optimal strategy. Since the required halo star counts are already available for many nearby galaxies, the most challenging part of this test is likely to be the optical surface photometry, which requires several nights of exposure time on a 4-8 m telescope, and the near-IR surface photometry, which is most readily carried out using the upcoming James Webb Space Telescope.

We present a HST multiband analysis of the large population of star clusters in the blue compact galaxy (BCG) Mrk 930. We have reconstructed the spectral energy distributions of the star clusters and estimated age, mass, and extinction for a representative sample. We observe a very young cluster population with 70% of the systems formed less than 10 Myr ago. The peak in the star cluster age distribution at 4 Myr is corroborated by the presence of Wolf-Rayet spectral features, and by the observed optical and IR lines ratios [OIII]/H_beta and [NeIII]/[NeII]. The recovered extinction in these very young clusters shows large variations, with a decrease at older ages. It is likely that our analysis is limited to the optically brightest objects (i.e. systems only partially embedded in their natal cocoons). We map the extinction across the galaxy using low-resolution spectra and the H_alpha/H_beta ratio, as obtained from ground-based narrow band imaging. We find that the mean optical extinction derived in the starburst regions is close to the averaged value observed in the clusters, but locally, do not trace the more extinguished clusters. Previous HST studies of BCGs have revealed a population of young and extremely red super star clusters. We detect a considerable fraction of clusters affected by a red excess also in Mrk 930. The nature of the red excess, which turns up at near-IR wavelengths remains unknown. We compare the cluster and the star formation history, the latter derived from the fit of spectral population synthesis models to the spectra. We find a general agreement between the two independently estimated quantities. Using the cluster properties we perform a study of the host environmental properties. We find that the cluster formation efficiency is significantly higher, suggesting a key role of the environment for the formation of these massive objects.[Abridged]

Deep surface photometry of disk galaxies at optical and near-IR wavelengths have revealed faint halos with colours much too red to be reconciled with the resolved stellar populations detected in the halos of the Milky Way or M31. Alternative scenarios involving high metallicities, nebular emission or large amounts of dust in these halos are also disfavoured. Here, we argue that extinction of the optical extragalactic background light in the halos of these galaxies may possibly be responsible for the reported colour anomalies. We also discuss how an independent measurement of the optical extragalactic background light might be accomplished by combining direct star counts with surface photometry for a single nearby galaxy.

We have used multi-band imaging to investigate the nature of the extreme starburst environment in Haro 11 galaxy. The central starburst region has been observed in 8 HST wavebands and at 2.16 micron at the ESO-VLT. We constructed integrated spectral energy distributions (SEDs) for about 200 star clusters and compared them with single stellar population models in order to derive ages, masses and extinctions of thestar clusters. The present starburst has lasted for 40 Myr, and shows a peak of cluster formation only 3.5 Myr old. With such an extremely young cluster population, Haro 11 represents a unique opportunity to investigate the youngest phase of the cluster formation process and evolution in starburst systems. Extinction tends to diminish as function of the cluster age, but the spread is large and for clusters in partial embedded phases (< 5 Myr). A fraction of low-mass (> 10^4 Msun), very young (1-3 Myr) clusters is missing, either because they are embedded inthe parental molecular cloud and heavily extinguished, or because of blending. Almost half of the cluster sample is affected by flux excesses at wavelengths 8000 \AA which cannot be explained by simple stellar evolutionary models. Fitting SED models over all wavebands leads to systematic overestimates of cluster ages and incorrect masses for the stellar population supplying the light in these clusters. We show that the red excess affects also the HST F814W filter, which is typically used to constrain cluster physical properties. The clusters which show the red excess are younger than 40 Myr; we propose possible physical explanations for the phenomenon. Finally, we estimate that Haro 11 hasproduced bound clusters at a rate almost a factor of 10 higher than the massive and regular spirals, like the Milky Way. (Abriged)

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.

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.

The faint regions of galaxies, groups and clusters hold important clues about how these objects formed, and surface photometry at optical and near-infrared wavelengths represents a powerful tool for studying such structures. Here, we identify a hitherto unrecognized problem with this technique, related to how the night sky flux is typically measured and subtracted from astronomical images. While most of the sky flux comes from regions between the observer and the target object, a small fraction - the extragalactic background light (EBL) - comes from behind. We argue that since this part of the sky flux can be subject to extinction by dust present in the galaxy/group/cluster studied, standard reduction procedures may lead to a systematic oversubtraction of the EBL. Even very small amounts of extinction can lead to spurious features in radial surface surface brightness profiles and colour maps of extended objects. We assess the likely impact of this effect on a number of topics in extragalactic astronomy where very deep surface photometry is currently attempted, including studies of stellar halos, starburst host galaxies, disc truncations and diffuse intragroup/intracluster light. We argue that EBL extinction may provide at least a partial explanation for the anomalously red colours reported for the halos of disc galaxies and the hosts of local starburst galaxies. EBL extinction effects also mimic truncations in discs with unusually high dust opacities, but are unlikely to be the cause of such features in general. Failure to account for EBL extinction can also give rise to a non-negligible underestimate of intragroup and intracluster light at the faintest surface brightness levels currently probed. (Abridged)

Red halos are faint, extended and extremely red structures that have been reported around various types of galaxies since the mid-1990s. The colours of these halos are too red to be reconciled with any hitherto known type of stellar population, and instead indicative of a very bottom-heavy stellar initial mass function (IMF). Due to the large mass-to-light ratios of such stellar halos, they could contribute substantially to the baryonic masses of galaxies while adding very little to their overall luminosities. The red halos of galaxies therefore constitute potential reservoirs for some of the baryons still missing from inventories in the low-redshift Universe. While most studies of red halos have focused on disk galaxies, a red excess has also been reported in the faint outskirts of blue compact galaxies (BCGs). A bottom-heavy IMF can explain the colours of these structures as well, but due to model degeneracies, stellar populations with standard IMFs and abnormally high metallicities have also been demonstrated to fit the data. Here, we show that due to recent developments in the field of spectral synthesis, the metallicities required in this alternative scenario may be less extreme than previously thought. This suggests that the red excess seen in the outskirts of BCGs may stem from a normal, intermediate-metallicity host galaxy rather than a red halo of the type seen around disk galaxies. The inferred host metallicity does, however, still require the host to be more metal-rich than the gas in the central starburst of BCGs, in contradiction with current simulations of how BCGs form.

Deep optical/near-IR surface photometry of galaxies outside the Local Group have revealed faint and very red halos around objects as diverse as disk galaxies and starbursting dwarf galaxies. The colours of these structures are too extreme to be reconciled with stellar populations similar to those seen in the stellar halos of the Milky Way or M31, and alternative explanations like dust reddening, high metallicities or nebular emission are also disfavoured. A stellar population obeying an extremely bottom-heavy initial mass function (IMF), is on the other hand consistent with all available data. Because of its high mass-to-light ratio, such a population would effectively behave as baryonic dark matter and could account for some of the baryons still missing in the low-redshift Universe. Here, we give an overview of current red halo detections, alternative explanations for the origin of the red colours and ongoing searches for red halos around types of galaxies for which this phenomenon has not yet been reported. A number of potential tests of the bottom-heavy IMF hypothesis are also discussed.

Optical and near-IR observations of the halos of disk galaxies and blue compact galaxies have revealed a very red spectral energy distribution, which cannot easily be reconciled with a normal, metal-poor stellar population like that in the stellar halo of the Milky Way. Here, spectral evolutionary models are used to explore the consequences of these observations. We demonstrate that a stellar population of low to intermediate metallicity, but with an extremely bottom-heavy initial mass function, can explain the red halos around both types of objects. Other previously suggested explanations, like nebular emission or very metal-rich stars, are shown to fail in this respect. This indicates that, if the reported halo colours are correct, halo populations dominated by low-mass stars may be a phenomenon common to galaxies of very different Hubble types. Potential tests of this hypothesis are discussed, along with its implications for the baryonic dark matter content of galaxies.