-
ODIN: Identifying Protoclusters and Cosmic Filaments Traced by Ly$α$-emitting Galaxies
Authors:
Vandana Ramakrishnan,
Kyoung-Soo Lee,
Maria Celeste Artale,
Eric Gawiser. Yujin Yang,
Changbom Park,
Robin Ciardullo,
Lucia Guaita,
Sang Hyeok Im,
Seongjae Kim,
Ankit Kumar,
Jaehyun Lee,
Seong-Kook Lee,
Byeongha Moon,
Nelson Padilla,
Alexandra Pope,
Roxana Popescu,
Hyunmi Song,
Paulina Troncoso,
Francisco Valdes,
Ann Zabludoff
Abstract:
To understand the formation and evolution of massive cosmic structures, studying them at high redshift, in the epoch when they formed the majority of their mass is essential. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey is undertaking the widest-area narrowband program to date, to use Ly$α$-emitting galaxies (LAEs) to trace the large-scale structure (LSS) of the Universe at t…
▽ More
To understand the formation and evolution of massive cosmic structures, studying them at high redshift, in the epoch when they formed the majority of their mass is essential. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey is undertaking the widest-area narrowband program to date, to use Ly$α$-emitting galaxies (LAEs) to trace the large-scale structure (LSS) of the Universe at three cosmic epochs. In this work, we present results at $z$ = 3.1 based on early ODIN data in the COSMOS field. We identify and characterize protoclusters and cosmic filaments using multiple methods and discuss their strengths and weaknesses. We then compare our observations against the IllustrisTNG suite of cosmological hydrodynamical simulations. The two are in excellent agreement, with a similar number and angular size of structures identified above a specified density threshold. We are able to recover the simulated protoclusters with $\log$(M$_{z=0}$/$M_\odot$) $\gtrsim$ 14.4 in $\sim$ 60\% of the cases. With these objects we show that the descendant masses of the protoclusters in our sample can be estimated purely based on our 2D measurements, finding a median $z$ = 0 mass of $\sim10^{14.5}$M$_\odot$. The lack of information on the radial extent of each protocluster introduces a $\sim$0.4~dex uncertainty in its descendant mass. Finally, we show that the recovery of the cosmic web in the vicinity of protoclusters is both efficient and accurate. The similarity of our observations and the simulations imply that our structure selection is likewise robust and efficient, demonstrating that LAEs are reliable tracers of the LSS.
△ Less
Submitted 12 June, 2024;
originally announced June 2024.
-
The New Worlds Simulations: Large-scale Simulations across Three Cosmologies
Authors:
Katrin Heitmann,
Thomas Uram,
Nicholas Frontiere,
Salman Habib,
Adrian Pope,
Silvio Rizzi,
Joe Insley
Abstract:
In this paper we describe the set of ``New Worlds Simulations'', three very large cosmology simulations, Qo'noS, Vulcan, and Ferenginar, that were carried out on the Summit supercomputer with the Hardware/Hybrid Cosmology Code, HACC. The gravity-only simulations follow the evolution of structure in the Universe by each employing 12,288^3 particles in (3 Gpc/h)^3 volumes, leading to a mass resoluti…
▽ More
In this paper we describe the set of ``New Worlds Simulations'', three very large cosmology simulations, Qo'noS, Vulcan, and Ferenginar, that were carried out on the Summit supercomputer with the Hardware/Hybrid Cosmology Code, HACC. The gravity-only simulations follow the evolution of structure in the Universe by each employing 12,288^3 particles in (3 Gpc/h)^3 volumes, leading to a mass resolution of m_p~10^9 Msun/h. The simulations cover three different cosmologies, one LambdaCDM model, consistent with measurements from Planck, one simulation with massive neutrinos, and one simulation with a varying dark energy equation of state. All simulations have the same phases to allow a detailed comparison of the results and the investigation of the impact of different cosmological parameters. We present measurements of some basic statistics, such as matter power spectra, correlation function, halo mass function and concentration-mass relation and investigate the differences due to the varying cosmologies. Given the large volume and high resolution, these simulations provide excellent bases for creating synthetic skies. A subset of the data is made publicly available as part of this paper.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Halfway to the Peak: The JWST MIRI 5.6 micron number counts and source population
Authors:
Leonid Sajkov,
Anna Sajina,
Alexandra Pope,
Stacey Alberts,
Lee Armus,
Duncan Farrah,
Jamie Lin,
Danilo Marchesini,
Jed McKinney,
Sylvain Veilleux,
Lin Yan,
Jason Young
Abstract:
We present an analysis of 8 JWST Mid-Infrared Instrument 5.6 micron images with 5sigma depths of ~0.1uJy. We detect 2854 sources within our combined area of 18.4 sq.arcmin -- a >4x increase in source density over earlier IRAC channel 3 data. We compute the MIRI 5.6um number counts including an analysis of the field-to-field variation. Relative to earlier published MIRI 5.6micron counts, our counts…
▽ More
We present an analysis of 8 JWST Mid-Infrared Instrument 5.6 micron images with 5sigma depths of ~0.1uJy. We detect 2854 sources within our combined area of 18.4 sq.arcmin -- a >4x increase in source density over earlier IRAC channel 3 data. We compute the MIRI 5.6um number counts including an analysis of the field-to-field variation. Relative to earlier published MIRI 5.6micron counts, our counts have a more pronounced knee, at roughly 2\,$μ$Jy. The location and amplitude of the counts at the knee are consistent with the Cowley et al. (2018) model predictions, although these models tend to overpredict the counts below the knee. In areas of overlap, 84% of the MIRI sources have a counterpart in the COSMOS2020 catalog. These MIRI sources have redshifts that are mostly in the z~0.5-2, with a tail out to z~5. They are predominantly moderate to low stellar masses $10^8-10^{10}$M$_{\odot}$) main sequence star-forming galaxies suggesting that with $\approx$2hr exposures, MIRI can reach well below $M^*$ at cosmic noon and reach higher mass systems out to z~5. Nearly 70% of the COSMOS2020 sources in areas of overlap now have a data point at 5.6micron (rest-frame near-IR at cosmic noon) which allows for more accurate stellar population parameter estimates. Finally, we discover 31 MIRI-bright sources not in COSMOS2020. A cross-match with IRAC channel 1 suggests that 10-20% of these are likely lower mass (M$_*\approx10^9$M$_{\odot}$), $z\sim1$ dusty galaxies. The rest (80--90%) are consistent with more massive, but still very dusty galaxies at z>3.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
First Constraints on the ISM Conditions of a Low Mass, Highly Obscured z=4.27 Main Sequence Galaxy
Authors:
Andrew Mizener,
Alexandra Pope,
Jed McKinney,
Patrick Kamieneski,
Katherine E. Whitaker,
Andrew Battisti,
Eric Murphy
Abstract:
We present the molecular gas content and ISM conditions of MACSJ0717 Az9, a strong gravitationally lensed $z=4.273$, $M_{*} \simeq 2\times10^9M_{\odot}$ star-forming galaxy with an unusually high ($\sim 80\%$) obscured star formation fraction. We detect CO(4-3) in two independent lensed images, as well as [N II]205$μ$m, with ALMA. We derive a molecular gas mass of log…
▽ More
We present the molecular gas content and ISM conditions of MACSJ0717 Az9, a strong gravitationally lensed $z=4.273$, $M_{*} \simeq 2\times10^9M_{\odot}$ star-forming galaxy with an unusually high ($\sim 80\%$) obscured star formation fraction. We detect CO(4-3) in two independent lensed images, as well as [N II]205$μ$m, with ALMA. We derive a molecular gas mass of log$_{10}[M_{H_{2}} (M_{\odot})] = 9.77$ making it moderately deficient in molecular gas compared to the lower redshift gas fraction scaling relation. Leveraging photodissociation region (PDR) models, we combine our CO(4-3) measurements with existing measurements of the [C II] 158$μ$m line and total infrared luminosity to model the PDR conditions. We find PDR conditions similar to local star-forming galaxies, with a mean hydrogen density log$_{10}$[$n_H$ $cm^{-3}$] = $4.80\pm0.39$ and a mean radiation field strength log$_{10}$[G$_0$ Habing] = $2.83\pm0.26$. Based on Band 3 continuum data, we derive an upper limit on the intrinsic dust mass of log$_{10}[M_{\rm dust} (M_{\odot})] < 7.73$, consistent with existing estimates. We use the 3D tilted-ring model fitting code 3D-Barolo to determine the kinematic properties of the CO(4-3) emitting gas. We find that it is rotationally dominated, with a $V/σ=4.6 \pm 1.7$, consistent with the kinematics of the [C II]. With PDR conditions remarkably similar to normal dusty star-forming galaxies at z ~ 0.2 and a stable molecular disk, our observations of Az9 suggest that the dust-obscured phase for a low-mass galaxy at z$\sim$4 is relatively long. Thus, Az9 may be representative of a more widespread population that has been missed due to insufficiently deep existing millimeter surveys.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Disentangling the co-evolution of galaxies and supermassive black holes with PRIMA
Authors:
L. Bisigello,
C. Gruppioni,
A. Bolatto,
L. Ciesla,
A. Pope,
L. Armus,
L.,
J. D. Smith,
R. Somerville,
L. Y. A. Yung,
R. J. Wright,
C. M. Bradford,
J. Glenn,
A. Feltre
Abstract:
The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-I…
▽ More
The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-IR provide unmatched diagnostic power to determine the properties of gas and dust, measure gas-phase metallicities and map cold galactic outflows in even the most obscured galaxies. By combining multi-band photometric surveys with low and high-resolution far-IR spectroscopy, the PRobe far-Infrared Mission for Astrophysics (PRIMA), a concept for a far-IR, 1.8m-diameter, cryogenically cooled observatory, will revolutionise the field of galaxy evolution by taking advantage of this IR toolkit to find and study dusty galaxies across galactic time. In this work, we make use of the phenomenological simulation SPRITZ and the Santa Cruz semi-analytical model to describe how a moderately deep multi-band PRIMA photometric survey can easily reach beyond previous IR missions to detect and study galaxies down to $10^{11}\,L_{\odot}$ beyond cosmic noon and at least up to z=4, even in the absence of gravitational lensing. By decomposing the spectral energy distribution (SED) of these photometrically selected galaxies, we show that PRIMA can be used to accurately measure the relative AGN power, the mass fraction contributed by polycyclic aromatic hydrocarbon (PAH) and the total IR luminosity. At the same time, spectroscopic follow up with PRIMA will allow to trace both the star formation and black hole accretion rates (SFR, BHAR), the gas phase metallicities and the mass outflow rates of cold gas in hundreds to thousands of individual galaxies to z=2.
△ Less
Submitted 3 June, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Overcoming Confusion Noise with Hyperspectral Imaging from PRIMAger
Authors:
James M. S. Donnellan,
Seb J. Oliver,
Matthieu Bethermin,
Longji Bing,
Alberto Bolatto,
Charles M. Bradford,
Denis Burgarella,
Laure Ciesla,
Jason Glenn,
Alexandra Pope,
Stephen Serjeant,
Raphael Shirley,
JD T. Smith,
Chris Sorrell
Abstract:
The PRobe far-Infrared Mission for Astrophysics (PRIMA) concept aims to perform mapping with spectral coverage and sensitivities inaccessible to previous FIR space telescopes. PRIMA's imaging instrument, PRIMAger, provides unique hyperspectral imaging simultaneously covering 25-235 $μ$m. We synthesise images representing a deep, 1500 hr deg$^{-2}$ PRIMAger survey, with realistic instrumental and c…
▽ More
The PRobe far-Infrared Mission for Astrophysics (PRIMA) concept aims to perform mapping with spectral coverage and sensitivities inaccessible to previous FIR space telescopes. PRIMA's imaging instrument, PRIMAger, provides unique hyperspectral imaging simultaneously covering 25-235 $μ$m. We synthesise images representing a deep, 1500 hr deg$^{-2}$ PRIMAger survey, with realistic instrumental and confusion noise. We demonstrate that we can construct catalogues of galaxies with a high purity ($>95$ per cent) at a source density of 42k deg$^{-2}$ using PRIMAger data alone. Using the XID+ deblending tool we show that we measure fluxes with an accuracy better than 20 per cent to flux levels of 0.16, 0.80, 9.7 and 15 mJy at 47.4, 79.7, 172, 235 $μ$m respectively. These are a factor of $\sim$2 and $\sim$3 fainter than the classical confusion limits for 72-96 $μ$m and 126-235 $μ$m, respectively. At $1.5 \leq z \leq 2$, we detect and accurately measure fluxes in 8-10 of the 10 channels covering 47-235 $μ$m for sources with $2 \leq$ log(SFR) $\leq 2.5$, a 0.5 dex improvement on what might be expected from the classical confusion limit. Recognising that PRIMager will operate in a context where high quality data will be available at other wavelengths, we investigate the benefits of introducing additional prior information. We show that by introducing even weak prior flux information when employing a higher source density catalogue (more than one source per beam) we can obtain accurate fluxes an order of magnitude below the classical confusion limit for 96-235 $μ$m.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Confusion of extragalactic sources in the far infrared: a baseline assessment of the performance of PRIMAger in intensity and polarization
Authors:
Matthieu Béthermin,
Alberto D. Bolatto,
François Boulanger,
Charles M. Bradford,
Denis Burgarella,
Laure Ciesla,
James Donnellan,
Brandon S. Hensley,
Jason Glenn,
Guilaine Lagache,
Enrique Lopez-Rodriguez,
Seb Oliver,
Alexandra Pope,
Marc Sauvage
Abstract:
Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in in…
▽ More
Because of their limited angular resolution, far-infrared telescopes are usually affected by confusion phenomenon. Since several galaxies can be located in the same instrumental beam, only the brightest objects emerge from the fluctuations caused by fainter sources. The probe far-infrared mission for astrophysics imager (PRIMAger) will observe the mid- and far-infrared (25-235 $μ$m) sky both in intensity and polarization. We aim to provide predictions of the confusion level and its consequences for future surveys. We produced simulated PRIMAger maps affected only by the confusion noise using the simulated infrared extragalactic sky (SIDES) semi-empirical simulation. We then estimated the confusion limit in these maps and extracted the sources using a basic blind extractor. By comparing the input galaxy catalog and the extracted source catalog, we derived various performance metrics as completeness, purity, and the accuracy of various measurements. In intensity, we predict that the confusion limit increases rapidly with increasing wavelength. The confusion limit in polarization is more than 100x lower. The measured flux density is dominated by the brightest galaxy in the beam, but other objects also contribute at longer wavelength (~30% at 235 $μ$m). We also show that galaxy clustering has a mild impact on confusion in intensity (up to 25%), while it is negligible in polarization. In intensity, a basic blind extraction will be sufficient to detect galaxies at the knee of the luminosity function up to z~3 and 10$^{11}$ M$_\odot$ main-sequence galaxies up to z~5. In polarization for a conservative sensitivity, we expect ~8 000 detections up to z=2.5 opening a totally new window on the high-z dust polarization. Finally, we show that intensity surveys at short wavelength and polarization surveys at long wavelength tend to reach confusion at similar depth. There is thus a strong synergy.
△ Less
Submitted 28 April, 2024; v1 submitted 5 April, 2024;
originally announced April 2024.
-
A Performance-Portable SYCL Implementation of CRK-HACC for Exascale
Authors:
Esteban M. Rangel,
S. John Pennycook,
Adrian Pope,
Nicholas Frontiere,
Zhiqiang Ma,
Varsha Madananth
Abstract:
The first generation of exascale systems will include a variety of machine architectures, featuring GPUs from multiple vendors. As a result, many developers are interested in adopting portable programming models to avoid maintaining multiple versions of their code. It is necessary to document experiences with such programming models to assist developers in understanding the advantages and disadvan…
▽ More
The first generation of exascale systems will include a variety of machine architectures, featuring GPUs from multiple vendors. As a result, many developers are interested in adopting portable programming models to avoid maintaining multiple versions of their code. It is necessary to document experiences with such programming models to assist developers in understanding the advantages and disadvantages of different approaches.
To this end, this paper evaluates the performance portability of a SYCL implementation of a large-scale cosmology application (CRK-HACC) running on GPUs from three different vendors: AMD, Intel, and NVIDIA. We detail the process of migrating the original code from CUDA to SYCL and show that specializing kernels for specific targets can greatly improve performance portability without significantly impacting programmer productivity. The SYCL version of CRK-HACC achieves a performance portability of 0.96 with a code divergence of almost 0, demonstrating that SYCL is a viable programming model for performance-portable applications.
△ Less
Submitted 24 October, 2023;
originally announced October 2023.
-
Halfway to the peak: Spatially resolved star formation and kinematics in a z=0.54 dusty galaxy with JWST/MIRI
Authors:
Jason Young,
Alexandra Pope,
Anna Sajina,
Lin Yan,
Thiago S Goncalves,
Miriam Eleazer,
Stacey Alberts,
Lee Armus,
Matteo Bonato,
Daniel A. Dale,
Duncan Farrah,
Carl Ferkinhoff,
Christopher C. Hayward,
Jed McKinney,
Eric J. Murphy,
Nicole Nesvadba,
Patrick Ogle,
Leonid Sajkov,
Sylvain Veilleux
Abstract:
We present JWST/MIRI/MRS observations of an infrared luminous disk galaxy, FLS1, at z=0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar II]6.99 micron, [Ar III]8.99 micron, [Ne II]12.81 micron, and [Ne III]15.55 micron, polycyclic aromatic hydroc…
▽ More
We present JWST/MIRI/MRS observations of an infrared luminous disk galaxy, FLS1, at z=0.54. With a lookback time of 5 Gyr, FLS1 is chronologically at the midpoint between the peak epoch of star formation and the present day. The MRS data provide maps of the atomic fine structure lines [Ar II]6.99 micron, [Ar III]8.99 micron, [Ne II]12.81 micron, and [Ne III]15.55 micron, polycyclic aromatic hydrocarbon (PAH) features at 3.3 micron, 6.2 micron, and 11.3 micron, and the warm molecular gas indicators H2S(5) and H2S(3); all these emission features are spatially resolved. We find that the PAH emission is more extended along the Northern side of the galaxy when compared to the well-studied star-formation tracer [Ne II]. The H2 rotational lines, which are shock indicators, are strongest and most extended on the Southern side of the galaxy. [Ar II] is the second brightest fine structure line detected in FLS1 and we show that it is a useful kinematic probe which can be detected with JWST out to z=3. Velocity maps of [Ar II] show a rotating disk with signs of turbulence. Our results provide an example of how spatially resolved mid-infrared spectroscopy can allow us to better understand the star formation and ISM conditions in a galaxy halfway back to the peak epoch of galaxy evolution.
△ Less
Submitted 12 October, 2023; v1 submitted 10 October, 2023;
originally announced October 2023.
-
The One-hundred-deg^2 DECam Imaging in Narrowbands (ODIN): Survey Design and Science Goals
Authors:
Kyoung-Soo Lee,
Eric Gawiser,
Changbom Park,
Yujin Yang,
Francisco Valdes,
Dustin Lang,
Vandana Ramakrishnan,
Byeongha Moon,
Nicole Firestone,
Stephen Appleby,
Maria Celeste Artale,
Moira Andrews,
Franz E. Bauer,
Barbara Benda,
Adam Broussard,
Yi-Kuan Chiang,
Robin Ciardullo,
Arjun Dey,
Rameen Farooq,
Caryl Gronwall,
Lucia Guaita,
Yun Huang,
Ho Seong Hwang,
Sanghyeok Im,
Woong-Seob Jeong
, et al. (17 additional authors not shown)
Abstract:
We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.…
▽ More
We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.2, 7.4, and 9.8 nm, corresponding to Lya at z=2.4, 3.1, and 4.5 and cosmic times of 2.8, 2.1, and 1.4 Gyr, respectively. When combined with even deeper, public broad-band data from Hyper Suprime-Cam, DECam, and in the future, LSST, the ODIN narrow-band images will enable the selection of over 100,000 Lya-emitting (LAE) galaxies at these epochs. ODIN-selected LAEs will identify protoclusters as galaxy overdensities, and the deep narrow-band images enable detection of highly extended Lya blobs (LABs). Primary science goals include measuring the clustering strength and dark matter halo connection of LAEs, LABs, and protoclusters, and their respective relationship to filaments in the cosmic web. The three epochs allow the redshift evolution of these properties to be determined during the period known as Cosmic Noon, where star formation was at its peak. The two narrow-band filter wavelengths are designed to enable interloper rejection and further scientific studies by revealing [O II] and [O III] at z=0.34, Lya and He II 1640 at z=3.1, and Lyman continuum plus Lya at z=4.5. Ancillary science includes similar studies of the lower-redshift emission-line galaxy samples and investigations of nearby star-forming galaxies resolved into numerous [O III] and [S II] emitting regions.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
Tracing the Total Stellar Mass and Star Formation of High-Redshift Protoclusters
Authors:
Roxana Popescu,
Alexandra Pope,
Kyoung-Soo Lee,
Stacey Alberts,
Yi-Kuan Chiang,
Sowon Lee,
Mark Brodwin,
Jed McKinney,
Vandana Ramakrishnan
Abstract:
As the progenitors of present-day galaxy clusters, protoclusters are excellent laboratories to study galaxy evolution. Since existing observations of protoclusters are limited to the detected constituent galaxies at UV and/or infrared wavelengths, the details of how typical galaxies grow in these young, pre-virialized structures remain uncertain. We measure the total stellar mass and star formatio…
▽ More
As the progenitors of present-day galaxy clusters, protoclusters are excellent laboratories to study galaxy evolution. Since existing observations of protoclusters are limited to the detected constituent galaxies at UV and/or infrared wavelengths, the details of how typical galaxies grow in these young, pre-virialized structures remain uncertain. We measure the total stellar mass and star formation within protoclusters, including the contribution from faint undetected members by performing a stacking analysis of 211 $z=2-4$ protoclusters selected as Planck cold sources. We stack WISE and Herschel/SPIRE images to measure the angular size and the spectral energy distribution of the integrated light from the protoclusters. The fluxes of protoclusters selected as Planck cold sources can be contaminated by line of sight interlopers. Using the WebSky simulation, we estimate that a single protocluster contributes $33\pm15$% of the flux of a Planck cold source on average. After this correction, we obtain a total star formation rate of $7.3\pm3.2 \times 10^3\ M_{\odot} {\rm yr}^{-1}$ and a total stellar mass of $4.9\pm 2.2\times 10^{12}\ M_{\odot}$. Our results indicate that protoclusters have, on average, 2x more star formation and 4x more stellar mass than the total contribution from individually-detected galaxies in spectroscopically-confirmed protoclusters. This suggests that much of the total flux within $z=2-4$ protoclusters comes from galaxies with luminosities lower than the detection limit of SPIRE ($L_{IR} < 3 \times 10^{12} L_{\odot}$). Lastly, we find that protoclusters subtend a half-light radius of 2.8' (4.2-5.8 cMpc) which is consistent with simulations.
△ Less
Submitted 28 December, 2023; v1 submitted 1 August, 2023;
originally announced August 2023.
-
The IR Compactness of Dusty Galaxies Set Star-formation and Dust Properties at z~0-2
Authors:
Jed McKinney,
Alexandra Pope,
Allison Kirkpatrick,
Lee Armus,
Tanio Diaz-Santos,
Carlos Gomez-Guijarro,
Maximilien Franco,
David Elbaz,
Christopher C. Hayward,
Hanae Inami,
Gergo Popping,
Mengyuan Xiao
Abstract:
Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to…
▽ More
Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to measure the extent of cold dust emission in 15 $z\sim2$ IR selected galaxies selected on the basis of having available mid-IR spectroscopy from Spitzer. We use the mid-IR spectra to constrain the relative balance between dust heating from star-formation and active galactic nuclei (AGN), and to measure emission from Polycylic Aromatic Hydrocarbons (PAHs) -- small dust grains that play a key role in the photoelectric heating of gas. In general, we find that dust-obscured star-formation at high IR surface densities exhibits similar properties at low- and high-redshift, namely: local luminous IR galaxies have comparable PAH luminosity to total dust mass ratios as high-$z$ galaxies, and star-formation at $z\sim0-2$ is more efficient at high IR surface densities despite the fact that our sample of high$-z$ galaxies are closer to the main-sequence than local luminous IR galaxies. High star-formation efficiencies are coincident with a decline in the PAH/IR luminosity ratio reminiscent of the deficit observed in far-infrared fine-structure lines. Changes in the gas and dust conditions arising from high star-formation surface densities might help drive the star-formation efficiency up. This could help explain high efficiencies needed to reconcile star-formation and gas volume densities in dusty galaxies at cosmic noon.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
ALMA reveals a stable rotating gas disk in a paradoxical low-mass, ultra-dusty galaxy at z = 4.274
Authors:
Alexandra Pope,
Jed McKinney,
Patrick Kamieneski,
Andrew Battisti,
Itziar Aretxaga,
Gabriel Brammer,
Jose M. Diego,
David H. Hughes,
Erica Keller,
Danilo Marchesini,
Andrew Mizener,
Alfredo Montana,
Eric Murphy,
Katherine E. Whitaker,
Grant Wilson,
Min Yun
Abstract:
We report ALMA detections of [CII] and dust continuum in Az9, a multiply-imaged galaxy behind the Frontier Field cluster MACSJ0717.5+3745. The bright [CII] emission line provides a spectroscopic redshift of z = 4.274. This strongly lensed (mu = 7 +/- 1) galaxy has an intrinsic stellar mass of only 2e9 Msun and a total star formation rate of 26 Msun/yr (~80% of which is dust obscured). Using public…
▽ More
We report ALMA detections of [CII] and dust continuum in Az9, a multiply-imaged galaxy behind the Frontier Field cluster MACSJ0717.5+3745. The bright [CII] emission line provides a spectroscopic redshift of z = 4.274. This strongly lensed (mu = 7 +/- 1) galaxy has an intrinsic stellar mass of only 2e9 Msun and a total star formation rate of 26 Msun/yr (~80% of which is dust obscured). Using public magnification maps, we reconstruct the [CII] emission in the source plane to reveal a stable, rotation-dominated disk with V/sigma = 5.3, which is > 2x higher than predicted from simulations for similarly high-redshift, low-mass galaxies. In the source plane, the [CII] disk has a half-light radius of 1.8 kpc and, along with the dust, is spatially offset from the peak of the stellar light by 1.4 kpc. Az9 is not deficient in [CII]; L[CII]/LIR = 0.0027 consistent with local and high redshift normal star forming galaxies. While dust-obscured star formation is expected to dominate in higher mass galaxies, such a large reservoir of dust and gas in a lower mass disk galaxy 1.4 Gyr after the Big Bang challenges our picture of early galaxy evolution. Furthermore, the prevalence of such low-mass dusty galaxies has important implications for the selection of the highest redshift dropout galaxies with JWST. As one of the lowest stellar mass galaxies at z > 4 to be detected in dust continuum and [CII], Az9 is an excellent laboratory in which to study early dust enrichment in the interstellar medium.
△ Less
Submitted 17 June, 2023;
originally announced June 2023.
-
Galaxy Clustering in the Mira-Titan Universe I: Emulators for the redshift space galaxy correlation function and galaxy-galaxy lensing
Authors:
Juliana Kwan,
Shun Saito,
Alexie Leauthaud,
Katrin Heitmann,
Salman Habib,
Nicholas Frontiere,
Hong Guo,
Song Huang,
Adrian Pope,
Sergio Rodríguez-Torres
Abstract:
We construct accurate emulators for the projected and redshift space galaxy correlation functions and excess surface density as measured by galaxy-galaxy lensing, based on Halo Occupation Distribution (HOD) modeling. Using the complete Mira-Titan suite of 111 $N$-body simulations, our emulators vary over eight cosmological parameters and include the effects of neutrino mass and dynamical dark ener…
▽ More
We construct accurate emulators for the projected and redshift space galaxy correlation functions and excess surface density as measured by galaxy-galaxy lensing, based on Halo Occupation Distribution (HOD) modeling. Using the complete Mira-Titan suite of 111 $N$-body simulations, our emulators vary over eight cosmological parameters and include the effects of neutrino mass and dynamical dark energy. We demonstrate that our emulators are sufficiently accurate for the analysis of the BOSS DR12 CMASS galaxy sample over the range 0.5 < r < 50 Mpc/h. Furthermore, we show that our emulators are capable of recovering unbiased cosmological constraints from realistic mock catalogs over the same range. Our mock catalog tests show the efficacy of combining small scale galaxy-galaxy lensing with redshift space clustering and that we can constrain the growth rate and σ_8 to 7% and 4.5% respectively for a CMASS-like sample using only the measurements covered by our emulator. With the inclusion of a CMB prior on H_0, this reduces to a 2% measurement on the growth rate.
△ Less
Submitted 6 June, 2023; v1 submitted 23 February, 2023;
originally announced February 2023.
-
UV-FIR SED modeling of AGN in IR-luminous galaxies up to z~2.5: Understanding the effects of torus models
Authors:
Alyssa D. Sokol,
M. Yun,
A. Pope,
A. Kirkpatrick,
K. Cooke
Abstract:
UV-FIR SED modeling is an effective way to disentangle emission between star formation (SF) and active galactic nuclei (AGN) in galaxies; however, this approach becomes uncertain for composite AGN/SF galaxies that comprise 50-70% of IR-samples. Cosmic X-ray background (XRB) models require a large fraction of obscured AGN to reproduce the observed XRB peak, motivating reliable SED analyses in objec…
▽ More
UV-FIR SED modeling is an effective way to disentangle emission between star formation (SF) and active galactic nuclei (AGN) in galaxies; however, this approach becomes uncertain for composite AGN/SF galaxies that comprise 50-70% of IR-samples. Cosmic X-ray background (XRB) models require a large fraction of obscured AGN to reproduce the observed XRB peak, motivating reliable SED analyses in objects where the AGN may be ``buried" in the galaxy and in the mid-IR to far-IR SED. In this paper, we study a 24$μ$m-selected ($S_{24}$ > 100$μ$Jy) sample of 95 galaxies with $0 \% < f_{MIR,AGN} < 100 \%$, 0.4 < z < 2.7, and $10^{11}$L$_{\odot}$ < L$_{IR}$ < $10^{13}$L$_{\odot}$. We test the performance of AGN models ranging in torus optical depth via SED fitting, comparing results with Spitzer MIR spectroscopy and X-ray observations. Best-fit torus optical depth can shed light on whether these galaxies host a luminous obscured AGN population. We find that permitting a broader AGN SED parameter space results in improved fit quality with higher optical depths, higher FIR AGN contributions, and higher $L_{Bol}$, impacting the bright-end of the $L_{Bol}$ luminosity function. Our results suggest there may be a population of dust-obscured composites that are bolometrically significant but have their AGN mostly hidden in the mid-IR SED. If so, literature applications of SED fitting that often simplify AGN models or omit optically thick tori may largely underestimate AGN contribution from composite sources, as these sources are both numerous and have solutions sensitive to the assumed range of AGN models.
△ Less
Submitted 20 February, 2023;
originally announced February 2023.
-
Modeling the Galaxy Distribution in Clusters using Halo Cores
Authors:
Danila Korytov,
Esteban Rangel,
Lindsey Bleem,
Nicholas Frontiere,
Salman Habib,
Katrin Heitmann,
Joseph Hollowed,
Adrian Pope
Abstract:
The galaxy distribution in dark matter-dominated halos is expected to approximately trace the details of the underlying dark matter substructure. In this paper we introduce halo `core-tracking' as a way to efficiently follow the small-scale substructure in cosmological simulations and apply the technique to model the galaxy distribution in observed clusters. The method relies on explicitly trackin…
▽ More
The galaxy distribution in dark matter-dominated halos is expected to approximately trace the details of the underlying dark matter substructure. In this paper we introduce halo `core-tracking' as a way to efficiently follow the small-scale substructure in cosmological simulations and apply the technique to model the galaxy distribution in observed clusters. The method relies on explicitly tracking the set of particles identified as belonging to a halo's central density core, once a halo has attained a certain threshold mass. The halo cores are then followed throughout the entire evolution of the simulation. The aim of core-tracking is to simplify substructure analysis tasks by avoiding the use of subhalos and, at the same time, to more easily account for the so-called ``orphan'' galaxies, which have lost substantial dark mass due to tidal stripping. We show that simple models based on halo cores can reproduce the number and spatial distribution of galaxies found in optically-selected clusters in the Sloan Digital Sky Survey. We also discuss future applications of the core-tracking methodology in studying the galaxy-halo connection.
△ Less
Submitted 11 July, 2023; v1 submitted 8 February, 2023;
originally announced February 2023.
-
Broad emission lines in optical spectra of hot dust-obscured galaxies can contribute significantly to JWST/NIRCam photometry
Authors:
Jed McKinney,
Luke Finnerty,
Caitlin Casey,
Maximilien Franco,
Arianna Long,
Seiji Fujimoto,
Jorge Zavala,
Olivia Cooper,
Hollis Akins,
Alexandra Pope,
Lee Armus,
B. T. Soifer,
Kirsten Larson,
Keith Matthews,
Jason Melbourne,
Michael Cushing
Abstract:
Selecting the first galaxies at z>7-10 from JWST surveys is complicated by z<6 contaminants with degenerate photometry. For example, strong optical nebular emission lines at z<6 may mimic JWST/NIRCam photometry of z>7-10 Lyman Break Galaxies (LBGs). Dust-obscured 3<z<6 galaxies in particular are potentially important contaminants, and their faint rest-optical spectra have been historically difficu…
▽ More
Selecting the first galaxies at z>7-10 from JWST surveys is complicated by z<6 contaminants with degenerate photometry. For example, strong optical nebular emission lines at z<6 may mimic JWST/NIRCam photometry of z>7-10 Lyman Break Galaxies (LBGs). Dust-obscured 3<z<6 galaxies in particular are potentially important contaminants, and their faint rest-optical spectra have been historically difficult to observe. A lack of optical emission line and continuum measures for 3<z<6 dusty galaxies now makes it difficult to test their expected JWST/NIRCam photometry for degenerate solutions with NIRCam dropouts. Towards this end, we quantify the contribution by strong emission lines to NIRCam photometry in a physically motivated manner by stacking 21 Keck II/NIRES spectra of hot, dust-obscured, massive ($\log\mathrm{M_*/M_\odot}\gtrsim10-11$) and infrared (IR) luminous galaxies at z~1-4. We derive an average spectrum and measure strong narrow (broad) [OIII]5007 and H$α$ features with equivalent widths of $130\pm20$ A ($150\pm50$ A) and $220\pm30$ A ($540\pm80$ A) respectively. These features can increase broadband NIRCam fluxes by factors of 1.2-1.7 (0.2-0.6 mag). Due to significant dust-attenuation ($A_V\sim6$), we find H$α$+[NII] to be significantly brighter than [OIII]+H$β$, and therefore find that emission-line dominated contaminants of high-z galaxy searches can only reproduce moderately blue perceived UV continua of $S_λ\proptoλ^β$ with $β>-1.5$ and z>4. While there are some redshifts (z~3.75) where our stack is more degenerate with the photometry of z>10 LBGs between $λ_{rest}\sim0.3-0.8\,μ$m, redder filter coverage beyond $λ_{obs}>3.5\,μ$m and far-IR/sub-mm follow-up may be useful for breaking the degeneracy and making a crucial separation between two fairly unconstrained populations, dust-obscured galaxies at z~3-6 and LBGs at z>10.
△ Less
Submitted 30 December, 2022;
originally announced January 2023.
-
The hidden side of cosmic star formation at z > 3: Bridging optically-dark and Lyman break galaxies with GOODS-ALMA
Authors:
Mengyuan Xiao,
David Elbaz,
Carlos Gómez-Guijarro,
Lucas Leroy,
Longji Bing,
Emanuele Daddi,
Benjamin Magnelli,
Maximilien Franco,
Luwenjia Zhou,
Mark Dickinson,
Tao Wang,
Wiphu Rujopakarn,
Georgios E. Magdis,
Ezequiel Treister,
Hanae Inami,
Ricardo Demarco,
Mark T. Sargent,
Xinwen Shu,
Jeyhan S. Kartaltepe,
David M. Alexander,
Matthieu Béthermin,
Frederic Bournaud,
Laure Ciesla,
Henry C. Ferguson,
Steven L. Finkelstein
, et al. (15 additional authors not shown)
Abstract:
Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order…
▽ More
Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift.
△ Less
Submitted 10 February, 2023; v1 submitted 6 October, 2022;
originally announced October 2022.
-
The Past and Future of Mid-Infrared Studies of AGN
Authors:
Anna Sajina,
Mark Lacy,
Alexandra Pope
Abstract:
Observational studies of AGN in the mid-infrared regime are crucial to our understanding of AGN and their role in the evolution of galaxies. Mid-IR-based selection of AGN is complementary to more traditional techniques allowing for a more complete census of AGN activity across cosmic time. Mid-IR observations including time variability and spatially resolved imaging have given us unique insights i…
▽ More
Observational studies of AGN in the mid-infrared regime are crucial to our understanding of AGN and their role in the evolution of galaxies. Mid-IR-based selection of AGN is complementary to more traditional techniques allowing for a more complete census of AGN activity across cosmic time. Mid-IR observations including time variability and spatially resolved imaging have given us unique insights into the nature of the obscuring structures around AGN. The wealth of fine structure, molecular, and dust features in the mid-IR allow us to simultaneously probe multiple components of the ISM allowing us to explore in detail the impact on the host galaxy by the presence of an AGN -- a crucial step toward understanding galaxy-SMBH co-evolution. This review gives a broad overview of this wide range of studies. It also aims to show the evolution of this field starting with its nascency in the 1960s, through major advances thanks to several generations of space-based and ground-based facilities, as well as the promise of upcoming facilities such as the {\sl James Webb Space Telescope (JWST)}.
△ Less
Submitted 5 October, 2022;
originally announced October 2022.
-
Measuring Star Formation and Black Hole Accretion Rates in Tandem using Mid-Infrared Spectra of Local Infrared-Luminous Galaxies
Authors:
Meredith Stone,
Alexandra Pope,
Jed McKinney,
Lee Armus,
Tanio Díaz-Santos,
Hanae Inami,
Allison Kirkpatrick,
Sabrina Stierwalt
Abstract:
We present the results of a stacking analysis performed on Spitzer/Infrared Spectrograph high-resolution mid-infrared spectra of luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). By binning on mid-infrared active galactic nucleus (AGN) fraction and stacking spectra, we detect bright emission lines [Ne II] and [Ne III], which trace star formation, and fainte…
▽ More
We present the results of a stacking analysis performed on Spitzer/Infrared Spectrograph high-resolution mid-infrared spectra of luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). By binning on mid-infrared active galactic nucleus (AGN) fraction and stacking spectra, we detect bright emission lines [Ne II] and [Ne III], which trace star formation, and fainter emission lines [Ne V] and [O IV], which trace AGN activity, throughout the sample. We find the [Ne II] luminosity is fairly constant across all AGN fraction bins, while the [O IV] and [Ne V] luminosities increase by over an order of magnitude. Our measured average line ratios, [Ne V]/[Ne II] and [O IV]/[Ne II], at low AGN fraction are similar to H II galaxies while the line ratios at high AGN fraction are similar to LINERs and Seyferts. We decompose the [O IV] luminosity into star-formation and AGN components by fitting the [O IV] luminosity as a function of the [Ne II] luminosity and the mid-infrared AGN fraction. The [O IV] luminosity in LIRGs is dominated by star formation for mid-infrared AGN fractions $\lesssim0.3$. With the corrected [O IV] luminosity, we calculate black hole accretion rates ranging from $10^{-5}$ M$_{\odot}$/yr at low AGN fractions to 0.2 M$_{\odot}$/yr at the highest AGN fractions. We find that using the [O IV] luminosity, without correcting for star formation, can lead to an overestimate of the BHAR by up to a factor of 30 in starburst dominated LIRGs. Finally, we show the BHAR/SFR ratio increases by more than three orders of magnitude as a function of mid-infrared AGN fraction in LIRGs.
△ Less
Submitted 10 June, 2022;
originally announced June 2022.
-
The dust-to-gas mass ratio of luminous galaxies as a function of their metallicity at cosmic noon
Authors:
Gergö Popping,
Irene Shivaei,
Ryan L. Sanders,
Tucker Jones,
Alexandra Pope,
Naveen A. Reddy,
Alice E. Shapley,
Alison L. Coil,
Mariska Kriek
Abstract:
We aim to quantify the relation between the dust-to-gas mass ratio (DTG) and gas-phase metallicity of $z=$2.1-2.5 luminous galaxies and contrast this high-redshift relation against analogous constraints at z$=$0. We present a sample of ten star-forming main-sequence galaxies in the redshift range $2.1<z<2.5$ with rest-optical emission-line information available from the MOSDEF survey and with ALMA…
▽ More
We aim to quantify the relation between the dust-to-gas mass ratio (DTG) and gas-phase metallicity of $z=$2.1-2.5 luminous galaxies and contrast this high-redshift relation against analogous constraints at z$=$0. We present a sample of ten star-forming main-sequence galaxies in the redshift range $2.1<z<2.5$ with rest-optical emission-line information available from the MOSDEF survey and with ALMA 1.2 millimetre and CO J$=$3-2 follow-up observations. The galaxies have stellar masses ranging from $10^{10.3}$ to $10^{10.6}\,\rm{M}_\odot$ and cover a range in star-formation rate from 35 to 145 $\rm{M}_\odot\,\rm{yr}^{-1}$. We calculated the gas-phase oxygen abundance of these galaxies from rest-optical nebular emission lines (8.4 < $12 + \log{(\rm{O/H})} < 8.8$, corresponding to 0.5 - 1.25 Z$_\odot$). We estimated the dust and H$_2$ masses of the galaxies (using a metallicity-dependent CO-to-H$_2$ conversion factor) from the 1.2~mm and CO J$=$3-2 observations, respectively, from which we estimated a DTG. We find that the galaxies in this sample follow the trends already observed between CO line luminosity and dust-continuum luminosity from $z=0$ to $z=3$, extending such trends to fainter galaxies at $2.1<z<2.5$ than observed to date. We find no second-order metallicity dependence in the CO - dust-continuum luminosity relation for the galaxies presented in this work. The DTGs of main-sequence galaxies at $2.1<z<2.5$ are consistent with an increase in the DTG with gas-phase metallicity. The metallicity dependence of the DTG is driven by the metallicity dependence of the CO-to-H$_2$ conversion factor. Galaxies at z$=$2.1-2.5 are furthermore consistent with the DTG-metallicity relation found at z$=$0 (i.e. with no significant evolution), providing relevant constraints for galaxy formation models.
△ Less
Submitted 19 December, 2022; v1 submitted 18 April, 2022;
originally announced April 2022.
-
Measuring the Total Ultraviolet Light from Galaxy Clusters at z=0.5-1.6: The Balance of Obscured and Unobscured Star-Formation
Authors:
Jed McKinney,
Vandana Ramakrishnan,
Kyoung-Soo Lee,
Alexandra Pope,
Stacey Alberts,
Yi-Kuan Chiang,
Roxana Popescu
Abstract:
Combined observations from UV to IR wavelengths are necessary to fully account for the star-formation in galaxy clusters. Low mass (log M/Msun<10) galaxies are typically not individualy detected, particularly at higher redshifts (z~1-2) where galaxy clusters are undergoing rapid transitions from hosting mostly active, dust-obscured star-forming galaxies to quiescent, passive galaxies. To account f…
▽ More
Combined observations from UV to IR wavelengths are necessary to fully account for the star-formation in galaxy clusters. Low mass (log M/Msun<10) galaxies are typically not individualy detected, particularly at higher redshifts (z~1-2) where galaxy clusters are undergoing rapid transitions from hosting mostly active, dust-obscured star-forming galaxies to quiescent, passive galaxies. To account for these undetected galaxies, we measure the total light emerging from GALEX/NUV stacks of galaxy clusters between z=0.5-1.6. Combined with existing measurements from Spitzer, WISE, and Herschel, we study the average UV through far-infrared (IR) spectral energy distribution (SED) of clusters. From the SEDs, we measure the total stellar mass and amount of dust-obscured and unobscured star-formation arising from all cluster-member galaxies, including the low mass population. The relative fraction of unobscured star-formation we observe in the UV is consistent with what is observed in field galaxies. There is tentative evidence for lower than expected unobscured star-formation at z~0.5, which may arise from rapid redshift evolution in the low mass quenching efficiency in clusters reported by other studies. Finally, the GALEX data places strong constraints on derived stellar-to-halo mass ratios at z<1 which anti-correlate with the total halo mass, consistent with trends found from local X-ray observations of clusters. The data exhibit steeper slopes than implementations of the cluster star-formation efficiency in semi-analytical models.
△ Less
Submitted 4 February, 2022;
originally announced February 2022.
-
IR SED and Dust Masses of Sub-solar Metallicity Galaxies at z~2.3
Authors:
Irene Shivaei,
Gergö Popping,
George Rieke,
Naveen Reddy,
Alexandra Pope,
Robert Kennicutt,
Bahram Mobasher,
Alison Coil,
Yoshinobu Fudamoto,
Mariska Kriek,
Jianwei Lyu,
Pascal Oesch,
Ryan Sanders,
Alice Shapley,
Brian Siana
Abstract:
We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$μ$m), we examine the IR SED of high-redshift sub…
▽ More
We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$μ$m), we examine the IR SED of high-redshift subsolar metallicity (~0.5 $Z_{\odot}$) LIRGs. We find that the data agree well with an average SED template of higher luminosity local low-metallicity dwarf galaxies (reduced $χ^2$ of 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ULIRGs, even in the most favorable case (with reduced $χ^2$ of 2.8), the templates are rejected at >98% confidence level. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar metallicity galaxies may result from different grain properties, a clumpy dust geometry, or a harder/more intense ionizing radiation field that heats the dust to higher temperatures. The obscured SFR indicated by the FIR emission of the subsolar metallicity galaxies is only ~ 60% of the total SFR, which is considerably lower than that of the local LIRGs with ~ 96-97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data can overestimate the Rayleigh-Jeans tail measurements at z~2 by a factor of 2-20, and these templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4dex. At a given stellar mass or metallicity, dust masses at z~2.3 are an order of magnitude higher than those at z~0. Given the predicted molecular gas mass fractions, the observed z~2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies at the same metallicity.
△ Less
Submitted 21 February, 2022; v1 submitted 11 January, 2022;
originally announced January 2022.
-
GOODS-ALMA 2.0: Starbursts in the main sequence reveal compact star formation regulating galaxy evolution prequenching
Authors:
C. Gómez-Guijarro,
D. Elbaz,
M. Xiao,
V. I. Kokorev,
G. E. Magdis,
B. Magnelli,
E. Daddi,
F. Valentino,
M. T. Sargent,
M. Dickinson,
M. Béthermin,
M. Franco,
A. Pope,
B. S. Kalita,
L. Ciesla,
R. Demarco,
H. Inami,
W. Rujopakarn,
X. Shu,
T. Wang,
L. Zhou,
D. M. Alexander,
F. Bournaud,
R. Chary,
H. C. Ferguson
, et al. (16 additional authors not shown)
Abstract:
Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived…
▽ More
Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence.
△ Less
Submitted 7 January, 2022;
originally announced January 2022.
-
Significant Molecular Gas Deficiencies in Star-forming Cluster Galaxies at $z\sim1.4$
Authors:
Stacey Alberts,
Jéa Adams,
Benjamin Gregg,
Alexandra Pope,
Christina C. Williams,
Peter R. M. Eisenhardt
Abstract:
We present the average gas properties derived from ALMA Band 6 dust continuum imaging of 126 massive (log $M_{\star} / M_{\odot} \gtrsim 10.5$), star-forming cluster galaxies across 11 galaxy clusters at $z=1-1.75$. Using stacking analysis on the ALMA images, combined with UV-far-infrared data, we quantify the average infrared SEDs and gas properties (molecular gas masses, $M_{\rm mol}$; gas deple…
▽ More
We present the average gas properties derived from ALMA Band 6 dust continuum imaging of 126 massive (log $M_{\star} / M_{\odot} \gtrsim 10.5$), star-forming cluster galaxies across 11 galaxy clusters at $z=1-1.75$. Using stacking analysis on the ALMA images, combined with UV-far-infrared data, we quantify the average infrared SEDs and gas properties (molecular gas masses, $M_{\rm mol}$; gas depletion timescales, $τ_{\rm depl}$; and gas fractions, f$_{\rm gas}$) as a function of cluster-centric radius and properties including stellar mass and distance from the Main Sequence. We find a significant dearth in the ALMA fluxes relative to that expected in the field $-$ with correspondingly low $M_{\rm mol}$ and f$_{\rm gas}$ and short $τ_{\rm depl}$ $-$ with weak or no dependence on cluster-centric radius out to twice the virial radius. The Herschel+ALMA SEDs indicate warmer dust temperatures ($\sim36-38$ K) than coeval field galaxies ($\sim30$ K). We perform a thorough comparison of the cluster galaxy gas properties to field galaxies, finding deficits of 2-3x, 3-4x, and 2-4x in $M_{\rm mol}$, $τ_{\rm depl}$, and f$_{\rm gas}$ compared to coeval field stacks and larger deficits compared to field scaling relations built primarily on detections. The cluster gas properties derived here are comparable with stacking analyses in (proto-)clusters in the literature and at odds with findings of field-like $τ_{\rm depl}$ and enhanced f$_{\rm gas}$ reported using CO and dust continuum detections. Our analysis suggests that environment has considerable impact on gas properties out to large radii, in good agreement with cosmological simulations which project gas depletion begins beyond the virial radius and largely completes by first passage of the cluster core.
△ Less
Submitted 4 January, 2022;
originally announced January 2022.
-
High Molecular-Gas to Dust Mass Ratios Predicted in Most Quiescent Galaxies
Authors:
Katherine E. Whitaker,
Desika Narayanan,
Christina C. Williams,
Qi Li,
Justin S. Spilker,
Romeel Davé,
Mohammad Akhshik,
Hollis B. Akins,
Rachel Bezanson,
Neal Katz,
Joel Leja,
Georgios E. Magdis,
Lamiya Mowla,
Erica J. Nelson,
Alexandra Pope,
George C. Privon,
Sune Toft,
Francesco Valentino
Abstract:
Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $δ_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum obse…
▽ More
Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $δ_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust has been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme $δ_{\mathrm{GDR}}$ ratios. We show that simulated massive quiescent galaxies at $0 < z < 3$ in the \textsc{simba} cosmological simulations have $δ_{\mathrm{GDR}}$ values that extend $>$4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of $>$3 orders of magnitude in dust mass to return to normal values of $δ_{\mathrm{GDR}}$ despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly pre-select exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities.
△ Less
Submitted 10 November, 2021; v1 submitted 9 November, 2021;
originally announced November 2021.
-
Quenching of star formation from a lack of inflowing gas to galaxies
Authors:
Katherine E. Whitaker,
Christina C. Williams,
Lamiya Mowla,
Justin S. Spilker,
Sune Toft,
Desika Narayanan,
Alexandra Pope,
Georgios E. Magdis,
Pieter G. van Dokkum,
Mohammad Akhshik,
Rachel Bezanson,
Gabriel B. Brammer,
Joel Leja,
Allison Man,
Erica J. Nelson,
Johan Richard,
Camilla Pacifici,
Keren Sharon,
Francesco Valentino
Abstract:
Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion…
▽ More
Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion have been based upon upper limits that are insufficiently sensitive to determine this robustly, or stacked emission with its problems of averaging. Here we report 1.3mm observations of dust emission from six strongly lensed galaxies where star formation has been quenched, with magnifications of up to a factor of 30. Four of the six galaxies are undetected in dust emission, with an estimated upper limit on the dust mass of 0.0001 times the stellar mass, and by proxy (assuming a Milky Way molecular gas-to-dust ratio) 0.01 times the stellar mass in molecular gas. This is two orders of magnitude less molecular gas per unit stellar mass than seen in star forming galaxies at similar redshifts. It remains difficult to extrapolate from these small samples, but these observations establish that gas depletion is responsible for a cessation of star formation in some fraction of high-redshift galaxies.
△ Less
Submitted 21 September, 2021;
originally announced September 2021.
-
Farpoint: A High-Resolution Cosmology Simulation at the Gigaparsec Scale
Authors:
Nicholas Frontiere,
Katrin Heitmann,
Esteban Rangel,
Patricia Larsen,
Adrian Pope,
Imran Sultan,
Thomas Uram,
Salman Habib,
Silvio Rizzi,
Joe Insley
Abstract:
In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000$h^{-1}$Mpc)$^3$ and evolves close to two trillion particles, corresponding to a mass resolution of $m_p\sim 4.6\cdot 10^7 h^{-1}$M$_\odot$. These specifications enable comprehensive investigations of the g…
▽ More
In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000$h^{-1}$Mpc)$^3$ and evolves close to two trillion particles, corresponding to a mass resolution of $m_p\sim 4.6\cdot 10^7 h^{-1}$M$_\odot$. These specifications enable comprehensive investigations of the galaxy-halo connection, capturing halos down to small masses. Further, the large volume resolves scales typical of modern surveys with good statistical coverage of high mass halos. The simulation was carried out on the GPU-accelerated system Summit, one of the fastest supercomputers currently available. We provide specifics about the Farpoint run and present an initial set of results. The high mass resolution facilitates precise measurements of important global statistics, such as the halo concentration-mass relation and the correlation function down to small scales. Selected subsets of the simulation data products are publicly available via the HACC Simulation Data Portal.
△ Less
Submitted 28 February, 2022; v1 submitted 4 September, 2021;
originally announced September 2021.
-
GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies
Authors:
C. Gómez-Guijarro,
D. Elbaz,
M. Xiao,
M. Béthermin,
M. Franco,
B. Magnelli,
E. Daddi,
M. Dickinson,
R. Demarco,
H. Inami,
W. Rujopakarn,
G. E. Magdis,
X. Shu,
R. Chary,
L. Zhou,
D. M. Alexander,
F. Bournaud,
L. Ciesla,
H. C. Ferguson,
S. L. Finkelstein,
M. Giavalisco,
D. Iono,
S. Juneau,
J. S. Kartaltepe,
G. Lagache
, et al. (14 additional authors not shown)
Abstract:
Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas…
▽ More
Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $σ= 68.4μ$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks.
△ Less
Submitted 21 December, 2021; v1 submitted 24 June, 2021;
originally announced June 2021.
-
Early Science with the Large Millimeter Telescope: a 1.1 mm AzTEC Survey of Red-$Herschel$ dusty star-forming galaxies
Authors:
A. Montaña,
J. A. Zavala,
I. Aretxaga,
D. H. Hughes,
R. J. Ivison,
A. Pope,
D. Sánchez-Argüelles,
G. W. Wilson,
M. Yun,
O. A. Cantua,
M. McCrackan,
M. J. Michałowski,
E. Valiante,
V. Arumugam,
C. M. Casey,
R. Chávez,
E. Colín-Beltrán,
H. Dannerbauer,
J. S. Dunlop,
L. Dunne,
S. Eales,
D. Ferrusca,
V. Gómez-Rivera,
A. I. Gómez-Ruiz,
V. H. de la Luz
, et al. (10 additional authors not shown)
Abstract:
We present LMT/AzTEC 1.1mm observations of $\sim100$ luminous high-redshift dusty star-forming galaxy candidates from the $\sim600\,$sq.deg $Herschel$-ATLAS survey, selected on the basis of their SPIRE red far-infrared colours and with $S_{500μ\rm m}=35-80$ mJy. With an effective $θ_{\rm FWHM}\approx9.5\,$ arcsec angular resolution, our observations reveal that at least 9 per cent of the targets b…
▽ More
We present LMT/AzTEC 1.1mm observations of $\sim100$ luminous high-redshift dusty star-forming galaxy candidates from the $\sim600\,$sq.deg $Herschel$-ATLAS survey, selected on the basis of their SPIRE red far-infrared colours and with $S_{500μ\rm m}=35-80$ mJy. With an effective $θ_{\rm FWHM}\approx9.5\,$ arcsec angular resolution, our observations reveal that at least 9 per cent of the targets break into multiple systems with SNR $\geq 4$ members. The fraction of multiple systems increases to $\sim23\,$ per cent (or more) if some non-detected targets are considered multiples, as suggested by the data. Combining the new AzTEC and deblended $Herschel$ photometry we derive photometric redshifts, IR luminosities, and star formation rates. While the median redshifts of the multiple and single systems are similar $(z_{\rm med}\approx3.6)$, the redshift distribution of the latter is skewed towards higher redshifts. Of the AzTEC sources $\sim85\,$ per cent lie at $z_{\rm phot}>3$ while $\sim33\,$ per cent are at $z_{\rm phot}>4$. This corresponds to a lower limit on the space density of ultra-red sources at $4<z<6$ of $\sim3\times10^{-7}\, \textrm{Mpc}^{-3}$ with a contribution to the obscured star-formation of $\gtrsim 8\times10^{-4}\, \textrm{M}_\odot \textrm{yr}^{-1} \textrm{Mpc}^{-3}$. Some of the multiple systems have members with photometric redshifts consistent among them suggesting possible physical associations. Given their angular separations, these systems are most likely galaxy over-densities and/or early-stage pre-coalescence mergers. Finally, we present 3mm LMT/RSR spectroscopic redshifts of six red-$Herschel$ galaxies at $z_{\rm spec}=3.85-6.03$, two of them (at $z \sim 4.7$) representing new redshift confirmations. Here we release the AzTEC and deblended $Herschel$ photometry as well as catalogues of the most promising interacting systems and $z>4$ galaxies.
△ Less
Submitted 6 June, 2021;
originally announced June 2021.
-
Dust-Enshrouded AGN can Dominate Host-Galaxy-Scale Cold-Dust Emission
Authors:
Jed McKinney,
Christopher C. Hayward,
Lee J. Rosenthal,
Juan Rafael Martinez-Galarza,
Alexandra Pope,
Anna Sajina,
Howard A. Smith
Abstract:
It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T~20-40K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at $λ\ga100$ \micron. However, these mode…
▽ More
It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T~20-40K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at $λ\ga100$ \micron. However, these models do not include AGN heating of host-galaxy-scale diffuse dust, which may have far-IR (FIR) colors comparable to cold diffuse dust heated by stars. To quantify the contribution of AGN heating to host-galaxy-scale cold dust emission at $λ\ga100$ \micron, we perform dust RT calculations on a simulated galaxy merger both including and excluding the bright AGN that it hosts. By differencing the SEDs yielded by RT calculations with and without AGN that are otherwise identical, we quantify the FIR cold dust emission arising solely from re-processed AGN photons. In extreme cases, AGN-heated host-galaxy-scale dust can increase galaxy-integrated FIR flux densities by factors of 2-4; star formation rates calculated from the FIR luminosity assuming no AGN contribution can overestimate the true value by comparable factors. Because the FIR colors of such systems are similar to those of purely star-forming galaxies and redder than torus models, broadband SED decomposition may be insufficient for disentangling the contributions of stars and heavily dust-enshrouded AGN in the most IR-luminous galaxies. We demonstrate how kpc-scale resolved observations can be used to identify deeply dust-enshrouded AGN with cool FIR colors when spectroscopic and/or X-ray detection methods are unavailable.
△ Less
Submitted 26 July, 2021; v1 submitted 23 March, 2021;
originally announced March 2021.
-
Cosmic evolution of the H2 mass density and the epoch of molecular gas
Authors:
T. K. Garratt,
K. E. K. Coppin,
J. E. Geach,
O. Almaini,
W. G. Hartley,
D. T. Maltby,
C. J. Simpson,
A. Wilkinson,
C. J. Conselice,
M. Franco,
R. J. Ivison,
M. P. Koprowski,
C. C. Lovell,
A. Pope,
D. Scott,
P. van der Werf
Abstract:
We present new empirical constraints on the evolution of $ρ_{\rm H_2}$, the cosmological mass density of molecular hydrogen, back to $z\approx2.5$. We employ a statistical approach measuring the average observed $850μ{\rm m}$ flux density of near-infrared selected galaxies as a function of redshift. The redshift range considered corresponds to a span where the $850μ{\rm m}$ band probes the Rayleig…
▽ More
We present new empirical constraints on the evolution of $ρ_{\rm H_2}$, the cosmological mass density of molecular hydrogen, back to $z\approx2.5$. We employ a statistical approach measuring the average observed $850μ{\rm m}$ flux density of near-infrared selected galaxies as a function of redshift. The redshift range considered corresponds to a span where the $850μ{\rm m}$ band probes the Rayleigh-Jeans tail of thermal dust emission in the rest-frame, and can therefore be used as an estimate of the mass of the interstellar medium (ISM). Our sample comprises of ${\approx}150,000$ galaxies in the UKIDSS-UDS field with near-infrared magnitudes $K_{\rm AB}\leq25$ mag and photometric redshifts with corresponding probability distribution functions derived from deep 12-band photometry. With a sample approximately 2 orders of magnitude larger than in previous works we significantly reduce statistical uncertainties on $ρ_{\rm H_2}$ to $z\approx2.5$. Our measurements are in broad agreement with recent direct estimates from blank field molecular gas surveys, finding that the epoch of molecular gas coincides with the peak epoch of star formation with $ρ_{\rm H_2}\approx2\times10^7\,{\rm M_\odot}\,{\rm Mpc^{-3}}$ at $z\approx2$. We demonstrate that $ρ_{\rm H_2}$ can be broadly modelled by inverting the star-formation rate density with a fixed or weakly evolving star-formation efficiency. This 'constant efficiency' model shows a similar evolution to our statistically derived $ρ_{\rm H_2}$, indicating that the dominant factor driving the peak star formation history at $z\approx2$ is a larger supply of molecular gas in galaxies rather than a significant evolution of the star-formation rate efficiency within individual galaxies.
△ Less
Submitted 15 March, 2021;
originally announced March 2021.
-
Early Science with the Large Millimeter Telescope: Constraining the Gas Fraction of a Compact Quiescent Galaxy at z=1.883
Authors:
Joyce N. Caliendo,
Katherine E. Whitaker,
Mohammad Akhshik,
Grant Wilson,
Christina C. Williams,
Justin S. Spilker,
Guillaume Mahler,
Alexandra Pope,
Keren Sharon,
Emmaly Aguilar,
Rachel Bezanson,
Miguel Chavez Dagastino,
Arturo I. Gómez-Ruiz,
Alfredo Montaña,
Sune Toft,
Miguel Velazquez De La Rosa,
Milagros Zeballos
Abstract:
We present constraints on the dust continuum flux and inferred gas content of a gravitationally lensed massive quiescent galaxy at $z$=1.883$\pm$0.001 using AzTEC 1.1mm imaging with the Large Millimeter Telescope. MRG-S0851 appears to be a prototypical massive compact quiescent galaxy, but has evidence that it experienced a centrally concentrated rejuvenation event in the last 100 Myr (see Akhshik…
▽ More
We present constraints on the dust continuum flux and inferred gas content of a gravitationally lensed massive quiescent galaxy at $z$=1.883$\pm$0.001 using AzTEC 1.1mm imaging with the Large Millimeter Telescope. MRG-S0851 appears to be a prototypical massive compact quiescent galaxy, but has evidence that it experienced a centrally concentrated rejuvenation event in the last 100 Myr (see Akhshik et al. 2020). This galaxy is undetected in the AzTEC image but we calculate an upper limit on the millimeter flux and use this to estimate the H$_2$ mass limit via an empirically calibrated relation that assumes a constant molecular gas-to-dust ratio of 150. We constrain the 3$σ$ upper limit of the H$_2$ fraction from the dust continuum in MRG-S0851 to be ${M_{H_2}/M_{\star}}$ $\leq$ 6.8%. MRG-S0851 has a low gas fraction limit with a moderately low sSFR owing to the recent rejuvenation episode, which together results in a relatively short depletion time of $<$0.6 Gyr if no further H$_2$ gas is accreted. Empirical and analytical models both predict that we should have detected molecular gas in MRG-S0851, especially given the rejuvenation episode; this suggests that cold gas and/or dust is rapidly depleted in at least some early quiescent galaxies.
△ Less
Submitted 4 March, 2021;
originally announced March 2021.
-
DESC DC2 Data Release Note
Authors:
LSST Dark Energy Science Collaboration,
Bela Abolfathi,
Robert Armstrong,
Humna Awan,
Yadu N. Babuji,
Franz Erik Bauer,
George Beckett,
Rahul Biswas,
Joanne R. Bogart,
Dominique Boutigny,
Kyle Chard,
James Chiang,
Johann Cohen-Tanugi,
Andrew J. Connolly,
Scott F. Daniel,
Seth W. Digel,
Alex Drlica-Wagner,
Richard Dubois,
Eric Gawiser,
Thomas Glanzman,
Salman Habib,
Andrew P. Hearin,
Katrin Heitmann,
Fabio Hernandez,
Renée Hložek
, et al. (32 additional authors not shown)
Abstract:
In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with th…
▽ More
In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with the LSST Science Pipelines (19.0.0). In this Note, we describe the public data release of the resulting object catalogs for the coadded images of five years of simulated observations along with associated truth catalogs. We include a brief description of the major features of the available data sets. To enable convenient access to the data products, we have developed a web portal connected to Globus data services. We describe how to access the data and provide example Jupyter Notebooks in Python to aid first interactions with the data. We welcome feedback and questions about the data release via a GitHub repository.
△ Less
Submitted 13 June, 2022; v1 submitted 12 January, 2021;
originally announced January 2021.
-
Regulating Star Formation in Nearby Dusty Galaxies: Low Photoelectric Efficiencies in the Most Compact Systems
Authors:
Jed McKinney,
Lee Armus,
Alexandra Pope,
Tanio Diaz-Santos,
Vassilis Charmandaris,
Hanae Inami,
Yiqing Song,
Aaron Evans
Abstract:
Star formation in galaxies is regulated by the heating and cooling in the interstellar medium. In particular, the processing of molecular gas into stars will depend strongly on the ratio of gas heating to gas cooling in the neutral gas around sites of recent star-formation. In this work, we combine mid-infrared (mid-IR) observations of Polycyclic Aromatic Hydrocarbons (PAHs), the dominant heating…
▽ More
Star formation in galaxies is regulated by the heating and cooling in the interstellar medium. In particular, the processing of molecular gas into stars will depend strongly on the ratio of gas heating to gas cooling in the neutral gas around sites of recent star-formation. In this work, we combine mid-infrared (mid-IR) observations of Polycyclic Aromatic Hydrocarbons (PAHs), the dominant heating mechanism of gas in the interstellar medium (ISM), with [C II], [O I], and [Si II] fine-structure emission, the strongest cooling channels in dense, neutral gas. The ratio of IR cooling line emission to PAH emission measures the photoelectric efficiency, a property of the ISM which dictates how much energy carried by ultraviolet photons gets transferred into the gas. We find that star-forming, IR luminous galaxies in the Great Observatories All-Sky LIRG Survey (GOALS) with high IR surface densities have low photoelectric efficiencies. These systems also have, on average, higher ratios of radiation field strength to gas densities, and larger average dust grain size distributions. The data support a scenario in which the most compact galaxies have more young star-forming regions per unit area, which exhibit less efficient gas heating. These conditions may be more common at high-z, and may help explain the higher star-formation rates at cosmic noon. We make predictions on how this can be investigated with JWST.
△ Less
Submitted 6 January, 2021; v1 submitted 4 January, 2021;
originally announced January 2021.
-
The LSST DESC DC2 Simulated Sky Survey
Authors:
LSST Dark Energy Science Collaboration,
Bela Abolfathi,
David Alonso,
Robert Armstrong,
Éric Aubourg,
Humna Awan,
Yadu N. Babuji,
Franz Erik Bauer,
Rachel Bean,
George Beckett,
Rahul Biswas,
Joanne R. Bogart,
Dominique Boutigny,
Kyle Chard,
James Chiang,
Chuck F. Claver,
Johann Cohen-Tanugi,
Céline Combet,
Andrew J. Connolly,
Scott F. Daniel,
Seth W. Digel,
Alex Drlica-Wagner,
Richard Dubois,
Emmanuel Gangler,
Eric Gawiser
, et al. (55 additional authors not shown)
Abstract:
We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses…
▽ More
We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep (WFD) area of approximately 300 deg^2 as well as a deep drilling field (DDF) of approximately 1 deg^2. We simulate 5 years of the planned 10-year survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the dataset to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic testbed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time-domain cosmology.
△ Less
Submitted 26 January, 2021; v1 submitted 12 October, 2020;
originally announced October 2020.
-
GOODS-ALMA: Optically dark ALMA galaxies shed light on a cluster in formation at z = 3.5
Authors:
L. Zhou,
D. Elbaz,
M. Franco,
B. Magnelli,
C. Schreiber,
T. Wang,
L. Ciesla,
E. Daddi,
M. Dickinson,
N. Nagar,
G. Magdis,
D. M. Alexander,
M. Béthermin,
R. Demarco,
J. Mullaney,
F. Bournaud,
H. Ferguson,
S. L. Finkelstein,
M. Giavalisco,
H. Inami,
D. Iono,
S. Juneau,
G. Lagache,
H. Messias,
K. Motohara
, et al. (8 additional authors not shown)
Abstract:
In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 2…
▽ More
In this paper we study the properties of the six optically dark galaxies detected in the 69 arcmin^2 GOODS-ALMA 1.1mm continuum survey. While none of them are listed in the deepest H-band based CANDELS catalog in the GOODS-South field down to H=28.16AB, we were able to de-blend two of them from their bright neighbor and measure an $H$-band flux for them. We note that AGS4 and AGS15 have H=25.23, 27.11AB respectively. Their extreme proximity (0.50", 0.27") to a bright optical source and their extreme faintness prevented them from being included in the H-band catalog. We present the spectroscopic scan follow-up of five of the six sources with ALMA band 4. All are detected in the 2mm continuum with signal-to-noise ratios higher than eight. One emission line is detected in AGS4 (ν_{obs} =151.44GHz with a S/N=8.58) and AGS17 (ν_{obs} =154.78GHz with a S/N=10.23), which we interpret in both cases as being due to the CO(6-5) line at z^{AGS4}_{spec}=3.556 and z^{AGS4}_{spec}=3.467, respectively. These redshifts match both the probability distribution of the photometric redshifts derived from the UV to near-infrared spectral energy distributions (SEDs) and the far-infrared SEDs for typical dust temperatures of galaxies at these redshifts. We present evidence that nearly 70% (4/6 of galaxies) of the optically dark galaxies belong to the same overdensity of galaxies at z~3.5. overdensity The most massive one, AGS24 (M_{\star} = 10^{11.32^{+0.02}_{-0.19}} M_{\odot}), is the most massive galaxy without an active galactic nucleus (AGN) at $z$\,>\,3 in the GOODS-ALMA field. It falls in the very center of the peak of the galaxy surface density, which suggests that the surrounding overdensity is a proto-cluster in the process of virialization and that AGS24 is the candidate progenitor of the future brightest cluster galaxy (BCG).
△ Less
Submitted 19 August, 2020;
originally announced August 2020.
-
Measuring the total infrared light from galaxy clusters at z=0.5-1.6: connecting stellar populations to dusty star formation
Authors:
Stacey Alberts,
Kyoung-Soo Lee,
Alexandra Pope,
Mark Brodwin,
Yi-Kuan Chiang,
Jed McKinney,
Rui Xue,
Yun Huang,
Michael Brown,
Arjun Dey,
Peter R. M. Eisenhardt,
Buell T. Jannuzi,
Roxana Popescu,
Vandana Ramakrishnan,
Spencer A. Stanford,
Benjamin J. Weiner
Abstract:
Massive galaxy clusters undergo strong evolution from z~1.6 to z~0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) spectrum can directly trace this transformation; however, such studies have largely been limited to the massive galaxy e…
▽ More
Massive galaxy clusters undergo strong evolution from z~1.6 to z~0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) spectrum can directly trace this transformation; however, such studies have largely been limited to the massive galaxy end of cluster populations. In this work, we present ``total light" stacking techniques spanning 3.4-500μm aimed at revealing the total cluster IR emission, including low mass members and potential intracluster dust. We detail our procedures for WISE, Spitzer, and Herschel imaging, including corrections to recover the total stacked emission in the case of high fractions of detected galaxies. We apply our stacking techniques to 232 well-studied massive (log M200/Msun~13.8) clusters across multiple z bins, recovering extended cluster emission at all wavelengths, typically at >5sigma. We measure the averaged near- to far-IR radial profiles and SEDs, quantifying the total stellar and dust content. The near-IR radial profiles are well described by an NFW model with a high (c~7) concentration parameter. Dust emission is similarly concentrated, albeit suppressed at small radii (r<0.2Mpc). The measured SEDs lack warm dust, consistent with the colder SEDs expected for low mass galaxies. We derive total stellar masses consistent with the theoretical Mhalo-M_star relation and specific-star formation rates that evolve strongly with redshift, echoing that of massive (log Mstar/Msun>10) cluster galaxies. Separating out the massive galaxy population reveals that the majority of cluster far-IR emission (~70-80%) is provided by the low mass constituents, which differs from field galaxies. This effect may be a combination of mass-dependent quenching and excess dust in low mass cluster galaxies.
△ Less
Submitted 3 July, 2020;
originally announced July 2020.
-
The Massive and Distant Clusters of WISE Survey IX: High Radio Activity in a Merging Cluster
Authors:
Emily Moravec,
Anthony Gonzalez,
Simon Dicker,
Stacey Alberts,
Mark Brodwin,
Tracy Clarke,
Thomas Connor,
Bandon Decker,
Mark Devlin,
Peter Eisenhardt,
Brian Mason,
Wenli Mo,
Tony Mroczkowski,
Alexandra Pope,
Charles Romero,
Craig Sarazin,
Jonathan Sievers,
Spencer Stanford,
Daniel Stern,
Dominika Wylezalek,
Fernando Zago
Abstract:
We present a multi-wavelength investigation of the radio galaxy population in the galaxy cluster MOO J1506+5137 at $z$=1.09$\pm$0.03, which in previous work we identified as having multiple complex radio sources. The combined dataset used in this work includes data from the Low-Frequency Array Two-metre Sky Survey (LoTSS), NSF's Karl G. Jansky Very Large Array (VLA), the Robert C. Byrd Green Bank…
▽ More
We present a multi-wavelength investigation of the radio galaxy population in the galaxy cluster MOO J1506+5137 at $z$=1.09$\pm$0.03, which in previous work we identified as having multiple complex radio sources. The combined dataset used in this work includes data from the Low-Frequency Array Two-metre Sky Survey (LoTSS), NSF's Karl G. Jansky Very Large Array (VLA), the Robert C. Byrd Green Bank Telescope (GBT), the Spitzer Space Telescope, and the Dark Energy Camera Legacy Survey (DECaLS). We find that there are five radio sources which are all located within 500 kpc ($\sim$1$^{\prime}$) of the cluster center and have radio luminosities $P_{\mathrm{1.4GHz}}$ > 1.6$\times$10$^{24}$ W Hz$^{-1}$. The typical host galaxies are among the highest stellar mass galaxies in the cluster. The exceptional radio activity among the massive galaxy population appears to be linked to the dynamical state of the cluster. The galaxy distribution suggests an ongoing merger, with a subgroup found to the northwest of the main cluster. Further, two of the five sources are classified as bent-tail sources with one being a potential wide-angle tail (WAT)/hybrid morphology radio source (HyMoRS) indicating a dynamic environment. The cluster also lies in a region of the mass-richness plane occupied by other merging clusters in the Massive and Distant Clusters of WISE Survey (MaDCoWS). The data suggest that during the merger phase radio activity can be dramatically enhanced, which would contribute to the observed trend of increased radio activity in clusters with increasing redshift.
△ Less
Submitted 20 January, 2021; v1 submitted 26 June, 2020;
originally announced June 2020.
-
The Last Journey. I. An Extreme-Scale Simulation on the Mira Supercomputer
Authors:
Katrin Heitmann,
Nicholas Frontiere,
Esteban Rangel,
Patricia Larsen,
Adrian Pope,
Imran Sultan,
Thomas Uram,
Salman Habib,
Hal Finkel,
Danila Korytov,
Eve Kovacs,
Silvio Rizzi,
Joe Insley
Abstract:
The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have b…
▽ More
The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have been run in situ to enable a diverse set of science projects, and at the same time, to keep the resulting data amount manageable. Analysis outputs have been generated starting at redshift z~10 to allow for construction of synthetic galaxy catalogs using a semi-analytic modeling approach in post-processing. As part of our in situ analysis pipeline we employ a new method for tracking halo sub-structures, introducing the concept of subhalo cores. The production of multi-wavelength synthetic sky maps is facilitated by generating particle lightcones in situ, also beginning at z~10. We provide an overview of the simulation set-up and the generated data products; a first set of analysis results is presented. A subset of the data is publicly available.
△ Less
Submitted 8 January, 2021; v1 submitted 2 June, 2020;
originally announced June 2020.
-
GOODS-ALMA: The slow downfall of star-formation in $z$ = 2-3 massive galaxies
Authors:
M. Franco,
D. Elbaz,
L. Zhou,
B. Magnelli,
C. Schreiber,
L. Ciesla,
M. Dickinson,
N. Nagar,
G. Magdis,
D. M. Alexander,
M. Béthermin,
R. Demarco,
E. Daddi,
T. Wang,
J. Mullaney,
M. Sargent,
H. Inami,
X. Shu,
F. Bournaud,
R. Chary,
R. T. Coogan,
H. Ferguson,
S. L. Finkelstein,
M. Giavalisco,
C. Gómez-Guijarro
, et al. (12 additional authors not shown)
Abstract:
We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA a…
▽ More
We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2.
△ Less
Submitted 7 October, 2020; v1 submitted 6 May, 2020;
originally announced May 2020.
-
GOODS-ALMA: Using IRAC and VLA to probe fainter millimeter galaxies
Authors:
M. Franco,
D. Elbaz,
L. Zhou,
B. Magnelli,
C. Schreiber,
L. Ciesla,
M. Dickinson,
N. Nagar,
G. Magdis,
D. M. Alexander,
M. Béthermin,
R. Demarco,
E. Daddi,
T. Wang,
J. Mullaney,
H. Inami,
X. Shu,
F. Bournaud,
R. Chary,
R. T. Coogan,
H. Ferguson,
S. L. Finkelstein,
M. Giavalisco,
C. Gómez-Guijarro,
D. Iono
, et al. (11 additional authors not shown)
Abstract:
In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $σ$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $μ$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$σ$ at which the numbe…
▽ More
In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $σ$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $μ$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$σ$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$σ$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues.
△ Less
Submitted 7 October, 2020; v1 submitted 6 May, 2020;
originally announced May 2020.
-
The Mira-Titan Universe. III. Emulation of the Halo Mass Function
Authors:
Sebastian Bocquet,
Katrin Heitmann,
Salman Habib,
Earl Lawrence,
Thomas Uram,
Nicholas Frontiere,
Adrian Pope,
Hal Finkel
Abstract:
We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological $N$-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs…
▽ More
We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological $N$-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs in the redshift range $z=[0.0, 2.0]$ and for masses $M_{200\mathrm{c}}\geq 10^{13}M_\odot/h$. The emulator covers an 8-dimensional hypercube spanned by {$Ω_\mathrm{m}h^2$, $Ω_\mathrm{b}h^2$, $Ω_νh^2$, $σ_8$, $h$, $n_s$, $w_0$, $w_a$}; spatial flatness is assumed. We obtain smooth halo mass functions by fitting piecewise second-order polynomials to the halo catalogs and employ Gaussian process regression to construct the emulator while keeping track of the statistical noise in the input halo catalogs and uncertainties in the regression process. For redshifts $z\lesssim1$, the typical emulator precision is better than $2\%$ for $10^{13}-10^{14} M_\odot/h$ and $<10\%$ for $M\simeq 10^{15}M_\odot/h$. For comparison, fitting functions using the traditional universal form for the halo mass function can be biased at up to 30\% at $M\simeq 10^{14}M_\odot/h$ for $z=0$. Our emulator is publicly available at \url{https://github.com/SebastianBocquet/MiraTitanHMFemulator}.
△ Less
Submitted 5 August, 2020; v1 submitted 26 March, 2020;
originally announced March 2020.
-
Measuring the Heating and Cooling of the Interstellar Medium at High redshift: PAH and [C II] Observations of the Same Star Forming Galaxies at z~2
Authors:
Jed McKinney,
Alexandra Pope,
Lee Armus,
Ranga-Ram Chary,
Tanio Diaz-Santos,
Mark E. Dickinson,
Allison Kirkpatrick
Abstract:
Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critic…
▽ More
Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C II] 158μm fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C II] in six z ~ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C II] emission, and constrain their relationship as a function of total infrared luminosity (LIR) and IR surface density. [C II] emission is detected in one galaxy at high signal-to-noise (34σ), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C II] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C II]/LIR and PAH/LIR observed in the literature. However, the ratio of [C II] to PAH emission at z ~ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O I] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished.
△ Less
Submitted 19 February, 2020;
originally announced February 2020.
-
An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: Dust attenuation in high-redshift Lyman break Galaxies
Authors:
M. P. Koprowski,
K. E. K. Coppin,
J. E. Geach,
U. Dudzeviciute,
Ian Smail,
O. Almaini,
Fangxia An,
A. W. Blain,
S. C. Chapman,
Chian-Chou Chen,
C. J. Conselice,
J. S. Dunlop,
D. Farrah,
B. Gullberg,
W. Hartley,
R. J. Ivison,
A. Karska,
D. Maltby,
M. J. Michałowski,
A. Pope,
S. Salim,
D. Scott,
C. J. Simpson,
J. M. Simpson,
A. M. Swinbank
, et al. (4 additional authors not shown)
Abstract:
We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV…
▽ More
We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.
△ Less
Submitted 2 January, 2020;
originally announced January 2020.
-
Origins Space Telescope Mission Concept Study Report
Authors:
M. Meixner,
A. Cooray,
D. Leisawitz,
J. Staguhn,
L. Armus,
C. Battersby,
J. Bauer,
E. Bergin,
C. M. Bradford,
K. Ennico-Smith,
J. Fortney,
T. Kataria,
G. Melnick,
S. Milam,
D. Narayanan,
D. Padgett,
K. Pontoppidan,
A. Pope,
T. Roellig,
K. Sandstrom,
K. Stevenson,
K. Su,
J. Vieira,
E. Wright,
J. Zmuidzinas
, et al. (44 additional authors not shown)
Abstract:
The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and lar…
▽ More
The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments.
Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.
△ Less
Submitted 23 December, 2019; v1 submitted 12 December, 2019;
originally announced December 2019.
-
The Massive and Distant Clusters of WISE Survey VII: The Environments and Properties of Radio Galaxies in Clusters at z~1
Authors:
Emily Moravec,
Anthony H. Gonzalez,
Daniel Stern,
Tracy Clarke,
Mark Brodwin,
Bandon Decker,
Peter R. M. Eisenhardt,
Wenli Mo,
Alexandra Pope,
Spencer A. Stanford,
Dominika Wylezalek
Abstract:
We present the results from a study with NSF's Karl G. Jansky Very Large Array (VLA) to determine the radio morphologies of extended radio sources and the properties of their host galaxies in 50 massive galaxy clusters at z~1. We find a majority of the radio morphologies to be Fanaroff-Riley (FR) type IIs. By analyzing the infrared counterparts of the radio sources, we find that ~40% of the host g…
▽ More
We present the results from a study with NSF's Karl G. Jansky Very Large Array (VLA) to determine the radio morphologies of extended radio sources and the properties of their host galaxies in 50 massive galaxy clusters at z~1. We find a majority of the radio morphologies to be Fanaroff-Riley (FR) type IIs. By analyzing the infrared counterparts of the radio sources, we find that ~40% of the host galaxies are the candidate brightest cluster galaxy (BCG) and ~83% are consistent with being one of the top six most massive galaxies in the cluster. We investigate the role of environmental factors on the radio-loud AGN population by examining correlations between environmental and radio-galaxy properties. We find that the highest stellar mass hosts ($M_{*} \gtrsim$ 4$\times 10^{11} M_{\odot}$) are confined to the cluster center and host compact jets. There is evidence for an increase in the size of the jets with cluster-centric radius, which may be attributed to the decreased ICM pressure confinement with increasing radius. Besides this correlation, there are no other significant correlations between the properties of the radio-AGN (luminosity, morphology, or size) and environmental properties (cluster richness and location within the cluster). The fact that there are more AGN in the cluster environment than the field at this epoch, combined with the lack of strong correlation between galaxy and environmental properties, argues that the cluster environment fosters radio activity but does not solely drive the evolution of these sources at this redshift.
△ Less
Submitted 21 November, 2019;
originally announced November 2019.
-
CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST
Authors:
Danila Korytov,
Andrew Hearin,
Eve Kovacs,
Patricia Larsen,
Esteban Rangel,
Joseph Hollowed,
Andrew J. Benson,
Katrin Heitmann,
Yao-Yuan Mao,
Anita Bahmanyar,
Chihway Chang,
Duncan Campbell,
Joseph Derose,
Hal Finkel,
Nicholas Frontiere,
Eric Gawiser,
Salman Habib,
Benjamin Joachimi,
François Lanusse,
Nan Li,
Rachel Mandelbaum,
Christopher Morrison,
Jeffrey A. Newman,
Adrian Pope,
Eli Rykoff
, et al. (5 additional authors not shown)
Abstract:
This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `…
▽ More
This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.
△ Less
Submitted 27 July, 2019; v1 submitted 15 July, 2019;
originally announced July 2019.
-
Simulations Find Our Accounting of Dust-Obscured Star Formation May Be Incomplete
Authors:
Eric Roebuck,
Anna Sajina,
Christopher C. Hayward,
Nicholas Martis,
Danilo Marchesini,
Nicholas Krefting,
Alexandra Pope
Abstract:
The bulk of the star-formation rate density peak at cosmic noon was obscured by dust. How accurately we can assess the role of dust obscured star-formation is affected by inherent biases in our empirical methods -- both those that rely on direct dust emission and those that rely on the inferred dust attenuation of starlight. We use a library of hydrodynamic simulations with radiative transfer to e…
▽ More
The bulk of the star-formation rate density peak at cosmic noon was obscured by dust. How accurately we can assess the role of dust obscured star-formation is affected by inherent biases in our empirical methods -- both those that rely on direct dust emission and those that rely on the inferred dust attenuation of starlight. We use a library of hydrodynamic simulations with radiative transfer to explore these biases. We find that for IR luminous galaxies that are in rapidly quenching systems (e.g. post-coalescence) standard luminosity-to-SFR relations can strongly overestimate the true SFRs. We propose using the $L_{IR}/L_{1.6}$ color to both help identify such systems and provide more accurate SFRs. Conversely, we find that the diagnostic UVJ plot misidentifies a subset of dusty star-forming galaxies. This is due to variability in the effective attenuation curves including being much grayer in the optical-to-near-IR regime than the Calzetti starburst law. This is in agreement with recent observations of IR-selected galaxies at cosmic noon. Our results support the view that we need a panchromatic approach from the rest-frame UV through the IR and SED modeling that includes realistic SFHs and allows for variable attenuation curves if we want to fully account for dust obscured star-formation across the epochs of greatest galaxy build-up.
△ Less
Submitted 7 July, 2019;
originally announced July 2019.
-
CO Emission in Infrared-Selected Active Galactic Nuclei
Authors:
Allison Kirkpatrick,
Chelsea Sharon,
Erica Keller,
Alexandra Pope
Abstract:
In order to better understand how active galactic nuclei (AGN) effect the interstellar media of their host galaxies, we perform a meta-analysis of the CO emission for a sample of $z=0.01-4$ galaxies from the literature with existing CO detections and well-constrained AGN contributions to the infrared (67 galaxies). Using either Spitzer/IRS mid-IR spectroscopy or Spitzer+Herschel colors we determin…
▽ More
In order to better understand how active galactic nuclei (AGN) effect the interstellar media of their host galaxies, we perform a meta-analysis of the CO emission for a sample of $z=0.01-4$ galaxies from the literature with existing CO detections and well-constrained AGN contributions to the infrared (67 galaxies). Using either Spitzer/IRS mid-IR spectroscopy or Spitzer+Herschel colors we determine the fraction of the infrared luminosity in each galaxy that can be attributed to heating by the AGN or stars. We calculate new average CO spectral line ratios (primarily from Carilli \& Walter 2013) to uniformly scale the higher-$J$ CO detections to the ground state and accurately determine our sample's molecular gas masses. We do not find significant differences in the gas depletion timescales/star formation efficiencies (SFEs) as a function of the mid-infrared AGN strength ($f_{\rm AGN}({\rm MIR})$ or $L_{\rm IR} ({\rm AGN})$), which indicates that the presence of an IR-bright AGN is not a sufficient sign-post of galaxy quenching. We also find that the dust-to-gas ratio is consistent for all sources, regardless of AGN emission, redshift, or $L_{\rm IR}$, indicating that dust is likely a reliable tracer of gas mass for massive dusty galaxies (albeit with a large degree of scatter). Lastly, if we classify galaxies as either AGN or star formation dominated, we do not find a robust statistically significant difference between their CO excitation.
△ Less
Submitted 16 May, 2019;
originally announced May 2019.