Search

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 $099.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 $, Comment: 27 pages, 10 figures. Accepted for publication in the Astrophysical Journal on March 23, 2022

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors., Contribution to Snowmass 2021. Note some text overlap with CMB-HD Astro2020 APC and RFI (arXiv:1906.10134, arXiv:2002.12714). Science case further broadened and updated

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}$-$λ_{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 $λ_{peak}=λ_{t}({L_{IR}}/{L_{t}})^η$ where $η=-0.09\pm0.01$, L$_{t}=10^{12}$ L$_{\odot}$, and $λ_{t}=92\pm2μ$m. We find no evidence for redshift evolution of galaxies' temperatures, or $λ_{peak}$, at fixed L$_{IR}$ from $099.99% confidence. Our finding does not preclude evolution in dust temperatures at fixed stellar mass, which is expected from a non-evolving L$_{IR}$-$λ_{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}$-$λ_{peak}$ evolution toward higher redshift ($z\sim5-6$) requires better sampling of galaxies' dust SEDs near their peaks (observed $\sim$200-600$μ$m) with $, 27 pages, 10 figures. Accepted for publication in the Astrophysical Journal on March 23, 2022

We present Ly$��$ and ultraviolet-continuum (UV) luminosity functions (LFs) of galaxies and active galactic nuclei (AGN) at $z=2.0-3.5$ determined by the un-targetted optical spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in $11.4$ deg$^2$ of fiber-spectra sky coverage, obtaining $18320$ galaxies spectroscopically identified with Ly$��$ emission, $2126$ of which host type 1 AGN showing broad (FWHM$~>1000$ km s$^{-1}$) Ly$��$ emission lines. We derive the Ly$��$ (UV) LF over 2 orders of magnitude covering bright galaxies and AGN in $\log L_\mathrm{Ly��}/\mathrm{[erg~s^{-1}]}=43.3-45.5$ ($-27-21$) type 1 AGN increases from $z\sim2$ to $z\sim0$ as opposed to the evolution of bright ($M_\mathrm{UV}, 27 pages, 19 figures, accepted for publication in ApJ

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$ $��$Jy and find their photometric redshifts span $07 ��$ in the density map and we confirm the photometric redshifts are consistent with previously measured spectra from a few galaxies at the cluster centers., 14 pages, 13 figures, accepted for publication in MNRAS on March 3 2020

We present cosmological ($30-400$ Mpc) distributions of neutral hydrogen (HI) in the inter-galactic medium (IGM) traced by Ly$��$ Emitters (LAEs) and QSOs at $z=2.1-2.5$, selected with the data of the on-going Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the eBOSS survey. We investigate spatial correlations of LAEs and QSOs with HI tomography maps reconstructed from HI Ly$��$ forest absorption in the spectra of background galaxies and QSOs obtained by the CLAMATO survey and this study, respectively. In the cosmological volume far from QSOs, we find that LAEs reside in regions of strong HI absorption, i.e. HI rich, which is consistent with results of previous galaxy-background QSO pair studies. Moreover, there is an anisotropy in the HI-distribution plot of transverse and line-of-sight distances; on average the HI absorption peak is blueshifted by $\sim 200$ km s$^{-1}$ from the LAE Ly$��$ redshift, reproducing the known average velocity offset between the Ly$��$ emission redshift and the galaxy systemic redshift. We have identified a $\sim$ 40-Mpc scale volume of HI underdensity that is a candidate for a giant HII bubble, where six QSOs and an LAE overdensity exist at $\left < z \right > =2.16$. The coincidence of the QSO and LAE overdensities with the HI underdensity indicates that the ionizing photon radiation of the QSOs has created a highly ionized volume of multiple proximity zones in a matter overdensity. Our results suggest an evolutionary picture where HI gas in an overdensity of galaxies becomes highly photoionized when QSOs emerge in the galaxies., 26 pages, 19 figures. Accepted for Publication in ApJ

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), Comment: Concept Paper submitted to Astro2020

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), Concept Paper submitted to Astro2020

We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of halos up to M_halo~10^12 M_sun (M_star~10^10.5 M_sun) at z>=5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust extinction, dust emission, and dust temperature within these simulated z>=5 galaxies. Our sample forms a tight correlation between infrared excess (IRX=F_IR/F_UV) and ultraviolet (UV)-continuum slope (beta_UV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX-beta_UV relation is mainly determined by the shape of the extinction curve and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (L_IR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given L_IR, the peak wavelength of the dust spectral energy distributions for z>=5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z=0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (LFs) at z=5-10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IRLFs and UV luminosity densities at z=5-10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z>=5 galaxies., 21 pages, 17 figures, key results shown in Figs. 5, 8, 10, 11, MNRAS accepted. Data products including mock galaxy SEDs and images are available upon request

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure., Comment: Science White paper submitted to the Astro2020 Decadal Survey

We present near-infrared integral-field spectroscopic observations targeting H$\alpha$ in eight sub-millimeter galaxies (SMGs) at $z$=1.3-2.5 using VLT/SINFONI, obtaining significant detections for six of them. The star formation rates derived from the H$\alpha$ emission are $\sim$100 M$_\odot$yr$^{-1}$, which account for only $\sim$ 20-30\% of the infrared-derived values, thus suggesting that these systems are very dusty. Two of these systems present [NII]/H$\alpha$ ratios indicative of the presence of an Active Galactic Nucleus (AGN). We mapped the spatial distribution and kinematics of the star forming regions in these galaxies on kpc-scales. In general, the H$\alpha$ morphologies tend to be highly irregular and/or clumpy, showing spatial extents of $\sim$3-11~kpc. We find evidence for significant spatial offsets, of $\sim$0.1-0.4$"$ or 1.2-3.4 kpc, between the H$\alpha$ and the continuum emission in three of the sources. Performing a kinemetry analysis we conclude that the majority of the sample is not consistent with disk-like rotation-dominated kinematics. Instead, they tend to show irregular and/or clumpy and turbulent velocity and velocity dispersion fields. This can be interpreted as evidence for scenario in which these extreme star formation episodes are triggered by galaxy-galaxy interactions and major mergers. In contrast to recent results for SMGs, these sources appear to follow the same relations between gas and star forming rate densities as less luminous and/or normal star forming galaxies., Comment: Astrophysical Journal, accepted. 17 page in emulateapj format, 6 figures

We present the results from a large 850 micron survey of the sigma Orionis cluster using the SCUBA-2 camera on the James Clerk Maxwell Telescope. The 0.5-degree diameter circular region we surveyed contains 297 young stellar objects with an age estimated at about 3Myr. We detect 9 of these objects, 8 of which have infrared excesses from an inner disc. We also serendipitously detect 3 non-stellar sources at > 5sigma that are likely background submillimetre galaxies. The 9 detected stars have inferred disc masses ranging from 5 to about 17MJup, assuming similar dust properties as Taurus discs and an ISM gas-to-dust ratio of 100. There is a net positive signal toward the positions of the individually undetected infrared excess sources indicating a mean disc mass of 0.5 MJup . Stacking the emission toward those stars without infrared excesses constrains their mean disc mass to less than 0.3MJup, or an equivalent Earth mass in dust. The submillimetre luminosity distribution is significantly different from that in the younger Taurus region, indicating disc mass evolution as star forming regions age and the infrared excess fraction decreases. Submillimeter Array observations reveal CO emission toward 4 sources demonstrating that some, but probably not much, molecular gas remains in these relatively evolved discs. These observations provide new constraints on the dust and gas mass of protoplanetary discs during the giant planet building phase and provide a reference level for future studies of disc evolution., 10 pages, 6 figures, accepted for publication in MNRAS

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model., 14 pages, 11 figures, ApJ accepted

We present near infrared integral field spectroscopic observations targeting Ha in eight submillimeter galaxies (SMGs) at z = 1.3 2.5 using the Very Large Telescope/Spectrograph for Integral Field Observations in the Near Infrared obtaining significant detections for six of them. The star formation rates derived from the Ha emission are ˜100 M ? yr 1 which account for only ˜20 30 of the infrared derived values thus suggesting that these systems are very dusty. Two of these systems present [N ii]/Ha ratios indicative of the presence of an active galactic nucleus. We mapped the spatial distribution and kinematics of the star forming regions in these galaxies on kiloparsec scales. In general the Ha morphologies tend to be highly irregular and/or clumpy showing spatial extents of ˜3 11 kpc. We find evidence for significant spatial offsets of ˜0.?1 0.?4 or 1.2 3.4 kpc between the Ha and the continuum emission in three of the sources. Performing a kinemetry analysis we conclude that the majority of the sample is not consistent with disk like rotation dominated kinematics. Instead they tend to show irregular and/or clumpy and turbulent velocity and velocity dispersion fields. This can be interpreted as evidence for a scenario in which these extreme star formation episodes are triggered by galaxy galaxy interactions and major mergers. In contrast to recent results for SMGs these sources appear to follow the same relations between gas and star forming rate densities as less luminous and/or normal star forming galaxies.

Spectral Energy Distribution (SED) fitting in the far-infrared (FIR) is greatly limited by a dearth of data and an excess of free parameters - from galaxies' dust composition, temperature, mass, orientation, opacity, to heating from AGN. This paper presents a simple FIR SED fitting technique joining a modified, single dust temperature greybody, representing the reprocessed starburst emission in the whole galaxy, to a mid-infrared powerlaw, which approximates hot-dust emission from AGN heating or clumpy, hot starbursting regions. This FIR SED can be used to measure infrared luminosities, dust temperatures and dust masses for both local and high-z galaxies with 3 to 10+ FIR photometric measurements. This fitting method is compared to infrared template SEDs in the literature using photometric data on 65 local luminous and ultraluminous infrared galaxies, (U)LIRGs. Despite relying only on 2-4 free parameters, the coupled greybody/powerlaw SED fitting described here produces better fits to photometric measurements than best-fit literature template SEDs (with residuals a factor of ~2 lower). A mean emissivity index of beta=1.60+-0.38 and mid-infrared powerlaw slope of alpha=2.0+-0.5 is measured; the former agrees with the widely presumed emissivity index of beta=1.5 and the latter is indicative of an optically-thin dust medium with a shallow radial density profile, ~r^-0.5. Adopting characteristic dust temperature as the inverse wavelength where the SED peaks, dust temperatures ~25-45K are measured for local (U)LIRGs, ~5-15K colder than previous estimates using only simple greybodies. This comparative study highlights the impact of SED fitting assumptions on the measurement of physical properties such as infrared luminosity (and thereby infrared-based star formation rate), dust temperature and dust mass, for both local and high-redshift galaxies. [abridged], 12 pages, 7 figures, MNRAS accepted

We present the results of SCUBA2 observations at 450 micron and 850 micron of the field lensed by the massive cluster A370. With a total survey area > 100 arcmin2 and 1 sigma sensitivities of 3.92 and 0.82 mJy/beam at 450 and 850 micron respectively, we find a secure sample of 20 sources at 450 micron and 26 sources at 850 micron with a signal-to-noise ratio > 4. Using the latest lensing model of A370 and Monte Carlo simulations, we derive the number counts at both wavelengths. The 450 micron number counts probe a factor of four deeper than the counts recently obtained from the Herschel Space Telescope at similar wavelengths, and we estimate that ~47-61% of the 450 micron extragalactic background light (EBL) resolved into individual sources with 450 micron fluxes greater than 4.5 mJy. The faint 450 micron sources in the 4 sigma sample have positional accuracies of 3 arcseconds, while brighter sources (signal-to-noise > 6 sigma) are good to 1.4 arcseconds. Using the deep radio map (1 sigma ~ 6 uJy) we find that the percentage of submillimeter sources having secure radio counterparts is 85% for 450 micron sources with intrinsic fluxes > 6 mJy and 67% for 850 micron sources with intrinsic fluxes > 4 mJy. We also find that 67% of the > 4 sigma 450 micron sources are detected at 850 micron, while the recovery rate at 450 micron of > 4 sigma 850 micron sources is 54%. Combined with the source redshifts estimated using millimetric flux ratios, the recovered rate is consistent with the scenario where both 450 micron and 20 cm emission preferentially select lower redshift dusty sources, while 850 micron emission traces a higher fraction of dusty sources at higher redshifts. [Abridge], 22 pages, 14 figures, ApJ in press. Minor changes to the submitted version

The relationship between galaxy star formation rates (SFRs) and stellar masses (M ) is reexamined using a mass selected sample of 62000 star forming galaxies at z 1010 M ? have a much lower average specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST., Comment: 52 pages, 13 figures, NRAO Next Generation Very Large Array Memos Series: http://library.nrao.edu/ngvla.shtml