Some recent literature has claimed there to be an evolution in galaxies' dust temperatures towards warmer (or colder) spectral energy distributions (SEDs) between low and high redshift. These conclusions are driven by both theoretical models and empirical measurement. Such claims sometimes contradict one another and are prone to biases in samples or SED fitting techniques. What has made direct comparisons difficult is that there is no uniform approach to fitting galaxies' infrared/millimeter SEDs. Here we aim to standardize the measurement of galaxies' dust temperatures with a python-based SED fitting procedure, MCIRSED. We draw on reference datasets observed by IRAS, Herschel, and Scuba-2 to test for redshift evolution out to $z\sim2$. We anchor our work to the L$_{IR}$-$\lambda_{peak}$ plane, where there is an empirically observed anti-correlation between IR luminosity and rest-frame peak wavelength (an observational proxy for luminosity-weighted dust temperature) such that $\lambda_{peak}=\lambda_{t}({L_{IR}}/{L_{t}})^{\eta}$ where $\eta=-0.09\pm0.01$, L$_{t}=10^{12}$ L$_{\odot}$, and $\lambda_{t}=92\pm2\mu$m. We find no evidence for redshift evolution of galaxies' temperatures, or $\lambda_{peak}$, at fixed L$_{IR}$ from $0<z<2$ with >99.99% confidence. Our finding does not preclude evolution in dust temperatures at fixed stellar mass, which is expected from a non-evolving L$_{IR}$-$\lambda_{peak}$ relation and a strongly evolving SFR-M$_\star$ relation. The breadth of dust temperatures at a given L$_{IR}$ is likely driven by variation in galaxies' dust geometries and sizes and does not evolve. Testing for L$_{IR}$-$\lambda_{peak}$ evolution toward higher redshift ($z\sim5-6$) requires better sampling of galaxies' dust SEDs near their peaks (observed $\sim$200-600$\mu$m) with $<$1 mJy sensitivity. This poses a significant challenge to current instrumentation.
Sinclaire M. Manning, Caitlin M. Casey, Jorge A. Zavala, Georgios E. Magdis, Patrick M. Drew, Jaclyn B. Champagne, Manuel Aravena, Matthieu Béthermin, David L. Clements, Steven L. Finkelstein, Seiji Fujimoto, Christopher C. Hayward, Jacqueline A. Hodge, Olivier Ilbert, Jeyhan S. Kartaltepe, Kirsten K. Knudsen, Anton M. Koekemoer, Allison W. S. Man, David B. Sanders, Kartik Sheth, et al (5) The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3<z<4$ massive quiescent galaxies.
Caitlin M. Casey, Jorge A. Zavala, Sinclaire M. Manning, Manuel Aravena, Matthieu Béthermin, Karina I. Caputi, Jaclyn B. Champagne, David L. Clements, Patrick Drew, Steven L. Finkelstein, Seiji Fujimoto, Christopher C. Hayward, Anton M. Koekemoer, Vasily Kokorev, Claudia del P. Lagos, Arianna S. Long, Georgios E. Magdis, Allison W.S. Man, Ikki Mitsuhashi, Gergö Popping, et al (7) We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$\sigma$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $3<z<6$. The volume density of 2mm-selected DSFGs is consistent with predictions from some cosmological simulations and is similar to the volume density of their hypothesized descendants: massive, quiescent galaxies at $z>2$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langle\beta\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe's first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe's first dusty galaxies.
J. A. Zavala, C. M. Casey, S. M. Manning, M. Aravena, M. Bethermin, K. I. Caputi, D. L. Clements, E. da Cunha, P. Drew, S. L. Finkelstein, S. Fujimoto, C. Hayward, J. Hodge, J. S. Kartaltepe, K. Knudsen, A. M. Koekemoer, A. S. Long, G. E. Magdis, A. W. S. Man, G. Popping, et al (8) We present the first results from the 2mm Mapping Obscuration to Reionization (MORA) survey, the largest ALMA contiguous blank-field survey to-date with a total area of 184 sq. arcmin and the only at 2mm to search for dusty star-forming galaxies (DSFGs). We use the 13 sources detected above 5sigma to estimate the first ALMA galaxy number counts at this wavelength. These number counts are then combined with the state-of-the-art galaxy number counts at 1.2mm and 3mm and with a backward evolution model to place constraints on the evolution of the IR luminosity function and dust-obscured star formation in the last 13 billion years. Our results suggest a steep redshift evolution on the space density of DSFGs and confirm the flattening of the IR luminosity function at faint luminosities, with a slope of $\alpha_{LF} = -0.42^{+0.02}_{-0.04}$. We conclude that the dust-obscured component, which peaks at z=2-2.5, has dominated the cosmic history of star formation for the past ~12 billion years, back to z~4. At z=5, the dust-obscured star formation is estimated to be ~35% of the total star formation rate density and decreases to 25%-20% at z=6-7, implying a minor contribution of dust-enshrouded star formation in the first billion years of the Universe. With the dust-obscured star formation history constrained up to the end of the epoch of reionization, our results provide a benchmark to test galaxy formation models, to study the galaxy mass assembly history, and to understand the dust and metal enrichment of the Universe at early times.
Sinclaire M. Manning, Caitlin M. Casey, Chao-Ling Hung, Richard Battye, Michael L. Brown, Neal Jackson, Filipe Abdalla, Scott Chapman, Constantinos Demetroullas, Patrick Drew, Christopher A. Hales, Ian Harrison, Christopher J. Riseley, David B. Sanders, Robert A. Watson We present optical and near-infrared imaging covering a $\sim$1.53 deg$^2$ region in the Super-Cluster Assisted Shear Survey (SuperCLASS) field, which aims to make the first robust weak lensing measurement at radio wavelengths. We derive photometric redshifts for $\approx$176,000 sources down to $i^\prime_{\rm AB}\sim24$ and present photometric redshifts for 1.4 GHz $e$-MERLIN and VLA detected radio sources found in the central 0.26 deg$^{2}$. We compile an initial catalog of 149 radio sources brighter than S$_{1.4}>75$ $\mu$Jy and find their photometric redshifts span $0<z_{\rm phot}<4$ with radio luminosities between $10^{21}-10^{25}$ W Hz$^{-1}$, with medians of $\langle z \rangle =0.55$ and $\langle L_{1.4}\rangle =1.9\times10^{23}$ W Hz$^{-1}$ respectively. We find 95% of the \uJy radio source sample (141/149) have SEDs best fit by star-forming templates while 5% (8/149) are better fit by AGN. Spectral indices are calculated for sources with radio observations from VLA and GMRT at 325 MHz, with an average spectral slope of $\alpha=0.59\pm0.04$. Using the full photometric redshift catalog we construct a density map at the redshift of the known galaxy clusters, $z=0.20\pm0.08$. Four of the five clusters are prominently detected at $>7 \sigma$ in the density map and we confirm the photometric redshifts are consistent with previously measured spectra from a few galaxies at the cluster centers.
The main sequence of galaxies, a correlation between the star formation rates and stellar masses of galaxies, has been observed out to $z\sim4$. Galaxies within the scatter of the correlation are typically interpreted to be secularly evolving while galaxies with star formation rates elevated above the main sequence are interpreted to be undergoing interactions or to be Toomre-unstable disks with starbursting clumps. In this paper we investigate the recent merger histories of three dusty star forming galaxies, identified by their bright submillimeter emission at $z\sim1.5$. We analyze rest-frame optical and UV imaging, rest-frame optical emission line kinematics using slit spectra obtained with MOSFIRE on Keck I, and calculate Gini and M$_{20}$ statistics for each galaxy and conclude two are merger-driven while the third is an isolated disk galaxy. The disk galaxy lies $\sim$4$\times$ above the main sequence, one merger lies within the scatter of the main sequence, and one merger lies $\sim$4$\times$ below the main sequence. This hints that the location of a galaxy with respect to the main sequence may not be a useful discriminator of the recent star formation history of high-M$_{\star}$ galaxies at $z\sim1$.
Jorge A. Zavala, C. M. Casey, J. B. Champagne, Y. Chiang, H. Dannerbauer, P. Drew, H. Fu, J. Spilker, L. Spitler, K. V. Tran, E. Treister, S. Toft We present ALMA Band 6 (nu=233GHz, lambda=1.3mm) continuum observations towards 68 'normal' star-forming galaxies within two Coma-like progenitor structures at z=2.10 and 2.47, from which ISM masses are derived, providing the largest census of molecular gas mass in overdense environments at these redshifts. Our sample comprises galaxies with a stellar mass range of 1x10^9M_sun - 4x10^11M_sun with a mean M_*~6x10^10M_sun. Combining these measurements with multiwavelength observations and SED modeling, we characterize the gas mass fraction and the star formation efficiency, and infer the impact of the environment on galaxies' evolution. Most of our detected galaxies (~70%) have star formation efficiencies and gas fractions similar to those found for coeval field galaxies and in agreement with the field scaling relations. However, we do find that the proto-clusters contain an increased fraction of massive, gas-poor galaxies, with low gas fractions (f_gas<6-10%) and red rest-frame ultraviolet/optical colors typical of post-starburst and passive galaxies. The relatively high abundance of passive galaxies suggests an accelerated evolution of massive galaxies in proto-cluster environments. The large fraction of quenched galaxies in these overdense structures also implies that environmental quenching takes place during the early phases of cluster assembly, even before virialization. From our data, we derive a quenching efficiency of E_q~0.45 and an upper limit on the quenching timescale of T_q<1Gyr.
Caitlin M. Casey, Jorge A. Zavala, Manuel Aravena, Matthieu Bethermin, Karina I. Caputi, Jaclyn B. Champagne, David L. Clements, Elisabete Da Cunha, Patrick Drew, Steven L. Finkelstein, Christopher C. Hayward, Jeyhan S. Kartaltepe, Kirsten Knudsen, Anton M. Koekemoer, Georgios E. Magdis, Allison Man, Sinclaire M. Manning, Nick Z. Scoville, Kartik Sheth, Justin Spilker, et al (6) We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).
Caitlin M. Casey, Peter Capak, Johannes Staguhn, Lee Armus, Andrew Blain, Matthieu Bethermin, Jaclyn Champagne, Asantha Cooray, Kristen Coppin, Patrick Drew, Eli Dwek, Steven Finkelstein, Maximilien Franco, James Geach, Jacqueline Hodge, Maciej Koprowski, Claudia Lagos, Desika Narayanan, Alexandra Pope, David Sanders, et al (7) Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.
Observations of the baryon to dark matter fraction in galaxies through cosmic time are a fundamental test for galaxy formation models. Recent observational studies have suggested that some disk galaxies at $z>1$ host declining rotation curves, in contrast with observations of low redshift disk galaxies where stellar or HI rotation curves flatten at large radii. We present an observational counterexample, a galaxy named DSFG850.95 at $z=1.555$ (4.1 Gyr after the big bang) that hosts a flat rotation curve between radii of $\sim$6--14 kpc (1.2--2.8 disk scale lengths) and has a dark matter fraction of $0.44\pm0.08$ at the H-band half light radius, similar to the Milky Way. We create position-velocity and position-dispersion diagrams using Keck/MOSFIRE spectroscopic observations of H$\alpha$ and [NII] emission features, which reveal a flat rotation velocity of $V_{\rm flat}=285\pm12$ km/s and an ionized gas velocity dispersion of $\sigma_{0}=48\pm4$ km/s. This galaxy has a rotation-dominated velocity field with $V_{\rm flat}/\sigma_{0}\sim6$. Ground-based H-band imaging reveals a disk with Sérsic index of $1.29\pm0.03$, an edge-on inclination angle of $87\pm2^{\circ}$, and an H-band half light radius of $8.4\pm0.1$ kpc. Our results point to DSFG850.95 being a massive, rotationally-supported disk galaxy with a high dark-matter-to-baryon fraction in the outer galaxy, similar to disk galaxies at low redshift.
We carry out a blind search of 3mm continuum sources using the ALMA Science Archive to derive the first galaxy number counts at this wavelength. The analyzed data are drawn from observations towards three extragalactic legacy fields: COSMOS, CDF-S, and the UDS comprising more than 130 individual ALMA Band 3 pointings and an effective survey area of ~200 sq. arcmin with a continuum sensitivity that allows for the direct detection of unlensed Dusty Star-Forming Galaxies (DSFGs) dust emission beyond the epoch of reionization. We present a catalog of 16 sources detected at >5sigma with flux densities S_3mm = 60 - 600 uJy from which number counts are derived. These number counts are then used to place constraints on the volume density of DSFGs with an empirical backward evolution model. Our measured 3mm number counts indicate that the contribution of DSFGs to the cosmic star formation rate density at z >~ 4 is non-negligible. This is contrary to the generally adopted assumption of a sharply decreasing contribution of obscured galaxies at z > 4 as inferred by optical and near-infrared surveys. This work demonstrates the power of ALMA 3mm observations which can reach outstanding continuum sensitivities during typical spectral line science programs. Further constraints on 3mm-selected galaxies will be essential to refine models of galaxy formation and evolution as well as models of early Universe dust production mechanisms.
Deep, pencil-beam surveys from ALMA at 1.1-1.3mm have uncovered an apparent absence of high-redshift dusty galaxies, with existing redshift distributions peaking around $z\sim1.5-2.5$. This has led to a perceived dearth of dusty systems at $z>4$, and the conclusion, according to some models, that the early Universe was relatively dust-poor. In this paper, we extend the backward evolution galaxy model described by Casey et al. (2018) to the ALMA regime (in depth and area) and determine that the measured number counts and redshift distributions from ALMA deep field surveys are fully consistent with constraints of the infrared luminosity function (IRLF) at $z<2.5$ determined by single-dish submillimeter and millimeter surveys conducted on much larger angular scales ($\sim1-10$deg$^{2}$). We find that measured 1.1-1.3mm number counts are most constraining for the measurement of the faint-end slope of the IRLF at $z<2.5$ instead of the prevalence of dusty galaxies at $z>4$. Recent studies have suggested that UV-selected galaxies at $z>4$ may be particularly dust-poor, but we find their millimeter-wave emission cannot rule out consistency with the Calzetti dust attenuation law even by assuming relatively typical, cold-dust ($T_{\rm dust}\approx30\,$K) SEDs. Our models suggest that the design of ALMA deep fields requires substantial revision to constrain the prevalence of $z>4$ early Universe obscured starbursts. The most promising avenue for detection and characterization of such early dusty galaxies will come from future ALMA 2mm blank field surveys covering a few hundred arcmin$^{2}$ and the combination of existing and future dual-purpose 3mm datasets.
Though half of cosmic starlight is absorbed by dust and reradiated at long wavelengths (3$\mu$m-3mm), constraints on the infrared through millimeter galaxy luminosity function (the `IRLF') are poor in comparison to the rest-frame ultraviolet and optical galaxy luminosity function, particularly at z>2.5. Here we present a backward evolution model for interpreting number counts, redshift distributions, and cross-band flux density correlations in the infrared and submillimeter sky, from 70$\mu$m-2mm, using a model for the IRLF out to the epoch of reionization. Mock submillimeter maps are generated by injecting sources according to the prescribed IRLF and flux densities drawn from model spectral energy distributions that mirror the distribution of SEDs observed in $0<z<5$ dusty star-forming galaxies (DSFGs). We explore two extreme hypothetical case-studies: a dust-poor early Universe model, where DSFGs contribute negligibly ($<$10%) to the integrated star-formation rate density at $z>4$, and an alternate dust-rich early Universe model, where DSFGs dominate $\sim$90% of $z>4$ star-formation. We find that current submm/mm datasets do not clearly rule out either of these extreme models. We suggest that future surveys at 2mm will be crucial to measuring the IRLF beyond $z\sim4$. The model framework developed in this paper serves as a unique tool for the interpretation of multiwavelength IR/submm extragalactic datasets and will enable more refined constraints on the IRLF than can be made from direct measurements of individual galaxies' integrated dust emission.
The age and evolutionary status of MWC349A, the unique emission line star with maser and laser radiation in hydrogen recombination lines, remain unknown because the spectrum of the star is veiled by bright emission from the ionized disk and wind. The major argument for this massive ($>$10 M_sun) star being evolved is its association with a close-by (2.4 arcsec) companion, MWC349B, whose B0III spectrum implies an age of a few Myrs. However, newly obtained high-resolution spectra of MWC349B reveal a difference $\approx$ 35 km/s in the radial velocities of the two stars, which makes their being gravitationally bound highly improbable. An estimate of the relative proper motion of the two stars seems to confirm this conclusion. This reopens the previously suggested possibility that MWC349A is a young massive star in a region of active star formation close to the Cyg OB2 association. MWC349B, which moves with a speed $\geq$ 35 km/s relative to Cyg OB2, may be a runaway star from this association.