An Arc of Young Stars in the Halo of M82

Early-type galaxies (ETGs) have been characterized as objects dominated by old stellar populations, containing little or no cold gas and dust, and thus, non-existent star formation. However, there are indications in the literature that some ETGs deviate from this: some have significant amounts of gas and dust, are forming stars, and/or display stellar substructures (tidal features, disks or shells, e.g., Kormendy et al. 1997, Rix, Carollo & Freeman 1999). A better understanding of the evolution of ETGs and the details of their “peculiarities” is critical to properly constrain models of galaxy formation. We present preliminary results on a photometric analysis of substructures in local ETGs, based on 3.6μm IRAC images from the Spitzer Survey of Stellar Structure in Galaxies (S4G; Sheth et al. 2010), which comprises one of the largest mid-IR photometric surveys of the local Universe. Relatively unhindered by extinction and dominated by the low-mass stellar populations that dominate a galaxy's stellar mass budget, the IR is the ideal waveband to trace the details of stellar structures in galaxies. Based on 2D GALFIT (Peng et al. 2002) decomposition, we find tidal features in 17% of 146 ETGs from S4G. For both the GALFIT model and the galaxy residual images, we calculate the total counts inside an annular region centered on the galaxy, where the inner radius is the effective radius of the galaxy. Assuming that a tidal feature and its host galaxy have the same mass-to-luminosity ratio (M/L), the ratio of the residual counts over model counts translates into the ratio of their stellar masses. We find that the tidal features in the majority of peculiar ETGs in our sample account for no more than 11% of the galaxy's total stellar mass. Considering that simulations (Canalizo et al. 2007) suggest an upper limit in relative stellar mass of 25% for shells resulting from a past major merger, the values we find support a merger origin. We are in the process of applying the decomposition method to GALEX UV images and optical SDSS images of these peculiar ETGs in order to characterize the underlying substructure and provide constraints on astrophysical properties such as star formation rates and stellar masses associated to these tidal features, based on broad-band SED template fitting techniques.

The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type galaxies. These discoveries motivate the hypothesis that the IMF may vary (possibly very slightly) across galaxies and across components of individual galaxies (e.g. bulges vs discs). In this paper we use a sample of 19 late-type strong gravitational lenses from the SWELLS survey to investigate the IMFs of the bulges and discs in late-type galaxies. We perform a joint analysis of the galaxies' total masses (constrained by strong gravitational lensing) and stellar masses (constrained by optical and near-infrared colours in the context of a stellar population synthesis [SPS] model, up to an IMF normalisation parameter). Using minimal assumptions apart from the physical constraint that the total stellar mass within any aperture must be less than the total mass within the aperture, we find that the bulges of the galaxies cannot have IMFs heavier (i.e. implying high mass per unit luminosity) than Salpeter, while the disc IMFs are not well constrained by this data set. We also discuss the necessity for hierarchical modelling when combining incomplete information about multiple astronomical objects. This modelling approach allows us to place upper limits on the size of any departures from universality. More data, including spatially resolved kinematics (as in paper V) and stellar population diagnostics over a range of bulge and disc masses, are needed to robustly quantify how the IMF varies within galaxies.

We study the gas and stellar mass content of galaxy groups and clusters in the fable suite of cosmological hydrodynamical simulations, including the evolution of their central brightest cluster galaxies (BCGs), satellite galaxies, and intracluster light (ICL). The total gas and stellar mass of fable clusters are in good agreement with observations and show negligible redshift evolution at fixed halo mass for $M_{500} \gtrsim 3 \times 10^{14} \, \mathrm{M}_{\odot }$ at z ≲ 1, in line with recent findings from Sunyaev–Zel’dovich (SZ)-selected cluster samples. Importantly, the simulations predict significant redshift evolution in these quantities in the low-mass ($M_{500} \sim 10^{14} \, \mathrm{M}_{\odot }$) regime, which will be testable with upcoming SZ surveys such as SPT-3G. Whilst the stellar masses of fable BCGs are in reasonable agreement with observations, the total stellar mass in satellite galaxies is lower than observed and the total mass in ICL is somewhat higher. This may be caused by enhanced tidal stripping of satellite galaxies due to their large sizes. BCGs are characterized by moderate stellar mass growth at z < 1 coincident with a late-time development of the ICL. The level of BCG mass growth is in good agreement with recent observations; however, we caution that the inferred growth depends sensitively on the mass definition. We further show that in situ star formation contributes more than half the mass of a BCG over its lifetime, the bulk of which is gained at z > 1 where star formation rates are highest. The stellar mass profiles of the BCG+ICL component are similar to observed profiles out to ∼100 kpc at z ≈ 0 and follow a close to power law shape out to several hundred kpc. We further demonstrate that the inferred size growth of BCGs can be severely biased by the choice of parametric model and the outer radius of the fit.

We present the surface brightness profile of M31's stellar halo out to a projected radius of 175 kpc. The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric M31 membership diagnostics is used to identify and reject foreground and background contaminants. This enables us to trace the stellar halo of M31 to larger projected distances and fainter surface brightnesses than previous photometric studies. The surface brightness profile of M31's halo follows a power-law with index -2.2 +/- 0.2 and extends to a projected distance of at least ~175 kpc (~ 2/3 of M31's virial radius), with no evidence of a downward break at large radii. The best-fit elliptical isophotes have b/a=0.94 with the major axis of the halo aligned along the minor axis of M31's disk, consistent with a prolate halo, although the data are also consistent with M31's halo having spherical symmetry. The fact that tidal debris features are kinematically cold is used to identify substructure in the spectroscopic fields out to projected radii of 90 kpc, and investigate the effect of this substructure on the surface brightness profile. The scatter in the surface brightness profile is reduced when kinematically identified tidal debris features in M31 are statistically subtracted; the remaining profile indicates a comparatively diffuse stellar component to M31's stellar halo exists to large distances. Beyond 90 kpc, kinematically cold tidal debris features can not be identified due to small number statistics; nevertheless, the significant field-to-field variation in surface brightness beyond 90 kpc suggests that the outermost region of M31's halo is also comprised to a significant degree of stars stripped from accreted objects.

We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.

We present a first measurement of the stellar mass component of galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect, using 3.6 um and 4.5 um photometry from the Spitzer Space Telescope. Our sample consists of 14 clusters detected by the Atacama Cosmology Telescope (ACT), which span the redshift range 0.27 < z < 1.07 (median z = 0.50), and have dynamical mass measurements, accurate to about 30 per cent, with median M500 = 6.9 x 10^{14} MSun. We measure the 3.6 um and 4.5 um galaxy luminosity functions, finding the characteristic magnitude (m*) and faint-end slope (alpha) to be similar to those for IR-selected cluster samples. We perform the first measurements of the scaling of SZ-observables (Y500 and y0) with both brightest cluster galaxy (BCG) stellar mass and total cluster stellar mass (M500star). We find a significant correlation between BCG stellar mass and Y500 (E(z)^{-2/3} DA^2 Y500 ~ M*^{1.2 +/- 0.6}), although we are not able to obtain a strong constraint on the slope of the relation due to the small sample size. Additionally, we obtain E(z)^{-2/3} DA^2 Y500 ~ M500star^{1.0 +/- 0.6} for the scaling with total stellar mass. The mass fraction in stars spans the range 0.006-0.034, with the second ranked cluster in terms of dynamical mass (ACT-CL J0237-4939) having an unusually low total stellar mass and the lowest stellar mass fraction. For the five clusters with gas mass measurements available in the literature, we see no evidence for a shortfall of baryons relative to the cosmic mean value.

We present a study on the stellar mass growth of the progenitors of local massive galaxies with a variety of number density selections with $n\le1\times10^{-4}\,\rm{Mpc^{-3}}$ (corresponding to $M_*=10^{11.24}\rm{M_{\odot}}$ at z=0.3) in the redshift range $0.3<z<3.0$. We select the progenitors of massive galaxies using a constant number density selection, and one which is adjusted to account for major mergers. We find that the progenitors of massive galaxies grow by a factor of four in total stellar mass over this redshift range. On average the stellar mass added via the processes of star formation, major and minor mergers account for $24\pm8\%$, $17\pm15\%$ and $34\pm14\%$, respectively, of the total galaxy stellar mass at $z=0.3$. Therefore $51\pm20\%$ of the total stellar mass in massive galaxies at $z=0.3$ is created externally to their z=3 progenitors. We explore the implication of these results on the cold gas accretion rate and size evolution of the progenitors of most massive galaxies over the same redshift range. We find an average gas accretion rate of $\sim66\pm32\,\rm{M_{\odot}yr^{-1}}$ over the redshift range of $1.5<z<3.0$. We find that the size evolution of a galaxy sample selected this way is on average lower than the findings of other investigations.

We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars beyond 100 kpc. While a wide range of metallicities is seen throughout the halo, the metal-rich peak of the metallicity distribution function becomes significantly less prominent with increasing radius. The metallicity profile of M31's stellar halo shows a continuous gradient from 9 to ~100 kpc, with a magnitude of -0.01 dex/kpc. The stellar velocity distributions in each field are used to identify stars that are likely associated with tidal debris features. The removal of tidal debris features does not significantly alter the metallicity gradient in M31's halo: a gradient is maintained in fields spanning 10 to 90 kpc. We analyze the halo metallicity profile, as well as the relative metallicities of stars associated with tidal debris features and the underlying halo population, in the context of current simulations of stellar halo formation. We argue that the large scale gradient in M31's halo implies M31 accreted at least one relatively massive progenitor in the past, while the field to field variation seen in the metallicity profile indicates that multiple smaller progenitors are likely to have contributed substantially to M31's outer halo.

Gaseous inflows are necessary suppliers of galaxies’ star-forming fuel, but are difficult to characterize at the survey scale. We use integral-field spectroscopic measurements of gas-phase metallicity and single-dish radio measurements of total atomic gas mass to estimate the magnitude and frequency of gaseous inflows incident on star-forming galaxies. We reveal a mutual correlation between steep oxygen abundance profiles between 0.25 and 1.5 R e , increased variability of metallicity between 1.25 and 1.75 R e , and elevated H i content at fixed total galaxy stellar mass. Employing a simple but intuitive inflow model, we find that galaxies with total stellar mass less than 1010.1 M ⊙ have local oxygen abundance profiles consistent with reinvigoration by inflows. Approximately 10%–25% of low-mass galaxies possess signatures of recent accretion, with estimated typical enhancements of approximately 10%–90% in local gas mass surface density. Higher-mass galaxies have limited evidence for such inflows. The large diversity of H i mass implies that inflow-associated gas ought to reside far from the star-forming disk. We therefore propose that a combination of high H i mass, steep metallicity profile between 0.25 and 1.5 R e , and wide metallicity distribution function between 1.25 and 1.75 R e be employed to target possible hosts of inflowing gas for high-resolution radio follow-up.

Correlations between the mass of a supermassive black hole (SMBH) and the properties of its host galaxy (e.g., total stellar mass M *, luminosity L host) suggest an evolutionary connection. A powerful test of a coevolution scenario is to measure the relations  BH –L host and  BH –M * at high redshift and compare with local estimates. For this purpose, we acquired Hubble Space Telescope (HST) imaging with WFC3 of 32 X-ray-selected broad-line (type 1) active galactic nuclei at 1.2 &lt; z &lt; 1.7 in deep survey fields. By applying state-of-the-art tools to decompose the HST images including available ACS data, we measured the host galaxy luminosity and stellar mass along with other properties through the two-dimensional model fitting. The black hole mass (  BH ) was determined using the broad Hα line, detected in the near-infrared with the Subaru Fiber Multi-Object Spectrograph, which potentially minimizes systematic effects using other indicators. We find that the observed ratio of  BH to total M * is 2.7× larger at z ∼ 1.5 than in the local universe, while the scatter is equivalent between the two epochs. A nonevolving mass ratio is consistent with the data at the 2σ–3σ confidence level when accounting for selection effects (estimated using two independent and complementary methods) and their uncertainties. The relationship between  BH and host galaxy total luminosity paints a similar picture. Therefore, our results cannot distinguish whether SMBHs and their total host stellar mass and luminosity proceed in lockstep or whether the growth of the former somewhat overshoots the latter, given the uncertainties. Based on a statistical estimate of the bulge-to-total mass fraction, the ratio  BH /M *,bulge is offset from the local value by a factor of ∼7, which is significant even accounting for selection effects. Taken together, these observations are consistent with a scenario in which stellar mass is subsequently transferred from an angular momentum–supported component of the galaxy to a pressure-supported one through secular processes or minor mergers at a faster rate than mass accretion onto the SMBH.

We present, for the first time, the relationship between local stellar mass surface density, Σ*, and N/O derived from SDSS-IV MaNGA data, using a sample of 792,765 high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modeling and partial correlation analysis, we find that Σ* alone is insufficient to predict the N/O in MaNGA spaxels and that there is an additional dependence on the local star formation rate surface density, ΣSFR. This effect is a factor of 3 stronger than the dependence of 12+log(O/H) on ΣSFR. Surprisingly, we find that the local N/O scaling relations also depend on the total galaxy stellar mass at fixed Σ* and the galaxy size at fixed stellar mass. We find that more compact galaxies are more nitrogen rich, even when Σ* and ΣSFR are controlled for. We show that ∼50% of the variance of N/O is explained by the total stellar mass and size. Thus, the evolution of nitrogen in galaxies is set by more than just local effects and does not simply track the buildup of oxygen in galaxies. The precise form of the N/O–O/H relation is therefore sensitive to the sample of galaxies from which it is derived. This result casts doubt on the universal applicability of nitrogen-based strong-line metallicity indicators derived in the local universe.

To investigate galaxy properties as a function of their total stellar mass, we obtained 21cm HI line observations at the 100-m class Nan\c{c}ay Radio Telescope of 2839 galaxies from the Sloan Digital Sky Survey (SDSS) in the Local Volume (900<cz<12,000 km/s), dubbed the Nan\c{c}ay Interstellar Baryons Legacy Extragalactic Survey (NIBLES) sample. They were selected evenly over their entire range of absolute SDSS z-band magnitudes (-13.5 to -24 mag), which were used as a proxy for their stellar masses. Here, a first, global presentation of the observations and basic results is given, their further analysis will be presented in other papers in this series. The galaxies were selected based on their properties, as listed in SDSS DR5. Comparing this photometry to their total HI masses, we noted that, for a few percent, the SDSS magnitudes appeared severely misunderestimated, as confirmed by our re-measurements for selected objects. Although using the later DR9 results eliminated this problem in most cases, 384 still required manual photometric source selection. Usable HI spectra were obtained for 2600 galaxies, of which 1733 (67%) were clearly detected and 174 (7%) marginally. The spectra for 241 other observed galaxies could not be used for further analysis owing to problems with either the HI or the SDSS data. We reached the target number of about 150 sources per half-magnitude bin over the Mz range -16.5 to -23 mag. Down to -21 mag the overall detection rate is rather constant at the ~75% level but it starts to decline steadily towards the 30% level at -23 mag. Making regression fits by comparing total HI and stellar masses for our sample, including our conservatively estimated HI upper limits for non-detections, we find the relationship log(M_HI/M*) = -0.59 log(M*) + 5.05, which lies significantly below the relationship found in the M_HI/M* - M* plane when only using HI detections.

We analysed the optical spectra of HII regions extracted from a sample of 350 galaxies of the CALIFA survey. We calculated total O/H abundances and, for the first time, N/O ratios using the semi-empirical routine HII-CHI-mistry, which, according to P\'erez-Montero (2014), is consistent with the direct method and reduces the uncertainty in the O/H derivation using [NII] lines owing to the dispersion in the O/H-N/O relation. Then we performed linear fittings to the abundances as a function of the de-projected galactocentric distances. The analysis of the radial distribution both for O/H and N/O in the non-interacting galaxies reveals that both average slopes are negative, but a non-negligible fraction of objects have a flat or even a positive gradient (at least 10\% for O/H and 4\% for N/O). The slopes normalised to the effective radius appear to have a slight dependence on the total stellar mass and the morphological type, as late low-mass objects tend to have flatter slopes. No clear relation is found, however, to explain the presence of inverted gradients in this sample, and there is no dependence between the average slopes and the presence of a bar. The relation between the resulting O/H and N/O linear fittings at the effective radius is much tighter (correlation coefficient $\rho_s$ = 0.80) than between O/H and N/O slopes ($\rho_s$ = 0.39) or for O/H and N/O in the individual \hii\ regions ($\rho_s$ = 0.37). These O/H and N/O values at the effective radius also correlate very tightly (less than 0.03 dex of dispersion) with total luminosity and stellar mass. The relation with other integrated properties, such as star formation rate, colour, or morphology, can be understood only in light of the found relation with mass.

We study galaxy pair samples selected from the Sloan Digital Sky Survey (SDSS-DR7) and we perform an analysis of minor and major mergers with the aim of investigating the dependence of galaxy properties on interactions. We build a galaxy pair catalog requiring rp < 25 kpc h-1 and Delta V < 350 km s-1 within redshift z<0.1. By visual inspection of SDSS images we removed false identifications and we classify the interactions into three categories: pairs undergoing merging, M; pairs with evident tidal features, T; and non disturbed, N. We also divide the pair sample into minor and major interactions according to the luminosity ratio of the galaxy members. We study star formation activity through colors and star formation rates. We find that 10% of the pairs are classified as M. These systems show an excess of young stellar populations as inferred from the Dn(4000) spectral index, colors, and star formation rates of the member galaxies, an effect which we argue, is directly related to the ongoing merging process. We find 30% of the pairs exhibiting tidal features (T pairs) with member galaxies showing evidence of old stellar populations. Regardless of the color distribution, we find a prominent blue peak in the strongest mergers, while pairs with tidal signs under a minor merger show a strong red peak. Therefore, our results show that galaxy interactions are important in driving the evolution of galaxy bimodality. By adding stellar masses and star formation rates of the two members of the pairs, we explore the global efficiency of star formation of the pairs as a whole. We find that, at a given total stellar mass, major mergers are significantly more efficient (a factor 2) in forming new stars, with respect to both minor mergers or a control sample of non-interacting galaxies.

We present the first results of a wide-field mapping survey of the M81 group conducted with Hyper Suprime-Cam on the Subaru Telescope. Our deep photometry reaches ∼2 mag below the tip of the red giant branch (RGB) and reveals the spatial distribution of both old and young stars over an area of kpc at the distance of M81. The young stars (∼30–160 Myr old) closely follow the neutral hydrogen distribution and can be found in a stellar stream between M81 and NGC 3077 and in numerous outlying stellar associations, including the known concentrations of Arp's Loop, Holmberg IX, an arc in the halo of M82, BK3N, and the Garland. Many of these groupings do not have counterparts in the RGB maps, suggesting they may be genuinely young systems. Our survey also reveals for the first time the very extended () halos of RGB stars around M81, M82, and NGC 3077, as well as faint tidal streams that link these systems. The halos of M82 and NGC 3077 exhibit highly disturbed morphologies, presumably a consequence of the recent gravitational encounter and their ongoing disruption. While the halos of M81 and NGC 3077 and the inner halo of M82 have similar colors, the outer halo of M82 is significantly bluer indicating it is more metal poor. Remarkably, our deep panoramic view of the M81 group demonstrates that the complexity long known to be present in HI is equally matched in the low surface brightness stellar component.