Showing 1–50 of 210 results for author: Pope, A

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.

Comments: 26 pages, 18 figures; submitted to ApJ

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.

Comments: 10 pages, 13 figures, simulation data portal: https://cosmology.alcf.anl.gov/transfer/newworlds

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.

Comments: 17 pages, 14 figures, submitted to ApJ

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.

Comments: 19 pages, 8 figures

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.

Comments: 16 pages, 17 figures, accepted for publication in A&A

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.

Comments: 14 pages, 11 figures

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.

Comments: 21 pages, 15 figures, 4 tables, submitted to A&A, minor update on the names of the sensitivity scenarii (for consistency with other papers)

Gemini 1.5: Unlocking multimodal understanding across millions of tokens of context

Authors: Gemini Team, Petko Georgiev, Ving Ian Lei, Ryan Burnell, Libin Bai, Anmol Gulati, Garrett Tanzer, Damien Vincent, Zhufeng Pan, Shibo Wang, Soroosh Mariooryad, Yifan Ding, Xinyang Geng, Fred Alcober, Roy Frostig, Mark Omernick, Lexi Walker, Cosmin Paduraru, Christina Sorokin, Andrea Tacchetti, Colin Gaffney, Samira Daruki, Olcan Sercinoglu, Zach Gleicher, Juliette Love , et al. (1092 additional authors not shown)

Abstract: In this report, we introduce the Gemini 1.5 family of models, representing the next generation of highly compute-efficient multimodal models capable of recalling and reasoning over fine-grained information from millions of tokens of context, including multiple long documents and hours of video and audio. The family includes two new models: (1) an updated Gemini 1.5 Pro, which exceeds the February… ▽ More In this report, we introduce the Gemini 1.5 family of models, representing the next generation of highly compute-efficient multimodal models capable of recalling and reasoning over fine-grained information from millions of tokens of context, including multiple long documents and hours of video and audio. The family includes two new models: (1) an updated Gemini 1.5 Pro, which exceeds the February version on the great majority of capabilities and benchmarks; (2) Gemini 1.5 Flash, a more lightweight variant designed for efficiency with minimal regression in quality. Gemini 1.5 models achieve near-perfect recall on long-context retrieval tasks across modalities, improve the state-of-the-art in long-document QA, long-video QA and long-context ASR, and match or surpass Gemini 1.0 Ultra's state-of-the-art performance across a broad set of benchmarks. Studying the limits of Gemini 1.5's long-context ability, we find continued improvement in next-token prediction and near-perfect retrieval (>99%) up to at least 10M tokens, a generational leap over existing models such as Claude 3.0 (200k) and GPT-4 Turbo (128k). Finally, we highlight real-world use cases, such as Gemini 1.5 collaborating with professionals on completing their tasks achieving 26 to 75% time savings across 10 different job categories, as well as surprising new capabilities of large language models at the frontier; when given a grammar manual for Kalamang, a language with fewer than 200 speakers worldwide, the model learns to translate English to Kalamang at a similar level to a person who learned from the same content. △ Less

Submitted 14 June, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Gemini: A Family of Highly Capable Multimodal Models

Authors: Gemini Team, Rohan Anil, Sebastian Borgeaud, Jean-Baptiste Alayrac, Jiahui Yu, Radu Soricut, Johan Schalkwyk, Andrew M. Dai, Anja Hauth, Katie Millican, David Silver, Melvin Johnson, Ioannis Antonoglou, Julian Schrittwieser, Amelia Glaese, Jilin Chen, Emily Pitler, Timothy Lillicrap, Angeliki Lazaridou, Orhan Firat, James Molloy, Michael Isard, Paul R. Barham, Tom Hennigan, Benjamin Lee , et al. (1325 additional authors not shown)

Abstract: This report introduces a new family of multimodal models, Gemini, that exhibit remarkable capabilities across image, audio, video, and text understanding. The Gemini family consists of Ultra, Pro, and Nano sizes, suitable for applications ranging from complex reasoning tasks to on-device memory-constrained use-cases. Evaluation on a broad range of benchmarks shows that our most-capable Gemini Ultr… ▽ More This report introduces a new family of multimodal models, Gemini, that exhibit remarkable capabilities across image, audio, video, and text understanding. The Gemini family consists of Ultra, Pro, and Nano sizes, suitable for applications ranging from complex reasoning tasks to on-device memory-constrained use-cases. Evaluation on a broad range of benchmarks shows that our most-capable Gemini Ultra model advances the state of the art in 30 of 32 of these benchmarks - notably being the first model to achieve human-expert performance on the well-studied exam benchmark MMLU, and improving the state of the art in every one of the 20 multimodal benchmarks we examined. We believe that the new capabilities of the Gemini family in cross-modal reasoning and language understanding will enable a wide variety of use cases. We discuss our approach toward post-training and deploying Gemini models responsibly to users through services including Gemini, Gemini Advanced, Google AI Studio, and Cloud Vertex AI. △ Less

Submitted 17 June, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

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.

Comments: 12 pages, 13 figures, 2023 International Workshop on Performance, Portability & Productivity in HPC

ACM Class: D.2.7; D.2.8; D.1.3; J.2

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.

Comments: 10 pages, 5 figures, 1 table

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.

Comments: 21 pages, 8 figures, submitted to Astrophysical Journal

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.

Comments: Published in ApJ

Journal ref: The Astrophysical Journal, 958, 12 (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.

Comments: 18 pages, 10 figures, accepted to ApJ

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.

Comments: Accepted for publication in ApJ Letters

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.

Comments: 31 pages, 13 figures, revised to match published version

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.

Comments: 20 pages, 10 figures

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.

Comments: 17 pages, 20 figures, 1 Appendix; version accepted by OJA

Journal ref: Open Journal of Astrophysics 6 (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.

Comments: 8 pages, 3 figures, 1 table, submitted to ApJL

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.

Comments: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

Journal ref: A&A 672, A18 (2023)

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.

Comments: 37 pages, 20 figures; Invited review (published) for the special issue "Recent Advances in Infrared Galaxies and AGN", edited by Anna Sajina and Asantha R. Cooray, in Universe

Journal ref: Universe 2022, 8(7), 356

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.

Comments: 17 pages, 15 figures. Accepted to ApJ

De novo design of protein target specific scaffold-based Inhibitors via Reinforcement Learning

Authors: Andrew D. McNaughton, Mridula S. Bontha, Carter R. Knutson, Jenna A. Pope, Neeraj Kumar

Abstract: Efficient design and discovery of target-driven molecules is a critical step in facilitating lead optimization in drug discovery. Current approaches to develop molecules for a target protein are intuition-driven, hampered by slow iterative design-test cycles due to computational challenges in utilizing 3D structural data, and ultimately limited by the expertise of the chemist - leading to bottlene… ▽ More Efficient design and discovery of target-driven molecules is a critical step in facilitating lead optimization in drug discovery. Current approaches to develop molecules for a target protein are intuition-driven, hampered by slow iterative design-test cycles due to computational challenges in utilizing 3D structural data, and ultimately limited by the expertise of the chemist - leading to bottlenecks in molecular design. In this contribution, we propose a novel framework, called 3D-MolGNN$_{RL}$, coupling reinforcement learning (RL) to a deep generative model based on 3D-Scaffold to generate target candidates specific to a protein building up atom by atom from the starting core scaffold. 3D-MolGNN$_{RL}$ provides an efficient way to optimize key features by multi-objective reward function within a protein pocket using parallel graph neural network models. The agent learns to build molecules in 3D space while optimizing the activity, binding affinity, potency, and synthetic accessibility of the candidates generated for infectious disease protein targets. Our approach can serve as an interpretable artificial intelligence (AI) tool for lead optimization with optimized activity, potency, and biophysical properties. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: Published at the MLDD workshop, ICLR 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.

Comments: Version accepted for publication in A&A

Journal ref: A&A 670, A138 (2023)

Soft Actor-Critic with Inhibitory Networks for Faster Retraining

Authors: Jaime S. Ide, Daria Mićović, Michael J. Guarino, Kevin Alcedo, David Rosenbluth, Adrian P. Pope

Abstract: Reusing previously trained models is critical in deep reinforcement learning to speed up training of new agents. However, it is unclear how to acquire new skills when objectives and constraints are in conflict with previously learned skills. Moreover, when retraining, there is an intrinsic conflict between exploiting what has already been learned and exploring new skills. In soft actor-critic (SAC… ▽ More Reusing previously trained models is critical in deep reinforcement learning to speed up training of new agents. However, it is unclear how to acquire new skills when objectives and constraints are in conflict with previously learned skills. Moreover, when retraining, there is an intrinsic conflict between exploiting what has already been learned and exploring new skills. In soft actor-critic (SAC) methods, a temperature parameter can be dynamically adjusted to weight the action entropy and balance the explore $\times$ exploit trade-off. However, controlling a single coefficient can be challenging within the context of retraining, even more so when goals are contradictory. In this work, inspired by neuroscience research, we propose a novel approach using inhibitory networks to allow separate and adaptive state value evaluations, as well as distinct automatic entropy tuning. Ultimately, our approach allows for controlling inhibition to handle conflict between exploiting less risky, acquired behaviors and exploring novel ones to overcome more challenging tasks. We validate our method through experiments in OpenAI Gym environments. △ Less

Submitted 7 February, 2022; v1 submitted 6 February, 2022; originally announced February 2022.

Comments: 16 pages including Appendix

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.

Comments: 13 pages, 5 figures, accepted to ApJ

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.

Comments: Accepted for publication in ApJ. 29 pages, 2 appendices, 10 figures, 3 tables

Journal ref: The Astrophysical Journal 2022, Volume 928, Issue 1, id.68, 26 pp

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.

Comments: Accepted by A&A. 26 pages, 13 figures

Journal ref: A&A 659, A196 (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.

Comments: Accepted for publication in ApJ

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.

Comments: 9 pages, 5 figures, version accepted for publication in the Astrophysical Journal Letters

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.

Comments: 17 pages, 3 figures. Authors' version. Published online by Nature on September 22, 2021

Enabling particle applications for exascale computing platforms

Authors: Susan M Mniszewski, James Belak, Jean-Luc Fattebert, Christian FA Negre, Stuart R Slattery, Adetokunbo A Adedoyin, Robert F Bird, Choongseok Chang, Guangye Chen, Stephane Ethier, Shane Fogerty, Salman Habib, Christoph Junghans, Damien Lebrun-Grandie, Jamaludin Mohd-Yusof, Stan G Moore, Daniel Osei-Kuffuor, Steven J Plimpton, Adrian Pope, Samuel Temple Reeve, Lee Ricketson, Aaron Scheinberg, Amil Y Sharma, Michael E Wall

Abstract: The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four sub-motifs: short-range particle-particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed partic… ▽ More The Exascale Computing Project (ECP) is invested in co-design to assure that key applications are ready for exascale computing. Within ECP, the Co-design Center for Particle Applications (CoPA) is addressing challenges faced by particle-based applications across four sub-motifs: short-range particle-particle interactions (e.g., those which often dominate molecular dynamics (MD) and smoothed particle hydrodynamics (SPH) methods), long-range particle-particle interactions (e.g., electrostatic MD and gravitational N-body), particle-in-cell (PIC) methods, and linear-scaling electronic structure and quantum molecular dynamics (QMD) algorithms. Our crosscutting co-designed technologies fall into two categories: proxy applications (or apps) and libraries. Proxy apps are vehicles used to evaluate the viability of incorporating various types of algorithms, data structures, and architecture-specific optimizations and the associated trade-offs; examples include ExaMiniMD, CabanaMD, CabanaPIC, and ExaSP2. Libraries are modular instantiations that multiple applications can utilize or be built upon; CoPA has developed the Cabana particle library, PROGRESS/BML libraries for QMD, and the SWFFT and fftMPI parallel FFT libraries. Success is measured by identifiable lessons learned that are translated either directly into parent production application codes or into libraries, with demonstrated performance and/or productivity improvement. The libraries and their use in CoPA's ECP application partner codes are also addressed. △ Less

Submitted 19 September, 2021; originally announced September 2021.

Comments: 26 pages, 17 figures

Report number: LA-UR-20-26599

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.

Journal ref: The Astrophysical Journal Supplement Series, 259(1):15, 2022

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.

Comments: Accepted by A&A. 30 pages, 18 figures

Journal ref: A&A 658, A43 (2022)

Launchpad: A Programming Model for Distributed Machine Learning Research

Authors: Fan Yang, Gabriel Barth-Maron, Piotr Stańczyk, Matthew Hoffman, Siqi Liu, Manuel Kroiss, Aedan Pope, Alban Rrustemi

Abstract: A major driver behind the success of modern machine learning algorithms has been their ability to process ever-larger amounts of data. As a result, the use of distributed systems in both research and production has become increasingly prevalent as a means to scale to this growing data. At the same time, however, distributing the learning process can drastically complicate the implementation of eve… ▽ More A major driver behind the success of modern machine learning algorithms has been their ability to process ever-larger amounts of data. As a result, the use of distributed systems in both research and production has become increasingly prevalent as a means to scale to this growing data. At the same time, however, distributing the learning process can drastically complicate the implementation of even simple algorithms. This is especially problematic as many machine learning practitioners are not well-versed in the design of distributed systems, let alone those that have complicated communication topologies. In this work we introduce Launchpad, a programming model that simplifies the process of defining and launching distributed systems that is specifically tailored towards a machine learning audience. We describe our framework, its design philosophy and implementation, and give a number of examples of common learning algorithms whose designs are greatly simplified by this approach. △ Less

Submitted 7 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.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 23 pages, 18 figures

The Petascale DTN Project: High Performance Data Transfer for HPC Facilities

Authors: Eli Dart, William Allcock, Wahid Bhimji, Tim Boerner, Ravinderjeet Cheema, Andrew Cherry, Brent Draney, Salman Habib, Damian Hazen, Jason Hill, Matt Kollross, Suzanne Parete-Koon, Daniel Pelfrey, Adrian Pope, Jeff Porter, David Wheeler

Abstract: The movement of large-scale (tens of Terabytes and larger) data sets between high performance computing (HPC) facilities is an important and increasingly critical capability. A growing number of scientific collaborations rely on HPC facilities for tasks which either require large-scale data sets as input or produce large-scale data sets as output. In order to enable the transfer of these data sets… ▽ More The movement of large-scale (tens of Terabytes and larger) data sets between high performance computing (HPC) facilities is an important and increasingly critical capability. A growing number of scientific collaborations rely on HPC facilities for tasks which either require large-scale data sets as input or produce large-scale data sets as output. In order to enable the transfer of these data sets as needed by the scientific community, HPC facilities must design and deploy the appropriate data transfer capabilities to allow users to do data placement at scale. This paper describes the Petascale DTN Project, an effort undertaken by four HPC facilities, which succeeded in achieving routine data transfer rates of over 1PB/week between the facilities. We describe the design and configuration of the Data Transfer Node (DTN) clusters used for large-scale data transfers at these facilities, the software tools used, and the performance tuning that enabled this capability. △ Less

Submitted 8 September, 2021; v1 submitted 26 May, 2021; originally announced May 2021.

Hierarchical Reinforcement Learning for Air-to-Air Combat

Authors: Adrian P. Pope, Jaime S. Ide, Daria Micovic, Henry Diaz, David Rosenbluth, Lee Ritholtz, Jason C. Twedt, Thayne T. Walker, Kevin Alcedo, Daniel Javorsek

Abstract: Artificial Intelligence (AI) is becoming a critical component in the defense industry, as recently demonstrated by DARPA`s AlphaDogfight Trials (ADT). ADT sought to vet the feasibility of AI algorithms capable of piloting an F-16 in simulated air-to-air combat. As a participant in ADT, Lockheed Martin`s (LM) approach combines a hierarchical architecture with maximum-entropy reinforcement learning… ▽ More Artificial Intelligence (AI) is becoming a critical component in the defense industry, as recently demonstrated by DARPA`s AlphaDogfight Trials (ADT). ADT sought to vet the feasibility of AI algorithms capable of piloting an F-16 in simulated air-to-air combat. As a participant in ADT, Lockheed Martin`s (LM) approach combines a hierarchical architecture with maximum-entropy reinforcement learning (RL), integrates expert knowledge through reward shaping, and supports modularity of policies. This approach achieved a $2^{nd}$ place finish in the final ADT event (among eight total competitors) and defeated a graduate of the US Air Force's (USAF) F-16 Weapons Instructor Course in match play. △ Less

Submitted 11 June, 2021; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: 10 pages, 10 figures, The 2021 International Conference on Unmanned Aircraft System (ICUAS 21), June 15-18, 2021, Athens, Greece

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.

Comments: 11 pages, 7 figures, Accepted to ApJ

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.

Comments: Accepted for publication in ApJ

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.

Comments: 8 pages, 4 figures, accepted for publication in Astrophysical Journal Letters in January 2021 (in press)

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.

Comments: 25 pages, 3 figures; 9 tables. A detailed changelog can be found in Appendix A. To obtain data, visit the DESC Data Portal at https://data.lsstdesc.org/

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.

Comments: 17 pages, 8 figures, 2 tables, accepted for publication in ApJ

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.

Comments: 39 pages, 19 figures, version accepted for publication in ApJS

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.

Comments: 17 pages, 13 figures, 3 tables, accepted by A&A

Journal ref: A&A 642, A155 (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.

Comments: 32 pages, 16 figures, 7 tables. Submitted to MNRAS, comments welcome!

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.

Comments: 17 pages and 8 figures. Accepted in ApJ for publication

Journal ref: 2020 ApJ 898 2

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.

Comments: 14 pages, 9 figures. Accepted for publication in ApJS. New visualization and new results for the matter correlation function added, minor edits. The Last Journey data products can be accessed here: https://cosmology.alcf.anl.gov/

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.

Comments: 18 pages, 8 figures, Accepted for publication in A&A

Journal ref: A&A 643, A30 (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.

Comments: 18 pages, 13 figures. Accepted for publication in A&A

Journal ref: A&A 643, A53 (2020)