Showing 1–50 of 215 results for author: Kauffmann, G

Simulated [CII] in high-z galaxies

Authors: N. Muñoz-Elgueta, F. Arrigoni Battaia, G. Kauffmann, R. Pakmor, S. Walch, A. Obreja, L. Buhlmann

Abstract: Extended [CII] emission on tens of kpc, also known as a [CII] halo, is being currently reported around z$\sim$4-6 star-forming galaxies, especially thanks to the statistics of the ALPINE survey. The [CII] emission is expected to trace dense cold gas in the inner CGM of these galaxies. The origin of this emission is still debated. In this paper, we present a post-processing model applied to Illustr… ▽ More Extended [CII] emission on tens of kpc, also known as a [CII] halo, is being currently reported around z$\sim$4-6 star-forming galaxies, especially thanks to the statistics of the ALPINE survey. The [CII] emission is expected to trace dense cold gas in the inner CGM of these galaxies. The origin of this emission is still debated. In this paper, we present a post-processing model applied to Illustris-TNG50 star-forming galaxies at $z\sim$4-6, and we compare our results with the ALPINE observations. By incorporating C$^{+}$ abundances derived from UV background and young stars as radiation sources, we generate mock observations, from which we extract surface-brightness (SB) profiles. We found that our model predicts similar [CII] emission values on galactic scales as the observations, providing validation for our approach. However, we find that the predicted [CII] emission in the inner circumgalactic medium (CGM) falls below the observed values by a factor of $\sim$10. We discuss several model limitations that may contribute to this discrepancy. We also find discrepancies with observations when comparing SB profiles of low and high-SFR galaxies. Unlike the observations, simulations exhibit no discernible difference in the extended [CII] emission between the two subsamples. This discrepancy may reflect shortcomings in feedback model of the simulation. Finally, our analysis suggests that the extended [CII] emission is likely a result of both gas from satellite galaxies and outflows from central galaxies, with satellites playing a dominant role within 0.6$<$R/R$_{\rm vir}<$1. A firm estimate of the importance of each contribution is beyond the scope of the current simulations. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: Accepted on July 23, A&A. 17 pages, 16 figures

FEASTS Combined with Interferometry (I): Overall Properties of Diffuse HI and Implications for Gas Accretion in Nearby Galaxies

Authors: Jing Wang, Xuchen Lin, Dong Yang, Lister Staveley-Smith, Fabian Walter, Q. Daniel Wang, Ran Wang, A. J. Battisti, Barbara Catinella, Hsiao-Wen Chen, Luca Cortese, D. B. Fisher, Luis C. Ho, Suoqing Ji, Peng Jiang, Guinevere Kauffmann, Xu Kong, Ziming Liu, Li Shao, Jie Wang, Lile Wang, Shun Wang

Abstract: We present a statistical study of the properties of diffuse HI in ten nearby galaxies, comparing the HI detected by the single-dish telescope FAST (FEASTS program) and the interferometer VLA (THINGS program), respectively. The THINGS' observation missed HI with a median of 23% due to the short-spacing problem of interferometry and limited sensitivity. We extract the diffuse HI by subtracting the d… ▽ More We present a statistical study of the properties of diffuse HI in ten nearby galaxies, comparing the HI detected by the single-dish telescope FAST (FEASTS program) and the interferometer VLA (THINGS program), respectively. The THINGS' observation missed HI with a median of 23% due to the short-spacing problem of interferometry and limited sensitivity. We extract the diffuse HI by subtracting the dense HI, which is obtained from the THINGS data with a uniform flux-density threshold, from the total HI detected by FAST. Among the sample, the median diffuse-HI fraction is 34%, and more diffuse HI is found in galaxies exhibiting more prominent tidal-interaction signatures. The diffuse HI we detected seems to be distributed in disk-like layers within a typical thickness of $1\,\text{kpc}$, different from the more halo-like diffuse HI detected around NGC 4631 in a previous study. Most of the diffuse HI is cospatial with the dense HI and has a typical column density of $10^{17.7}$-$10^{20.1}\,\text{cm}^{-2}$. The diffuse and dense HI exhibits a similar rotational motion, but the former lags by a median of 25% in at least the inner disks, and its velocity dispersions are typically twice as high. Based on a simplified estimation of circum-galactic medium properties and assuming pressure equilibrium, the volume density of diffuse HI appears to be constant within each individual galaxy, implying its role as a cooling interface. Comparing with existing models, these results are consistent with a possible link between tidal interactions, the formation of diffuse HI, and gas accretion. △ Less

Submitted 14 April, 2024; originally announced April 2024.

Comments: 45 pages, 23 figures. In press at ApJ. Data will be released at the FEASTS site upon publication

DIISC-III: Signatures of Stellar Disk Growth in Nearby Galaxies

Authors: Mansi Padave, Sanchayeeta Borthakur, Hansung B. Gim, David Thilker, Rolf A. Jansen, Jacqueline Monckiewicz, Robert C. Kennicutt, Guinevere Kauffmann, Andrew J. Fox, Emmanuel Momjian, Timothy Heckman

Abstract: We explore the growth of the stellar disks in 14 nearby spiral galaxies as part of the Deciphering the Interplay between the Interstellar medium, Stars, and the Circumgalactic medium (DIISC) survey. We study the radial distribution of specific star formation rates (sSFR) and investigate the ratio of the difference in the outer and inner sSFR ($Δ_{sSFR}~={\rm sSFR}_{out}-{\rm sSFR}_{in}$) of the di… ▽ More We explore the growth of the stellar disks in 14 nearby spiral galaxies as part of the Deciphering the Interplay between the Interstellar medium, Stars, and the Circumgalactic medium (DIISC) survey. We study the radial distribution of specific star formation rates (sSFR) and investigate the ratio of the difference in the outer and inner sSFR ($Δ_{sSFR}~={\rm sSFR}_{out}-{\rm sSFR}_{in}$) of the disk and the total sSFR, $Δ_{sSFR}$/sSFR to quantify disk growth. We find $Δ_{sSFR}$/sSFR and the HI gas fraction to show a mild correlation of Spearman's $ρ=0.30$, indicating that star formation and disk growth are likely to proceed outward in galactic disks with high HI gas fractions. The HI gas fractions and $Δ_{sSFR}$/sSFR of the galaxies also increase with the distance to the nearest L$_\star$ neighbor, suggesting that galaxies are likely to sustain their ISM cold gas and exhibit inside-out growth in isolated environments. However, the HI content in their circumgalactic medium, probed by the Ly$α$ equivalent width (W$_{Lyα}$) excess, is observed to be suppressed in isolated environments, apparent from the strong anti-correlation between the W$_{Lyα}$ excess and the distance to the 5$^{\rm th}$ nearest L$_\star$ neighbor (Spearman's $ρ=-0.62$). As expected, W$_{Lyα}$ is also found to be suppressed in cluster galaxies. We find no relation between the W$_{Lyα}$ excess of the detected CGM absorber and $Δ_{sSFR}$/sSFR implying that the enhancement and suppression of the circumgalactic HI gas does not affect the direction in which star formation proceeds in a galactic disk or vice-versa. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ

A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN-starburst connection?

Authors: Guinevere Kauffmann, Iker Millán-Irigoyen, Paul Crowther, Claudia Maraston

Abstract: We have selected galaxies with very high levels of H$α$ emission (EQW(H$α$)$>$700 Å.) in their central regions from the final data release of the MaNGA survey . Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range $9.5 < \log M_*/(M_{\odot}) < 11.5$. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populati… ▽ More We have selected galaxies with very high levels of H$α$ emission (EQW(H$α$)$>$700 Å.) in their central regions from the final data release of the MaNGA survey . Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range $9.5 < \log M_*/(M_{\odot}) < 11.5$. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2-dimensional space of extinction-corrected H$α$ EQW and [SIII]/[SII] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100Å$<$EQW(H$α$)$<$200Å, we obtain the following main results: 1) Clear signatures of excess Wolf Rayet stars are found in half of the H$α$ excess regions, 2) Galaxy regions containing excess Wolf Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [NeIII] is detected in 4 out of 8 of the WR regions and there are tentative [FeX] detections in 2 galaxies. 3) Regions of the galaxy with excess Wolf Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high wavelength-resolution theoretical atmosphere libraries for Wolf Rayet stars. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: accepted for publication in MNRAS

Circumgalactic Ly$α$ emission around submillimeter-bright galaxies with different quasar contributions

Authors: Vale González Lobos, Fabrizio Arrigoni Battaia, Seok-Jun Chang, Max Gronke, Guinevere Kauffmann, Chian-Chou Chen, Hai Fu, Aura Obreja, Emanuele P. Farina

Abstract: We present VLT/MUSE observations targeting the extended Lyman-$α$ (Ly$α$) emission of five high-redshift ($z\sim$3-4) submillimeter galaxies (SMGs) with increasing quasar (QSO) radiation: two SMGs, two SMGs hosting a QSO, and one SMG hosting a QSO with a SMG companion (QSO+SMG). These sources should be located in dark matter halos of comparable masses (average mass of… ▽ More We present VLT/MUSE observations targeting the extended Lyman-$α$ (Ly$α$) emission of five high-redshift ($z\sim$3-4) submillimeter galaxies (SMGs) with increasing quasar (QSO) radiation: two SMGs, two SMGs hosting a QSO, and one SMG hosting a QSO with a SMG companion (QSO+SMG). These sources should be located in dark matter halos of comparable masses (average mass of $M_{\rm DM}\sim10^{12.2}\,{\rm M}_\odot$). We quantify the luminosity and extent of the Ly$α$ emission, together with its kinematics, and examine four Ly$α$ powering mechanisms: photoionization from QSOs/star formation, shocks by galactic/QSO outflows, gravitational cooling radiation, and Ly$α$ photons resonant scattering. We find a variety of Ly$α$ luminosities and extents, with the QSO+SMG system displaying the most extended and bright nebula, followed by the SMGs hosting a QSO, and finally the undetected circumgalactic medium (CGM) of SMGs. This diversity implies that gravitational cooling is unlikely to be the main powering mechanism. We show that photoionization from the QSO and QSO outflows can contribute to power the emission for average densities $n_{\rm H}>0.5\,$cm$^{-3}$. Moreover, the observed Ly$α$ luminosities scale with the QSO's budget of Ly$α$ photons modulo the dust content in each galaxy, highlighting a possible contribution from resonant scattering of QSO's radiation in powering the nebulae. We find larger Ly$α$ linewidths (FWHM$\gtrsim 1200\,$km$\,$s$^{-1}$) than usually reported around radio-quiet systems, pointing to large-scale outflows. A statistical survey targeting similar high-redshift massive systems with known host properties is needed to confirm our findings. △ Less

Submitted 30 August, 2023; originally announced August 2023.

Comments: 17 pages, 9 figures, accepted for publication in ApJ

Journal ref: A&A 679, A41 (2023)

Formation of star clusters and enrichment by massive stars in simulations of low-metallicity galaxies with a fully sampled initial stellar mass function

Authors: Natalia Lahén, Thorsten Naab, Guinevere Kauffmann, Dorottya Szécsi, Jessica May Hislop, Antti Rantala, Alexandra Kozyreva, Stefanie Walch, Chia-Yu Hu

Abstract: We present new GRIFFIN project hydrodynamical simulations that model the formation of galactic star cluster populations in low-metallicity ($Z=0.00021$) dwarf galaxies, including radiation, supernova and stellar wind feedback of individual massive stars. In the simulations, stars are sampled from the stellar initial mass function (IMF) down to the hydrogen burning limit of $0.08$ M$_\odot$. Mass c… ▽ More We present new GRIFFIN project hydrodynamical simulations that model the formation of galactic star cluster populations in low-metallicity ($Z=0.00021$) dwarf galaxies, including radiation, supernova and stellar wind feedback of individual massive stars. In the simulations, stars are sampled from the stellar initial mass function (IMF) down to the hydrogen burning limit of $0.08$ M$_\odot$. Mass conservation is enforced within a radius of $1$ pc for the formation of massive stars. We find that massive stars are preferentially found in star clusters and follow a correlation set at birth between the highest initial stellar mass and the star cluster mass that differs from pure stochastic IMF sampling. With a fully sampled IMF, star clusters lose mass in the galactic tidal field according to mass-loss rates observed in nearby galaxies. Of the released stellar feedback, $60\%$ of the supernova material and up to $35\%$ of the wind material reside either in the hot interstellar medium (ISM) or in gaseous, metal enriched outflows. While stellar winds (instantaneously) and supernovae (delayed) start enriching the ISM right after the first massive stars form, the formation of supernova-enriched stars and star clusters is significantly delayed (by $>50$ Myr) compared to the formation of stars and star clusters enriched by stellar winds. Overall, supernova ejecta dominate the enrichment by mass, while the number of enriched stars is determined by continuous stellar winds. These results present a concept for the formation of chemically distinct populations of stars in bound star clusters, reminiscent of multiple populations in globular clusters. △ Less

Submitted 20 April, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 26 pages, 23 figures. Accepted for publication in MNRAS

Flows around galaxies. I. The dependency of galaxy connectivity on cosmic environments and effects on the star-formation rate

Authors: Daniela Galárraga-Espinosa, Enrico Garaldi, Guinevere Kauffmann

Abstract: With the aim of bringing substantial insight to the fundamental question of how galaxies acquire their material for star-formation, we present the first comprehensive characterisation of the galaxy connectivity (i.e. the number of small-scale filamentary streams connected to a galaxy) in relation with the cosmic environment, and a statistical exploration of the impact of connectivity on the star-f… ▽ More With the aim of bringing substantial insight to the fundamental question of how galaxies acquire their material for star-formation, we present the first comprehensive characterisation of the galaxy connectivity (i.e. the number of small-scale filamentary streams connected to a galaxy) in relation with the cosmic environment, and a statistical exploration of the impact of connectivity on the star-formation rate at z=2. We detect kpc-scale filaments directly connected to galaxies by applying the DisPerSE filament finder to the DM density around 2942 central galaxies ($M_* > 10^{8}$ $\mathrm{M}_\odot / h$) of the TNG50-1 simulation. Our results demonstrate that galaxy connectivity spans a broad range (from 0 to 9), with more than half of the galaxies connected to two or three streams. We examine a variety of factors that could influence the connectivity finding out that it increases with mass, decreases with local density for low mass galaxies, and does not depend on local environment, estimated by the Delaunay tessellation, for high mass galaxies. We further classify galaxies according to their location in different cosmic web environments, and we highlight the influence of the large-scale structure on the number of connected streams. Our results reflect the different strengths of the cosmic tides, which can prevent the formation of coherent streams feeding the galaxies, or even disconnect the galaxy from its local web. Finally, we show that, at fixed local density, the star-formation rate (SFR) of low mass galaxies is up to $5.9σ$ enhanced due to connectivity. This SFR boost is even more significant ($6.3σ$) for galaxies embedded in cosmic filaments, where the available matter reservoirs are large. A milder impact is found for high mass galaxies, hinting at different relative efficiencies of matter inflow via small-scale streams in galaxies of different masses. △ Less

Submitted 27 April, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: published in A&A

Journal ref: A&A 671, A160 (2023)

The Physical Origin of Galactic Conformity: From Theory to Observation

Authors: Mohammadreza Ayromlou, Guinevere Kauffmann, Abhijeet Anand, Simon D. M. White

Abstract: We employ several galaxy formation models, in particular, L-GALAXIES, IllustrisTNG, and EAGLE, as well as observational samples from SDSS and DESI, to investigate galactic conformity, the observed large-scale correlation between the star-formation properties of central (primary) galaxies and those of their neighbours. To analyse the models and observations uniformly, we introduce CENSAT, a new alg… ▽ More We employ several galaxy formation models, in particular, L-GALAXIES, IllustrisTNG, and EAGLE, as well as observational samples from SDSS and DESI, to investigate galactic conformity, the observed large-scale correlation between the star-formation properties of central (primary) galaxies and those of their neighbours. To analyse the models and observations uniformly, we introduce CENSAT, a new algorithm to define whether a galaxy is a central or a satellite system based on an isolation criterion. We find that the conformity signal is present, up to at least 5 Mpc from the centres of low- and intermediate-mass centrals in the latest version of L-GALAXIES (Ayromlou et al. 2021), IllustrisTNG, and EAGLE, as well as in SDSS and DESI observational samples. In comparison, the conformity signal is substantially weaker in an older version of L-GALAXIES (Henriques et al. 2020). One of the main differences between this older model and the other models is its neglect of ram-pressure stripping of the gas reservoirs of galaxies except within the boundaries of massive cluster haloes. Our observational comparisons demonstrate that this difference significantly affects the observed large-scale conformity signal. Furthermore, by examining the contribution of backsplash, fly-by, central, and satellite galaxies to the conformity signal, we show that much, but not all, of it arises from primary galaxies near massive systems. Remaining tensions between the models and observations may be solved by modifying the physical prescriptions for how feedback processes affect the distribution and kinematics of gas and the environment around galaxies out to scales of several Megaparsecs. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: Submitted to MNRAS

A study of 1000 galaxies with unusually young and massive stars in the SDSS: a search for hidden black holes

Authors: Guinevere Kauffmann, Claudia Maraston, Johan Comparat, Paul Crowther

Abstract: We select 1076 galaxies with extinction-corrected Halpha equivalent widths too large to be explained with a Kroupa (2001) IMF, and compare these with a control sample of galaxies that is matched in stellar mass, redshift and 4000 AA break strength, but with normal Halpha equivalent widths. Our goal is to study how processes such as black hole growth and energetic feedback processes from massive st… ▽ More We select 1076 galaxies with extinction-corrected Halpha equivalent widths too large to be explained with a Kroupa (2001) IMF, and compare these with a control sample of galaxies that is matched in stellar mass, redshift and 4000 AA break strength, but with normal Halpha equivalent widths. Our goal is to study how processes such as black hole growth and energetic feedback processes from massive stars differ between galaxies with extreme central Halpha emission and galaxies with normal young central stellar populations. The stellar mass distribution of Halpha excess galaxies is peaked at 3 \times 10^10 Msun and almost all fall well within the star-forming locus in the [OIII]/Hbeta versus [NII]/Halpha BPT disgram. Halpha excess galaxies are twice as likely to exhibit Halpha line asymmetries and 1.55 times more likely to be detected at 1 GHz in the VLA FIRST survey compared to control sample galaxies. The radio luminosity per unit stellar mass decreases with the stellar age of the system. Using stacked spectra, we demonstrate that [NeV] emission is not present in the very youngest of the radio-quiet Halpha excess galaxies with detectable Wolf-Rayet features, suggesting that black hole growth has not yet commenced in such systems. [NeV] emission is detected in Halpha excess galaxies with radio detections and the strength of the line correlates with the radio luminosity. This is the clearest indication for a possible population of black holes that may be forming in a subset of the Halpha excess population. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: 17 pages, 22 figures, accepted in MNRAS

SDSS-IV MaNGA: Exploring the local scaling relations for N/O

Authors: Adam L. Schaefer, Christy Tremonti, Guinevere Kauffmann, Brett H. Andrews, Matthew A. Bershady, Nicholas F. Boardman, Kevin Bundy, Niv Drory, José G. Fernández-Trincado, Holly P. Preece, Rogério Riffel, Rogemar A. Riffel, Sebastián F. Sánchez

Abstract: We present, for the first time, the relationship between local stellar mass surface density, $\mathrm{Σ_{*}}$, and N/O derived from SDSS-IV MaNGA data, using a sample of $792765$ high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modelling and partial correlation analysis, we find that $\mathrm{Σ_{*}}$ alone is insufficient to predict the N/O in MaNGA spaxels,… ▽ More We present, for the first time, the relationship between local stellar mass surface density, $\mathrm{Σ_{*}}$, and N/O derived from SDSS-IV MaNGA data, using a sample of $792765$ high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modelling and partial correlation analysis, we find that $\mathrm{Σ_{*}}$ 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, $\mathrm{Σ_{SFR}}$. This effect is a factor of $3$ stronger than the dependence of 12+log(O/H) on $\mathrm{Σ_{SFR}}$. Surprisingly, we find that the local N/O scaling relations also depend on the total galaxy stellar mass at fixed $Σ_{*}$ as well as the galaxy size at fixed stellar mass. We find that more compact galaxies are more nitrogen rich, even when $\mathrm{Σ_{*}}$ and $\mathrm{Σ_{SFR}}$ are controlled for. We show that $\sim50\%$ 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 build up 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. △ Less

Submitted 31 March, 2022; originally announced March 2022.

Comments: 22 pages, 14 figures. Accepted for publication in the Astrophysical Journal

APEX at the QSO MUSEUM: molecular gas reservoirs associated with $z\sim$3 quasars and their link to the extended Ly$α$ emission

Authors: N. Muñoz-Elgueta, F. Arrigoni Battaia, G. Kauffmann, C. De Breuck, C. García-Vergara, A. Zanella, E. P. Farina, R. Decarli

Abstract: Cool gas (T$\sim$10$^{4}$~K) traced by hydrogen Ly$α$ emission is now routinely detected around $z\sim3$ quasars, but little is known about their molecular gas reservoirs. Here, we present an APEX spectroscopic survey of the CO(6-5), CO(7-6) and [CI](2-1) emission lines for 9 quasars from the QSO MUSEUM survey which have similar UV luminosities, but very diverse Ly$α$ nebulae. These observations (… ▽ More Cool gas (T$\sim$10$^{4}$~K) traced by hydrogen Ly$α$ emission is now routinely detected around $z\sim3$ quasars, but little is known about their molecular gas reservoirs. Here, we present an APEX spectroscopic survey of the CO(6-5), CO(7-6) and [CI](2-1) emission lines for 9 quasars from the QSO MUSEUM survey which have similar UV luminosities, but very diverse Ly$α$ nebulae. These observations ($\langle~\rm rms~\rangle=2.6$~mJy in 300~km~s$^{-1}$) detected three CO(6-5) lines with 3.4$\leq I_{\rm CO(6-5)} \leq$5.1~Jy~km~s$^{-1}$, 620$\leq$FWHM$\leq$707~km~s$^{-1}$, and three [CI](2-1) lines with 2.3$\leq I_{\rm [CI](2-1)} \leq$15.7~Jy~km~s$^{-1}$, 329$\leq$FWHM$\leq$943~km~s$^{-1}$. For the CO and [CI] detected sources, we constrain the molecular gas reservoirs to be $\rm M_{H_{2}} = (0.4-6.9) \times 10^{11} M_{\odot}$, while the non-detections imply $\rm M_{H_{2}} < 1.1\times 10^{11} M_{\odot}$. We compare our observations with the extended Ly$α$ properties to understand the link between the cool and the molecular gas phases. We find large velocity shifts between the bulk of Ly$α$ and the molecular gas systemic redshift in five sources (from $\sim$-400 to $\sim+$1200~km~s$^{-1}$). The sources with the largest shifts have the largest Ly$α$ line widths in the sample, suggesting more turbulent gas conditions and/or large-scale inflows/outflows around these quasars. We also find that the brightest ($I_{\rm[CI](2-1)}=15.7\pm3.7~\rm Jy~km~s^{-1}$) and the widest (FWHM$\sim$900~km~s$^{-1}$) lines are detected for the smallest and dimmest Ly$α$ nebulae. From this, we speculate that host galaxy obscuration can play an important role in reducing the ionizing and Ly$α$ photons able to escape to halo scales, and/or that these systems are hosted by more massive halos. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: Accepted on January 3, MNRAS. 25 pages, 17 figures

Cool circumgalactic gas in galaxy clusters: connecting the DESI legacy imaging survey and SDSS DR16 MgII absorbers

Authors: Abhijeet Anand, Guinevere Kauffmann, Dylan Nelson

Abstract: We investigate the cool gas absorption in galaxy clusters by cross-correlating MgII absorbers detected in quasar spectra from Data Release 16 of the Sloan Digital Sky Survey (SDSS) with galaxy clusters identified in the Dark Energy Spectroscopic Instrument (DESI) survey. We find significant covering fractions ($1-5\, \%$ within $r_{500}$, depending on the chosen redshift interval), $\sim 4-5$ time… ▽ More We investigate the cool gas absorption in galaxy clusters by cross-correlating MgII absorbers detected in quasar spectra from Data Release 16 of the Sloan Digital Sky Survey (SDSS) with galaxy clusters identified in the Dark Energy Spectroscopic Instrument (DESI) survey. We find significant covering fractions ($1-5\, \%$ within $r_{500}$, depending on the chosen redshift interval), $\sim 4-5$ times higher than around random sightlines. While the covering fraction of cool gas in clusters decreases with increasing mass of the central galaxy, the total MgII mass within $r_{\rm 500}$ is nonetheless $\sim 10$ times higher than for SDSS luminous red galaxies (LRGs). The MgII covering fraction versus impact parameter is well described by a power law in the inner regions and a exponential function at larger distances. The characteristic scale of the transition between these two regimes is smaller for large equivalent width absorbers. Cross-correlating MgII absorption with photo-z selected cluster member galaxies from DESI reveals a statistically significant connection. The median projected distance between MgII absorbers and the nearest cluster member is $\sim200$ kpc, compared to $\sim500$ kpc in random mocks with the same galaxy density profiles. We do not find a correlation between MgII strength and the star formation rate of the closest cluster neighbour. This suggests that cool gas in clusters, as traced by MgII absorption, is: (i) associated with satellite galaxies, (ii) dominated by cold gas clouds in the intracluster medium, rather than by the interstellar medium of galaxies, and (iii) may originate in part from gas stripped from these cluster satellites in the past. △ Less

Submitted 4 April, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Comments: 19 pages, revised version accepted for publication in MNRAS

A panchromatic view of star cluster formation in a simulated dwarf galaxy starburst

Authors: Natalia Lahén, Thorsten Naab, Guinevere Kauffmann

Abstract: We present a photometric analysis of star and star cluster (SC) formation in a high-resolution simulation of a dwarf galaxy starburst that allows the formation of individual stars to be followed. Previous work demonstrated that the properties of the SCs formed in the simulation are in good agreement with observations. In this paper, we create mock spectral energy distributions and broad-band photo… ▽ More We present a photometric analysis of star and star cluster (SC) formation in a high-resolution simulation of a dwarf galaxy starburst that allows the formation of individual stars to be followed. Previous work demonstrated that the properties of the SCs formed in the simulation are in good agreement with observations. In this paper, we create mock spectral energy distributions and broad-band photometric images using the radiative transfer code SKIRT 9. We test several observational star formation rate (SFR) tracers and find that $24$ $μ$m, total infrared and H$α$ trace the underlying SFR during the (post)starburst phase, while UV tracers yield a more accurate picture of star formation during quiescent phases prior to and after the merger. We then place the simulated galaxy at distances of $10$ and $50$ Mpc and use aperture photometry at Hubble Space Telescope resolution to analyse the simulated SC population. During the starburst phase, a hierarchically forming set of SCs leads inaccurate source separation because of crowding. This results in estimated SC mass function slopes that are up to $\sim0.3$ shallower than the true slope of $\sim-1.9$ to $-2$ found for the bound clusters identified from the particle data in the simulation. The masses of the largest clusters are overestimated by a factor of up to $2.9$ due to unresolved clusters within the apertures. The aperture-based analysis also produces a relation between cluster formation efficiency and SFR surface density that is slightly flatter than that recovered from bound clusters. The differences are strongest in quiescent SF environments. △ Less

Submitted 28 June, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 21 pages, accepted for publication in MNRAS

DIISC-I: The Discovery of Kinematically Anomalous HI Clouds in M 100

Authors: Hansung B. Gim, Sanchayeeta Borthakur, Emmanuel Momjian, Mansi Padave, Rolf A. Jansen, Dylan Nelson, Timothy M. Heckman, Robert C. Kennicutt Jr., Andrew J. Fox, Jorge L. Pineda, David Thilker, Guinevere Kauffmann, Jason Tumlinson

Abstract: We report the discovery of two kinematically anomalous atomic hydrogen (HI) clouds in M 100 (NGC 4321), which was observed as part of the Deciphering the Interplay between the Interstellar medium, Stars, and the Circumgalactic medium (DIISC) survey in HI 21 cm at 3.3 km s$^{-1}$ spectroscopic and 44 arcsec$\times$30 arcsec spatial resolution using the Karl G. Jansky Very Large Array. These clouds… ▽ More We report the discovery of two kinematically anomalous atomic hydrogen (HI) clouds in M 100 (NGC 4321), which was observed as part of the Deciphering the Interplay between the Interstellar medium, Stars, and the Circumgalactic medium (DIISC) survey in HI 21 cm at 3.3 km s$^{-1}$ spectroscopic and 44 arcsec$\times$30 arcsec spatial resolution using the Karl G. Jansky Very Large Array. These clouds were identified as structures that show significant kinematic offsets from the rotating disk of M100. The velocity offsets of 40 km s$^{-1}$ observed in these clouds are comparable to the offsets seen in intermediate-velocity clouds (IVCs) in the circumgalactic medium (CGM) of the Milky Way and nearby galaxies. We find that one anomalous cloud in M 100 is associated with star-forming regions detected in H$α$ and far-ultraviolet imaging. Our investigation shows that anomalous clouds in M 100 may originate from multiple mechanisms, such as star formation feedback-driven outflows, ram-pressure stripping, and tidal interactions with satellite galaxies. Moreover, we do not detect any cool CGM at 38.8 kpc from the center of M 100, giving an upper limit of N(HI) $\le$ $1.7\times10^{13}$ cm$^{-2}$ (3$σ$). Since M 100 is in the Virgo cluster, the non-existence of neutral/cool CGM is a likely pathway for turning it into a red galaxy. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 19 pages, 3 tables, 7 figures, Accepted for publication at ApJ

A study of the central stellar populations of galaxies in SDSS-IV MaNGA: identification of a sub-sample with unusually young and massive stars

Authors: Guinevere Kauffmann

Abstract: This paper describes a search for galaxy centers with clear indications of unusual stellar populations with an initial mass function flatter than Salpeter at high stellar masses. Out of a sample of 668 face-on galaxies with stellar masses in the range 10^10- 10^11 M_sol, I identify 15 galaxies with young to intermediate age central stellar populations with unusual stellar population gradients in t… ▽ More This paper describes a search for galaxy centers with clear indications of unusual stellar populations with an initial mass function flatter than Salpeter at high stellar masses. Out of a sample of 668 face-on galaxies with stellar masses in the range 10^10- 10^11 M_sol, I identify 15 galaxies with young to intermediate age central stellar populations with unusual stellar population gradients in the inner regions of the galaxy. In these galaxies, the 4000 Angstrom break is either flat or rising towards the center of the galaxy, indicating that the central regions host evolved stars, but the H$α$ equivalent width also rises steeply in the central regions. The ionization parameter [OIII]/[OII] is typically low in these galactic centers, indicating that ionizing sources are stellar rather than AGN. Wolf Rayet features characteristic of hot young stars are often found in the spectra and these also get progressively stronger at smaller galactocentric radii. These outliers are compared to a control sample of galaxies of similar mass with young inner stellar populations, but where the gradients in Halpha equivalent width and 4000 Angstrom break follow each other more closely. The outliers exhibit central Wolf Rayet red bump excesses much more frequently, they have higher central stellar and ionized gas metallicities, and they are also more frequently detected at 20 cm radio wavelengths. I highlight one outlier where the ionized gas is clearly being strongly perturbed and blown out either by massive stars after they explode as supernovae, or by energy injection from matter falling onto a black hole. △ Less

Submitted 11 July, 2021; originally announced July 2021.

Comments: 16 pages, 16 figures, accepted for publication in MNRAS

A study of outer disk stellar populations of face-on star-forming galaxies in SDSS-IV MaNGA: causes of H$α$ deficiency

Authors: Guinevere Kauffmann

Abstract: Integral field unit (IFU) spectra of face-on star-forming galaxies from the MaNGA survey are stacked in radial bins so as to reach a S/N high enough to measure emission lines and Lick indices out to 2.5-3 R_e. Two thirds of galaxies have stellar populations in the outer disks that are older, more metal poor and less dusty than in the inner disks. Recent bursts of star formation have occurred more… ▽ More Integral field unit (IFU) spectra of face-on star-forming galaxies from the MaNGA survey are stacked in radial bins so as to reach a S/N high enough to measure emission lines and Lick indices out to 2.5-3 R_e. Two thirds of galaxies have stellar populations in the outer disks that are older, more metal poor and less dusty than in the inner disks. Recent bursts of star formation have occurred more frequently in the outer disk, but extinction-corrected Halpha equivalent widths are significantly lower at fixed D_n(4000) in these regions. I examine the properties of a subset of galaxies with the the most H$α$ deficient outer disks. These regions contain young stellar populations that must have formed within the last 0.5 Gyr, but extinction-corrected Halpha values well below the values predicted for a standard Kroupa IMF. The Halpha deficient galaxies have flat D_n(4000) and Hdelta_A profiles with little radial fluctuation, indicating that star formation has occurred extremely uniformly across the entire disk. The H$α$ line profiles indicate that the ionized gas kinematics is also very regular across the disk. The main clue to the origin of the Halpha deficiency is that it sets in at the same radius where the dust extinction abruptly decreases, suggesting a mode of star formation deficient in massive stars in quiescent, HI-dominated gas. Finally, I have carried out a search for galaxies with signatures of unusual Halpha kinematics and find that 15% of the sample exhibit evidence for significant ionized gas that is displaced from the systemic velocity of the disk. △ Less

Submitted 11 July, 2021; originally announced July 2021.

Comments: 16 pages, 15 figures, accepted for publication in MNRAS

The parametrization of gas flows in discs in the Auriga simulations

Authors: Periklis Okalidis, Robert J. J. Grand, Robert M. Yates, Guinevere Kauffmann

Abstract: We study the radial motions of cold, star-forming gas in the secular evolution phase of a set of 14 magnetohydrodynamical cosmological zoom-in simulations of Milky Way-mass galaxies. We study the radial transport of material within the disc plane in a series of concentric rings. For the gas in each ring at a given time we compute two quantities as a function of time and radius: 1) the radial bulk… ▽ More We study the radial motions of cold, star-forming gas in the secular evolution phase of a set of 14 magnetohydrodynamical cosmological zoom-in simulations of Milky Way-mass galaxies. We study the radial transport of material within the disc plane in a series of concentric rings. For the gas in each ring at a given time we compute two quantities as a function of time and radius: 1) the radial bulk flow of the gas; and 2) the radial spread of the gas relative to the bulk flow. Averaging the data from all the halos, we find that the radial spread increases with radius in the form of a power law with strong secondary dependencies on the fraction of accreted material and the local radial velocity dispersion of the gas. We find that the bulk motion of gas is well described in the inner disc regions by a radially-independent mean inward flow speed of $-$2.4 km s$^{-1}$. The spread around this value relates to the change in angular momentum of the gas and also the amount of accreted material. These scalings from fully cosmological, MHD simulations of galaxy formation can then be used in semi-analytic models to better parameterise the radial flow of gas in discs. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 17 pages, 13 figures, Accepted for publication in MNRAS

Characterizing the abundance, properties, and kinematics of the cool circumgalactic medium of galaxies in absorption with SDSS DR16

Authors: Abhijeet Anand, Dylan Nelson, Guinevere Kauffmann

Abstract: In order to study the circumgalactic medium (CGM) of galaxies we develop an automated pipeline to estimate the optical continuum of quasars and detect intervening metal absorption line systems with a matched kernel convolution technique and adaptive S/N criteria. We process $\sim$ one million quasars in the latest Data Release 16 (DR16) of the Sloan Digital Sky Survey (SDSS) and compile a large sa… ▽ More In order to study the circumgalactic medium (CGM) of galaxies we develop an automated pipeline to estimate the optical continuum of quasars and detect intervening metal absorption line systems with a matched kernel convolution technique and adaptive S/N criteria. We process $\sim$ one million quasars in the latest Data Release 16 (DR16) of the Sloan Digital Sky Survey (SDSS) and compile a large sample of $\sim$ 160,000 MgII absorbers, together with $\sim$ 70,000 FeII systems, in the redshift range $0.35<z_{abs}<2.3$. Combining these with the SDSS DR16 spectroscopy of $\sim1.1$ million luminous red galaxies (LRGs) and $\sim 200,000$ emission line galaxies (ELGs), we investigate the nature of cold gas absorption at $0.5<z<1$. These large samples allow us to characterize the scale dependence of MgII with greater accuracy than in previous work. We find that there is a strong enhancement of MgII absorption within $\sim 50$ kpc of ELGs, and the covering fraction within $0.5r_{\rm vir}$ of ELGs is 2-5 times higher than for LRGs. Beyond 50 kpc, there is a sharp decline in MgII for both kinds of galaxies, indicating a transition to the regime where the CGM is tightly linked with the dark matter halo. The MgII covering fraction correlates strongly with stellar mass for LRGs, but weakly for ELGs, where covering fractions increase with star formation rate. Our analysis implies that cool circumgalactic gas has a different physical origin for star forming versus quiescent galaxies. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 23 pages, 20 Figures, 3 Appendices, MNRAS accepted. The MgII absorber catalogue is publicly available at https://wwwmpa.mpa-garching.mpg.de/SDSS/MgII/

Galaxy formation with L-GALAXIES: Modelling the environmental dependency of galaxy evolution and comparing with observations

Authors: Mohammadreza Ayromlou, Guinevere Kauffmann, Robert M. Yates, Dylan Nelson, Simon D. M. White

Abstract: We present a variation of the recently updated Munich semi-analytical galaxy formation model, L-Galaxies, with a new gas stripping method. Extending earlier work, we directly measure the local environmental properties of galaxies to formulate a more accurate treatment of ram-pressure stripping for all galaxies. We fully re-calibrate the modified L-Galaxies model using a Markov Chain Monte Carlo (M… ▽ More We present a variation of the recently updated Munich semi-analytical galaxy formation model, L-Galaxies, with a new gas stripping method. Extending earlier work, we directly measure the local environmental properties of galaxies to formulate a more accurate treatment of ram-pressure stripping for all galaxies. We fully re-calibrate the modified L-Galaxies model using a Markov Chain Monte Carlo (MCMC) method with the stellar mass function and quenched fraction of galaxies at $0\leq z\leq2$ as constraints. Due to this re-calibration, global galaxy population relations, including the stellar mass function, quenched fractions versus galaxy mass and HI mass function are all largely unchanged and remain consistent with observations. By comparing to data on galaxy properties in different environments from the SDSS and HSC surveys, we demonstrate that our modified model improves the agreement with the quenched fractions and star formation rates of galaxies as a function of environment, stellar mass, and redshift. Overall, in the vicinity of haloes with total mass $10^{12}$ to $10^{15}\rm M_{\odot}$ at $z=0$, our new model produces higher quenched fractions and stronger environmental dependencies, better recovering observed trends with halocentric distance up to several virial radii. By analysing the actual amount of gas stripped from galaxies in our model, we show that those in the vicinity of massive haloes lose a large fraction of their hot halo gas before they become satellites. We demonstrate that this affects galaxy quenching both within and beyond the halo boundary. This is likely to influence the correlations between galaxies up to tens of megaparsecs. △ Less

Submitted 10 November, 2020; originally announced November 2020.

Comments: Submitted to MNRAS. The full model output and model description are publicly available at https://lgalaxiespublicrelease.github.io

L-GALAXIES 2020: The evolution of radial metallicity profiles and global metallicities in disc galaxies

Authors: Robert M. Yates, Bruno M. B. Henriques, Jian Fu, Guinevere Kauffmann, Peter A. Thomas, Qi Guo, Simon D. M. White, Patricia Schady

Abstract: We present a modified version of the L-GALAXIES 2020 semi-analytic model of galaxy evolution, which includes significantly increased direct metal enrichment of the circumgalactic medium (CGM) by supernovae (SNe). These more metal-rich outflows do not require increased mass-loading factors, in contrast to some other galaxy evolution models. This modified L-GALAXIES 2020 model is able to simultaneou… ▽ More We present a modified version of the L-GALAXIES 2020 semi-analytic model of galaxy evolution, which includes significantly increased direct metal enrichment of the circumgalactic medium (CGM) by supernovae (SNe). These more metal-rich outflows do not require increased mass-loading factors, in contrast to some other galaxy evolution models. This modified L-GALAXIES 2020 model is able to simultaneously reproduce the gas-phase metallicity $(Z_{\rm g})$ and stellar metallicity $(Z_{*})$ radial profiles observed in nearby disc galaxies by MaNGA and MUSE, as well as the observed mass - metallicity relations for gas and stars at $z=0$ and their evolution back to $z\sim{}2-3$. A direct CGM enrichment fraction of $\sim{}90\%$ for SNe-II is preferred. We find that massive disc galaxies have slightly flatter $Z_{\rm g}$ profiles than their lower-mass counterparts in L-GALAXIES 2020, due to more efficient enrichment of their outskirts via inside-out growth and metal-rich accretion. Such a weak, positive correlation between stellar mass and $Z_{\rm g}$ profile slope is also seen in our MaNGA-DR15 sample of 571 star-forming disc galaxies. Although, below ${\rm log}(M_{*}/{\rm M}_{\odot})\sim{}10.0$ this observational result is strongly dependent on the metallicity diagnostic and morphological selection chosen. In addition, a lowered maximum SN-II progenitor mass of $25{\rm M}_{\odot}$, reflecting recent theoretical and observational estimates, can also provide a good match to observed metallicity profiles at $z=0$ in L-GALAXIES 2020. However, this model version fails to reproduce an evolution in $Z_{\rm g}$ at fixed mass over cosmic time, or the magnesium abundances observed in the intracluster medium (ICM). △ Less

Submitted 3 June, 2021; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: 22 pages, 19 figures. Accepted by MNRAS. (Missing citation hyperlinks corrected)

A study of the HI gas fractions of galaxies at z ~ 1

Authors: Wei Zhang, Guinevere Kauffmann, Jing Wang, Yanmei Chen, Jian Fu, Hong Wu

Abstract: Due to the fact that HI mass measurements are not available for large galaxy samples at high redshifts, we apply a photometric estimator of the HI-to-stellar mass ratio (M_HI/M_*) calibrated using a local Universe sample of galaxies to a sample of galaxies at z ~ 1 in the DEEP2 survey. We use these HI mass estimates to calculate HI mass functions (HIMFs) and cosmic HI mass densities (Omega_HI), an… ▽ More Due to the fact that HI mass measurements are not available for large galaxy samples at high redshifts, we apply a photometric estimator of the HI-to-stellar mass ratio (M_HI/M_*) calibrated using a local Universe sample of galaxies to a sample of galaxies at z ~ 1 in the DEEP2 survey. We use these HI mass estimates to calculate HI mass functions (HIMFs) and cosmic HI mass densities (Omega_HI), and to examine the correlation between star formation rate and HI gas content, for galaxies at z ~ 1. We have estimated HI gas masses for ~ 7,000 galaxies in the DEEP2 survey with redshifts in the range 0.75 < z < 1.4 and stellar masses M_* > 10^{10} M_solar, using a combination of the rest-frame ultraviolet-optical colour (NUV - r) and stellar mass density (mu_*) as a way to estimate M_HI/M_*. It is found that the high mass end of high-z HI mass function (HIMF) is quite similar to that of the local HIMF. The lower limit of Omega_HI,limit = 2.1 * 10^{-4} h_70^{-1}, obtained by directly integrating the HI mass of galaxies with M_* > 10^{10} M_solar, confirms that massive star-forming galaxies do not dominate the neutral gas at z ~ 1. We study the evolution of the HI mass to stellar mass ratio from z ~ 1 to today and find a steeper relation between HI gas mass fraction and stellar mass at higher redshifts. Specifically, galaxies with M_* = 10^{11} M_solar at z ~ 1 are found to have 3 - 4 times higher neutral gas fractions than local galaxies, while the increase is as high as 4 - 12 times at M_* = 10^{10} M_solar. The quantity M_HI/SFR exhibits very large scatter, and the scatter increases from a factor of 5 - 7 at z = 0 to factors close to a hundred at z = 1. This implies that there is no relation between HI gas and star formation in high redshift galaxies. The HI gas must be linked to cosmological gas accretion processes at high redshifts. △ Less

Submitted 8 February, 2021; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: 10 pages, 13 figures, A&A accepted

Journal ref: A&A 648, A25 (2021)

Lyman-αabsorption beyond the disk of simulated spiral galaxies

Authors: Bernhard Röttgers, Thorsten Naab, Miha Cernetic, Romeel Davé, Guinevere Kauffmann, Sanchayeeta Borthakur, Horst Foidl

Abstract: We present an analysis of the origin and properties of the circum-galactic medium (CGM) in a suite of 11 cosmological zoom simulations resembling present day spiral galaxies. On average the galaxies retain about 50\% of the cosmic fraction in baryons, almost equally divided into disc (interstellar medium) gas, cool CGM gas and warm-hot CGM gas. At radii smaller than 50 kpc the CGM is dominated by… ▽ More We present an analysis of the origin and properties of the circum-galactic medium (CGM) in a suite of 11 cosmological zoom simulations resembling present day spiral galaxies. On average the galaxies retain about 50\% of the cosmic fraction in baryons, almost equally divided into disc (interstellar medium) gas, cool CGM gas and warm-hot CGM gas. At radii smaller than 50 kpc the CGM is dominated by recycled warm-hot gas injected from the central galaxy, while at larger radii it is dominated by cool gas accreted onto the halo. The recycled gas typically accounts for one-third of the CGM mass. We introduce the novel publicly available analysis tool \textsc{pygad} to compute ion abundances and mock absorption spectra. For Lyman-$α$ absorption we find good agreement of the simulated equivalent width (EW) distribution and observations out to large radii. Disc galaxies with quiescent assembly histories show significantly more absorption along the disc major axis. By comparing the EW and HI column densities we find that CGM Lyman-$α$ absorbers are best represented by an effective line-width $b\approx 50 - 70$ km s$^{-1}$ that increases mildly with halo mass, larger than typically assumed. △ Less

Submitted 18 May, 2020; originally announced May 2020.

Comments: 18 pages, 20 figures

Comparing galaxy formation in the L-GALAXIES semi-analytical model and the IllustrisTNG simulations

Authors: Mohammadreza Ayromlou, Dylan Nelson, Robert M. Yates, Guinevere Kauffmann, Malin Renneby, Simon D. M. White

Abstract: We perform a comparison, object-by-object and statistically, between the Munich semi-analytical model, L-Galaxies, and the IllustrisTNG hydrodynamical simulations. By running L-Galaxies on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate and gas content of galaxies, as well as the baryon… ▽ More We perform a comparison, object-by-object and statistically, between the Munich semi-analytical model, L-Galaxies, and the IllustrisTNG hydrodynamical simulations. By running L-Galaxies on the IllustrisTNG dark matter-only merger trees, we identify the same galaxies in the two models. This allows us to compare the stellar mass, star formation rate and gas content of galaxies, as well as the baryonic content of subhaloes and haloes in the two models. We find that both the stellar mass functions and the stellar masses of individual galaxies agree to better than $\sim0.2\,$dex. On the other hand, specific star formation rates and gas contents can differ more substantially. At $z=0$ the transition between low-mass star-forming galaxies and high-mass, quenched galaxies occurs at a stellar mass scale $\sim0.5\,$dex lower in IllustrisTNG than in L-Galaxies. IllustrisTNG also produces substantially more quenched galaxies at higher redshifts. Both models predict a halo baryon fraction close to the cosmic value for clusters, but IllustrisTNG predicts lower baryon fractions in group environments. These differences are due primarily to differences in modelling feedback from stars and supermassive black holes. The gas content and star formation rates of galaxies in and around clusters and groups differ substantially, with IllustrisTNG satellites less star-forming and less gas-rich. We show that environmental processes such as ram-pressure stripping are stronger and operate to larger distances and for a broader host mass range in IllustrisTNG. We suggest that the treatment of galaxy evolution in the semi-analytic model needs to be improved by prescriptions which capture local environmental effects more accurately. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: Submitted to MNRAS. The catalogue of L-Galaxies run on TNG is publicly available at https://tng-project.org/ayromlou20

L-GALAXIES 2020: Spatially resolved cold gas phases, star formation and chemical enrichment in galactic discs

Authors: Bruno Henriques, Rob Yates, Jian Fu, Qi Guo, Guinevere Kauffmann, Chaichalit Srisawat, Peter Thomas, Simon White

Abstract: We have updated the Munich galaxy formation model, L-Galaxies, to follow the radial distributions of stars and atomic and molecular gas in galaxy discs. We include an H2-based star-formation law, as well as a detailed chemical-enrichment model with explicit mass-dependent delay times for SN-II, SN-Ia and AGB stars. Information about the star formation, feedback and chemical-enrichment histories of… ▽ More We have updated the Munich galaxy formation model, L-Galaxies, to follow the radial distributions of stars and atomic and molecular gas in galaxy discs. We include an H2-based star-formation law, as well as a detailed chemical-enrichment model with explicit mass-dependent delay times for SN-II, SN-Ia and AGB stars. Information about the star formation, feedback and chemical-enrichment histories of discs is stored in 12 concentric rings. The new model retains the success of its predecessor in reproducing the observed evolution of the galaxy population, in particular, stellar mass functions and passive fractions over the redshift range 0<=z<=3 and mass range 8<=log(M_*/Msun)<=12, the black hole-bulge mass relation at z=0, galaxy morphology as a function of stellar mass and the mass-metallicity relations of both stellar and gas components. In addition, its detailed modelling of the radial structure of discs allows qualitatively new comparisons with observation, most notably with the relative sizes and masses of the stellar, atomic and molecular components in discs. Good agreement is found with recent data. Comparison of results obtained for simulations differing in mass resolution by more than two orders of magnitude shows that all important distributions are numerically well converged even for this more detailed model. An examination of metallicity and surface-density gradients in the stars and gas indicates that our new model, with star formation, chemical enrichment and feedback calculated self-consistently on local disc scales, reproduces some but not all of the trends seen in recent many-galaxy IFU surveys. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: published in MNRAS. A full description of the model is available as MNRAS online supplementary material and at https://lgalaxiespublicrelease.github.io/. Data from the model are publicly released in SQL-queryable form and in flat files at http://www.mpa-garching.mpg.de/millennium

Journal ref: 2020MNRAS.491.5795H

Gas accretion and galactic fountain flows in the Auriga cosmological simulations: angular momentum and metal re-distribution

Authors: Robert J. J. Grand, Freeke van de Voort, Jolanta Zjupa, Francesca Fragkoudi, Facundo A. Gómez, Guinevere Kauffmann, Federico Marinacci, Rüdiger Pakmor, Volker Springel, Simon D. M. White

Abstract: Using a set of 15 high-resolution magnetohydrodynamic cosmological simulations of Milky Way formation, we investigate the origin of the baryonic material found in stars at redshift zero. We find that roughly half of this material originates from subhalo/satellite systems and half is smoothly accreted from the Inter-Galactic Medium (IGM). About $90 \%$ of all material has been ejected and re-accret… ▽ More Using a set of 15 high-resolution magnetohydrodynamic cosmological simulations of Milky Way formation, we investigate the origin of the baryonic material found in stars at redshift zero. We find that roughly half of this material originates from subhalo/satellite systems and half is smoothly accreted from the Inter-Galactic Medium (IGM). About $90 \%$ of all material has been ejected and re-accreted in galactic winds at least once. The vast majority of smoothly accreted gas enters into a galactic fountain that extends to a median galactocentric distance of $\sim 20$ kpc with a median recycling timescale of $\sim 500$ Myr. We demonstrate that, in most cases, galactic fountains acquire angular momentum via mixing of low-angular momentum, wind-recycled gas with high-angular momentum gas in the Circum-Galactic Medium (CGM). Prograde mergers boost this activity by helping to align the disc and CGM rotation axes, whereas retrograde mergers cause the fountain to lose angular momentum. Fountain flows that promote angular momentum growth are conducive to smooth evolution on tracks quasi-parallel to the disc sequence of the stellar mass-specific angular momentum plane, whereas retrograde minor mergers, major mergers and bar-driven secular evolution move galaxies towards the bulge-sequence. Finally, we demonstrate that fountain flows act to flatten and narrow the radial metallicity gradient and metallicity dispersion of disc stars, respectively. Thus, the evolution of galactic fountains depends strongly on the cosmological merger history and is crucial for the chemo-dynamical evolution of Milky Way-sized disc galaxies. △ Less

Submitted 15 October, 2019; v1 submitted 9 September, 2019; originally announced September 2019.

Comments: MNRAS; accepted. Changes include: clarifications and additional discussion, particularly of the fountain flow statistics shown in Fig. 7; and a resolution study that shows the results are converged for two resolution levels (~10^5 and ~10^4 Msun per baryonic element)

Emission line properties of the most luminous AGN in massive galaxies at intermediate redshifts

Authors: Guinevere Kauffmann, Claudia Maraston

Abstract: We have analyzed the emission line properties of 6019 Type II AGN at redshifts between 0.4-0.8 with [OIII] luminosities greater than 3x10^8 Lsun, characteristic of the Type II quasars first identified in population studies by Zakamska et al (2003). The AGN are drawn from the CMASS sample of galaxies with stellar masses greater than 10^11 Msun that were studied as part of the Baryon Oscillation Spe… ▽ More We have analyzed the emission line properties of 6019 Type II AGN at redshifts between 0.4-0.8 with [OIII] luminosities greater than 3x10^8 Lsun, characteristic of the Type II quasars first identified in population studies by Zakamska et al (2003). The AGN are drawn from the CMASS sample of galaxies with stellar masses greater than 10^11 Msun that were studied as part of the Baryon Oscillation Spectroscopic Survey (BOSS) and comprise 0.5% of the total population of these galaxies. Individual spectra have low S/N, so the analysis is carried out on stacked spectra in bins of [OIII] luminosity and estimated stellar age. The emission line ratios of the stacks are well fit with simple uniform-density photoionization models with metallicities between solar and twice solar. In the stacks, a number of emission lines are found to have distinct broad components requiring a double Gaussian rather than a single Gaussian fit, indicative of outflowing ionized gas. These are: [OIII]4959, [OIII]5007, [OII]3727,3729 and Halpha$. Higher ionization lines such as [NeIII]3869 and [NeV]3345 are detected in the stacks, but are well fit by single Gaussians. The broad components typically contain a third of the total line flux and have widths of 600 km/s for the oxygen lines and 900 km/s for Halpha. The fraction of the flux in the broad component and its width are independent of [OIII] luminosity, stellar age, radio and mid-IR luminosity. The stellar mass of the galaxy is the only parameter we could identify that influences the width of the broad line component. △ Less

Submitted 14 August, 2019; originally announced August 2019.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

The morphology and kinematics of the gaseous circumgalactic medium of Milky Way mass galaxies -- II. comparison of IllustrisTNG and Illustris simulation results

Authors: Guinevere Kauffmann, Dylan Nelson, Sanchayeeta Borthakur, Timothy Heckman, Lars Hernquist, Federico Marinacci, Ruediger Pakmor, Annelisa Pillepich

Abstract: We have carried out a controlled comparison of the structural and kinematic properties of the circumgalactic medium (CGM) around Milky Way mass galaxies in the Illustris and IllustrisTNG simulations. Very striking differences are found. At z=0, gas column density and temperature profiles at large radii ($\sim 100$ kpc) correlate strongly with disk gas mass fraction in Illustris, but not in TNG. Th… ▽ More We have carried out a controlled comparison of the structural and kinematic properties of the circumgalactic medium (CGM) around Milky Way mass galaxies in the Illustris and IllustrisTNG simulations. Very striking differences are found. At z=0, gas column density and temperature profiles at large radii ($\sim 100$ kpc) correlate strongly with disk gas mass fraction in Illustris, but not in TNG. The neutral gas at large radii is preferentially aligned in the plane of the disk in TNG, whereas it is much more isotropic in Illustris. The vertical coherence scale of the rotationally supported gas in the CGM is linked to the gas mass fraction of the galaxy in Illustris, but not in TNG. A tracer particle analysis allows us to show how these differences can be understood as a consequence of the different sub-grid models of feedback in the two simulations. A study of spatially matched galaxies in the two simulations shows that in TNG, feedback by supernovae and AGN helps to create an extended smooth reservoir of hot gas at high redshifts, that then cools to form a thin, rotationally-supported disk at later times. In Illustris, AGN dump heat in the form of hot gas bubbles that push diffuse material at large radii out of the halo. The disk is formed by accretion of colder, recycled material, and this results in more vertically extended gas distributions above and below the Galactic plane. We conclude that variations in the structure of gas around Milky Way mass galaxies are a sensitive probe of feedback physics in simulations and are worthy of more observational consideration in future. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 16 pages, 14 figures, accepted for publication in MNRAS

A New Method to Quantify Environment and Model Ram-Pressure Stripping in N-Body Simulations

Authors: Mohammadreza Ayromlou, Dylan Nelson, Robert M. Yates, Guinevere Kauffmann, Simon D. M. White

Abstract: We introduce a Local Background Environment (LBE) estimator that can be measured in and around every galaxy or its dark matter subhalo in high-resolution cosmological simulations. The LBE is designed to capture the influence of environmental effects such as ram-pressure stripping on the formation and evolution of galaxies in semi-analytical models. We define the LBE directly from the particle data… ▽ More We introduce a Local Background Environment (LBE) estimator that can be measured in and around every galaxy or its dark matter subhalo in high-resolution cosmological simulations. The LBE is designed to capture the influence of environmental effects such as ram-pressure stripping on the formation and evolution of galaxies in semi-analytical models. We define the LBE directly from the particle data within an adaptive spherical shell, and devise a Gaussian mixture estimator to separate background particles from previously unidentified subhalo particles. Analyzing the LBE properties, we find that the LBE of satellite galaxies is not at rest with respect to their host halo, in contrast to typical assumptions. The orientations of the velocities of a subhalo and its LBE are well aligned in the outer infall regions of haloes, but decorrelated near halo center. Significantly, there is no abrupt change in LBE velocity or density at the halo virial radius. This suggests that stripping should also happen beyond this radius. Therefore, we use the time-evolving LBE of galaxies to develop a method to better account for ram-pressure stripping within the Munich semi-analytical model, L-Galaxies. Overall, our new approach results in a significant increase in gas stripping across cosmic time. Central galaxies, as well as satellites beyond the virial radius, can lose a significant fraction of their hot halo gas. As a result, the gas fractions and star formation rates of satellite galaxies are suppressed relative to the fiducial model, although the stellar masses and global stellar mass functions are largely unchanged. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: Submitted to MNRAS

A 1-dimensional hydrodynamic model for accretion, cooling and heating of gas in dark matter halos from $z=6$ to $z=0$

Authors: Prakriti Pal Choudhury, Guinevere Kauffmann, Prateek Sharma

Abstract: We study an idealized 1D model for the evolution of hot gas in dark matter halos for redshifts $z=[0,6]$. We introduce a numerical setup incorporating cosmological accretion of gas, along with the growth of the halo, based on the Van den Bosch model for the average growth of halos as a function of cosmic time. We evolve one-dimensional Lagrangian shells with radiative cooling of the gas and heatin… ▽ More We study an idealized 1D model for the evolution of hot gas in dark matter halos for redshifts $z=[0,6]$. We introduce a numerical setup incorporating cosmological accretion of gas, along with the growth of the halo, based on the Van den Bosch model for the average growth of halos as a function of cosmic time. We evolve one-dimensional Lagrangian shells with radiative cooling of the gas and heating due to feedback from the gas cooling and moving in toward the center. A simple Bondi accretion model on to a central black hole is used to include feedback heating. The setup captures some of the key characteristics of spherically symmetric accretion onto the halos: formation of virial shocks slightly outside $r_{200}$ and long-term thermal balance in the form of cooling and heating cycles. The gas density outside our initial halos at $z=6$ is constrained by requiring that the baryon fraction within the virial radius for non-radiative evolution be equal to the universal value at almost all times. The total mass in the cold phase (taken to be $\sim 10^4$ K) within $40$ kpc is tracked as a function of the halo mass and redshift. We compare the evolution of the cold gas mass to the observed stellar-mass versus halo mass relations, following which, we can constrain the feedback energy required for different halo masses and redshifts. We also compare and match the hot gas density and temperature profiles for our most massive halo to those of clusters observed upto redshift $2$. Our model is thus an improvement over the semi-analytic models in which isothermal condition and $ρ\propto r^{-2}$ are assumed. △ Less

Submitted 25 February, 2019; v1 submitted 15 August, 2018; originally announced August 2018.

Comments: 17 pages, 13 figures, 3 tables. Accepted in MNRAS. Comments are welcome

Properties of AGN selected by their mid-IR colours: evidence for a physically distinct mode of black hole growth

Authors: Guinevere Kauffmann

Abstract: We study the narrow emission line properties and stellar populations of a sample of 1385 AGN selected to have strong excess emission at mid-infrared wavelengths based on comparing Wide-field Infrared Survey Explorer W1-W2 band colours with optical stellar absorption line indicators. Our goal is to understand whether the physical conditions in the interstellar medium of these objects differ from th… ▽ More We study the narrow emission line properties and stellar populations of a sample of 1385 AGN selected to have strong excess emission at mid-infrared wavelengths based on comparing Wide-field Infrared Survey Explorer W1-W2 band colours with optical stellar absorption line indicators. Our goal is to understand whether the physical conditions in the interstellar medium of these objects differ from those of local AGN selected by their optical emission line ratios. To enable this comparison, we construct a control sample of 50,000 optically-selected AGN with the same redshifts that do not have strong mid-IR excess emission. The mid-IR excess and control samples differ strongly in [OIII] line luminosity, ionized gas excitation mechanism, ionization state and electron density. We show that the radio-detected, mid-IR excess AGN constitute the most luminous and highly ionized AGN in the local Universe and they contribute primarily to the growth of black holes in the most massive galaxies. At least half of this black hole growth is occurring in galaxies with recent starbursts. The morphologies of these systems indicate that the starbursts have probably been triggered by galaxy-galaxy mergers and interactions. The most luminous AGN in our mid-IR excess sample have properties that are similar to the Type II quasars studied at higher redshifts. In contrast, the control sample constitute a class of lower ionization, less luminous AGN in more quiescent galaxies that contribute primarily to the growth of low mass black holes. △ Less

Submitted 9 August, 2018; originally announced August 2018.

Comments: 15 pages, 17 figures, accepted for publication in MNRAS

Journal ref: 2018MNRAS.tmp.1935K

The Auriga Stellar Haloes: Connecting stellar population properties with accretion and merging history

Authors: Antonela Monachesi, Facundo A. Gómez, Robert J. J. Grand, Christine M. Simpson, Guinevere Kauffmann, Sebastián Bustamante, Federico Marinacci, Rüdiger Pakmor, Volker Springel, Carlos S. Frenk, Simon D. M. White, Patricia B. Tissera

Abstract: We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and… ▽ More We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass--metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies. △ Less

Submitted 21 February, 2019; v1 submitted 20 April, 2018; originally announced April 2018.

Comments: Accepted to MNRAS. 30 pages, 19 figures

Large scale correlations in gas traced by MgII absorbers around low mass galaxies

Authors: Guinevere Kauffmann

Abstract: The physical origin of the large-scale conformity in the colours and specific star formation rates of isolated low mass central galaxies and their neighbours on scales in excess of 1 Mpc is still under debate. One possible scenario is that gas is heated over large scales by feedback from active galactic nuclei (AGN), leading to coherent modulation of cooling and star formation between well-separat… ▽ More The physical origin of the large-scale conformity in the colours and specific star formation rates of isolated low mass central galaxies and their neighbours on scales in excess of 1 Mpc is still under debate. One possible scenario is that gas is heated over large scales by feedback from active galactic nuclei (AGN), leading to coherent modulation of cooling and star formation between well-separated galaxies. In this Letter, the metal line absorption catalogue of Zhu & Menard (2013) is used to probe gas out to large projected radii around a sample of a million galaxies with stellar masses ~10^{10} M_{sun} and photometric redshifts in the range 0.4<z<0.8 selected from Sloan Digital Sky Survey imaging data. This galaxy sample covers an effective volume of 2.2 Gpc^3. A statistically significant excess of MgII absorbers is present around the red low mass galaxies compared to their blue counterparts out to projected radii of 10 Mpc. In addition, the equivalent width distribution function of MgII absorbers around low mass galaxies is shown to be strongly affected by the presence of a nearby (R_p<2 Mpc) radio-loud AGN out to projected radii of 5 Mpc. △ Less

Submitted 17 December, 2017; originally announced December 2017.

Comments: 5 pages, accepted in MNRAS Letters

The abundance, distribution, and physical nature of highly ionized oxygen OVI, OVII, and OVIII in IllustrisTNG

Authors: Dylan Nelson, Guinevere Kauffmann, Annalisa Pillepich, Shy Genel, Volker Springel, Ruediger Pakmor, Lars Hernquist, Rainer Weinberger, Paul Torrey, Mark Vogelsberger, Federico Marinacci

Abstract: We explore the abundance, spatial distribution, and physical properties of the OVI, OVII, and OVIII ions of oxygen in circumgalactic and intergalactic media (the CGM, IGM, and WHIM). We use the TNG100 and TNG300 large volume cosmological magneto-hydrodynamical simulations. Modeling the ionization states of simulated oxygen, we find good agreement with observations of the low-redshift OVI column de… ▽ More We explore the abundance, spatial distribution, and physical properties of the OVI, OVII, and OVIII ions of oxygen in circumgalactic and intergalactic media (the CGM, IGM, and WHIM). We use the TNG100 and TNG300 large volume cosmological magneto-hydrodynamical simulations. Modeling the ionization states of simulated oxygen, we find good agreement with observations of the low-redshift OVI column density distribution function (CDDF), and present its evolution for all three ions from z=0 to z=4. Producing mock quasar absorption line spectral surveys, we show that the IllustrisTNG simulations are fully consistent with constraints on the OVI content of the CGM from COS-Halos and other low redshift observations, producing columns as high as observed. We measure the total amount of mass and average column densities of each ion using hundreds of thousands of simulated galaxies spanning 10^11 < Mhalo/Msun < 10^15 corresponding to 10^9 < M*/Msun < 10^12 in stellar mass. Stacked radial profiles of OVI are computed in 3D number density and 2D projected column, decomposing into the 1-halo and 2-halo terms. Relating halo OVI to properties of the central galaxy, we find a correlation between the (g-r) color of a galaxy and the total amount of OVI in its CGM. In comparison to the COS-Halos finding, this leads to a dichotomy of columns around star-forming versus passive galaxies at fixed stellar (or halo) mass. We demonstrate that this correlation is a direct result of blackhole feedback associated with quenching and represents a causal consequence of galactic-scale baryonic feedback impacting the physical state of the circumgalactic medium. △ Less

Submitted 15 March, 2018; v1 submitted 30 November, 2017; originally announced December 2017.

Comments: MNRAS. For an interactive OVI/OVIII visualization, see the 2D explorer at http://www.tng-project.org (fig 5 updated in v2)

The Physical Properties of Galaxies with Unusually Red Mid-Infrared Colours

Authors: Guinevere Kauffmann

Abstract: The goal of this paper is to investigate the physical nature of galaxies in the redshift range $0.02<z<0.15$ that have strong excess emission at mid-IR wavelengths and to determine whether they host a population of accreting black holes that cannot be identified using optical emission lines. We show that at fixed stellar mass $M_*$ and $D_n(4000)$, the distribution of [3.4]-[4.6] $μ$m (WISE W1-W2… ▽ More The goal of this paper is to investigate the physical nature of galaxies in the redshift range $0.02<z<0.15$ that have strong excess emission at mid-IR wavelengths and to determine whether they host a population of accreting black holes that cannot be identified using optical emission lines. We show that at fixed stellar mass $M_*$ and $D_n(4000)$, the distribution of [3.4]-[4.6] $μ$m (WISE W1-W2 band) colours is sharply peaked, with a long tail to much redder W1-W2 colours. We introduce a procedure to pull out the red outlier population based on a combination of three stellar population diagnostics. When compared with optically-selected AGN, red outliers are more likely to be found in massive galaxies, and they tend to have lower stellar mass densities, younger stellar ages and higher dust content than optically-selected AGN hosts. They are twice as likely to be detected at radio wavelengths. We examine W1-W2 colour profiles for a subset of the nearest, reddest outliers and find that most are not centrally peaked, indicating that the hot dust emission is spread throughout the galaxy. We find that radio luminosity is the quantity that is most predictive of a redder central W1-W2 colour. Radio-loud galaxies with centrally concentrated hot dust emission are almost always morphologically disturbed, with compact, unresolved emission at 1.4 Ghz. Eighty percent of such systems are identifiable as AGN using optical emission line diagnostics. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: 12 pages, 15 figures, accepted in MNRAS

xCOLD GASS: the complete IRAM-30m legacy survey of molecular gas for galaxy evolution studies

Authors: Amélie Saintonge, Barbara Catinella, Linda J. Tacconi, Guinevere Kauffmann, Reinhard Genzel, Luca Cortese, Romeel Davé, Thomas J. Fletcher, Javier Graciá-Carpio, Carsten Kramer, Timothy M. Heckman, Steven Janowiecki, Katharina Lutz, David Rosario, David Schiminovich, Karl Schuster, Jing Wang, Stijn Wuyts, Sanchayeeta Borthakur, Isabella Lamperti, Guido W. Roberts-Borsani

Abstract: We introduce xCOLD GASS, a legacy survey providing a census of molecular gas in the local Universe. Building upon the original COLD GASS survey, we present here the full sample of 532 galaxies with CO(1-0) measurements from the IRAM-30m telescope. The sample is mass-selected in the redshift interval $0.01<z<0.05$ from SDSS, and therefore representative of the local galaxy population with M… ▽ More We introduce xCOLD GASS, a legacy survey providing a census of molecular gas in the local Universe. Building upon the original COLD GASS survey, we present here the full sample of 532 galaxies with CO(1-0) measurements from the IRAM-30m telescope. The sample is mass-selected in the redshift interval $0.01<z<0.05$ from SDSS, and therefore representative of the local galaxy population with M$_{\ast}>10^9$M$_{\odot}$. The CO(1-0) flux measurements are complemented by observations of the CO(2-1) line with both the IRAM-30m and APEX telescopes, HI observations from Arecibo, and photometry from SDSS, WISE and GALEX. Combining the IRAM and APEX data, we find that the CO(2-1) to CO(1-0) luminosity ratio for integrated measurements is $r_{21}=0.79\pm0.03$, with no systematic variations across the sample. The CO(1-0) luminosity function is constructed and best fit with a Schechter function with parameters {$L_{\mathrm{CO}}^* = (7.77\pm2.11) \times 10^9\,\mathrm{K\,km\,s^{-1}\, pc^{2}}$, $φ^{*} = (9.84\pm5.41) \times 10^{-4} \, \mathrm{Mpc^{-3}}$ and $α= -1.19\pm0.05$}. With the sample now complete down to stellar masses of $10^9$M$_{\odot}$, we are able to extend our study of gas scaling relations and confirm that both molecular gas fraction and depletion timescale vary with specific star formation rate (or offset from the star-formation main sequence) much more strongly than they depend on stellar mass. Comparing the xCOLD GASS results with outputs from hydrodynamic and semi-analytic models, we highlight the constraining power of cold gas scaling relations on models of galaxy formation. △ Less

Submitted 5 October, 2017; originally announced October 2017.

Comments: 24 pages, 16 figures. Accepted for publication in ApJS. The full xCOLD GASS survey data products are available at http://www.star.ucl.ac.uk/xCOLDGASS/

Journal ref: ApJS, 233 (2017) 22

First results from the IllustrisTNG simulations: the galaxy color bimodality

Authors: Dylan Nelson, Annalisa Pillepich, Volker Springel, Rainer Weinberger, Lars Hernquist, Ruediger Pakmor, Shy Genel, Paul Torrey, Mark Vogelsberger, Guinevere Kauffmann, Federico Marinacci, Jill Naiman

Abstract: We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger. Here we present first results on the galaxy color bimodality at low reds… ▽ More We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger. Here we present first results on the galaxy color bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g-r) colors of 10^9 < M*/Msun < 10^12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey (SDSS). We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement in comparison to the observations, with a sharp transition in median color from blue to red at a characteristic M* ~ 10^10.5 Msun. Investigating the build-up of the color-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy color transition in the TNG model is supermassive blackhole feedback in its low-accretion state. Across the entire population we measure a median color transition timescale dt_green of ~1.6 Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, redder galaxies have lower SFRs, gas fractions, and gas metallicities; their stellar populations are also older and their large-scale interstellar magnetic fields weaker than in bluer galaxies. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* > 10^11 Msun which redden at z<1 accumulate on average ~25% of their final z=0 mass post-reddening; at the same time, ~18% of such massive galaxies acquire half or more of their final stellar mass while on the red sequence. △ Less

Submitted 31 January, 2018; v1 submitted 11 July, 2017; originally announced July 2017.

Comments: The IllustrisTNG project website is http://www.tng-project.org

The SDSS-IV MaNGA Sample: Design, Optimization, and Usage Considerations

Authors: David A. Wake, Kevin Bundy, Aleksandar M. Diamond-Stanic, Renbin Yan, Michael R. Blanton, Matthew A. Bershady, José R. Sánchez-Gallego, Niv Drory, Amy Jones, Guinevere Kauffmann, David R. Law, Cheng Li, Nicholas MacDonald, Karen Masters, Daniel Thomas, Jeremy Tinker, Anne-Marie Weijmans, Joel R. Brownstein

Abstract: We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey de… ▽ More We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey defined in terms of maximizing S/N, spatial resolution, and sample size. Our selection strategy makes use of redshift limits that only depend on i-band absolute magnitude ($M_i$), or, for a small subset of our sample, $M_i$ and color (NUV-i). Such a strategy ensures that all galaxies span the same range in angular size irrespective of luminosity and are therefore covered evenly by the adopted range of IFU sizes. We define three samples: the Primary and Secondary samples are selected to have a flat number density with respect to $M_i$ and are targeted to have spectroscopic coverage to 1.5 and 2.5 effective radii (Re), respectively. The Color-Enhanced supplement increases the number of galaxies in the low-density regions of color-magnitude space by extending the redshift limits of the Primary sample in the appropriate color bins. The samples cover the stellar mass range $5\times10^8 \leq M_* \leq 3\times10^{11} M_{\odot}$ and are sampled at median physical resolutions of 1.37 kpc and 2.5 kpc for the Primary and Secondary samples respectively. We provide weights that will statistically correct for our luminosity and color-dependent selection function and IFU allocation strategy, thus correcting the observed sample to a volume limited sample. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: Accepted for publication in AJ

Clustering of MgII absorption line systems around massive galaxies: an important constraint on feedback processes in galaxy formation

Authors: Guinevere Kauffmann, Dylan Nelson, Brice Menard, Guangtun Zhu

Abstract: We use the latest version of the metal line absorption catalogue of Zhu & Ménard (2013) to study the clustering of MgII absorbers around massive galaxies (~10^11.5 M_sun), quasars and radio-loud AGN with redshifts between 0.4 and 0.75. Clustering is evaluated in two dimensions, by binning absorbers both in projected radius and in velocity separation. Excess MgII is detected around massive galaxies… ▽ More We use the latest version of the metal line absorption catalogue of Zhu & Ménard (2013) to study the clustering of MgII absorbers around massive galaxies (~10^11.5 M_sun), quasars and radio-loud AGN with redshifts between 0.4 and 0.75. Clustering is evaluated in two dimensions, by binning absorbers both in projected radius and in velocity separation. Excess MgII is detected around massive galaxies out to R_p=20 Mpc. At projected radii less than 800 kpc, the excess extends out to velocity separations of 10,000 km/s. The extent of the high velocity tail within this radius is independent of the mean stellar age of the galaxy and whether or not it harbours an active galactic nucleus. We interpret our results using the publicly available Illustris and Millennium simulations. Models where the MgII absorbers trace the dark matter particle or subhalo distributions do not fit the data. They overpredict the clustering on small scales and do not reproduce the excess high velocity separation MgII absorbers seen within the virial radius of the halo. The Illustris simulations which include thermal, but not mechanical feedback from AGN, also do not provide an adequate fit to the properties of the cool halo gas within the virial radius. We propose that the large velocity separation MgII absorbers trace gas that has been pushed out of the dark matter halos, possibly by multiple episodes of AGN-driven mechanical feedback acting over long timescales. △ Less

Submitted 14 March, 2017; originally announced March 2017.

Comments: 10 pages, 11 figures, accepted in MNRAS

SDSS IV MaNGA: Deep observations of extra-planar, diffuse ionized gas around late-type galaxies from stacked IFU spectra

Authors: A. Jones, G. Kauffmann, R. D'Souza, D. Bizyaev, D. Law, L. Haffner, Y. Bahé, B. Andrews, M. Bershady, J. Brownstein, K. Bundy, B. Cherinka, A. Diamond-Stanic, N. Drory, R. A. Riffel, S. F. Sánchez D. Thomas, D. Wake, R. Yan, K. Zhang

Abstract: We have conducted a study of extra-planar diffuse ionized gas using the first year data from the MaNGA IFU survey. We have stacked spectra from 49 edge-on, late-type galaxies as a function of distance from the midplane of the galaxy. With this technique we can detect the bright emission lines Halpha, Hbeta, [OII]3726, 3729, [OIII]5007, [NII]6549, 6584, and [SII]6717, 6731 out to about 4 kpc above… ▽ More We have conducted a study of extra-planar diffuse ionized gas using the first year data from the MaNGA IFU survey. We have stacked spectra from 49 edge-on, late-type galaxies as a function of distance from the midplane of the galaxy. With this technique we can detect the bright emission lines Halpha, Hbeta, [OII]3726, 3729, [OIII]5007, [NII]6549, 6584, and [SII]6717, 6731 out to about 4 kpc above the midplane. With 16 galaxies we can extend this analysis out to about 9 kpc, i.e. a distance of ~2R_e, vertically from the midplane. In the halo, the surface brightnesses of the [OII] and Halpha emission lines are comparable, unlike in the disk where Halpha dominates. When we split the sample by specific star formation rate, concentration index, and stellar mass, each subsample's emission line surface brightness profiles and ratios differ, indicating that extra-planar gas properties can vary. The emission line surface brightnesses of the gas around high specific star formation rate galaxies are higher at all distances, and the line ratios are closer to ratios characteristic of HII regions compared with low specific star formation rate galaxies. The less concentrated and lower stellar mass samples exhibit line ratios that are more like HII regions at larger distances than their more concentrated and higher stellar mass counterparts. The largest difference between different subsamples occurs when the galaxies are split by stellar mass. We additionally infer that gas far from the midplane in more massive galaxies has the highest temperatures and steepest radial temperature gradients based on their [NII]/Halpha and [OII]/Halpha ratios between the disk and the halo. △ Less

Submitted 12 December, 2016; originally announced December 2016.

Comments: 19 pages, 14 figures, 2 appendices. Accepted in A&A (see https://drive.google.com/open?id=0B_MhxrO1dvVNUUUtb0FXdllqb0E for full version)

Journal ref: A&A 599, A141 (2017)

The Properties of the Circumgalactic Medium in Red and Blue Galaxies: Results from the COS-GASS+COS-Halos Surveys

Authors: Sanchayeeta Borthakur, Timothy Heckman, Jason Tumlinson, Rongmon Bordoloi, Guinevere Kauffmann, Barbara Catinella, David Schiminovich, Romeel Dave, Sean M. Moran, Amelie Saintonge

Abstract: We use the combined data from the COS-GASS and COS-Halos surveys to characterize the Circum-Galactic Medium (CGM) surrounding typical low-redshift galaxies in the mass range $\rm~M_*\sim~10^{9.5-11.5}~M_{\odot} $, and over a range of impact parameters extending to just beyond the halo virial radius ($\rm~R_{vir}$). We find the radial scale length of the distributions of the equivalent widths of th… ▽ More We use the combined data from the COS-GASS and COS-Halos surveys to characterize the Circum-Galactic Medium (CGM) surrounding typical low-redshift galaxies in the mass range $\rm~M_*\sim~10^{9.5-11.5}~M_{\odot} $, and over a range of impact parameters extending to just beyond the halo virial radius ($\rm~R_{vir}$). We find the radial scale length of the distributions of the equivalent widths of the Lyman~$α$ and Si III absorbers to be 0.9 and 0.4 $\rm~R_{vir}$, respectively. The radial distribution of equivalent widths is relatively uniform for the blue galaxies, but highly patchy (low covering fraction) for the red galaxies. We also find that the Lyman~$α$ and Si III equivalent widths show significant positive correlations with the specific star-formation rate (sSFR) of the galaxy. We find a surprising lack of correlations between the halo mass (virial velocity) and either the velocity dispersions or velocity offsets of the Lyman~$α$ lines. The ratio of the velocity offset to the velocity dispersion for the Lyman~$α$ absorbers has a mean value of $\sim$ 4, suggesting that a given the line-of-sight is intersecting a dynamically coherent structure in the CGM rather than a sea of orbiting clouds. The kinematic properties of the CGM are similar in the blue and red galaxies, although we find that a significantly larger fraction of the blue galaxies have large Lyman~$α$ velocity offsets (>200 km s$^{-1}$). We show that - if the CGM clouds represent future fuel for star-formation - our new results could imply a large drop in the specific star-formation rate across the galaxy mass-range we probe. △ Less

Submitted 20 September, 2016; originally announced September 2016.

Comments: 18 pages, 18 figures, 2 tables. Accepted for publication in ApJ

The morphology and kinematics of neutral hydrogen in the vicinity of z=0 galaxies with Milky Way masses -- a study with the Illustris simulation

Authors: Guinevere Kauffmann, Sanchayeeta Borthakur, Dylan Nelson

Abstract: We analyze the properties of the circumgalactic gas (CGM) around 120 galaxies with stellar and dark matter halo masses similar to that of the Milky Way. We focus on the morphology and kinematics of the neutral hydrogen and how this depends on f_g, the ratio of gas-to-stellar mass within the optical radius. In gas-rich galaxies with f_g > 0.1, gas temperatures rise slowly from center of the halo ou… ▽ More We analyze the properties of the circumgalactic gas (CGM) around 120 galaxies with stellar and dark matter halo masses similar to that of the Milky Way. We focus on the morphology and kinematics of the neutral hydrogen and how this depends on f_g, the ratio of gas-to-stellar mass within the optical radius. In gas-rich galaxies with f_g > 0.1, gas temperatures rise slowly from center of the halo out to the virial radius and average neutral gas column densities remain above 10^19 atoms cm^-2 out to radii of 50-70 kpc. In gas-poor galaxies with f_g < 0.1, gas temperatures rise quickly outside the edge of the disk to 10^6 K, and then remain fixed out to radii of 100 kpc. The column density of neutral gas quickly drops below 10^19 atoms cm^-2 at radii of 10 kpc. Neutral gas distributions are also more asymmetric in gas-poor galaxies. Most of the differences between gas-poor and gas-rich galaxies in the Illustris simulation can be attributed to the effects of "radio-mode" AGN feedback. In the Illustris simulation, the circumgalactic gas is found to rotate coherently about the center of the galaxy with a maximum rotational velocity of around 200 km/s. In gas-rich galaxies, the average coherence length of the rotating gas is 40 kpc, compared to 10 kpc in gas-poor galaxies. In the very most gas-rich systems, the CGM can rotate coherently over scales of 70-100 kpc. We discuss our results in the context of recent observations of the CGM in low mass galaxies via UV absorption-line spectroscopy and deep 21cm observations of edge-on spiral galaxies. △ Less

Submitted 15 August, 2016; v1 submitted 11 August, 2016; originally announced August 2016.

Comments: 16 pages, 16 figures, accepted for publication in MNRAS

The Data Reduction Pipeline for the SDSS-IV MaNGA IFU Galaxy Survey

Authors: David R. Law, Brian Cherinka, Renbin Yan, Brett H. Andrews, Matthew A. Bershady, Dmitry Bizyaev, Guillermo A. Blanc, Michael R. Blanton, Adam S. Bolton, Joel R. Brownstein, Kevin Bundy, Yanmei Chen, Niv Drory, Richard D'Souza, Hai Fu, Amy Jones, Guinevere Kauffmann, Nicholas MacDonald, Karen L. Masters, Jeffrey A. Newman, John K. Parejko, José R. Sánchez-Gallego, Sebastian F. Sánchez, David J. Schlegel, Daniel Thomas , et al. (4 additional authors not shown)

Abstract: Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit explora… ▽ More Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ~ 100 million raw-frame spectra and ~ 10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline (DRP) algorithms and centralized metadata framework that produces sky-subtracted, spectrophotometrically calibrated spectra and rectified 3-D data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13 (DR13), we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ~ 8500 Angstroms and reach a typical 10-sigma limiting continuum surface brightness mu = 23.5 AB/arcsec^2 in a five arcsec diameter aperture in the g band. The wavelength calibration of the MaNGA data is accurate to 5 km/s rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of sigma = 72 km/s. △ Less

Submitted 2 August, 2016; v1 submitted 28 July, 2016; originally announced July 2016.

Comments: 37 pages, 22 figures. Accepted for publication in AJ. v2 updates arXiv reference to Yan+16b, v3 fixes some embedded DR13 links. Version with full resolution figures is available at http://www.stsci.edu/~dlaw/Papers/MaNGA_Data.pdf

Resolution-independent modeling of environmental effects in semi-analytic models of galaxy formation that include ram-pressure stripping of both hot and cold gas

Authors: Yu Luo, Xi Kang, Guinevere Kauffmann, Jian Fu

Abstract: The quenching of star formation in satellite galaxies is observed over a wide range of dark matter halo masses and galaxy environments. In the recent Guo et al (2011) and Fu et al (2013) semi-analytic + N-body models, the gaseous environment of the satellite galaxy is governed by the properties of the dark matter subhalo in which it resides. This quantity depends of the resolution of the N-body si… ▽ More The quenching of star formation in satellite galaxies is observed over a wide range of dark matter halo masses and galaxy environments. In the recent Guo et al (2011) and Fu et al (2013) semi-analytic + N-body models, the gaseous environment of the satellite galaxy is governed by the properties of the dark matter subhalo in which it resides. This quantity depends of the resolution of the N-body simulation, leading to a divergent fraction of quenched satellites in high- and low-resolution simulations. Here, we incorporate an analytic model to trace the subhaloes below the resolution limit. We demonstrate that we then obtain better converged results between the Millennium I and II simulations, especially for the satellites in the massive haloes ($\rm log M_{halo}=[14,15]$). We also include a new physical model for the ram-pressure stripping of cold gas in satellite galaxies. However, we find very clear discrepancies with observed trends in quenched satellite galaxy fractions as a function of stellar mass at fixed halo mass. At fixed halo mass, the quenched fraction of satellites does not depend on stellar mass in the models, but increases strongly with mass in the data. In addition to the over-prediction of low-mass passive satellites, the models also predict too few quenched central galaxies with low stellar masses, so the problems in reproducing quenched fractions are not purely of environmental origin. Further improvements to the treatment of the gas-physical processes regulating the star formation histories of galaxies are clearly necessary to resolve these problems. △ Less

Submitted 29 January, 2016; originally announced February 2016.

Comments: 15 papes, 11 figures, accepted for publication in MNRAS

On the stellar halo metallicity profile of Milky Way-like galaxies in the Auriga simulations

Authors: Antonela Monachesi, Facundo A. Gomez, Robert J. J. Grand, Guinevere Kauffmann, Federico Marinacci, Rüdiger Pakmor, Volker Springel, Carlos S. Frenk

Abstract: A recent observational study of haloes of nearby Milky Way-like galaxies shows that only half (four out of eight) of the current sample exhibits strong negative metallicity ([Fe/H]) gradients. This is at odds with predictions from hydrodynamical simulations where such gradients are ubiquitous. In this Letter, we use high resolution cosmological hydrodynamical simulations to study the [Fe/H] distri… ▽ More A recent observational study of haloes of nearby Milky Way-like galaxies shows that only half (four out of eight) of the current sample exhibits strong negative metallicity ([Fe/H]) gradients. This is at odds with predictions from hydrodynamical simulations where such gradients are ubiquitous. In this Letter, we use high resolution cosmological hydrodynamical simulations to study the [Fe/H] distribution of galactic haloes. We find that kinematically selected stellar haloes, including both in-situ and accreted particles, have an oblate [Fe/H] distribution. Spherical [Fe/H] radial profiles show strong negative gradients within 100 kpc, in agreement with previous numerical results. However, the projected median [Fe/H] profiles along the galactic disc minor axis, typically obtained in observations, are significantly flatter. The median [Fe/H] values at a given radius are larger for the spherical profiles than for the minor axis profiles by as much as 0.4 dex within the inner 50 kpc. Similar results are obtained if only the accreted stellar component is considered indicating that the differences between spherical and minor axis profiles are not purely driven by heated disc star particles formed in situ. Our study highlights the importance of performing careful comparisons between models and observations of halo [Fe/H] distributions. △ Less

Submitted 23 March, 2016; v1 submitted 9 December, 2015; originally announced December 2015.

Comments: 6 pages, 4 figures. Accepted for publication in MNRAS Letters; very minor edits to match the accepted version

The distribution of atomic hydrogen in EAGLE galaxies: morphologies, profiles, and HI holes

Authors: Yannick M. Bahe, Robert A. Crain, Guinevere Kauffmann, Richard G. Bower, Joop Schaye, Michelle Furlong, Claudia Lagos, Matthieu Schaller, James W. Trayford, Claudio Dalla Vecchia, Tom Theuns

Abstract: We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical pr… ▽ More We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical prescription from the literature. The resulting neutral hydrogen fractions, M_(HI+H2)/M_star, agree with observations to better than 0.1 dex for galaxies with M_star between 10^10 and 10^11 M_Sun. Our fiducial, empirical H2 model based on gas pressure results in galactic HI mass fractions, M_HI/M_star, that agree with observations from the GASS survey to better than 0.3 dex, but the alternative theoretical H2 formula leads to a negative offset in M_HI/M_star of up to 0.5 dex. Visual inspection reveals that most HI disks in simulated HI-rich galaxies are vertically disturbed, plausibly due to recent accretion events. Many galaxies (up to 80 per cent) contain spuriously large HI holes, which are likely formed as a consequence of the feedback implementation in EAGLE. The HI mass-size relation of all simulated galaxies is close to (but 16 per cent steeper than) observed, and when only galaxies without large holes in the HI disc are considered, the agreement becomes excellent (better than 0.1 dex). The presence of large HI holes also makes the radial HI surface density profiles somewhat too low in the centre, at Σ_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (Σ_HI < 1 M_Sun pc^-2), the simulated profiles agree quantitatively with observations. Scaled by HI size, the simulated profiles of HI-rich (M_HI > 10^9.8 M_Sun) and control galaxies (10^9.1 M_Sun > M_HI > 10^9.8 M_Sun) follow each other closely, as observed. (Abridged) △ Less

Submitted 16 November, 2015; originally announced November 2015.

Comments: 23 pages, 17 figures (16 pages, 9 figures without appendices), accepted by MNRAS

The Massive End of the Stellar Mass Function

Authors: R. D'Souza, S. Vegetti, G. Kauffmann

Abstract: We derive average flux corrections to the \texttt{Model} magnitudes of the Sloan Digital Sky Survey (SDSS) galaxies by stacking together mosaics of similar galaxies in bins of stellar mass and concentration. Extra flux is detected in the outer low surface brightness part of the galaxies, leading to corrections ranging from 0.05 to 0.32 mag for the highest stellar mass galaxies. We apply these corr… ▽ More We derive average flux corrections to the \texttt{Model} magnitudes of the Sloan Digital Sky Survey (SDSS) galaxies by stacking together mosaics of similar galaxies in bins of stellar mass and concentration. Extra flux is detected in the outer low surface brightness part of the galaxies, leading to corrections ranging from 0.05 to 0.32 mag for the highest stellar mass galaxies. We apply these corrections to the MPA-JHU (Max-Planck Institute for Astrophysics - John Hopkins University) stellar masses for a complete sample of half a million galaxies from the SDSS survey to derive a corrected galaxy stellar mass function at $z=0.1$ in the stellar mass range $9.5<\log(M_\ast/M_\odot)<12.0$. We find that the flux corrections and the use of the MPA-JHU stellar masses have a significant impact on the massive end of the stellar mass function, making the slope significantly shallower than that estimated by Li \& White (2009), but steeper than derived by Bernardi et al. (2013). This corresponds to a mean comoving stellar mass density of galaxies with stellar masses $\log(M_\ast/M_\odot) \ge 11.0$ that is a factor of 3.36 larger than the estimate by Li \& White (2009), but is 43\% smaller than reported by Bernardi et al. (2013). △ Less

Submitted 24 September, 2015; originally announced September 2015.

Comments: 11 pages, 8 figures, Accepted to MNRAS

Physical origin of the large-scale conformity in the specific star formation rates of galaxies

Authors: Guinevere Kauffmann

Abstract: Two explanations have been put forward to explain the observed conformity between the colours and specific star formation rates (SFR/$M_*$) of galaxies on large scales: 1) the formation times of their surrounding dark matter halos are correlated (commonly referred to as "assembly bias"), 2) gas is heated over large scales at early times, leading to coherent modulation of cooling and star formation… ▽ More Two explanations have been put forward to explain the observed conformity between the colours and specific star formation rates (SFR/$M_*$) of galaxies on large scales: 1) the formation times of their surrounding dark matter halos are correlated (commonly referred to as "assembly bias"), 2) gas is heated over large scales at early times, leading to coherent modulation of cooling and star formation between well-separated galaxies (commonly referred to as "pre-heating") . To distinguish between the pre-heating and assembly bias scenarios, we search for relics of energetic feedback events in the neighbourhood of central galaxies with different specific star formation rates. We find a significant excess of very high mass ($\log M_* > 11.3$) galaxies out to a distance of 2.5 Mpc around low SFR/$M_*$ central galaxies compared to control samples of higher SFR/$M_*$ central galaxies with the same stellar mass and redshift. We also find that very massive galaxies in the neighbourhood of low SFR/$M_*$ galaxies have much higher probability of hosting radio loud active galactic nuclei. The radio-loud AGN fraction in neighbours with $\log M_* > 11.3$ is four times higher around passive, non star-forming centrals at projected distances of 1 Mpc and two times higher at projected distances of 4 Mpc. Finally, we carry out an investigation of conformity effects in the recently publicly-released Illustris cosmological hydrodynamical simulation, which includes energetic input both from quasars and from radio mode accretion onto black holes. We do not find conformity effects of comparable amplitude on large scales in the simulations and we propose that gas needs to be pushed out of dark matter halos more efficiently at high redshifts. △ Less

Submitted 10 August, 2015; originally announced August 2015.

Comments: 8 pages, 8 figures, submitted to MNRAS

Connection Between the Circumgalactic Medium and the Interstellar Medium of Galaxies: Results from the COS-GASS Survey

Authors: Sanchayeeta Borthakur, Timothy M. Heckman, Jason Tumlinson, Rongmon Bordoloi, Christopher Thom, Barbara Catinella, David Schiminovich, Romeel Dave, Guinevere Kauffmann, Sean M. Moran, Amelie Saintonge

Abstract: We present a study exploring the nature and properties of the Circum-Galactic Medium (CGM) and its connection to the atomic gas content in the interstellar medium (ISM) of galaxies as traced by the HI 21cm line. Our sample includes 45 low-z (0.026-0.049) galaxies from the GALEX Arecibo SDSS Survey. Their CGM was probed via absorption in the spectra of background Quasi-Stellar Objects at impact par… ▽ More We present a study exploring the nature and properties of the Circum-Galactic Medium (CGM) and its connection to the atomic gas content in the interstellar medium (ISM) of galaxies as traced by the HI 21cm line. Our sample includes 45 low-z (0.026-0.049) galaxies from the GALEX Arecibo SDSS Survey. Their CGM was probed via absorption in the spectra of background Quasi-Stellar Objects at impact parameters of 63 to 231kpc. The spectra were obtained with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. We detected neutral hydrogen (Ly$α$ absorption-lines) in the CGM of 92% of the galaxies. We find the radial profile of the CGM as traced by the Ly$α$ equivalent width can be fit as an exponential with a scale length of roughly the virial radius of the dark matter halo. We found no correlation between the orientation of sightline relative to the galaxy major axis and the Ly$α$ equivalent width. The velocity spread of the circumgalactic gas is consistent with that seen in the atomic gas in the interstellar medium. We find a strong correlation (99.8% confidence) between the gas fraction (M(HI)/M*) and the impact-parameter-corrected Ly$α$ equivalent width. This is stronger than the analogous correlation between corrected Ly$α$ equivalent width and SFR/M* (97.5% confidence). These results imply a physical connection between the HI disk and the CGM, which is on scales an order-of-magnitude larger. This is consistent with the picture in which the HI disk is nourished by accretion of gas from the CGM. △ Less

Submitted 15 September, 2015; v1 submitted 6 April, 2015; originally announced April 2015.

Comments: 13 pages, 9 figures, and 2 tables. Submitted to ApJ

A systematic study of the inner rotation curves of galaxies observed as part of the GASS and COLD GASS surveys

Authors: Guinevere Kauffmann, Mei-ling Huang, Sean Moran, Timothy Heckman

Abstract: We present a systematic analysis of the rotation curves of 187 galaxies with masses greater than 10^10 M_sol, with atomic gas masses from the GALEX Arecibo Sloan Survey (GASS), and with follow-up long-slit spectroscopy from the MMT. Our analysis focuses on stellar rotation curves derived by fitting stellar template spectra to the galaxy spectra binned along the slit. In this way, we are able to ob… ▽ More We present a systematic analysis of the rotation curves of 187 galaxies with masses greater than 10^10 M_sol, with atomic gas masses from the GALEX Arecibo Sloan Survey (GASS), and with follow-up long-slit spectroscopy from the MMT. Our analysis focuses on stellar rotation curves derived by fitting stellar template spectra to the galaxy spectra binned along the slit. In this way, we are able to obtain accurate rotation velocity measurements for a factor of 2 more galaxies than possible with the Halpha line. Galaxies with high atomic gas mass fractions are the most dark-matter dominated galaxies in our sample and have dark matter halo density profiles that are well fit by Navarro, Frenk & White profiles with an average concentration parameter of 10. The inner slopes and of the rotation curves correlate more strongly with stellar population age than with galaxy mass or structural parameters. At fixed stellar mass, the rotation curves of more actively star-forming galaxies have steeper inner slopes than less actively star-forming galaxies. The ratio between the galaxy specific angular momentum and the total specific angular momentum of its dark matter halo, R_j, correlates strongly with galaxy mass, structure and gas content. Low mass, disk-dominated galaxies with atomic gas mass fractions greater than 20% have median values of R_j of around 1, but massive, bulge-dominated galaxies have R_j=0.2-0.3. We argue that these trends can be understood in a picture where gas inflows triggered by disk instabilities lead to the formation of passive, bulge-dominated galaxies with low specific angular momentum. △ Less

Submitted 19 March, 2015; originally announced March 2015.

Comments: 10 pages, 12 figures, submitted to MNRAS

HI scaling relations of galaxies in the environment of HI-rich and control galaxies observed by the Bluedisk project

Authors: Enci Wang, Jing Wang, Guinevere Kauffmann, Gyula I. G. Jozsa, Cheng Li

Abstract: Our work is based on the "Bluedisk" project, a program to map the neutral gas in a sample of 25 HI-rich spirals and a similar number of control galaxies with the Westerbork Synthesis Radio Telescope (WSRT). In this paper we focus on the HI properties of the galaxies in the environment of our targeted galaxies. In total, we extract 65 galaxies from the WSRT cubes with stellar masses between… ▽ More Our work is based on the "Bluedisk" project, a program to map the neutral gas in a sample of 25 HI-rich spirals and a similar number of control galaxies with the Westerbork Synthesis Radio Telescope (WSRT). In this paper we focus on the HI properties of the galaxies in the environment of our targeted galaxies. In total, we extract 65 galaxies from the WSRT cubes with stellar masses between $10^8M_{\odot}$ and $10^{11}M_{\odot}$. Most of these galaxies are located on the same HI mass-size relation and "HI-plane" as normal spiral galaxies. We find that companions around HI-rich galaxies tend to be HI-rich as well and to have larger R90,HI/R50,HI. This suggests a scenario of "HI conformity", similar to the colour conformity found by Weinmann et al. (2006): galaxies tend to adopt the HI properties of their neighbours. We visually inspect the outliers from the HI mass-size relation and galaxies which are offset from the HI plane and find that they show morphological and kinematical signatures of recent interactions with their environment. We speculate that these outliers have been disturbed by tidal or ram-pressure stripping processes, or in a few cases, by accretion events. △ Less

Submitted 10 March, 2015; originally announced March 2015.

Comments: 16 pages, 12 figures; accepted for publication in MNRAS