Search

We study the stellar mass function (SMF) and the co-evolution with dark matter halos via abundance matching in the largest redshift range to date $0.25$, we find increased abundances of massive (log$\, M_{\star}/M_{\odot}>10.5$) implying integrated star formation efficiencies (SFE) $\epsilon_{\star}\equiv M_{\star}\, f_{\rm b}^{-1} M_{\rm halo}^{-1} \gtrsim 0.5$. We find a flattening of the SMF at the high-mass end that is better described by a double power law at $z>5.5$. At $z \lesssim 5.5$ it transitions to a Schechter law which coincides with the emergence of the first massive quiescent galaxies in the Universe. We trace the cosmic stellar mass density (SMD) and infer the star formation rate density (SFRD), which at $z>7.5$ agrees remarkably with recent \JWST{} UV luminosity function-derived estimates. However, at $z \lesssim 3.5$, we find significant tension ($\sim 0.3$ dex) with the cosmic star formation (SF) history from instantaneous SF measures, the causes of which remain poorly understood. We infer the stellar-to-halo mass relation (SHMR) and the SFE from abundance matching out to $z=12$, finding a non-monotonic evolution. The SFE has the characteristic strong dependence with mass in the range of $0.02 - 0.2$, and mildly decreases at the low mass end out to $z\sim3.5$. At $z\sim3.5$ the SFE increases sharply from $\sim 0.1$ to approach high SFE of $0.8-1$ by $z\sim 10$ for log$(M_{\rm h}/M_{\odot})\approx11.5$, albeit with large uncertainties. Finally, we use the SHMR to track the SFE and stellar mass growth throughout the halo history and find that they do not grow at the same rate -- from the earliest times up until $z\sim3.5$ the halo growth rate outpaces galaxy assembly, but at $z>3.5$ halo growth stagnates and accumulated gas reservoirs keep the SF going and galaxies outpace halos.

We study of the role of galaxy-galaxy interactions and disk instabilities in producing starburst activity in galaxies out to z=4. For this, we use a sample of 387 galaxies with robust total star formation rate measurements from Herschel, gas masses from ALMA, stellar masses and redshifts from multi-band photometry, and JWST/NIRCam rest-frame optical imaging. Using mass-controlled samples, we find an increased fraction of interacting galaxies in the starburst regime at all redshifts out to z=4. This increase correlates with star formation efficiency (SFE), but not with gas fraction. However, the correlation is weak (and only significant out to z=2), which could be explained by the short duration of SFE increase during interaction. In addition, we find that isolated disk galaxies make up a significant fraction of the starburst population. The fraction of such galaxies with star-forming clumps ("clumpy disks") is significantly increased compared to the main-sequence disk population. Furthermore, this fraction directly correlates with SFE. This is direct observational evidence for a long-term increase of SFE maintained due to disk instabilities, contributing to the majority of starburst galaxies in our sample and hence to substantial mass growth in these systems. This result could also be of importance for explaining the growth of the most massive galaxies at z>6., Comment: 14 pages, 5 figures. Submitted to ApJ. Comments welcome

Measurements of the size of dust continuum emission are an important tool for constraining the spatial extent of star formation and hence the build-up of stellar mass. Compact dust emission has generally been observed at Cosmic Noon (z~2-3). However, at earlier epochs, toward the end of the Reionization (z~4-6), only the sizes of a handful of IR-bright galaxies have been measured. In this work, we derive the dust emission sizes of main-sequence galaxies at z~5 from the ALPINE survey. We measure the dust effective radius r_e,FIR in the uv-plane in Band 7 of ALMA for seven ALPINE galaxies with resolved emission and we compare it with rest-frame UV and [CII]158$\mu$m measurements. We study the r_e,FIR-L_IR scaling relation by considering our dust size measurements and all the data in literature at z~4-6. Finally, we compare our size measurements with predictions from simulations. The dust emission in the selected ALPINE galaxies is rather extended (r_e,FIR~1.5-3 kpc), similar to [CII]158 um but a factor of ~2 larger than the rest-frame UV emission. Putting together all the measurements at z~5, spanning 2 decades in luminosity from L_IR ~ 10^11 L_sun to L_IR ~ 10^13 L_sun, the data highlight a steeply increasing trend of the r_e,FIR-L_IR relation at L_IR< 10^12 L_sun, followed by a downturn and a decreasing trend at brighter luminosities. Finally, simulations that extend up to the stellar masses of the ALPINE galaxies considered in the present work predict a sub-set of galaxies (~25% at 10^10 M_sun < M_star < 10^11 M_sun) with sizes as large as those measured., Comment: 8 pages, Accepted by A&A

Context: From redshift 6 to redshift $\approx$ 4 galaxies grow rapidly from low mass galaxies towards the more mature massive galaxies we see at the cosmic noon. Growth via gas accretion and mergers undoubtedly shape this evolution - however, there currently exists much uncertainty over the contribution of each of these processes to the overall evolution of galaxies. Furthermore, previous characterisations of the morphology of galaxies in the molecular gas phase has been limited by the coarse resolution of previous observations. Aims: The goal of this paper is to derive the morpho-kinematic properties of 3 main-sequence systems at $z\sim4.5$, drawn from the ALPINE survey, using brand new high-resolution ALMA data in band 7. The objects were previously characterised as one merger with three components, and and two dispersion-dominated galaxies. Methods: We use intensity and velocity maps, position-velocity diagrams and radial profiles of [CII], in combination with dust continuum maps, to analyse the morphology and kinematics of the 3 systems.} Results: In general, we find that the high-resolution ALMA data reveal more complex morpho-kinematic properties. We identify in one galaxy interaction-induced clumps, showing the profound effect that mergers have on the molecular gas in galaxies, consistent with what is suggested in recent simulations. A galaxy that was previously classified as dispersion dominated turned out to show two bright [CII] emission regions, that could either be merging galaxies or massive star-forming regions within the galaxy itself. The high resolution data for the other dispersion dominated object also revealed clumps of [CII] that were not previously identified. Within the sample, we might also detect star-formation powered outflows (or outflows from Active Galactic Nuclei) which appear to be fuelling diffuse gas regions and enriching the circumgalactic medium., Comment: Submitted to A&A, 13 pages

The Kennicutt-Schmidt (KS) relation between the gas and the star formation rate (SFR) surface density ($\Sigma_{\rm gas}$-$\Sigma_{\rm SFR}$) is essential to understand star formation processes in galaxies. So far, it has been measured up to z~2.5 in main-sequence galaxies. In this letter, we aim to put constraints at z~4.5 using a sample of four massive main-sequence galaxies observed by ALMA at high resolution. We obtained ~0.3"-resolution [CII] and continuum maps of our objects, which we then converted into gas and obscured SFR surface density maps. In addition, we produced unobscured SFR surface density maps by convolving Hubble ancillary data in the rest-frame UV. We then derived the average $\Sigma_{\rm SFR}$ in various $\Sigma_{\rm gas}$ bins, and estimated the uncertainties using a Monte Carlo sampling. Our galaxy sample follows the KS relation measured in main-sequence galaxies at lower redshift and is slightly lower than predictions from simulations. Our data points probe the high end both in terms of $\Sigma_{\rm gas}$ and $\Sigma_{\rm gas}$, and gas depletion timescales (285-843 Myr) remain similar to z~2 objects. However, three of our objects are clearly morphologically disturbed, and we could have expected shorter gas depletion timescales (~100 Myr) similar to merger-driven starbursts at lower redshifts. This suggests that the mechanisms triggering starbursts at high redshift may be different than in the low- and intermediate-z Universe., Comment: 9 pages, 7 figures, 2 tables, accepted by A&A (letter)

We investigate how the recovery of galaxy star formation rates (SFRs) using energy-balance spectral energy distribution (SED) fitting codes depends on their recent star formation histories (SFHs). We use the Magphys and Prospector codes to fit 6,706 synthetic spectral energy distributions of simulated massive galaxies at $1 < z < 8$ from the Feedback in Realistic Environments (FIRE) project. We identify a previously-unknown systematic error in the Magphys results due to bursty star formation: the derived SFRs can differ from the truth by as much as 1 dex, at large statistical significance ($>5\sigma$), depending on the details of their recent SFH. SFRs inferred using Prospector with non-parametric SFHs do not exhibit this trend. We show that using parametric SFHs (pSFHs) causes SFR uncertainties to be underestimated by a factor of up to $5\times$. Although this undoubtedly contributes to the significance of the systematic, it cannot explain the largest biases in the SFRs of the starbursting galaxies, which could be caused by details of the stochastic prior sampling or the burst implementation in the Magphys libraries. We advise against using pSFHs and urge careful consideration of starbursts when SED modelling galaxies where the SFR may have changed significantly over the last ~100 Myr, such as recently quenched galaxies, or those experiencing a burst. This concern is especially relevant, e.g. when fitting JWST observations of very high-redshift galaxies., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

Galaxies that are invisible in deep optical-NIR imaging but detected at longer wavelengths have been the focus of several recent observational studies, with speculation that they could constitute a substantial missing population and even dominate the cosmic star formation rate density at $z\gtrsim4$. The depths now achievable with JWST at the longest wavelengths probed by HST, coupled with the transformative resolution at longer wavelengths, are already enabling detailed, spatially-resolved characterisation of sources that were invisible to HST, often known as `HST-dark' galaxies. However, until now, there has been little theoretical work to compare against. We present the first simulation-based study of this population, using highly-resolved galaxies from the Feedback in Realistic Environments (FIRE) project, with multi-wavelength images along several lines of sight forward-modelled using radiative transfer. We naturally recover a population of modelled sources that meet commonly-used selection criteria ($H_{\rm{AB}}>27\,\rm{mag}$ and $H_{\rm{AB}}-\rm{F444W}>2.3$). These simulated HST-dark galaxies lie at high redshifts ($z=4-7$), have high levels of dust attenuation ($A_{V}=2-4$), and display compact recent star formation ($R_{1/2,\,\rm{4.4\,\mu\rm{m}}}\lesssim1\,\rm{kpc}$). Orientation is very important: for all but one of the 17 simulated galaxy snapshots with HST-dark sightlines, there exist other sightlines that do not meet the criteria. This result has important implications for comparisons between observations and models that do not resolve the detailed star-dust geometry, such as semi-analytic models or coarsely-resolved hydrodynamical simulations. Critically, we demonstrate that HST-dark sources are not an unexpected or exotic population, but a subset of high-redshift, highly-dust-attenuated sources viewed along certain lines of sight., Comment: 12 pages, 8 figures. Accepted for publication in ApJ

Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED fitting codes. Using JWST/NIRCam images, we also produce two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. The lens is a very massive and quiescent (sSFR < 10^(-13) yr-1) elliptical galaxy at z = 2.02 \pm 0.02 with a total mass Mtot(

Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially `decoupled'). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such galaxies' true properties. In this work, we test how well the SED fitting code Magphys can recover the stellar mass, star formation rate (SFR), specific SFR, dust mass and luminosity by fitting 6,706 synthetic SEDs generated from four zoom-in simulations of dusty, high-redshift galaxies from the FIRE project via dust continuum radiative transfer. Comparing our panchromatic results (using wavelengths 0.4-500$\mu$m, and spanning $1

Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multi-wavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ~1 kpc to ~0.1 kpc in <40 Myr when AGN winds are not included, but increases to ~2 kpc when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over 35 Myr, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations., Comment: 14 pages, 10 figures. Accepted for publication in MNRAS

How galaxies form, assemble, and cease their star-formation is a central question within the modern landscape of galaxy evolution studies. These processes are indelibly imprinted on the galaxy stellar mass function (SMF). We present constraints on the shape and evolution of the SMF, the quiescent galaxy fraction, and the cosmic stellar mass density across 90% of the history of the Universe from $z=7.5\rightarrow0.2$ via the COSMOS survey. Now with deeper and more homogeneous near-infrared coverage exploited by the COSMOS2020 catalog, we leverage the large 1.27 deg$^{2}$ effective area to improve sample statistics and understand cosmic variance particularly for rare, massive galaxies and push to higher redshifts with greater confidence and mass completeness than previous studies. We divide the total stellar mass function into star-forming and quiescent sub-samples through $NUVrJ$ color-color selection. Measurements are then fitted with Schechter functions to infer the intrinsic SMF, the evolution of its key parameters, and the cosmic stellar mass density out to $z=7.5$. We find a smooth, monotonic evolution in the galaxy SMF since $z=7.5$, in agreement with previous studies. The number density of star-forming systems seems to have undergone remarkably consistent growth spanning four decades in stellar mass from $z=7.5\rightarrow2$ whereupon high-mass systems become predominantly quiescent (i.e. downsizing). An excess of massive systems at $z\sim2.5-5.5$ with strikingly red colors, some newly identified, increase the observed number densities to the point where the SMF cannot be reconciled with a Schechter function. Systematics including cosmic variance and/or AGN contamination are unlikely to fully explain this excess, and so we speculate that there may be contributions from dust-obscured objects similar to those found in FIR surveys. (abridged), Comment: 39 pages, 24 figures, accepted for publication in A&A. Data files containing key measurements are available for download: https://doi.org/10.5281/zenodo.7808832

Recent years have seen growing interest in post-processing cosmological simulations with radiative transfer codes to predict observable fluxes for simulated galaxies. However, this can be slow, and requires a number of assumptions in cases where simulations do not resolve the ISM. Zoom-in simulations better resolve the detailed structure of the ISM and the geometry of stars and gas, however statistics are limited due to the computational cost of simulating even a single halo. In this paper, we make use of a set of high resolution, cosmological zoom-in simulations of massive M_star>10^10.5M_sol at z=2), star-forming galaxies from the FIRE suite. We run the SKIRT radiative transfer code on hundreds of snapshots in the redshift range 1.5

Recent observations indicate that galactic outflows are ubiquitous in high redshift galaxies, including normal star forming galaxies, quasar hosts, and dusty star forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star formation histories (SFH) and galactic outflow properties of galaxies in massive haloes ($10^{12}M_{\odot}4$ DSFGs, with SFRs of $\sim 1000\ M_{\odot}\rm yr^{-1}$ and molecular gas masses of $M_{\rm mol}\sim 10^{10}\ M_{\odot}$. However, the simulated galaxies are characterised by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of $\sim 0.1$, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium (ISM). We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of $30\%$. Finally, compared to FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high redshift sample. We argue that both effects originate from the higher total and gas surface densities that characterise high$-z$ massive systems., Comment: 17 pages, 12 figures. Comments welcome

Single flux density measurements at observed-frame sub-millimeter and millimeter wavelengths are commonly used to probe dust and gas masses in galaxies. In this Letter, we explore the robustness of this method to infer dust mass, focusing on quiescent galaxies, using a series of controlled experiments on four massive haloes from the Feedback in Realistic Environments (FIRE) project. Our starting point is four star-forming, central galaxies at seven redshifts between z=1.5 and z=4.5. We generate modified quiescent galaxies that have been quenched for 100Myr, 500Myr, or 1Gyr prior to each of the studied redshifts by re-assigning stellar ages. We derive spectral energy distributions for each fiducial and modified galaxy using radiative transfer. We demonstrate that the dust mass inferred is highly dependent on the assumed dust temperature, T_dust, which is often unconstrained observationally. Motivated by recent work on quiescent galaxies that assumed T_dust~25K, we show that the ratio between dust mass and 1.3mm flux density can be higher than inferred by up to an order of magnitude, due to the considerably lower dust temperatures seen in non star-forming galaxies. This can lead to an underestimation of dust mass (and, when sub-mm flux density is used as a proxy for molecular gas content, gas mass). This underestimation is most severe at higher redshifts, where the observed-frame 1.3mm flux density probes rest-frame wavelengths far from the Rayleigh-Jeans regime, and hence depends super-linearly on dust temperature. We fit relations between ratios of rest-frame far-infrared flux densities and mass-weighted dust temperature that can be used to constrain dust temperatures from observations and hence derive more reliable dust and molecular gas masses., Comment: 10 pages, 4 figures; accepted for publication in ApJL

We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope SZ survey (SPT-SZ) and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array (ALMA). This work builds on the previous average rest-frame spectrum, given in Spilker et al. (2014) for the first 22 sources, and is comprised of a total of 78 sources, normalized by their respective apparent dust masses. The spectrum spans $1.9$$<$z$<$$6.9$ and covers rest-frame frequencies of 240$-$800 GHz. Combining this data with low-J CO observations from the Australia Telescope Compact Array (ATCA), we detect multiple bright line features from $^{12}$CO, $[$CI$]$, and H$_2$O, as well as fainter molecular transitions from $^{13}$CO, HCN, HCO$^+$, HNC, CN, H$_2$O$^+$, and CH. We use these detections, along with limits from other molecules, to characterize the typical properties of the interstellar medium (ISM) for these high redshift DSFGs. We are able to divide the large sample into subsets in order to explore how the average spectrum changes with various galaxy properties, such as effective dust temperature. We find that systems with hotter dust temperatures exhibit differences in the bright $^{12}$CO emission lines, and contain either warmer and more excited dense gas tracers, or larger dense gas reservoirs. These observations will serve as a reference point to studies of the ISM in distant luminous DSFGs (L$_{\mathrm{IR}}$$>$$10^{12}$L$_\odot$), and will inform studies of chemical evolution before the peak epoch of star formation at $z=2-3$., Comment: 26 pages, 10 figures, accepted for publication in ApJ

High-redshift dusty star-forming galaxies with very high star formation rates (500 -- 3000 M$_{\odot}$ yr$^{-1}$) are key to understanding the formation of the most extreme galaxies in the early Universe. Characterising the gas reservoir of these systems can reveal the driving factor behind the high star formation. Using molecular gas tracers like high-J CO lines, neutral carbon lines and the dust continuum, we can estimate the gas density and radiation field intensity in their interstellar media. In this paper, we present high resolution ($\sim$0.4$^{\prime\prime}$) observations of CO(7-6), [CI](2-1) and dust continuum of 3 lensed galaxies from the SPT-SMG sample at $z\sim$ 3 with the Atacama Large Millimeter/submillimeter Array. Our sources have high intrinsic star-formation rates ($>$850 M$_{\odot}$yr$^{-1}$) and rather short depletion timescales ($<$100 Myr). Based on the L$_{[\rm CI](2-1)}$/L$_{\rm CO(7-6)}$ and L$_{[\rm CI](2-1)}$/L$_{\rm IR}$ ratios, our galaxy sample has similar radiation field intensities and gas densities compared to other submillimetre galaxies. We perform visibility-based lens modelling on these objects to reconstruct the kinematics in the source plane. We find that the cold gas masses of the sources are compatible with simple dynamical mass estimates using ULIRG-like values of the CO-H$_2$ conversion factor $\alpha_{\rm CO}$ but not Milky Way-like values. We find diverse source kinematics in our sample: SPT0103-45 and SPT2147-50 are likely rotating disks while SPT2357-51 is possibly a major merger. The analysis presented in the paper could be extended to a larger sample to determine better statistics of morphologies and interstellar medium properties of high-$z$ dusty star-forming galaxies., Comment: 29 pages, 26 figures, Accepted for publication in A&A. Figure 16 and the related conclusions have been revised

We present Gemini-S and {\it Spitzer}-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this proto-cluster core. Three of the four LBGs are new systems, while one appears to be a counterpart of one of the nine observed SMGs. We identify a candidate brightest cluster galaxy (BCG) with a stellar mass of $(3.2^{+2.5}_{-1.4})\times10^{11}\,{\rm M}_{\odot}$. The stellar masses of the eight other SMGs place them on, above, and below the main sequence of star formation at $z\approx4.5$. The cumulative stellar mass for the SPT2349-56 core is at least $(11.5\pm2.9)\times10^{11}\,{\rm M}_{\odot}$, a sizeable fraction of the stellar mass in local BCGs, and close to the universal baryon fraction (0.16) relative to the virial mass of the core ($10^{13}\,{\rm M}_{\odot}$). As all 14 of these SMGs are destined to quickly merge, we conclude that the proto-cluster core has already developed a significant stellar mass at this early stage, comparable to $z=1$ BCGs. Importantly, we also find that the SPT2349-56 core structure would be difficult to uncover in optical surveys, with none of the ALMA sources being easily identifiable or constrained through $g,r,$ and $i$ colour-selection in deep optical surveys and only a modest overdensity of LBGs over the extended core structure. SPT2349-56 therefore represents a truly dust-obscured phase of a massive cluster core under formation.

The South Pole Telescope (SPT) has systematically identified 81 high-redshift, strongly gravitationally lensed, dusty star-forming galaxies (DSFGs) in a 2500 square degree cosmological mm-wave survey. We present the final spectroscopic redshift survey of this flux-limited ($S_{870\, \mathrm{\mu m}} > 25\, \mathrm{mJy}$) sample, initially selected at $1.4$ mm. The redshift survey was conducted with the Atacama Large Millimeter/submillimeter Array across the $3$ mm spectral window, targeting carbon monoxide line emission. By combining these measurements with ancillary data, the SPT sample is now spectroscopically complete, with redshifts spanning $1.9$$<$$z$$<$$6.9$ and a median of $z=3.9 \pm 0.2$. We present the mm through far-infrared photometry and spectral energy density fits for all sources, along with their inferred intrinsic properties. Comparing the properties of the SPT sources to the unlensed DSFG population, we demonstrate that the SPT-selected DSFGs represent the most extreme infrared-luminous galaxies, even after accounting for strong gravitational lensing. The SPT sources have a median star formation rate of $2.3(2)\times 10^3\, \mathrm{M_\odot yr^{-1}}$ and a median dust mass of $1.4(1)\times10^9\, \mathrm{M_\odot}$. However, the inferred gas depletion timescales of the SPT sources are comparable to those of unlensed DSFGs, once redshift is taken into account. This SPT sample contains roughly half of the known spectroscopically confirmed DSFGs at $z$$>$$5$, making this the largest sample of high-redshift DSFGs to-date, and enabling the "high-redshift tail" of extremely luminous DSFGs to be measured. Though galaxy formation models struggle to account for the SPT redshift distribution, the larger sample statistics from this complete and well-defined survey will help inform future theoretical efforts., Comment: 56 pages, 25 figures, accepted for publication in ApJ

We present the first detailed study of the spatially-resolved dust continuum emission of simulated galaxies at 1

We present the implementation and the first results of cosmic ray (CR) feedback in the Feedback In Realistic Environments (FIRE) simulations. We investigate CR feedback in non-cosmological simulations of dwarf, sub-$L\star$ starburst, and $L\star$ galaxies with different propagation models, including advection, isotropic and anisotropic diffusion, and streaming along field lines with different transport coefficients. We simulate CR diffusion and streaming simultaneously in galaxies with high resolution, using a two moment method. We forward-model and compare to observations of $\gamma$-ray emission from nearby and starburst galaxies. We reproduce the $\gamma$-ray observations of dwarf and $L\star$ galaxies with constant isotropic diffusion coefficient $\kappa \sim 3\times 10^{29}\,{\rm cm^{2}\,s^{-1}}$. Advection-only and streaming-only models produce order-of-magnitude too large $\gamma$-ray luminosities in dwarf and $L\star$ galaxies. We show that in models that match the $\gamma$-ray observations, most CRs escape low-gas-density galaxies (e.g.\ dwarfs) before significant collisional losses, while starburst galaxies are CR proton calorimeters. While adiabatic losses can be significant, they occur only after CRs escape galaxies, so they are only of secondary importance for $\gamma$-ray emissivities. Models where CRs are ``trapped'' in the star-forming disk have lower star formation efficiency, but these models are ruled out by $\gamma$-ray observations. For models with constant $\kappa$ that match the $\gamma$-ray observations, CRs form extended halos with scale heights of several kpc to several tens of kpc., Comment: 31 pages, 26 figures, accepted for publication in MNRAS

The origin of the high SFR observed in high-z dusty star-forming galaxies is still unknown. Large fractions of dense molecular gas might provide part of the explanation, but there are few observational constraints on the amount of dense gas in high-z systems dominated by star formation. We present the results of our ALMA program targeting dense-gas tracers (HCN(5-4), HCO+(5-4), and HNC(5-4)) in 5 strongly lensed galaxies from the SPT SMG sample. We detected two of these lines (SNR>5) in SPT-125-47 at z=2.51 and tentatively detected all three (SNR~3) in SPT0551-50 at z=3.16. Since a significant fraction of our target lines is not detected, we developed a statistical method to derive unbiased mean properties taking into account both detections and non-detections. On average, the HCN(5-4) and HCO+(5-4) luminosities of our sources are a factor of ~1.7 fainter than expected, based on the local L'HCN(5-4)-LIR relation, but this offset corresponds to only ~2 sigma. We find that both the HCO+/HCN and HNC/HCN flux ratios are compatible with unity. The first ratio is expected for PDRs while the second is consistent with PDRs or XDRs and/or mid-IR pumping of HNC. Our sources are at the high end of the local relation between the star formation efficiency, determined using the LIR/[CI] and LIR/CO ratios, and the dense gas fraction, estimated using the HCN/[CI] and HCN/CO ratios. In SPT0125-47, we found that the velocity profiles of the lines tracing dense (HCN, HCO+) and lower-density (CO, [CI]) gas are similar. In addition to these lines, we obtained one robust and one tentative detection of 13CO(4-3) and found an average I12CO(4-3)/I13CO(4-3) flux ratio of 26.1$_{-3.5}^{+4.5}$, indicating a young but not pristine interstellar medium. We argue that the combination of large and slightly enriched gas reservoirs and high dense-gas fractions could explain the prodigious star formation in these systems., Comment: 19 pages, 12 figures, 6 tables, Accepted by A&A

Galaxies grow inefficiently, with only a few percent of the available gas converted into stars each free-fall time. Feedback processes, such as outflowing winds driven by radiation pressure, supernovae or supermassive black hole accretion, can act to halt star formation if they heat or expel the gas supply. We report a molecular outflow launched from a dust-rich star-forming galaxy at redshift 5.3, one billion years after the Big Bang. The outflow reaches velocities up to 800 km/s relative to the galaxy, is resolved into multiple clumps, and carries mass at a rate within a factor of two of the star formation rate. Our results show that molecular outflows can remove a large fraction of the gas available for star formation from galaxies at high redshift., Comment: Author's version of paper to appear in Science on 7 September 2018

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today., Comment: To appear in April 26 issue of Nature

According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch., Comment: Nature, published online on 06 December 2017

Sub-milllimetre galaxies (SMGs) are some of the most luminous star-forming galaxies in the Universe, however their properties remain hard to determine due to the difficulty of identifying their optical\slash near-infrared counterparts. One of the key steps to determining the nature of SMGs is measuring a redshift distribution representative of the whole population. We do this by applying statistical techniques to a sample of 761 850$\mu$m sources from the SCUBA-2 Cosmology Legacy Survey observations of the UKIDSS Ultra-Deep Survey (UDS) Field. We detect excess galaxies around $> 98.4$ per cent of the 850$\mu$m positions in the deep UDS catalogue, giving us the first 850$\mu$m selected sample to have virtually complete optical\slash near-infrared redshift information. Under the reasonable assumption that the redshifts of the excess galaxies are representative of the SMGs themselves, we derive a median SMG redshift of $z = 2.05 \pm 0.03$, with 68 per cent of SMGs residing between $1.07 < z < 3.06$. We find an average of $1.52\pm 0.09$ excess $K$-band galaxies within 12 arc sec of an 850$\mu$m position, with an average stellar mass of $2.2\pm 0.1 \times 10^{10}$ M$_\odot$. While the vast majority of excess galaxies are star-forming, $8.0 \pm 2.1$ per cent have passive rest-frame colours, and are therefore unlikely to be detected at sub-millimetre wavelengths even in deep interferometry. We show that brighter SMGs lie at higher redshifts, and use our SMG redshift distribution -- along with the assumption of a universal far-infrared SED -- to estimate that SMGs contribute around 30 per cent of the cosmic star formation rate density between $0.5 < z < 5.0$., Comment: 10 pages, 8 figures, accepted for publication in MNRAS

We report the discovery and constrain the physical conditions of the interstellar medium of the highest-redshift millimeter-selected dusty star-forming galaxy (DSFG) to date, SPT-S J031132-5823.4 (hereafter SPT0311-58), at $z=6.900 +/- 0.002$. SPT0311-58 was discovered via its 1.4mm thermal dust continuum emission in the South Pole Telescope (SPT)-SZ survey. The spectroscopic redshift was determined through an ALMA 3mm frequency scan that detected CO(6-5), CO(7-6) and [CI](2-1), and subsequently confirmed by detections of CO(3-2) with ATCA and [CII] with APEX. We constrain the properties of the ISM in SPT0311-58 with a radiative transfer analysis of the dust continuum photometry and the CO and [CI] line emission. This allows us to determine the gas content without ad hoc assumptions about gas mass scaling factors. SPT0311-58 is extremely massive, with an intrinsic gas mass of $M_{\rm gas} = 3.3 \pm 1.9 \times10^{11}\,M_{\odot}$. Its large mass and intense star formation is very rare for a source well into the Epoch of Reionization., Comment: Accepted by ApJL, Fixed compilation problem in v.2

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge., Comment: 11 pages, 4 figures, accepted for publication in ApJ

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts ( = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4., Comment: 26 pages, 15 figures, updated to match version accepted for publication in MNRAS

The Herschel Multi-tiered Extragalactic Survey (HerMES) has identified large numbers of dusty star-forming galaxies (DSFGs) over a wide range in redshift. A detailed understanding of these DSFGs is hampered by the limited spatial resolution of Herschel. We present 870um 0.45" resolution imaging from the Atacama Large Millimeter/submillimeter Array (ALMA) of 29 HerMES DSFGs with far-infrared (FIR) flux densities in between the brightest of sources found by Herschel and fainter DSFGs found in ground-based sub-millimeter (sub-mm) surveys. We identify 62 sources down to the 5-sigma point-source sensitivity limit in our ALMA sample (sigma~0.2mJy), of which 6 are strongly lensed (showing multiple images) and 36 experience significant amplification (mu>1.1). To characterize the properties of the ALMA sources, we introduce and make use of uvmcmcfit, a publicly available Markov chain Monte Carlo analysis tool for interferometric observations of lensed galaxies. Our lens models tentatively favor intrinsic number counts for DSFGs with a steep fall off above 8mJy at 880um. Nearly 70% of the Herschel sources comprise multiple ALMA counterparts, consistent with previous research indicating that the multiplicity rate is high in bright sub-mm sources. Our ALMA sources are located significantly closer to each other than expected based on results from theoretical models as well as fainter DSFGs identified in the LABOCA ECDFS Submillimeter Survey. The high multiplicity rate and low projected separations argue in favor of interactions and mergers driving the prodigious emission from the brightest DSFGs as well as the sharp downturn above S_880=8mJy., Comment: Accepted to ApJ. 26 pages, 8 figures, 5 tables

We present a merger-driven evolutionary model for the production of luminous (LIRGs) and ultraluminous infrared galaxies (ULIRGs) with warm IR colours. Our results show that simulations of gas--rich major mergers including star formation, black hole growth, and feedback can produce warm (U)LIRGs. We also find that while the warm evolutionary phase is associated with increased AGN activity, star formation alone may be sufficient to produce warm IR colours. However, the transition can be suppressed entirely - even when there is a significant AGN contribution - when we assume a single-phase ISM, which maximizes the attenuation. Finally, our evolutionary models are consistent with the 25-to-60 micron flux density ratio vs. L_HX/L_IR relation for local LIRGs and ULIRGs, and predict the observed scatter in IR color at fixed L_HX/L_IR. Therefore, our models suggest a cautionary note in the interpretation of warm IR colours: while associated with periods of active black hole growth, they are probably produced by a complex mix of star formation and AGN activity intermediate between the cold star formation dominated phase and the birth of a bright, unobscured quasar., Comment: 5 pages, accepted for publication in MNRAS Letters

We have searched for Type Ia supernovae (SNe Ia) in the local (d < 60 Mpc) Universe using Northern Sky Variability Survey (NSVS) data collected from the nightly optical surveys of the Robotic Optical Transient Search Experiment} (ROTSE) Telescope. It was hoped that SNe Ia would provide a means to find previously-unknown low surface brightness (LSB) galaxies or displaced stars that would otherwise be very difficult to detect. The ROTSE data allowed us to survey 19,000 square degrees at declinations north of 0 degrees, but we did not find a single SN Ia in a period of time covering roughly one year. Using known SNe Ia rates in bright galaxies, we set an upper limit on the optical luminosity density, L_B, of LSBs in the local Universe. Using mean LSB baryonic and dynamical mass-to-light ratios, we find 95% upper limits for LSBs of L_B \le 2.53 x 10^8 L_{B, solar} Mpc^{-3}, Omega_b \le 0.0040, and Omega_m \le 0.036. We conclude that LSBs and displaced stars are not a major constituent of matter in the local Universe., Comment: 17 pages, 4 figures, to be published in The Astrophysical Journal