Search

We present Mantis Shrimp, a multi-survey deep learning model for photometric redshift estimation that fuses ultra-violet (GALEX), optical (PanSTARRS), and infrared (UnWISE) imagery. Machine learning is now an established approach for photometric redshift estimation, with generally acknowledged higher performance in areas with a high density of spectroscopically identified galaxies over template-based methods. Multiple works have shown that image-based convolutional neural networks can outperform tabular-based color/magnitude models. In comparison to tabular models, image models have additional design complexities: it is largely unknown how to fuse inputs from different instruments which have different resolutions or noise properties. The Mantis Shrimp model estimates the conditional density estimate of redshift using cutout images. The density estimates are well calibrated and the point estimates perform well in the distribution of available spectroscopically confirmed galaxies with (bias = 1e-2), scatter (NMAD = 2.44e-2) and catastrophic outlier rate ($\eta$=17.53$\%$). We find that early fusion approaches (e.g., resampling and stacking images from different instruments) match the performance of late fusion approaches (e.g., concatenating latent space representations), so that the design choice ultimately is left to the user. Finally, we study how the models learn to use information across bands, finding evidence that our models successfully incorporates information from all surveys. The applicability of our model to the analysis of large populations of galaxies is limited by the speed of downloading cutouts from external servers; however, our model could be useful in smaller studies such as generating priors over redshift for stellar population synthesis.

We present ORACLE, the first hierarchical deep-learning model for real-time, context-aware classification of transient and variable astrophysical phenomena. ORACLE is a recurrent neural network with Gated Recurrent Units (GRUs), and has been trained using a custom hierarchical cross-entropy loss function to provide high-confidence classifications along an observationally-driven taxonomy with as little as a single photometric observation. Contextual information for each object, including host galaxy photometric redshift, offset, ellipticity and brightness, is concatenated to the light curve embedding and used to make a final prediction. Training on $\sim$0.5M events from the Extended LSST Astronomical Time-Series Classification Challenge, we achieve a top-level (Transient vs Variable) macro-averaged precision of 0.96 using only 1 day of photometric observations after the first detection in addition to contextual information, for each event; this increases to $>$0.99 once 64 days of the light curve has been obtained, and 0.83 at 1024 days after first detection for 19-way classification (including supernova sub-types, active galactic nuclei, variable stars, microlensing events, and kilonovae). We also compare ORACLE with other state-of-the-art classifiers and report comparable performance for the 19-way classification task, in addition to delivering accurate top-level classifications much earlier. The code and model weights used in this work are publicly available at our associated GitHub repository (https://github.com/uiucsn/ELAsTiCC-Classification)., Comment: 29 pages, 19 figures, 9 tables. Submitted to ApJ

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards ($16.5 < V < 19.5$) alongside the three CALSPEC standards, from 912 \r{A} to 32 $\mu$m. The framework is the first of its kind to jointly infer photometric zeropoints and WD parameters ($\log g$, $T_{\text{eff}}$, $A_V$, $R_V$) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic HST/WFC3 UVIS and IR fluxes, HST/STIS UV spectroscopy and ground-based optical spectroscopy to sub-percent precision. Photometric residuals for the sample are the lowest yet yielding $<0.004$ mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred SEDs from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory's Legacy Survey of Space and Time, and the Nancy Grace Roman Space Telescope., Comment: 32 pages, 24 figures, 5 tables, submitted to MNRAS

Fainter standard stars are essential for the calibration of larger telescopes. This work adds to the CALSPEC (calibration spectra) database 19 faint white dwarfs (WDs) with all-sky coverage and V magnitudes between 16.5 and 18.7. Included for these stars is new UV (ultraviolet) HST (Hubble Space Telescope) STIS (Space Telescope Imaging Spectrometer) spectrophotometry between 1150 and 3000~\AA\ with a resolution of $\sim$500. Pure hydrogen WD models are fit to these UV spectra and to six-band HST/WFC3 (Wide Field Camera 3) photometry at 0.28 to 1.6~\micron\ to construct predicted model SEDs (spectral energy distributions) covering wavelengths from 900~\AA\ to the JWST (James Webb Space Telescope) limit of 30~\micron\ using well-established CALSPEC procedures for producing flux standards with the goal of 1\% accuracy.

The discovery of a second electromagnetic counterpart to a gravitational wave event represents a critical goal in the field of multi-messenger astronomy. In order to determine the optimal strategy for achieving this goal, we perform comprehensive simulations comparing two potential paths forward: continuing the current LIGO-Virgo-KAGRA (LVK) observing run, O4, versus temporarily shutting down the detectors for upgrades before beginning the next observing run, O5. Our simulations incorporate current O4 instrument sensitivities and duty cycles, as well as projected configurations for O5, while accounting for variables such as binary neutron star merger rates, system properties, viewing angles, dust extinction, and kilonova (KN) observables. Our results indicate that a KN discovery would occur $125^{+253}_{-125}$~days (middle 50\% interval) sooner in O5 compared to O4, suggesting that extending O4 would lead to faster discovery if the shutdown period between runs is $>$4~months. Moreover, for 88\% of our simulations, continuing O4 results in earlier KN discovery when compared to the expected two-year shutdown between O4 and O5. Given these findings and the critical importance of avoiding a $>$10 year gap between first and second electromagnetic counterpart discoveries, we suggest LVK consider extending O4 operations for as long as feasible prior to shutting down for critical upgrades., Comment: 10 pages, 10 figures, 3 tables. Submitted to PASP

We present the Early Data Release of the Multi-Cycle JWST-NEXUS Treasury program (2024-2028), which includes NIRCam imaging and WFSS observations from the first (partial) NEXUS-Wide epoch covering the central 100 ${\rm arcmin^2}$ of the NEXUS field, located near the North Ecliptic Pole and within the Euclid Ultra-Deep Field. We release reduced NIRCam mosaics (F090W, F115W, F150W, F200W, F356W, F444W), photometric source catalogs, as well as preliminary WFSS spectra (in F322W2 and F444W) for the subset of bright sources (F356W$<$21 mag or F444W$<$21 mag). These observations fully cover the NEXUS-Deep area, and anchor the long-term baseline of the program. These data will be used for initial target selection for the NIRSpec/MSA spectroscopy starting from June 2025. The NIRCam imaging reaches depths of 27.4--28.2 (AB) mags in F090W--F444W. Upcoming NEXUS-Wide epochs will expand the area to the full $\sim 400\,{\rm arcmin^2}$, and improve the NIRCam exposure depths in the Wide tier by a factor of three. In addition, this central region will be repeatedly covered by the NEXUS-Deep observations (NIRCam imaging and NIRSpec/MSA PRISM spectroscopy) over 18 epochs with a $\sim 2$-month cadence. We demonstrate the data quality of the first NEXUS observations, and showcase some example science cases enabled by these data., Comment: 14 pages, 9 figures, 2 tables. Data products are publicly accessible at https://ariel.astro.illinois.edu/nexus/edr/ Online interactive map for quick visualization of released images and WFSS spectra can be found at https://ariel.astro.illinois.edu/nexus/map/

The Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory is planned to begin in the Fall of 2025. The LSST survey cadence has been designed via a community-driven process regulated by the Survey Cadence Optimization Committee (SCOC), which recommended up to 3% of the observing time to carry out Target of Opportunity (ToO) observations. Experts from the scientific community, Rubin Observatory personnel, and members of the SCOC were brought together to deliver a recommendation for the implementation of the ToO program during a workshop held in March 2024. Four main science cases were identified: gravitational wave multi-messenger astronomy, high energy neutrinos, Galactic supernovae, and small potentially hazardous asteroids possible impactors. Additional science cases were identified and briefly addressed in the documents, including lensed or poorly localized gamma-ray bursts and twilight discoveries. Trigger prioritization, automated response, and detailed strategies were discussed for each science case. This document represents the outcome of the Rubin ToO 2024 workshop, with additional contributions from members of the Rubin Science Collaborations. The implementation of the selection criteria and strategies presented in this document has been endorsed in the SCOC Phase 3 Recommendations document (PSTN-056). Although the ToO program is still to be finalized, this document serves as a baseline plan for ToO observations with the Rubin Observatory.

NEXUS is a JWST Multi-Cycle (Cycles 3-5; 368 primary hrs) GO Treasury imaging and spectroscopic survey around the North Ecliptic Pole. It contains two overlapping tiers. The Wide tier ($\sim 400~{\rm arcmin}^2$) performs NIRCam/WFSS 2.4-5 micron grism spectroscopy with three epochs over 3 years (final continuum ${\rm S/N/pixel>3}$ at F444W$<22.2$). The Deep tier ($\sim 50~{\rm arcmin}^2$) performs high-multiplexing NIRSpec 0.6-5.3 micron MOS/PRISM spectroscopy for $\sim 10,000$ targets, over 18 epochs with a 2-month cadence (epoch/final continuum ${\rm S/N/pixel>3}$ at F200W$\lesssim 27/29$). All epochs have simultaneous multi-band NIRCam and MIRI imaging ($5\sigma$ final depths of $\sim 28-29$ in NIRCam and $\sim 25$ in MIRI). The field is within the continuous viewing zone of JWST, and is fully covered by the Euclid Ultra-Deep Field, with 0.9-2 micron deep Euclid spectroscopy and cadenced photometry. NEXUS has three science pillars. First, with its massive and nearly complete (flux-limited) spectroscopic samples and deep photometry, it will perform efficient classification and physical characterization of galaxies and AGNs from $z\sim 1$ to Cosmic Dawn. With the large contiguous area coverage, it will measure the spatial clustering and demography of the first galaxies and SMBHs at $z>6$. Second, multi-epoch observations enable systematic time-domain investigations, focusing on $z>3$ transients and low-mass AGN reverberation mapping. Third, the comprehensive data set will enable knowledge transfer to other legacy fields, create data challenges, and initiate benchmark work for future space missions. With rapid public releases of processed data and an open invitation for collaboration, NEXUS aims for broad and swift community engagement, to become a powerhouse to drive transformative advancements in multiple key science areas of astronomy., Comment: 19 pages

The parametrized post-Einsteinian (ppE) framework and its variants are widely used to probe gravity through gravitational-wave tests that apply to a large class of theories beyond general relativity. However, the ppE framework is not truly theory-agnostic as it only captures certain types of deviations from general relativity: those that admit a post-Newtonian series representation in the inspiral of coalescencing compact objects. Moreover, each type of deviation in the ppE framework has to be tested separately, making the whole process computationally inefficient and expensive, possibly obscuring the theoretical interpretation of potential deviations that could be detected in the future. We here present the neural post-Einsteinian (npE) framework, an extension of the ppE formalism that overcomes the above weaknesses using deep-learning neural networks. The core of the npE framework is a variantional autoencoder that maps the discrete ppE theories into a continuous latent space in a well-organized manner. This design enables the npE framework to test many theories simultaneously and to select the theory that best describes the observation in a single parameter estimation run. The smooth extension of the ppE parametrization also allows for more general types of deviations to be searched for with the npE model. We showcase the application of the new npE framework to future tests of general relativity with the fifth observing run of the LIGO-Virgo-KAGRA collaboration. In particular, the npE framework is demonstrated to efficiently explore modifications to general relativity beyond what can be mapped by the ppE framework, including modifications coming from higher-order curvature corrections to the Einstein-Hilbert action at high post-Newtonian order, and dark-photon interactions in possibly hidden sectors of matter that do not admit a post-Newtonian representation., Comment: 30 pages, 15 figures, submitted to PRD

The availability of large, public, multi-modal astronomical datasets presents an opportunity to execute novel research that straddles the line between science of AI and science of astronomy. Photometric redshift estimation is a well-established subfield of astronomy. Prior works show that computer vision models typically outperform catalog-based models, but these models face additional complexities when incorporating images from more than one instrument or sensor. In this report, we detail our progress creating Mantis Shrimp, a multi-survey computer vision model for photometric redshift estimation that fuses ultra-violet (GALEX), optical (PanSTARRS), and infrared (UnWISE) imagery. We use deep learning interpretability diagnostics to measure how the model leverages information from the different inputs. We reason about the behavior of the CNNs from the interpretability metrics, specifically framing the result in terms of physically-grounded knowledge of galaxy properties., Comment: 4 pages, 1 figure, 1 table. Submitted to AI4Differential Equations in Science Workshop at ICLR24. Public repository unavailable while under institutional review

The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and production of both astronomical data and software. This report is the culmination of a workshop held in February 2023 on the Future of Astronomical Data Infrastructure. Attended by 70 scientists and software professionals from ground-based and space-based missions and archives spanning the entire spectrum of astronomical research, the group deliberated on the prevailing state of software and data infrastructure in astronomy, identified pressing issues, and explored potential solutions. In this report, we describe the ecosystem of astronomical data, its existing flaws, and the many gaps, duplication, inconsistencies, barriers to access, drags on productivity, missed opportunities, and risks to the long-term integrity of essential data sets. We also highlight the successes and failures in a set of deep dives into several different illustrative components of the ecosystem, included as an appendix., Comment: 59 pages; please send comments and/or questions to foadi@googlegroups.com

We present a comprehensive, configurable open-source framework for estimating the rate of electromagnetic detection of kilonovae (KNe) associated with gravitational wave detections of binary neutron star (BNS) mergers. We simulate the current LIGO-Virgo-KAGRA (LVK) observing run (O4) using up-to-date sensitivity and up-time values as well as the next observing run (O5) using predicted sensitivities. We find the number of discoverable kilonovae during LVK O4 to be ${ 1}_{- 1}^{+ 4}$ or ${ 2 }_{- 2 }^{+ 3 }$, (at 90% confidence) depending on the distribution of NS masses in coalescing binaries, with the number increasing by an order of magnitude during O5 to ${ 19 }_{- 11 }^{+ 24 }$. Regardless of mass model, we predict at most five detectable KNe (at 95% confidence) in O4. We also produce optical and near-infrared light curves that correspond to the physical properties of each merging system. We have collated important information for allocating observing resources and directing search and follow-up observations including distributions of peak magnitudes in several broad bands and timescales for which specific facilities can detect each KN. The framework is easily adaptable, and new simulations can quickly be produced as input information such as merger rates and NS mass distributions are refined. Finally, we compare our suite of simulations to the thus-far completed portion of O4 (as of October 14, 2023), finding a median number of discoverable KNe of 0 and a 95-percentile upper limit of 2, consistent with no detection so far in O4., Comment: 16 pages, 13 figures, submitted to MNRAS

Modern large-scale physics experiments create datasets with sizes and streaming rates that can exceed those from industry leaders such as Google Cloud and Netflix. Fully processing these datasets requires both sufficient compute power and efficient workflows. Recent advances in Machine Learning (ML) and Artificial Intelligence (AI) can either improve or replace existing domain-specific algorithms to increase workflow efficiency. Not only can these algorithms improve the physics performance of current algorithms, but they can often be executed more quickly, especially when run on coprocessors such as GPUs or FPGAs. In the winter of 2023, MIT hosted the Accelerating Physics with ML at MIT workshop, which brought together researchers from gravitational-wave physics, multi-messenger astrophysics, and particle physics to discuss and share current efforts to integrate ML tools into their workflows. The following white paper highlights examples of algorithms and computing frameworks discussed during this workshop and summarizes the expected computing needs for the immediate future of the involved fields., Comment: Whitepaper resulting from Accelerating Physics with ML@MIT workshop in Jan/Feb 2023

Hot DA white dwarfs have fully radiative pure hydrogen atmospheres that are the least complicated to model. Pulsationally stable, they are fully characterized by their effective temperature Teff, and surface gravity log g, which can be deduced from their optical spectra and used in model atmospheres to predict their spectral energy distribution (SED). Based on this, three bright DAWDs have defined the spectrophotometric flux scale of the CALSPEC system of HST. In this paper we add 32 new fainter (16.5 < V < 19.5) DAWDs spread over the whole sky and within the dynamic range of large telescopes. Using ground based spectra and panchromatic photometry with HST/WFC3, a new hierarchical analysis process demonstrates consistency between model and observed fluxes above the terrestrial atmosphere to < 0.004 mag rms from 2700 {\AA} to 7750 {\AA} and to 0.008 mag rms at 1.6{\mu}m for the total set of 35 DAWDs. These DAWDs are thus established as spectrophotometric standards with unprecedented accuracy from the near ultraviolet to the near-infrared, suitable for both ground and space based observatories. They are embedded in existing surveys like SDSS, PanSTARRS and GAIA, and will be naturally included in the LSST survey by Rubin Observatory. With additional data and analysis to extend the validity of their SEDs further into the IR, these spectrophotometric standard stars could be used for JWST, as well as for the Roman and Euclid observatories., Comment: Accepted for publication in Astrophysical Journal. Corrected error in Table 10 and associated Fig 7 in which RP and BP values for the 3 CALSPEC standards had been transposed

This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5 to 6 images per night per filter in $g$, $r$, $i$, and/or $z$, have become publicly available (the proprietary period for this program is waived). We describe the real-time difference-image pipeline and how alerts are distributed to brokers via the same distribution system as the Zwicky Transient Facility (ZTF). In this paper, we focus on the two extragalactic deep fields (COSMOS and ELAIS-S1), characterizing the detected sources and demonstrating that the survey design is effective for probing the discovery space of faint and fast variable and transient sources. We describe and make publicly available 4413 calibrated light curves based on difference-image detection photometry of transients and variables in the extragalactic fields. We also present preliminary scientific analysis regarding Solar System small bodies, stellar flares and variables, Galactic anomaly detection, fast-rising transients and variables, supernovae, and active galactic nuclei., Comment: 22 pages, 17 figures, 2 tables. Accepted to MNRAS

The local distance ladder estimate of the Hubble constant ($H_0$) is important in cosmology, given the recent tension with the early universe inference. We estimate $H_0$ from the Type Ia supernova (SN~Ia) distance ladder, inferring SN~Ia distances with the hierarchical Bayesian SED model, BayeSN. This method has a notable advantage of being able to continuously model the optical and near-infrared (NIR) SN~Ia light curves simultaneously. We use two independent distance indicators, Cepheids or the tip of the red giant branch (TRGB), to calibrate a Hubble-flow sample of 67 SNe~Ia with optical and NIR data. We estimate $H_0 = 74.82 \pm 0.97$ (stat) $\pm\, 0.84$ (sys) km\,s$^{-1}$\,Mpc$^{-1}$ when using the calibration with Cepheid distances to 37 host galaxies of 41 SNe~Ia, and $70.92 \pm 1.14$ (stat) $\pm\,1.49$ (sys) km\,s$^{-1}$\,Mpc$^{-1}$ when using the calibration with TRGB distances to 15 host galaxies of 18 SNe~Ia. For both methods, we find a low intrinsic scatter $\sigma_{\rm int} \lesssim 0.1$ mag. We test various selection criteria and do not find significant shifts in the estimate of $H_0$. Simultaneous modelling of the optical and NIR yields up to $\sim$15\% reduction in $H_0$ uncertainty compared to the equivalent optical-only cases. With improvements expected in other rungs of the distance ladder, leveraging joint optical-NIR SN~Ia data can be critical to reducing the $H_0$ error budget., Comment: 10 pages, 8 figures, published in MNRAS

We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian framework for analyzing SN Ia siblings. At its core is the cosmology-independent relative intrinsic scatter parameter, $\sigma_{Rel}$: the dispersion of siblings distance estimates relative to one another within a galaxy. It quantifies the contribution towards the total intrinsic scatter, $\sigma_0$, from within-galaxy variations about the siblings' common properties. It also affects the combined-distance uncertainty. We present analytic formulae for computing a $\sigma_{Rel}$-posterior from individual siblings distances (estimated using any SN-model). Applying a newly trained BayeSN model, we fit the light curves of each sibling in NGC 3147 individually, to yield consistent distance estimates. However, the wide $\sigma_{Rel}$-posterior means $\sigma_{Rel}\approx\sigma_0$ is not ruled out. We thus combine the distances by marginalizing over $\sigma_{Rel}$ with an informative prior: $\sigma_{Rel}\sim U(0,\sigma_0)$. Simultaneously fitting the trio's light curves improves constraints on distance, and each sibling's individual dust parameters, compared to individual fits. Higher correlation also tightens dust parameter constraints. Therefore, $\sigma_{Rel}$-marginalization yields robust estimates of siblings distances for cosmology, and dust parameters for siblings-host correlation studies. Incorporating NGC 3147's Cepheid-distance yields $H_0=78.4\pm 6.5\,$km/s/Mpc. Our work motivates analyses of homogeneous siblings samples, to constrain $\sigma_{Rel}$, and its SN-model dependence., Comment: Accepted for publication in ApJ. 22 pages, 11 figures, 8 tables. Single-galaxy siblings analysis

We present multi-band ATLAS photometry for SN 2019tsf, a stripped-envelope Type Ib supernova (SESN). The SN shows a triple-peaked light curve and a late (re-)brightening, making it unique among stripped-envelope systems. The re-brightening observations represent the latest photometric measurements of a multi-peaked Type Ib SN to date. As late-time photometry and spectroscopy suggest no hydrogen, the potential circumstellar material (CSM) must be H-poor. Moreover, late (>150 days) spectra show no signs of narrow emission lines, further disfavouring CSM interaction. On the contrary, an extended CSM structure is seen through a follow-up radio campaign with Karl G. Jansky Very Large Array (VLA), indicating a source of bright optically thick radio emission at late times, which is highly unusual among H-poor SESNe. We attribute this phenomenology to an interaction of the supernova ejecta with spherically-asymmetric CSM, potentially disk-like, and we present several models that can potentially explain the origin of this rare Type Ib supernova. The warped disc model paints a novel picture, where the tertiary companion perturbs the progenitors CSM, that can explain the multi-peaked light curves of SNe, and here we apply it to SN 2019tsf. This SN 2019tsf is likely a member of a new sub-class of Type Ib SNe and among the recently discovered class of SNe that undergo mass transfer at the moment of explosion, Comment: 23 pages, Comments are welcome, Submitted to ApJ

As we observe a rapidly growing number of astrophysical transients, we learn more about the diverse host galaxy environments in which they occur. Host galaxy information can be used to purify samples of cosmological Type Ia supernovae, uncover the progenitor systems of individual classes, and facilitate low-latency follow-up of rare and peculiar explosions. In this work, we develop a novel data-driven methodology to simulate the time-domain sky that includes detailed modeling of the probability density function for multiple transient classes conditioned on host galaxy magnitudes, colours, star formation rates, and masses. We have designed these simulations to optimize photometric classification and analysis in upcoming large synoptic surveys. We integrate host galaxy information into the SNANA simulation framework to construct the Simulated Catalogue of Optical Transients and Correlated Hosts (SCOTCH), a publicly-available catalogue of 5 million idealized transient light curves in LSST passbands and their host galaxy properties over the redshift range $0

This invited Snowmass 2021 White Paper highlights the power of joint-analysis of astronomical transients in advancing HEP Science and presents research activities that can realize the opportunities that come with current and upcoming projects. Transients of interest include gravitational wave events, neutrino events, strongly-lensed quasars and supernovae, and Type~Ia supernovae specifically. These transients can serve as probes of cosmological distances in the Universe and as cosmic laboratories of extreme strong-gravity, high-energy physics. Joint analysis refers to work that requires significant coordination from multiple experiments or facilities so encompasses Multi-Messenger Astronomy and optical transient discovery and distributed follow-up programs., Comment: Minor updates to align with the feedback from the Snowmass Community Summer Study Workshop

Methods based on machine learning have recently made substantial inroads in many corners of cosmology. Through this process, new computational tools, new perspectives on data collection, model development, analysis, and discovery, as well as new communities and educational pathways have emerged. Despite rapid progress, substantial potential at the intersection of cosmology and machine learning remains untapped. In this white paper, we summarize current and ongoing developments relating to the application of machine learning within cosmology and provide a set of recommendations aimed at maximizing the scientific impact of these burgeoning tools over the coming decade through both technical development as well as the fostering of emerging communities., Comment: Contribution to Snowmass 2021. 32 pages