Showing 1–50 of 93 results for author: Peek, J E G

Designing an Evaluation Framework for Large Language Models in Astronomy Research

Authors: John F. Wu, Alina Hyk, Kiera McCormick, Christine Ye, Simone Astarita, Elina Baral, Jo Ciuca, Jesse Cranney, Anjalie Field, Kartheik Iyer, Philipp Koehn, Jenn Kotler, Sandor Kruk, Michelle Ntampaka, Charles O'Neill, Joshua E. G. Peek, Sanjib Sharma, Mikaeel Yunus

Abstract: Large Language Models (LLMs) are shifting how scientific research is done. It is imperative to understand how researchers interact with these models and how scientific sub-communities like astronomy might benefit from them. However, there is currently no standard for evaluating the use of LLMs in astronomy. Therefore, we present the experimental design for an evaluation study on how astronomy rese… ▽ More Large Language Models (LLMs) are shifting how scientific research is done. It is imperative to understand how researchers interact with these models and how scientific sub-communities like astronomy might benefit from them. However, there is currently no standard for evaluating the use of LLMs in astronomy. Therefore, we present the experimental design for an evaluation study on how astronomy researchers interact with LLMs. We deploy a Slack chatbot that can answer queries from users via Retrieval-Augmented Generation (RAG); these responses are grounded in astronomy papers from arXiv. We record and anonymize user questions and chatbot answers, user upvotes and downvotes to LLM responses, user feedback to the LLM, and retrieved documents and similarity scores with the query. Our data collection method will enable future dynamic evaluations of LLM tools for astronomy. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 7 pages, 3 figures. Code available at https://github.com/jsalt2024-evaluating-llms-for-astronomy/astro-arxiv-bot

PHANGS-ML: dissecting multiphase gas and dust in nearby galaxies using machine learning

Authors: Dalya Baron, Karin M. Sandstrom, Erik Rosolowsky, Oleg V. Egorov, Ralf S. Klessen, Adam K. Leroy, Médéric Boquien, Eva Schinnerer, Francesco Belfiore, Brent Groves, Jérémy Chastenet, Daniel A. Dale, Guillermo A. Blanc, José E. Méndez-Delgado, Eric W. Koch, Kathryn Grasha, Mélanie Chevance, David A. Thilker, Dario Colombo, Thomas G. Williams, Debosmita Pathak, Jessica Sutter, Toby Brown, John F. Wu, J. E. G. Peek , et al. (3 additional authors not shown)

Abstract: The PHANGS survey uses ALMA, HST, VLT, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multi-wavelength dataset makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here we use unsupervised machine learning algorit… ▽ More The PHANGS survey uses ALMA, HST, VLT, and JWST to obtain an unprecedented high-resolution view of nearby galaxies, covering millions of spatially independent regions. The high dimensionality of such a diverse multi-wavelength dataset makes it challenging to identify new trends, particularly when they connect observables from different wavelengths. Here we use unsupervised machine learning algorithms to mine this information-rich dataset to identify novel patterns. We focus on three of the PHANGS-JWST galaxies, for which we extract properties pertaining to their stellar populations; warm ionized and cold molecular gas; and Polycyclic Aromatic Hydrocarbons (PAHs), as measured over 150 pc-scale regions. We show that we can divide the regions into groups with distinct multiphase gas and PAH properties. In the process, we identify previously-unknown galaxy-wide correlations between PAH band and optical line ratios and use our identified groups to interpret them. The correlations we measure can be naturally explained in a scenario where the PAHs and the ionized gas are exposed to different parts of the same radiation field that varies spatially across the galaxies. This scenario has several implications for nearby galaxies: (i) The uniform PAH ionized fraction on 150 pc scales suggests significant self-regulation in the ISM, (ii) the PAH 11.3/7.7 \mic~ band ratio may be used to constrain the shape of the non-ionizing far-ultraviolet to optical part of the radiation field, and (iii) the varying radiation field affects line ratios that are commonly used as PAH size diagnostics. Neglecting this effect leads to incorrect or biased PAH sizes. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Main results in figures 6 and 12. Submitted to ApJ, and comments are welcome!

A Galactic Eclipse: The Small Magellanic Cloud is Forming Stars in Two, Superimposed Systems

Authors: Claire E. Murray, Sten Hasselquist, Joshua E. G. Peek, Christina Willecke Lindberg, Andres Almeida, Yumi Choi, Jessica E. M. Craig, Helga Denes, John M. Dickey, Enrico M. Di Teodoro, Christoph Federrath, Isabella A. Gerrard, Steven J. Gibson, Denis Leahy, Min-Young Lee, Callum Lynn, Yik Ki Ma, Antoine Marchal, N. M. McClure-Griffiths, David Nidever, Hiep Nguyen, Nickolas M. Pingel, Elizabeth Tarantino, Lucero Uscanga, Jacco Th. van Loon

Abstract: The structure and dynamics of the star-forming disk of the Small Magellanic Cloud (SMC) have long confounded us. The SMC is widely used as a prototype for galactic physics at low metallicity, and yet we fundamentally lack an understanding of the structure of its interstellar medium (ISM). In this work, we present a new model for the SMC by comparing the kinematics of young, massive stars with the… ▽ More The structure and dynamics of the star-forming disk of the Small Magellanic Cloud (SMC) have long confounded us. The SMC is widely used as a prototype for galactic physics at low metallicity, and yet we fundamentally lack an understanding of the structure of its interstellar medium (ISM). In this work, we present a new model for the SMC by comparing the kinematics of young, massive stars with the structure of the ISM traced by high-resolution observations of neutral atomic hydrogen (HI) from the Galactic Australian Square Kilometer Array Pathfinder survey (GASKAP-HI). Specifically, we identify thousands of young, massive stars with precise radial velocity constraints from the Gaia and APOGEE surveys and match these stars to the ISM structures in which they likely formed. By comparing the average dust extinction towards these stars, we find evidence that the SMC is composed of two structures with distinct stellar and gaseous chemical compositions. We construct a simple model that successfully reproduces the observations and shows that the ISM of the SMC is arranged into two, superimposed, star-forming systems with similar gas mass separated by ~5 kpc along the line of sight. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: ApJ accepted. 20 pages, 18 figures

Constructing Impactful Machine Learning Research for Astronomy: Best Practices for Researchers and Reviewers

Authors: D. Huppenkothen, M. Ntampaka, M. Ho, M. Fouesneau, B. Nord, J. E. G. Peek, M. Walmsley, J. F. Wu, C. Avestruz, T. Buck, M. Brescia, D. P. Finkbeiner, A. D. Goulding, T. Kacprzak, P. Melchior, M. Pasquato, N. Ramachandra, Y. -S. Ting, G. van de Ven, S. Villar, V. A. Villar, E. Zinger

Abstract: Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best pr… ▽ More Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best practices, challenges, and drawbacks, which, at present, are often reported on incompletely in the astrophysical literature. With this paper, we aim to provide a primer to the astronomical community, including authors, reviewers, and editors, on how to implement machine learning models and report their results in a way that ensures the accuracy of the results, reproducibility of the findings, and usefulness of the method. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 14 pages, 3 figures; submitted to the Bulletin of the American Astronomical Society

A panoptic view of the Taurus molecular cloud I. The cloud dynamics revealed by gas emission and 3D dust

Authors: J. D. Soler, C. Zucker, J. E. G. Peek, M. Heyer, P. F. Goldsmith, S. C. O. Glover, S. Molinari, R. S. Klessen, P. Hennebelle, L. Testi, T. Colman, M. Benedettini, D. Elia, C. Mininni, S. Pezzuto, E. Schisano, A. Traficante

Abstract: We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) metho… ▽ More We present a study of the three-dimensional (3D) distribution of interstellar dust derived from stellar extinction observations toward the Taurus molecular cloud (MC) and its relation with the neutral atomic hydrogen (HI) emission at 21 cm wavelength and the carbon monoxide $^{12}$CO and $^{13}$CO emission in the $J=1\rightarrow0$ transition. We used the histogram of oriented gradients (HOG) method to match the morphology in a 3D reconstruction of the dust density (3D dust) and the distribution of the gas tracers' emission. The result of the HOG analysis is a map of the relationship between the distances and radial velocities. The HOG comparison between the 3D dust and the HI emission indicates a morphological match at the distance of Taurus but an anti-correlation between the dust density and the HI emission, which uncovers a significant amount of cold HI within the Taurus MC. The HOG between the 3D dust and $^{12}$CO reveals a pattern in radial velocities and distances that is consistent with converging motions of the gas in the Taurus MC, with the near side of the cloud moving at higher velocities and the far side moving at lower velocities. This convergence of flows is likely triggered by the large-scale gas compression caused by the interaction of the Local Bubble and the Per-Tau shell, with Taurus lying at the intersection of the two bubble surfaces. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 18 pages, 21 figures. Accepted for publication in Astronomy & Astrophysics (22MAY2023)

Journal ref: A&A 675, A206 (2023)

Identification of galaxy shreds in large photometric catalogs using Convolutional Neural Networks

Authors: Enrico M. Di Teodoro, Josh E. G. Peek, John F. Wu

Abstract: Contamination from galaxy fragments, identified as sources, is a major issue in large photometric galaxy catalogs. In this paper, we prove that this problem can be easily addressed with computer vision techniques. We use image cutouts to train a convolutional neural network (CNN) to identify catalogued sources that are in reality just star formation regions and/or shreds of larger galaxies. The CN… ▽ More Contamination from galaxy fragments, identified as sources, is a major issue in large photometric galaxy catalogs. In this paper, we prove that this problem can be easily addressed with computer vision techniques. We use image cutouts to train a convolutional neural network (CNN) to identify catalogued sources that are in reality just star formation regions and/or shreds of larger galaxies. The CNN reaches an accuracy ~98% on our testing datasets. We apply this CNN to galaxy catalogs from three amongst the largest surveys available today: the Sloan Digital Sky Survey (SDSS), the DESI Legacy Imaging Surveys and the Panoramic Survey Telescope and Rapid Response System Survey (Pan-STARSS). We find that, even when strict selection criteria are used, all catalogs still show a ~5% level of contamination from galaxy shreds. Our CNN gives a simple yet effective solution to clean galaxy catalogs from these contaminants. △ Less

Submitted 23 January, 2023; originally announced January 2023.

Comments: Accepted for publication in AJ

Target Selection and Sample Characterization for the DESI LOW-Z Secondary Target Program

Authors: Elise Darragh-Ford, John F. Wu, Yao-Yuan Mao, Risa H. Wechsler, Marla Geha, Jaime E. Forero-Romero, ChangHoon Hahn, Nitya Kallivayalil, John Moustakas, Ethan O. Nadler, Marta Nowotka, J. E. G. Peek, Erik J. Tollerud, Benjamin Weiner, J. Aguilar, S. Ahlen, D. Brooks, A. P. Cooper, A. de la Macorra, A. Dey, K. Fanning, A. Font-Ribera, S. Gontcho A Gontcho, K. Honscheid, T. Kisner , et al. (17 additional authors not shown)

Abstract: We introduce the DESI LOW-Z Secondary Target Survey, which combines the wide-area capabilities of the Dark Energy Spectroscopic Instrument (DESI) with an efficient, low-redshift target selection method. Our selection consists of a set of color and surface brightness cuts, combined with modern machine learning methods, to target low-redshift dwarf galaxies ($z$ < 0.03) between $19 < r < 21$ with hi… ▽ More We introduce the DESI LOW-Z Secondary Target Survey, which combines the wide-area capabilities of the Dark Energy Spectroscopic Instrument (DESI) with an efficient, low-redshift target selection method. Our selection consists of a set of color and surface brightness cuts, combined with modern machine learning methods, to target low-redshift dwarf galaxies ($z$ < 0.03) between $19 < r < 21$ with high completeness. We employ a convolutional neural network (CNN) to select high-priority targets. The LOW-Z survey has already obtained over 22,000 redshifts of dwarf galaxies (M$_* < 10^9$ M$_\odot$), comparable to the number of dwarf galaxies discovered in SDSS-DR8 and GAMA. As a spare fiber survey, LOW-Z currently receives fiber allocation for just ~50% of its targets. However, we estimate that our selection is highly complete: for galaxies at $z < 0.03$ within our magnitude limits, we achieve better than 95% completeness with ~1% efficiency using catalog-level photometric cuts. We also demonstrate that our CNN selections $z<0.03$ galaxies from the photometric cuts subsample at least ten times more efficiently while maintaining high completeness. The full five-year DESI program will expand the LOW-Z sample, densely mapping the low-redshift Universe, providing an unprecedented sample of dwarf galaxies, and providing critical information about how to pursue effective and efficient low-redshift surveys. △ Less

Submitted 12 June, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 24 pages, 14 figures, data to reproduce figures: https://zenodo.org/record/7422591

Measuring the 8621 Å Diffuse Interstellar Band in Gaia DR3 RVS Spectra: Obtaining a Clean Catalog by Marginalizing over Stellar Types

Authors: Andrew K. Saydjari, Ana Sofía M. Uzsoy, Catherine Zucker, J. E. G. Peek, Douglas P. Finkbeiner

Abstract: Diffuse interstellar bands (DIBs) are broad absorption features associated with interstellar dust and can serve as chemical and kinematic tracers. Conventional measurements of DIBs in stellar spectra are complicated by residuals between observations and best-fit stellar models. To overcome this, we simultaneously model the spectrum as a combination of stellar, dust, and residual components, with f… ▽ More Diffuse interstellar bands (DIBs) are broad absorption features associated with interstellar dust and can serve as chemical and kinematic tracers. Conventional measurements of DIBs in stellar spectra are complicated by residuals between observations and best-fit stellar models. To overcome this, we simultaneously model the spectrum as a combination of stellar, dust, and residual components, with full posteriors on the joint distribution of the components. This decomposition is obtained by modeling each component as a draw from a high-dimensional Gaussian distribution in the data-space (the observed spectrum) -- a method we call "Marginalized Analytic Data-space Gaussian Inference for Component Separation" (MADGICS). We use a data-driven prior for the stellar component, which avoids missing stellar features not well-modeled by synthetic spectra. This technique provides statistically rigorous uncertainties and detection thresholds, which are required to work in the low signal-to-noise regime that is commonplace for dusty lines of sight. We reprocess all public Gaia DR3 RVS spectra and present an improved 8621 Å DIB catalog, free of detectable stellar line contamination. We constrain the rest-frame wavelength to $8623.14 \pm 0.087$ Å (vacuum), find no significant evidence for DIBs in the Local Bubble from the $1/6^{\rm{th}}$ of RVS spectra that are public, and show unprecedented correlation with kinematic substructure in Galactic CO maps. We validate the catalog, its reported uncertainties, and biases using synthetic injection tests. We believe MADGICS provides a viable path forward for large-scale spectral line measurements in the presence of complex spectral contamination. △ Less

Submitted 30 August, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 28 pages, 27 figures, accepted at ApJ

Disconnecting the Dots: Re-Examining the Nature of Stellar "Strings" in the Milky Way

Authors: Catherine Zucker, J. E. G. Peek, Sarah R. Loebman

Abstract: Recent analyses of Gaia data have resulted in the identification of new stellar structures, including a new class of extended stellar filaments called stellar "strings", first proposed by Kounkel and Covey 2019. We explore the spatial, kinematic, and chemical composition of strings to demonstrate that these newfound structures are largely inconsistent with being physical objects whose members shar… ▽ More Recent analyses of Gaia data have resulted in the identification of new stellar structures, including a new class of extended stellar filaments called stellar "strings", first proposed by Kounkel and Covey 2019. We explore the spatial, kinematic, and chemical composition of strings to demonstrate that these newfound structures are largely inconsistent with being physical objects whose members share a common origin. Examining the 3D spatial distribution of string members, we find that the spatial dispersion around the claimed string spine does not improve in the latest Gaia DR3 data release -- despite tangible gains in the signal-to-noise of the parallax measurements -- counter to expectations of a bona fide structure. Using the radial velocity dispersion of the strings (averaging $\rmσ_{V_r} = 16\;km\;s^{-1}$) to estimate their virial masses, we find that all strings are gravitationally unbound. Given the finding that the strings are dispersing, the reported stellar ages of the strings are typically $120\times$ larger than their measured dispersal times. Finally, we validate prior work that stellar strings are more chemically homogeneous than their local field stars, but show it is possible to obtain the same signatures of chemical homogeneity by drawing random samples of stars from spatially, temporally, and kinematically unrelated open clusters. Our results show that while some strings may be composed of real sub-structures, there is no consistent evidence for larger string-like connections over the sample. These results underline the need for caution in over-interpreting the significance of these strings and their role in understanding the star formation history of the Milky Way. △ Less

Submitted 12 August, 2022; v1 submitted 27 May, 2022; originally announced May 2022.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: ApJ 936 160 (2022)

The emptiness inside: Finding gaps, valleys, and lacunae with geometric data analysis

Authors: Gabriella Contardo, David W. Hogg, Jason A. S. Hunt, Joshua E. G. Peek, Yen-Chi Chen

Abstract: Discoveries of gaps in data have been important in astrophysics. For example, there are kinematic gaps opened by resonances in dynamical systems, or exoplanets of a certain radius that are empirically rare. A gap in a data set is a kind of anomaly, but in an unusual sense: Instead of being a single outlier data point, situated far from other data points, it is a region of the space, or a set of po… ▽ More Discoveries of gaps in data have been important in astrophysics. For example, there are kinematic gaps opened by resonances in dynamical systems, or exoplanets of a certain radius that are empirically rare. A gap in a data set is a kind of anomaly, but in an unusual sense: Instead of being a single outlier data point, situated far from other data points, it is a region of the space, or a set of points, that is anomalous compared to its surroundings. Gaps are both interesting and hard to find and characterize, especially when they have non-trivial shapes. We present in this paper a statistic that can be used to estimate the (local) "gappiness" of a point in the data space. It uses the gradient and Hessian of the density estimate (and thus requires a twice-differentiable density estimator). This statistic can be computed at (almost) any point in the space and does not rely on optimization; it allows to highlight under-dense regions of any dimensionality and shape in a general and efficient way. We illustrate our method on the velocity distribution of nearby stars in the Milky Way disk plane, which exhibits gaps that could originate from different processes. Identifying and characterizing those gaps could help determine their origins. We provide in an Appendix implementation notes and additional considerations for finding under-densities in data, using critical points and the properties of the Hessian of the density. △ Less

Submitted 5 September, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 17 pages, 10 figures. Submitted to AJ. Comments welcomed. Revision: added 3D gridding + restructured outline: implementation notes (Quadratic Kernel) and methods for approx critical points and 1d-valley now in Annex

Burton's Curse: The Impact of Bulk Flows on the Galactic Longitude-Velocity Diagram and the Illusion of a Continuous Perseus Arm

Authors: J. E. G. Peek, Kirill Tchernyshyov, Marc-Antoine Miville-Deschenes

Abstract: In this work we demonstrate that the Perseus Arm is not a continuous structure of molecular gas in the second quadrant. We first show that the observed, distanced-resolved velocity structure of the Galaxy in the outer disk is capable of creating illusory spiral arms, as was first proposed by Burton (1971). Second, we measure the distances to a collection of CO clouds at velocities consistent with… ▽ More In this work we demonstrate that the Perseus Arm is not a continuous structure of molecular gas in the second quadrant. We first show that the observed, distanced-resolved velocity structure of the Galaxy in the outer disk is capable of creating illusory spiral arms, as was first proposed by Burton (1971). Second, we measure the distances to a collection of CO clouds at velocities consistent with the Perseus arm with $135^\circ < l < 160^\circ$. We find these distances using 3D dust maps from Green et al. (2019). We determine that these molecular cloud do not preferentially lie at the distance of a purported Perseus arm, but rather extend over 3 kpc in distance, with some evidence for a closer, high pitch angle structure between 1 and 1.5 kpc away. Finally, we demonstrate that velocity perturbations of the amplitude found near the Perseus arm can wreak havoc on our interpretation of the longitude-velocity diagram for more than half of the Milky Way disk. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 14 pages, 7 figures. Accepted to the Astrophysical Journal

Extending the SAGA Survey (xSAGA) I: Satellite Radial Profiles as a Function of Host Galaxy Properties

Authors: John F. Wu, J. E. G. Peek, Erik J. Tollerud, Yao-Yuan Mao, Ethan O. Nadler, Marla Geha, Risa H. Wechsler, Nitya Kallivayalil, Benjamin J. Weiner

Abstract: We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-$z$ galaxies on the basis of optical imaging, and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify $z < 0.03$ gal… ▽ More We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-$z$ galaxies on the basis of optical imaging, and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify $z < 0.03$ galaxies from more distant objects using image cutouts from the DESI Legacy Imaging Surveys. From the sample of $> 100,000$ CNN-selected low-$z$ galaxies, we identify $>20,000$ probable satellites located between 36-300 projected kpc from NASA-Sloan Atlas central galaxies in the stellar mass range $9.5 < \log(M_\star/M_\odot) < 11$. We characterize the incompleteness and contamination for CNN-selected samples, and apply corrections in order to estimate the true number of satellites as a function of projected radial distance from their hosts. Satellite richness depends strongly on host stellar mass, such that more massive host galaxies have more satellites, and on host morphology, such that elliptical hosts have more satellites than disky hosts with comparable stellar masses. We also find a strong inverse correlation between satellite richness and the magnitude gap between a host and its brightest satellite. The normalized satellite radial distribution between 36-300 kpc does not depend strongly on host stellar mass, morphology, or magnitude gap. The satellite abundances and radial distributions we measure are in reasonable agreement with predictions from hydrodynamic simulations. Our results deliver unprecedented statistical power for studying satellite galaxy populations, and highlight the promise of using machine learning for extending galaxy samples of wide-area surveys. △ Less

Submitted 25 January, 2022; v1 submitted 2 December, 2021; originally announced December 2021.

Comments: 27 pages, 18 figures, 1 table. Accepted to ApJ. Code available at https://github.com/jwuphysics/xsaga

Journal ref: ApJ (2022), 927, 121

Radial motions and radial gas flows in local spiral galaxies

Authors: Enrico M. Di Teodoro, Josh E. G. Peek

Abstract: We determine radial velocities and mass flow rates in a sample of 54 local spiral galaxies by modelling high-resolution and high-sensitivity data of the atomic hydrogen emission line. We found that, although radial inflow motions seem to be slightly preferred over outflow motions, their magnitude is generally small. Most galaxies show radial flows of only a few km/s throughout their HI disks, eith… ▽ More We determine radial velocities and mass flow rates in a sample of 54 local spiral galaxies by modelling high-resolution and high-sensitivity data of the atomic hydrogen emission line. We found that, although radial inflow motions seem to be slightly preferred over outflow motions, their magnitude is generally small. Most galaxies show radial flows of only a few km/s throughout their HI disks, either inwards or outwards, without any clear increase in magnitude in the outermost regions, as we would expect for continuous radial accretion. Gas mass flow rates for most galaxies are less than 1 M$_\odot$/yr. Over the entire sample, we estimated an average inflow rate of 0.3 M$_\odot$/yr outside the optical disk and of 0.1 M$_\odot$/yr in the outskirts of the HI disks. These inflow rates are about 5-10 times smaller than the average star formation rate of 1.4 M$_\odot$/yr. Our study suggests that there is no clear evidence for systematic radial accretion inflows that alone could feed and sustain the star formation process in the inner regions of local spiral galaxies at its current rate. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: Accepted for publications in ApJ. 12 pages, 6 figures + appendixes

On the Three-Dimensional Structure of Local Molecular Clouds

Authors: Catherine Zucker, Alyssa Goodman, João Alves, Shmuel Bialy, Eric W. Koch, Joshua S. Speagle, Michael M. Foley, Douglas Finkbeiner, Reimar Leike, Torsten Enßlin, Joshua E. G. Peek, Gordian Edenhofer

Abstract: We leverage the 1 pc spatial resolution of the Leike et al. 2020 3D dust map to characterize the three-dimensional structure of nearby molecular clouds ($d \lesssim 400$ pc). We start by "skeletonizing" the clouds in 3D volume density space to determine their "spines," which we project on the sky to constrain cloud distances with $\approx 1\%$ uncertainty. For each cloud, we determine an average r… ▽ More We leverage the 1 pc spatial resolution of the Leike et al. 2020 3D dust map to characterize the three-dimensional structure of nearby molecular clouds ($d \lesssim 400$ pc). We start by "skeletonizing" the clouds in 3D volume density space to determine their "spines," which we project on the sky to constrain cloud distances with $\approx 1\%$ uncertainty. For each cloud, we determine an average radial volume density profile around its 3D spine and fit the profiles using Gaussian and Plummer functions. The radial volume density profiles are well-described by a two-component Gaussian function, consistent with clouds having broad, lower-density outer envelopes and narrow, higher-density inner layers. The ratio of the outer to inner envelope widths is $\approx 3:1$. We hypothesize that these two components may be tracing a transition between atomic and diffuse molecular gas or between the unstable and cold neutral medium. Plummer-like models can also provide a good fit, with molecular clouds exhibiting shallow power-law wings with density, $n$, falling off like $n^{-2}$ at large radii. Using Bayesian model selection, we find that parameterizing the clouds' profiles using a single Gaussian is disfavored. We compare our results with 2D dust extinction maps, finding that the 3D dust recovers the total cloud mass from integrated approaches with fidelity, deviating only at higher levels of extinction ($A_V \gtrsim 2 - 3$ mag). The 3D cloud structure described here will enable comparisons with synthetic clouds generated in simulations, offering unprecedented insight into the origins and fates of molecular clouds in the interstellar medium. △ Less

Submitted 27 September, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: Published in ApJ. For data behind the figures see https://dataverse.harvard.edu/dataverse/3D_cloud_structure . For interactive gallery of local clouds see https://faun.rc.fas.harvard.edu/czucker/Paper_Figures/3D_Cloud_Topologies/gallery.html . Companion paper to Bialy et al. 2021 (arXiv:2109.09763). For press see https://www.cfa.harvard.edu/news/gigantic-cavity-space-sheds-new-light-how-stars-form

Journal ref: ApJ 919 35 (2021)

Best Practices for Data Publication in the Astronomical Literature

Authors: Tracy X. Chen, Marion Schmitz, Joseph M. Mazzarella, Xiuqin Wu, Julian C. van Eyken, Alberto Accomazzi, Rachel L. Akeson, Mark Allen, Rachael Beaton, G. Bruce Berriman, Andrew W. Boyle, Marianne Brouty, Ben Chan, Jessie L. Christiansen, David R. Ciardi, David Cook, Raffaele D'Abrusco, Rick Ebert, Cren Frayer, Benjamin J. Fulton, Christopher Gelino, George Helou, Calen B. Henderson, Justin Howell, Joyce Kim , et al. (20 additional authors not shown)

Abstract: We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to… ▽ More We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to streamline the extraction, preservation, integration and cross-linking of valuable data from astrophysics literature into major astronomical databases, and consequently facilitate new modes of science discovery that will better exploit the vast quantities of panchromatic and multi-dimensional data associated with the literature. We encourage authors, journal editors, referees, and publishers to implement the best practices reviewed here, as well as related recommendations from international astronomical organizations such as the International Astronomical Union (IAU) for publication of nomenclature, data, and metadata. A convenient Checklist of Recommendations for Publishing Data in the Literature is included for authors to consult before the submission of the final version of their journal articles and associated data files. We recommend that publishers of journals in astronomy and astrophysics incorporate a link to this document in their Instructions to Authors. △ Less

Submitted 16 April, 2022; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 19 pages, 1 figure, 3 tables, accepted for publication in ApJS

The Future Of The Arecibo Observatory: The Next Generation Arecibo Telescope

Authors: D. Anish Roshi, N. Aponte, E. Araya, H. Arce, L. A. Baker, W. Baan, T. M. Becker, J. K. Breakall, R. G. Brown, C. G. M. Brum, M. Busch, D. B. Campbell, T. Cohen, F. Cordova, J. S. Deneva, M. Devogele, T. Dolch, F. O. Fernandez-Rodriguez, T. Ghosh, P. F. Goldsmith, L. I. Gurvits, M. Haynes, C. Heiles, J. W. T. Hessel, D. Hickson , et al. (49 additional authors not shown)

Abstract: The Arecibo Observatory (AO) is a multidisciplinary research and education facility that is recognized worldwide as a leading facility in astronomy, planetary, and atmospheric and space sciences. AO's cornerstone research instrument was the 305-m William E. Gordon telescope. On December 1, 2020, the 305-m telescope collapsed and was irreparably damaged. In the three weeks following the collapse, A… ▽ More The Arecibo Observatory (AO) is a multidisciplinary research and education facility that is recognized worldwide as a leading facility in astronomy, planetary, and atmospheric and space sciences. AO's cornerstone research instrument was the 305-m William E. Gordon telescope. On December 1, 2020, the 305-m telescope collapsed and was irreparably damaged. In the three weeks following the collapse, AO's scientific and engineering staff and the AO users community initiated extensive discussions on the future of the observatory. The community is in overwhelming agreement that there is a need to build an enhanced, next-generation radar-radio telescope at the AO site. From these discussions, we established the set of science requirements the new facility should enable. These requirements can be summarized briefly as: 5 MW of continuous wave transmitter power at 2 - 6 GHz, 10 MW of peak transmitter power at 430 MHz (also at 220MHz under consideration), zenith angle coverage 0 to 48 deg, frequency coverage 0.2 to 30 GHz and increased Field-of-View. These requirements determine the unique specifications of the proposed new instrument. The telescope design concept we suggest consists of a compact array of fixed dishes on a tiltable, plate-like structure with a collecting area equivalent to a 300m dish. This concept, referred to as the Next Generation Arecibo Telescope (NGAT), meets all of the desired specifications and provides significant new science capabilities to all three research groups at AO. This whitepaper presents a sample of the wide variety of the science that can be achieved with the NGAT, the details of the telescope design concept and the need for the new telescope to be located at the AO site. We also discuss other AO science activities that interlock with the NGAT in the white paper. △ Less

Submitted 1 April, 2021; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 82 pages (executive summary 10 pages), 21 figures, Arecibo observatory white paper (Updated with the complete author list and minor edits)

The Gas Content and Stripping of Local Group Dwarf Galaxies

Authors: Mary E. Putman, Yong Zheng, Adrian M. Price-Whelan, Jana Grcevich, Amalya C. Johnson, Erik Tollerud, Joshua E. G. Peek

Abstract: The gas content of the complete compilation of Local Group dwarf galaxies (119 within 2 Mpc) is presented using HI survey data. Within the virial radius of the Milky Way (224 kpc here), 53 of 55 dwarf galaxies are devoid of gas to limits of M$_{\rm HI}<10^4$ M$_\odot$. Within the virial radius of M31 (266 kpc), 27 of 30 dwarf galaxies are devoid of gas (with limits typically $<10^5$ M$_\odot$). Be… ▽ More The gas content of the complete compilation of Local Group dwarf galaxies (119 within 2 Mpc) is presented using HI survey data. Within the virial radius of the Milky Way (224 kpc here), 53 of 55 dwarf galaxies are devoid of gas to limits of M$_{\rm HI}<10^4$ M$_\odot$. Within the virial radius of M31 (266 kpc), 27 of 30 dwarf galaxies are devoid of gas (with limits typically $<10^5$ M$_\odot$). Beyond the virial radii of the Milky Way and M31, the majority of the dwarf galaxies have detected HI gas and have HI masses higher than the limits. When the relationship between gas content and distance is investigated using a Local Group virial radius, more of the non-detected dwarf galaxies are within this radius (85$\pm1$ of the 93 non-detected dwarf galaxies) than within the virial radii of the Milky Way and M31. Using the Gaia proper motion measurements available for 38 dwarf galaxies, the minimum gas density required to completely strip them of gas is calculated. Halo densities between $10^{-5}$ and $5 \times 10^{-4}$ cm$^{-3}$ are typically required for instantaneous stripping at perigalacticon. When compared to halo density with radius expectations from simulations and observations, 80% of the dwarf galaxies with proper motions are consistent with being stripped by ram pressure at Milky Way pericenter. The results suggest a diffuse gaseous galactic halo medium is important in quenching dwarf galaxies, and that a Local Group medium also potentially plays a role. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 26 pages, 7 figures, 3 tables, Accepted by ApJ

Characterizing the Circumgalactic Medium of the Lowest-Mass Galaxies: A Case Study of IC 1613

Authors: Yong Zheng, Andrew Emerick, Mary E. Putman, Jessica K. Werk, Evan N. Kirby, Joshua E. G. Peek

Abstract: Using 10 sightlines observed with the Hubble Space Telescope/Cosmic Origins Spectrograph, we study the circumgalactic medium (CGM) and outflows of IC1613, which is a low-mass ($M_*\sim10^8~M_\odot$), dwarf irregular galaxy on the outskirts of the Local Group. Among the sightlines, 4 are pointed towards UV-bright stars in IC1613, and the other 6 sightlines are background QSOs at impact parameters f… ▽ More Using 10 sightlines observed with the Hubble Space Telescope/Cosmic Origins Spectrograph, we study the circumgalactic medium (CGM) and outflows of IC1613, which is a low-mass ($M_*\sim10^8~M_\odot$), dwarf irregular galaxy on the outskirts of the Local Group. Among the sightlines, 4 are pointed towards UV-bright stars in IC1613, and the other 6 sightlines are background QSOs at impact parameters from 6 kpc ($<0.1R_{200}$) to 61 kpc ($0.6R_{200}$). We detect a number of Si II, Si III, Si IV, C II, and C IV absorbers, most of which have velocities less than the escape velocity of IC1613 and thus are gravitationally bound. The line strengths of these ion absorbers are consistent with the CGM absorbers detected in dwarf galaxies at low redshifts. Assuming that Si II, Si III, and Si IV comprise nearly 100% of the total silicon, we find 3% ($\sim$8$\times$10$^3~{\rm M_\odot}$), 2% ($\sim$7$\times$10$^3~{\rm M_\odot}$), and 32--42% [$\sim$(1.0--1.3)$\times$10$^5~{\rm M_\odot}$] of the silicon mass in the stars, interstellar medium, and within $0.6R_{200}$ of the CGM of IC1613. We also estimate the metal outflow rate to be ${\rm \dot{M}_{out, Z}\geq1.1\times10^{-5}~M_\odot~yr^{-1}}$ and the instantaneous metal mass loading factor to be $η_{\rm Z}\geq0.004$, which are in broad agreement with available observation and simulation values. This work is the first time a dwarf galaxy of such low mass is probed by a number of both QSO and stellar sightlines, and it shows that the CGM of low-mass gas-rich galaxies can be a large reservoir enriched with metals from past and ongoing outflows. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: Resubmitted to ApJ with the referee's comments addressed. Minor changes were made to the original draft

Predicting galaxy spectra from images with hybrid convolutional neural networks

Authors: John F. Wu, J. E. G. Peek

Abstract: Galaxies can be described by features of their optical spectra such as oxygen emission lines, or morphological features such as spiral arms. Although spectroscopy provides a rich description of the physical processes that govern galaxy evolution, spectroscopic data are observationally expensive to obtain. For the first time, we are able to robustly predict galaxy spectra directly from broad-band i… ▽ More Galaxies can be described by features of their optical spectra such as oxygen emission lines, or morphological features such as spiral arms. Although spectroscopy provides a rich description of the physical processes that govern galaxy evolution, spectroscopic data are observationally expensive to obtain. For the first time, we are able to robustly predict galaxy spectra directly from broad-band imaging. We present a powerful new approach using a hybrid convolutional neural network with deconvolution instead of batch normalization; this hybrid CNN outperforms other models in our tests. The learned mapping between galaxy imaging and spectra will be transformative for future wide-field surveys, such as with the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope, by multiplying the scientific returns for spectroscopically-limited galaxy samples. △ Less

Submitted 30 November, 2020; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: 5 pages, 2 figures, accepted to the Machine Learning and the Physical Sciences workshop at NeurIPS 2020. Code available at https://github.com/jwuphysics/predicting-spectra-from-images/

The HI Structure of the Local Volume Dwarf Galaxy Pisces A

Authors: Luca Beale, Jennifer Donovan Meyer, Erik J. Tollerud, Mary E. Putman, J. E. G. Peek

Abstract: Dedicated HI surveys have recently led to a growing category of low-mass galaxies found in the Local Volume. We present synthesis imaging of one such galaxy, Pisces A, a low-mass dwarf originally confirmed via optical imaging and spectroscopy of neutral hydrogen (HI) sources in the Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey. Using HI observations taken with the Karl G. Jansky Very Lar… ▽ More Dedicated HI surveys have recently led to a growing category of low-mass galaxies found in the Local Volume. We present synthesis imaging of one such galaxy, Pisces A, a low-mass dwarf originally confirmed via optical imaging and spectroscopy of neutral hydrogen (HI) sources in the Galactic Arecibo L-band Feed Array HI (GALFA-HI) survey. Using HI observations taken with the Karl G. Jansky Very Large Array (JVLA), we characterize the kinematic structure of the gas and connect it to the galaxy's environment and evolutionary history. While the galaxy shows overall ordered rotation, a number of kinematic features indicate a disturbed gas morphology. These features are suggestive of a tumultuous recent history, and represent $\sim 3.5$% of the total baryonic mass. We find a total baryon fraction $f_{\rm bary} = 0.13$ if we include these features. We also quantify the cosmic environment of Pisces A, finding an apparent alignment of the disturbed gas with nearby, large scale filamentary structure at the edge of the Local Void. We consider several scenarios for the origin of the disturbed gas, including gas stripping via ram pressure or galaxy-galaxy interactions, as well as accretion and ram pressure compression. Though we cannot rule out a past interaction with a companion, our observations best support the suggestion that the neutral gas morphology and recent star formation in Pisces A is a direct result of its interactions with the IGM. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Comments: 16 pages, 12 figures, accepted to ApJ

Extracting the cold neutral medium from HI emission with deep learning: Implications for Galactic foregrounds at high latitude

Authors: Claire E. Murray, J. E. G. Peek, Chang-Goo Kim

Abstract: Resolving the phase structure of neutral hydrogen (HI) is crucial for understanding the life cycle of the interstellar medium (ISM). However, accurate measurements of HI temperature and density are limited by the availability of background continuum sources for measuring HI absorption. Here we test the use of deep learning for extracting HI properties over large areas without optical depth informa… ▽ More Resolving the phase structure of neutral hydrogen (HI) is crucial for understanding the life cycle of the interstellar medium (ISM). However, accurate measurements of HI temperature and density are limited by the availability of background continuum sources for measuring HI absorption. Here we test the use of deep learning for extracting HI properties over large areas without optical depth information. We train a 1D convolutional neural network using synthetic observations of 3D numerical simulations of the ISM to predict the fraction of cold neutral medium (f_CNM) and the correction to the optically-thin HI column density for optical depth (R_HI) from $21\rm\,cm$ emission alone. We restrict our analysis to high Galactic latitudes ($|b|>30 deg), where the complexity of spectral line profiles is minimized. We verify that the network accurately predicts f_CNM and R_HI by comparing the results with direct constraints from 21 cm absorption. By applying the network to the GALFA-HI survey, we generate large-area maps of f_CNM and R_HI. Although the overall contribution to the total HI column of cold neutral medium (CNM)-rich structures is small (~5%), we find that these structures are ubiquitous. Our results are consistent with the picture that small-scale structures observed in 21 cm emission aligned with the magnetic field are dominated by CNM. Finally, we demonstrate that the observed correlation between HI column density and dust reddening (E(B-V)) declines with increasing R_HI, indicating that future efforts to quantify foreground Galactic E(B-V) using HI, even at high latitudes, should increase fidelity by accounting for HI phase structure. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: ApJ accepted

DeepMerge: Classifying High-redshift Merging Galaxies with Deep Neural Networks

Authors: A. Ćiprijanović, G. F. Snyder, B. Nord, J. E. G. Peek

Abstract: We investigate and demonstrate the use of convolutional neural networks (CNNs) for the task of distinguishing between merging and non-merging galaxies in simulated images, and for the first time at high redshifts (i.e. $z=2$). We extract images of merging and non-merging galaxies from the Illustris-1 cosmological simulation and apply observational and experimental noise that mimics that from the H… ▽ More We investigate and demonstrate the use of convolutional neural networks (CNNs) for the task of distinguishing between merging and non-merging galaxies in simulated images, and for the first time at high redshifts (i.e. $z=2$). We extract images of merging and non-merging galaxies from the Illustris-1 cosmological simulation and apply observational and experimental noise that mimics that from the Hubble Space Telescope; the data without noise form a "pristine" data set and that with noise form a "noisy" data set. The test set classification accuracy of the CNN is $79\%$ for pristine and $76\%$ for noisy. The CNN outperforms a Random Forest classifier, which was shown to be superior to conventional one- or two-dimensional statistical methods (Concentration, Asymmetry, the Gini, $M_{20}$ statistics etc.), which are commonly used when classifying merging galaxies. We also investigate the selection effects of the classifier with respect to merger state and star formation rate, finding no bias. Finally, we extract Grad-CAMs (Gradient-weighted Class Activation Mapping) from the results to further assess and interrogate the fidelity of the classification model. △ Less

Submitted 24 April, 2020; originally announced April 2020.

Comments: 17 pages, 8 figures, submitted to Astronomy & Computing

Algorithms and Statistical Models for Scientific Discovery in the Petabyte Era

Authors: Brian Nord, Andrew J. Connolly, Jamie Kinney, Jeremy Kubica, Gautaum Narayan, Joshua E. G. Peek, Chad Schafer, Erik J. Tollerud, Camille Avestruz, G. Jogesh Babu, Simon Birrer, Douglas Burke, João Caldeira, Douglas A. Caldwell, Joleen K. Carlberg, Yen-Chi Chen, Chuanfei Dong, Eric D. Feigelson, V. Zach Golkhou, Vinay Kashyap, T. S. Li, Thomas Loredo, Luisa Lucie-Smith, Kaisey S. Mandel, J. R. Martínez-Galarza , et al. (13 additional authors not shown)

Abstract: The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our… ▽ More The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/). △ Less

Submitted 4 November, 2019; originally announced November 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1905.05116

Report number: FERMILAB-FN-1093-A-AE-SCD

The 3D Kinematics of Gas in the Small Magellanic Cloud

Authors: Claire E. Murray, J. E. G. Peek, Enrico M. Di Teodoro, N. M. McClure-Griffiths, John M. Dickey, Helga Denes

Abstract: We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities… ▽ More We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities match with dominant HI components. We compare the observed radial and tangential velocities of these stars with predictions from the state-of-the-art rotating disk model based on high-resolution 21 cm observations of the SMC from the Australian Square Kilometer Array Pathfinder telescope. We find that the observed kinematics of gas-tracing stars are inconsistent with disk rotation. We conclude that the kinematics of gas in the SMC are more complex than can be inferred from the integrated radial velocity field. As a result of violent tidal interactions with the LMC, non-rotational motions are prevalent throughout the SMC, and it is likely composed of distinct sub-structures overlapping along the line of sight. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 9 pages, 5 figures, 1 Appendix; ApJ accepted

Realizing the potential of astrostatistics and astroinformatics

Authors: Gwendolyn Eadie, Thomas J. Loredo, Ashish A. Mahabal, Aneta Siemiginowska, Eric Feigelson, Eric B. Ford, S. G. Djorgovski, Matthew Graham, Zeljko Ivezic, Kirk Borne, Jessi Cisewski-Kehe, J. E. G. Peek, Chad Schafer, Padma A. Yanamandra-Fisher, C. Alex Young

Abstract: This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy. This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Comments: 14 pages, 1 figure; submitted to the Decadal Survey on Astronomy and Astrophysics (Astro2020) on 10 July 2019; see https://sites.nationalacademies.org/DEPS/Astro2020/DEPS_192906

Small-Scale HI Channel Map Structure is Cold: Evidence from Na I Absorption at High Galactic Latitudes

Authors: J. E. G. Peek, S. E. Clark

Abstract: The spatial distribution of neutral hydrogen (H I) emission is a powerful probe of interstellar medium physics. The small-scale structure in H I channel maps is often assumed to probe the velocity field rather than real density structures. In this work we directly test this assumption, using high-resolution GALFA-H I observations and 50,985 quasar spectra from SDSS. We measure the equivalent width… ▽ More The spatial distribution of neutral hydrogen (H I) emission is a powerful probe of interstellar medium physics. The small-scale structure in H I channel maps is often assumed to probe the velocity field rather than real density structures. In this work we directly test this assumption, using high-resolution GALFA-H I observations and 50,985 quasar spectra from SDSS. We measure the equivalent widths of interstellar Na I D$_1$ and Na I D$_2$ absorption, and robustly conclude that they depend more strongly on the column density of small-scale structure in H I than on either the large-scale H I structure or the total H I column. This is inconsistent with the hypothesis that small-scale channel map structure is driven by velocity crowding. Instead, the data favor the interpretation that this emission structure predominantly originates in cold, dense interstellar material, consistent with a clumpy cold neutral medium. △ Less

Submitted 20 September, 2019; originally announced September 2019.

Comments: Submitted to the ApJL

Robust Archives Maximize Scientific Accessibility

Authors: J. E. G. Peek, Vandana Desai, Richard L. White, Raffaele D'Abrusco, Joseph M. Mazzarella, Carolyn Grant, Jenny L. Novacescu, Elena Scire, Sherry Winkelman

Abstract: We present a bibliographic analysis of Chandra, Hubble, and Spitzer publications. We find (a) archival data are used in >60% of the publication output and (b) archives for these missions enable a much broader set of institutions and countries to scientifically use data from these missions. Specifically, we find that authors from institutions that have published few papers from a given mission publ… ▽ More We present a bibliographic analysis of Chandra, Hubble, and Spitzer publications. We find (a) archival data are used in >60% of the publication output and (b) archives for these missions enable a much broader set of institutions and countries to scientifically use data from these missions. Specifically, we find that authors from institutions that have published few papers from a given mission publish 2/3 archival publications, while those with many publications typically have 1/3 archival publications. We also show that countries with lower GDP per capita overwhelmingly produce archival publications, while countries with higher GDP per capital produce guest observer and archival publications in equal amounts. We argue that robust archives are thus not only critical for the scientific productivity of mission data, but also the scientific accessibility of mission data. We argue that the astronomical community should support archives to maximize the overall scientific societal impact of astronomy, and represent an excellent investment in astronomy's future. △ Less

Submitted 14 July, 2019; originally announced July 2019.

Comments: White Paper submitted to the NAS call for Astro2020 Decadal Survey APC papers

Petabytes to Science

Authors: Amanda E. Bauer, Eric C. Bellm, Adam S. Bolton, Surajit Chaudhuri, A. J. Connolly, Kelle L. Cruz, Vandana Desai, Alex Drlica-Wagner, Frossie Economou, Niall Gaffney, J. Kavelaars, J. Kinney, Ting S. Li, B. Lundgren, R. Margutti, G. Narayan, B. Nord, Dara J. Norman, W. O'Mullane, S. Padhi, J. E. G. Peek, C. Schafer, Megan E. Schwamb, Arfon M. Smith, Erik J. Tollerud , et al. (2 additional authors not shown)

Abstract: A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauth… ▽ More A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey. △ Less

Submitted 17 November, 2019; v1 submitted 13 May, 2019; originally announced May 2019.

Comments: 70 pages 2 figures - this contains a few fixes

Do Androids Dream of Magnetic Fields? Using Neural Networks to Interpret the Turbulent Interstellar Medium

Authors: J. E. G. Peek, Blakesley Burkhart

Abstract: The interstellar medium (ISM) of galaxies is composed of a turbulent magnetized plasma. In order to quantitatively measure relevant turbulent parameters of the ISM, a wide variety of statistical techniques and metrics have been developed that are often tested using numerical simulations and analytic formalism. These metrics are typically based on the Fourier power spectrum, which does not capture… ▽ More The interstellar medium (ISM) of galaxies is composed of a turbulent magnetized plasma. In order to quantitatively measure relevant turbulent parameters of the ISM, a wide variety of statistical techniques and metrics have been developed that are often tested using numerical simulations and analytic formalism. These metrics are typically based on the Fourier power spectrum, which does not capture the Fourier phase information that carries the morphological characteristics of images. In this work we use density slices of magnetohydrodyanmic turbulence simulations to demonstrate that a modern tool, convolutional neural networks, can capture significant information encoded in the Fourier phases. We train the neural network to distinguish between two simulations with different levels of magnetization. We find that, even given a tiny slice of simulation data, a relatively simple network can distinguish sub-Alfvénic (strong magnetic field) and super-Alfvénic (weak magnetic field) turbulence >98% of the time, even when all spectral amplitude information is stripped from the images. In order to better understand how the neural network is picking out differences betweem the two classes of simulations we apply a neural network analysis method called "saliency maps". The saliency map analysis shows that sharp ridge-like features are a distinguishing morphological characteristic in such simulations. Our analysis provides a way forward for deeper understanding of the relationship between magnetohydrodyanmic turbulence and gas morphology and motivates further applications of neural networks for studies of turbulence. We make publicly available all data and software needed to reproduce our results. △ Less

Submitted 21 August, 2019; v1 submitted 2 May, 2019; originally announced May 2019.

Comments: Accepted to ApJL

The Panchromatic Circumgalactic Medium

Authors: Q. Daniel Wang, Joseph N. Burchett, Nicolas Lehner, John M. O'Meara, Molly S. Peeples, J. E. G. Peek, Marc Rafelski, Jason Tumlinson, Jessica Werk, Dennis Zaritsky

Abstract: Galaxies are surrounded by extended atmospheres, which are often called the circumgalactic medium (CGM) and are the least understood part of galactic ecosystems. The CGM serves as a reservoir of both diffuse, metal-poor gas accreted from the intergalactic medium, and metal-rich gas that is either ejected from galaxies by energetic feedback or stripped from infalling satellites. As such, the CGM is… ▽ More Galaxies are surrounded by extended atmospheres, which are often called the circumgalactic medium (CGM) and are the least understood part of galactic ecosystems. The CGM serves as a reservoir of both diffuse, metal-poor gas accreted from the intergalactic medium, and metal-rich gas that is either ejected from galaxies by energetic feedback or stripped from infalling satellites. As such, the CGM is empirically multi-phased and complex in dynamics. Significant progress has been made in the past decade or so in observing the cosmic-ray/B-field, as well as various phases of the CGM. But basic questions remain to be answered. First, what are the energy, mass, and metal contents of the CGM? More specifically, how are they spatially distributed and partitioned in the different components? Moreover, how are they linked to properties of host galaxies and their global clustering and intergalactic medium environments? Lastly, what are the origin, state, and life-cycle of the CGM? This question explores the dynamics of the CGM. Here we illustrate how these questions may be addressed with multi-wavelength observations of the CGM. △ Less

Submitted 25 March, 2019; originally announced March 2019.

Comments: Astro2020 Science White Paper

Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics

Authors: A. Siemiginowska, G. Eadie, I. Czekala, E. Feigelson, E. B. Ford, V. Kashyap, M. Kuhn, T. Loredo, M. Ntampaka, A. Stevens, A. Avelino, K. Borne, T. Budavari, B. Burkhart, J. Cisewski-Kehe, F. Civano, I. Chilingarian, D. A. van Dyk, G. Fabbiano, D. P. Finkbeiner, D. Foreman-Mackey, P. Freeman, A. Fruscione, A. A. Goodman, M. Graham , et al. (27 additional authors not shown)

Abstract: Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New met… ▽ More Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers

Understanding the circumgalactic medium is critical for understanding galaxy evolution

Authors: Molly S. Peeples, Peter Behroozi, Rongmon Bordoloi, Alyson Brooks, James S. Bullock, Joseph N. Burchett, Hsiao-Wen Chen, John Chisholm, Charlotte Christensen, Alison Coil, Lauren Corlies, Aleksandar Diamond-Stanic, Megan Donahue, Claude-André Faucher-Giguère, Henry Ferguson, Drummond Fielding, Andrew J. Fox, David M. French, Steven R. Furlanetto, Mario Gennaro, Karoline M. Gilbert, Erika Hamden, Nimish Hathi, Matthew Hayes, Alaina Henry , et al. (47 additional authors not shown)

Abstract: Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within… ▽ More Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Astro2020 Decadal Science White Paper

Journal ref: Astro2020: Decadal Survey on Astronomy and Astrophysics, science white papers, no. 368; Bulletin of the American Astronomical Society, Vol. 51, Issue 3, id. 368 (2019)

High-Dimensional Dust Mapping

Authors: G. Zasowski, D. P. Finkbeiner, G. M. Green, J. A. Kollmeier, D. M. Nataf, J. E. G. Peek, E. F. Schlafly, V. Silva Aguirre, J. S. Speagle, K. Tchernyshyov, J. D. Trujillo, C. Zucker

Abstract: Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM)… ▽ More Galactic interstellar dust has a profound impact not only on our observations of objects throughout the Universe, but also on the morphology, star formation, and chemical evolution of the Galaxy. The advent of massive imaging and spectroscopic surveys (particularly in the infrared) places us on the threshold of being able to map the properties and dynamics of dust and the interstellar medium (ISM) in three dimensions throughout the Milky Way disk and bulge. These developments will enable a fundamentally new understanding of dust properties, including how grains respond to their local environment and how those environments affect dust attenuation of background objects of interest. Distance-resolved maps of dust motion also hold great promise for tracing the flow of interstellar material throughout the Galaxy on a variety of scales, from bar-streaming motions to the collapse and dissolution of individual molecular clouds. These advances require optical and infrared imaging of stars throughout the Galactic midplane, stretching many kiloparsecs from the Sun, matched with very dense spectroscopic coverage to probe the ISM's fine-grained structure. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 call for science white papers. 5 pages plus references, 1 figure

The Role of Machine Learning in the Next Decade of Cosmology

Authors: Michelle Ntampaka, Camille Avestruz, Steven Boada, Joao Caldeira, Jessi Cisewski-Kehe, Rosanne Di Stefano, Cora Dvorkin, August E. Evrard, Arya Farahi, Doug Finkbeiner, Shy Genel, Alyssa Goodman, Andy Goulding, Shirley Ho, Arthur Kosowsky, Paul La Plante, Francois Lanusse, Michelle Lochner, Rachel Mandelbaum, Daisuke Nagai, Jeffrey A. Newman, Brian Nord, J. E. G. Peek, Austin Peel, Barnabas Poczos , et al. (5 additional authors not shown)

Abstract: In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster an… ▽ More In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors. △ Less

Submitted 14 January, 2021; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: Submitted to the Astro2020 call for science white papers

The Physical Nature of Neutral Hydrogen Intensity Structure

Authors: S. E. Clark, J. E. G. Peek, M. -A. Miville-Deschênes

Abstract: We investigate the physical properties of structures seen in channel map observations of 21-cm neutral hydrogen (HI) emission. HI intensity maps display prominent linear structures that are well aligned with the ambient magnetic field in the diffuse interstellar medium (ISM). Some literature hold that these structures are "velocity caustics", fluctuations imprinted by the turbulent velocity field,… ▽ More We investigate the physical properties of structures seen in channel map observations of 21-cm neutral hydrogen (HI) emission. HI intensity maps display prominent linear structures that are well aligned with the ambient magnetic field in the diffuse interstellar medium (ISM). Some literature hold that these structures are "velocity caustics", fluctuations imprinted by the turbulent velocity field, and are not three-dimensional density structures in the ISM. We test this hypothesis by stacking probes of the density field -- broadband far infrared (FIR) emission and the integrated HI column density (NHI) -- at the locations of linear HI intensity structures. We find that the HI intensity features are real density structures and not velocity caustics. We generalize the investigation to all small-scale structure in HI channel maps, and analyze this correlation as a function of velocity channel width, finding no measurable contribution from velocity caustics to the HI channel map emission. Further, we find that small-scale HI channel maps structures have elevated FIR/NHI, implying that this emission originates from a colder, denser phase of the ISM than the surrounding material. The data are consistent with a multi-phase diffuse ISM in which small-scale structures in narrow HI channel maps are preferentially cold neutral medium (CNM) that is anisotropically distributed and aligned with the local magnetic field. The shallow spatial power spectrum (SPS) of narrow HI channels is often attributed to velocity caustics. We conjecture instead that the small-scale structure and narrow linewidths typical of CNM explain the observed relationship between the SPS and channel width. △ Less

Submitted 9 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 16 pages, 8 figures, minor edits to reflect published version

Journal ref: The Astrophysical Journal, Volume 874, Issue 2, article id. 171, 12 pp. (2019)

LSST Observing Strategy White Paper: LSST Observations of WFIRST Deep Fields

Authors: R. J. Foley, A. M. Koekemoer, D. N. Spergel, F. B. Bianco, P. Capak, L. Dai, O. Dore, G. G. Fazio, H. Ferguson, A. V. Filippenko, B. Frye, L. Galbany, E. Gawiser, C. Gronwall, N. P. Hathi, C. Hirata, R. Hounsell, S. W. Jha, A. G. Kim, P. L. Kelly, J. W. Kruk, S. Malhotra, K. S. Mandel, R. Margutti, D. Marrone , et al. (16 additional authors not shown)

Abstract: The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. I… ▽ More The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: White Paper in response to LSST Call for Observing Strategy Input

Catching the Birth of a Dark Molecular Cloud for the First Time

Authors: Pei Zuo, Di Li, J. E. G. Peek, Qiang Chang, Xia Zhang, Nicholas chapman, Paul F. Goldsmith, Zhi-Yu Zhang

Abstract: The majority of hydrogen in the interstellar medium (ISM) is in atomic form. The transition from atoms to molecules and, in particular, the formation of the H$_2$ molecule, is a key step in cosmic structure formation en route to stars. Quantifying H$_2$ formation in space is difficult, due to the confusion in the emission of atomic hydrogen (HI) and the lack of a H$_2$ signal from the cold ISM. He… ▽ More The majority of hydrogen in the interstellar medium (ISM) is in atomic form. The transition from atoms to molecules and, in particular, the formation of the H$_2$ molecule, is a key step in cosmic structure formation en route to stars. Quantifying H$_2$ formation in space is difficult, due to the confusion in the emission of atomic hydrogen (HI) and the lack of a H$_2$ signal from the cold ISM. Here we present the discovery of a rare, isolated dark cloud currently undergoing H$_2$ formation, as evidenced by a prominent "ring" of HI self-absorption. Through a combined analysis of HI narrow self-absorption, CO emission, dust emission, and extinction, we directly measured, for the first time, the [HI]/[H$_2$] abundance varying from 2% to 0.2%, within one region. These measured HI abundances are orders of magnitude higher than usually assumed initial conditions for protoplanetary disk models. None of the fast cloud formation model could produce such low atomic hydrogen abundance. We derived a cloud formation timescale of 6$\times$10$^6$ years, consistent with the global Galactic star formation rate, and favoring the classical star formation picture over fast star formation models. Our measurements also help constrain the H$_2$ formation rate, under various ISM conditions. △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: 7 pages, 2 figures, accepted by ApJ

Kinetic Tomography. II. A second method for mapping the velocity field of the Milky Way Interstellar Medium and a comparison with spiral structure models

Authors: Kirill Tchernyshyov, J. E. G. Peek, Gail Zasowski

Abstract: In this work, we derive a spatially resolved map of the line-of-sight velocity of the interstellar medium and use it, along with a second map of line-of-sight velocity from Paper I of this series, to determine the nature of gaseous spiral structure in the Milky Way. This map is derived from measurements of the 1.527 $μ$m diffuse interstellar band (DIB) in stellar spectra from the APOGEE survey and… ▽ More In this work, we derive a spatially resolved map of the line-of-sight velocity of the interstellar medium and use it, along with a second map of line-of-sight velocity from Paper I of this series, to determine the nature of gaseous spiral structure in the Milky Way. This map is derived from measurements of the 1.527 $μ$m diffuse interstellar band (DIB) in stellar spectra from the APOGEE survey and covers the nearest 4-5 kpc of the Northern Galactic plane. We cross-check this new DIB-based line-of-sight velocity map with the map derived in Paper I and find that they agree. We then compare these maps with line-of-sight velocity maps derived from simulations of quasi-stationary density wave spiral structure and dynamic, or material, spiral structure in a Milky Way-like galaxy. While none of the maps derived from these simulations is an exact match to the measured velocity field of the Milky Way, the measurements are more consistent with simulations of dynamic spiral structure. In the dynamic spiral structure simulation that best matches the measurements, the Perseus spiral arm is being disrupted. △ Less

Submitted 3 August, 2018; originally announced August 2018.

Comments: Submitted to AAS journals

Multi-Messenger Astrophysics: Harnessing the Data Revolution

Authors: Gabrielle Allen, Warren Anderson, Erik Blaufuss, Joshua S. Bloom, Patrick Brady, Sarah Burke-Spolaor, S. Bradley Cenko, Andrew Connolly, Peter Couvares, Derek Fox, Avishay Gal-Yam, Suvi Gezari, Alyssa Goodman, Darren Grant, Paul Groot, James Guillochon, Chad Hanna, David W. Hogg, Kelly Holley-Bockelmann, D. Andrew Howell, David Kaplan, Erik Katsavounidis, Marek Kowalski, Luis Lehner, Daniel Muthukrishna , et al. (5 additional authors not shown)

Abstract: The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugura… ▽ More The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe. △ Less

Submitted 12 July, 2018; originally announced July 2018.

Comments: workshop summary, 20 pages, by the participants of the Cyberinfrastructure for Multi-Messenger Astrophysics NSF-funded Workshop Held on 23-24 May 2018 at the University of Maryland

Optically thick HI does not dominate dark gas in the local ISM

Authors: Claire E. Murray, J. E. G. Peek, Min-Young Lee, Snezana Stanimirovic

Abstract: The local interstellar medium (ISM) is suffused with "dark" gas, identified by excess infrared and gamma ray emission, yet undetected by standard ISM tracers such as neutral hydrogen (HI) or carbon monoxide emission. Based on observed dust properties from Planck, recent studies have argued that HI mixed with dust is strongly saturated and that dark gas is dominated by optically-thick HI. We test t… ▽ More The local interstellar medium (ISM) is suffused with "dark" gas, identified by excess infrared and gamma ray emission, yet undetected by standard ISM tracers such as neutral hydrogen (HI) or carbon monoxide emission. Based on observed dust properties from Planck, recent studies have argued that HI mixed with dust is strongly saturated and that dark gas is dominated by optically-thick HI. We test this hypothesis by reproducing this model using data from Planck and new 21 cm emission maps from GALFA-HI -- the first large-area 21cm emission survey with comparable angular resolution to Planck. We compare the results with those from a large sample of HI column densities based on direct observations of HI optical depth, and find that the inferred column density corrections are significantly lower than those inferred by the Planck-based model. Further, we rule out the hypothesis that the pencil-beam HI absorption sight lines preferentially miss opaque "blobs" with small covering fraction, as these structures require densities and pressures which are incompatible with ISM conditions. Our results support the picture that excess dust emission in the local ISM is not dominated by optically-thick HI, but is rather a combination of intrinsic changes in dust grain emissivities and H2 missed by CO observations. △ Less

Submitted 14 June, 2018; v1 submitted 4 June, 2018; originally announced June 2018.

Comments: 10 pages, 7 figures; ApJ accepted

A Model for Data Citation in Astronomical Research using Digital Object Identifiers (DOIs)

Authors: Jenny Novacescu, Joshua E. G. Peek, Sarah Weissman, Scott W. Fleming, Karen Levay, Elizabeth Fraser

Abstract: Standardizing and incentivizing the use of digital object identifiers (DOIs) to aggregate and identify both data analyzed and data generated by a research project will advance the field of astronomy to match best practices in other research fields like geosciences and medicine. Increase in the use of DOIs will prepare the discipline for changing expectations among funding agencies and publishers,… ▽ More Standardizing and incentivizing the use of digital object identifiers (DOIs) to aggregate and identify both data analyzed and data generated by a research project will advance the field of astronomy to match best practices in other research fields like geosciences and medicine. Increase in the use of DOIs will prepare the discipline for changing expectations among funding agencies and publishers, who increasingly expect accurate and thorough data citation to accompany scientific outputs. The use of DOIs ensures a robust, sustainable, and interoperable approach to data citation in which due credit is given to researchers and institutions who produce and maintain the primary data. We describe in this work the advantages of DOIs for data citation and best practices for integrating a DOI service in an astronomical archive. We report on a pilot project carried out in collaboration with AAS Journals. During the course of the 1.5 year pilot, over 75% of submitting authors opted to use the integrated DOI service to clearly identify data analyzed during their research project when prompted at the time of paper submission. △ Less

Submitted 28 December, 2017; originally announced January 2018.

Comments: 13 pages, 3 figures. Accepted on Dec 19, 2017 for publication in Astrophysical Journal Supplement Series

Where are All the Gas-Bearing Local Dwarf Galaxies? Quantifying Possible Impacts of Reionization

Authors: Erik J. Tollerud, J. E. G. Peek

Abstract: We present an approach for comparing the detections and non-detections of Local Group (LG) dwarf galaxies in large HI surveys to the predictions of a suite of n-body simulations of the LG. This approach depends primarily on a set of empirical scaling relations to connect the simulations to the observations, rather than making strong theoretical assumptions. We then apply this methodology to the Ga… ▽ More We present an approach for comparing the detections and non-detections of Local Group (LG) dwarf galaxies in large HI surveys to the predictions of a suite of n-body simulations of the LG. This approach depends primarily on a set of empirical scaling relations to connect the simulations to the observations, rather than making strong theoretical assumptions. We then apply this methodology to the Galactic Arecibo L-band Feed Array HI Compact Cloud Catalog (GALFA-HI CCC), and compare it to the ELVIS suite of simulations. This approach reveals a strong tension between the naïve results of the model and the observations: while there are no LG dwarfs in the GALFA-HI CCC, the simulations predict $\sim 10$. Applying a simple model of reionization can resolve this tension by preventing low-mass halos from forming gas. However, and if this effect operates as expected, the observations provide a constraint on the mass scale of dwarf galaxy that reionization impacts. Combined with the observed properties of Leo T, the halo virial mass scale at which reionization impacts dwarf galaxy gas content is constrained to be $\sim 10^{8.5} M_\odot$, independent of any assumptions about star formation. △ Less

Submitted 5 March, 2018; v1 submitted 1 November, 2017; originally announced November 2017.

Comments: 10 pages, 7 figures, accepted to ApJ

Revealing the Milky Way's Hidden Circumgalactic Medium with the Cosmic Origins Spectrograph Quasar Database for Galactic Absorption Lines

Authors: Y. Zheng, J. E. G. Peek, M. E. Putman, J. K. Werk

Abstract: Every quasar (QSO) spectrum contains absorption-line signatures from the interstellar medium, disk-halo interface, and circumgalactic medium (CGM) of the Milky Way (MW). We analyze Hubble Space Telescope/Cosmic Origins Spectrograph (COS) spectra of 132 QSOs to study the significance and origin of SiIV absorption at $|v_{\rm LSR}|\leq100$ km/s in the Galactic halo. The gas in the north predominantl… ▽ More Every quasar (QSO) spectrum contains absorption-line signatures from the interstellar medium, disk-halo interface, and circumgalactic medium (CGM) of the Milky Way (MW). We analyze Hubble Space Telescope/Cosmic Origins Spectrograph (COS) spectra of 132 QSOs to study the significance and origin of SiIV absorption at $|v_{\rm LSR}|\leq100$ km/s in the Galactic halo. The gas in the north predominantly falls in at $-50\lesssim v_{\rm LSR}\lesssim 0$ km/s, whereas in the south, no such pattern is observed. The SiIV column density has an average and a standard deviation of $\langle N_{\rm SiIV}\rangle=(3.8\pm1.4)\times10^{13}$ cm$^{-2}$. At $|b|\gtrsim 30$ degree, $N_{\rm SiIV}$ does not significantly correlate with $b$, which cannot be explained by a commonly adopted flat-slab geometry. We propose a two-component model to reconstruct the $N_{\rm SiIV}$-$b$ distribution: a plane-parallel component $N_{\rm DH}^{\perp}$ to account for the MW's disk-halo interface and a global component $N_{\rm G}$ to reproduce the weak dependence on $b$. We find $N_{\rm DH}^{\perp}=1.3^{+4.7}_{-0.7}\times10^{12}$ cm$^{-2}$ and $N_{\rm G}=(3.4\pm0.3)\times10^{13}$ cm$^{-2}$ on the basis of Bayesian analyses and block bootstrapping. The global component is most likely to have a Galactic origin, although its exact location is uncertain. If it were associated with the MW's CGM, we would find $M_{\rm gas, all}\gtrsim4.7\times10^9\ M_{\odot} (\frac{C_f}{1})(\frac{R}{75\ {\rm kpc}})^2 (\frac{f_{\rm SiIV}}{0.3})^{-1}(\frac{Z}{0.3\ Z_{\odot}})^{-1}$ for the cool gas at all velocities in the Galactic halo. Our analyses show that there is likely a considerable amount of gas at $|v_{\rm LSR}|\leq100$ km s$^{-1}$ hidden in the MW's CGM. Along with this work, we make our QSO dataset publicly available as the COS Quasar Database for Galactic Absorption Lines (COS-GAL). △ Less

Submitted 14 February, 2019; v1 submitted 29 October, 2017; originally announced October 2017.

Comments: Published on ApJ. https://iopscience.iop.org/article/10.3847/1538-4357/aaf6eb/meta

Journal ref: ApJ, 871, 35 (2019)

Calibration of ultraviolet, mid-infrared and radio star formation rate indicators

Authors: Michael J. I. Brown, John Moustakas, Robert C. Kennicutt, Nicolas J. Bonne, Huib T. Intema, Francesco de Gasperin, Mederic Boquien, T. H. Jarrett, Michelle E. Cluver, J. -D. T. Smith, Elisabete da Cunha, Masatoshi Imanishi, Lee Armus, Bernhard R. Brandl, J. E. G. Peek

Abstract: We present calibrations for star formation rate indicators in the ultraviolet, mid-infrared and radio continuum bands, including one of the first direct calibrations of 150 MHz as a star formation rate indicator. Our calibrations utilize 66 nearby star forming galaxies with Balmer decrement corrected H-alpha luminosities, which span 5 orders of magnitude in star formation rate and have absolute ma… ▽ More We present calibrations for star formation rate indicators in the ultraviolet, mid-infrared and radio continuum bands, including one of the first direct calibrations of 150 MHz as a star formation rate indicator. Our calibrations utilize 66 nearby star forming galaxies with Balmer decrement corrected H-alpha luminosities, which span 5 orders of magnitude in star formation rate and have absolute magnitudes of -24<M_r<-12. Most of our photometry and spectrophotometry is measured from the same region of each galaxy, and our spectrophotometry has been validated with SDSS photometry, so our random and systematic errors are small relative to the intrinsic scatter seen in star formation rate indicator calibrations. We find WISE W4 (22.8 micron), Spitzer 24 micron and 1.4 GHz have tight correlations with Balmer decrement corrected H-alpha luminosity, with scatter of only 0.2 dex. Our calibrations are comparable to those from the prior literature for L* galaxies, but for dwarf galaxies our calibrations can give star formation rates that are far greater than those derived from much of the prior literature. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: Accepted for publication in the Astrophysical Journal; 39 pages, 14 figures and 4 tables

The Discovery and Origin of A Very-High Velocity Cloud Toward M33

Authors: Y. Zheng, J. K. Werk, J. E. G. Peek, M. E. Putman

Abstract: We report the detection of a largely ionized very-high velocity cloud (VHVC; $v_{\rm LSR}\sim-350$ km/s) toward M33 with the Hubble Space Telescope/Cosmic Origin Spectrograph. The VHVC is detected in OI, CII, SiII, and SiIII absorption along five sightlines separated by ~0.06-0.4 degree. On sub-degree scales, the velocities and ionic column densities of the VHVC remain relatively smooth with stand… ▽ More We report the detection of a largely ionized very-high velocity cloud (VHVC; $v_{\rm LSR}\sim-350$ km/s) toward M33 with the Hubble Space Telescope/Cosmic Origin Spectrograph. The VHVC is detected in OI, CII, SiII, and SiIII absorption along five sightlines separated by ~0.06-0.4 degree. On sub-degree scales, the velocities and ionic column densities of the VHVC remain relatively smooth with standard deviations of +/-14 km/s and +/-0.15 dex between the sightlines, respectively. The VHVC has a metallicity of [OI/HI]=-0.56+/-0.17 dex (Z=0.28+/-0.11 Z$_{\odot}$). Despite the position-velocity proximity of the VHVC to the ionized Magellanic Stream, the VHVC's higher metallicity makes it unlikely to be associated with the Stream, highlighting the complex velocity structure of this region of sky. We investigate the VHVC's possible origin by revisiting its surrounding HI environment. We find that the VHVC may be: (1) a MW CGM cloud, (2) related to a nearby HI VHVC -- Wright's Cloud, or (3) connected to M33's northern warp. Furthermore, the VHVC could be a bridge connecting Wright's Cloud and M33's northern warp, which would make it a Magellanic-like structure in the halo of M33. △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: Accepted for publication in ApJ

Journal ref: 2017ApJ...840...65Z

Mapping the Extinction Curve in 3D: Structure on Kiloparsec Scales

Authors: E. F. Schlafly, J. E. G. Peek, D. P. Finkbeiner, G. M. Green

Abstract: Near-infrared spectroscopy from APOGEE and wide-field optical photometry from Pan-STARRS1 have recently made possible precise measurements of the shape of the extinction curve for tens of thousands of stars, parameterized by R(V). These measurements revealed structures in R(V) with large angular scales, which are challenging to explain in existing dust paradigms. In this work, we combine three-dim… ▽ More Near-infrared spectroscopy from APOGEE and wide-field optical photometry from Pan-STARRS1 have recently made possible precise measurements of the shape of the extinction curve for tens of thousands of stars, parameterized by R(V). These measurements revealed structures in R(V) with large angular scales, which are challenging to explain in existing dust paradigms. In this work, we combine three-dimensional maps of dust column density with R(V) measurements to constrain the three-dimensional distribution of R(V) in the Milky Way. We find that variations in R(V) are correlated on kiloparsec scales. In particular, most of the dust within one kiloparsec in the outer Galaxy, including many local molecular clouds (Orion, Taurus, Perseus, California, Cepheus), has a significantly lower R(V) than more distant dust in the Milky Way. These results provide new input to models of dust evolution and processing, and complicate application of locally derived extinction curves to more distant regions of the Milky Way and to other galaxies. △ Less

Submitted 23 March, 2017; v1 submitted 8 December, 2016; originally announced December 2016.

Comments: 12 pages, 9 figures, updating paper to published version

HST/COS Observations of Ionized Gas Accretion at the Disk-halo Interface of M33

Authors: Y. Zheng, J. E. G. Peek, J. K. Werk, M. E. Putman

Abstract: We report the detection of accreting ionized gas at the disk-halo interface of the nearby galaxy M33. We analyze HST/COS absorption-line spectra of seven ultraviolet-bright stars evenly distributed across the disk of M33. We find Si IV absorption components consistently redshifted relative to the bulk M33's ISM absorption along all the sightlines. The Si IV detection indicates an enriched, disk-wi… ▽ More We report the detection of accreting ionized gas at the disk-halo interface of the nearby galaxy M33. We analyze HST/COS absorption-line spectra of seven ultraviolet-bright stars evenly distributed across the disk of M33. We find Si IV absorption components consistently redshifted relative to the bulk M33's ISM absorption along all the sightlines. The Si IV detection indicates an enriched, disk-wide, ionized gas inflow toward the disk. This inflow is most likely multi-phase as the redshifted components can also be observed in ions with lower ionization states (e.g., S II, P II, Fe II, Si II). Kinematic modeling of the inflow is consistent with an accreting layer at the disk-halo interface of M33, which has an accretion velocity of 110$^{+15}_{-20}$ km s$^{-1}$ at a distance of 1.5$^{+1.0}_{-1.0}$ kiloparsec above the disk. The modeling indicates a total mass of $\sim3.9\times10^7$ M$_{\odot}$ for the accreting material at the disk-halo interface on the near side of the M33 disk , with an accretion rate of $\sim2.9$ M$_{\odot}$ yr$^{-1}$. The high accretion rate and the level of metal-enrichment suggest the inflow is likely to be the fall back of M33 gas from a galactic fountain and/or the gas pulled loosed during a close interaction between M31 and M33. Our study of M33 is the first to unambiguously reveal the existence of a disk-wide, ionized gas inflow beyond the Milky Way, providing a better understanding of gas accretion in the vicinity of a galaxy disk. △ Less

Submitted 29 November, 2016; originally announced November 2016.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: 2017ApJ...834..179Z

Kinetic Tomography I: A Method for Mapping the Milky Way's Interstellar Medium in Four Dimensions

Authors: Kirill Tchernyshyov, J. E. G. Peek

Abstract: We have developed a method for deriving the distribution of the Milky Way's interstellar medium as a function of longitude, latitude, distance and line-of-sight velocity. This method takes as input maps of reddening as a function of longitude, latitude, and distance and maps of line emission as a function of longitude, latitude, and line-of-sight velocity. We have applied this method to datasets c… ▽ More We have developed a method for deriving the distribution of the Milky Way's interstellar medium as a function of longitude, latitude, distance and line-of-sight velocity. This method takes as input maps of reddening as a function of longitude, latitude, and distance and maps of line emission as a function of longitude, latitude, and line-of-sight velocity. We have applied this method to datasets covering much of the Galactic plane. The output of this method correctly reproduces the line-of-sight velocities of high-mass star forming regions with known distances from Reid et al. (2014) and qualitatively agrees with results from the Milky Way kinematics literature. These maps will be useful for measuring flows of gas around the Milky Way's spiral arms and into and out of giant molecular clouds. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 11 pages, 5 figures, accepted for publication in The Astronomical Journal

A High-Velocity Cloud Impact Forming a Supershell in the Milky Way

Authors: Geumsook Park, Bon-Chul Koo, Ji-hyun Kang, Steven J. Gibson, J. E. G. Peek, Kevin A. Douglas, Eric J. Korpela, Carl E. Heiles

Abstract: Neutral atomic hydrogen (HI) gas in interstellar space is largely organized into filaments, loops, and shells, the most prominent of which are "supershells". These gigantic structures requiring $\gtrsim 3 \times 10^{52}$ erg to form are generally thought to be produced by either the explosion of multiple supernovae (SNe) in OB associations or alternatively by the impact of high-velocity clouds (HV… ▽ More Neutral atomic hydrogen (HI) gas in interstellar space is largely organized into filaments, loops, and shells, the most prominent of which are "supershells". These gigantic structures requiring $\gtrsim 3 \times 10^{52}$ erg to form are generally thought to be produced by either the explosion of multiple supernovae (SNe) in OB associations or alternatively by the impact of high-velocity clouds (HVCs) falling to the Galactic disk. Here we report the detection of a kiloparsec (kpc)-size supershell in the outskirts of the Milky Way with the compact HVC 040+01$-$282 (hereafter CHVC040) at its geometrical center using the "Inner-Galaxy Arecibo L-band Feed Array" HI 21-cm survey data. The morphological and physical properties of both objects suggest that CHVC040, which is either a fragment of a nearby disrupted galaxy or a cloud originated from an intergalactic accreting flow, collided with the disk $\sim 5$ Myrs ago to form the supershell. Our result shows that some compact HVCs can survive their trip through the Galactic halo and inject energy and momentum into the Milky Way disk. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 10 pages, 5 figures, Accepted for publication in the Astrophysical Journal Letters on July 26, 2016

The First Distance Constraint on the Renegade High Velocity Cloud Complex WD

Authors: J. E. G. Peek, Rongmon Bordoloi, Hugues Sana, Julia Roman-Duval, Jason Tumlinson, Yong Zheng

Abstract: We present medium-resolution, near-ultraviolet VLT/FLAMES observations of the star USNO-A0600-15865535. We adapt a standard method of stellar typing to our measurement of the shape of the Balmer epsilon absorption line to demonstrates that USNO-A0600-15865535 is a blue horizontal branch star, residing in the lower stellar halo at a distance of 4.4 kpc from the Sun. We measure the H & K lines of si… ▽ More We present medium-resolution, near-ultraviolet VLT/FLAMES observations of the star USNO-A0600-15865535. We adapt a standard method of stellar typing to our measurement of the shape of the Balmer epsilon absorption line to demonstrates that USNO-A0600-15865535 is a blue horizontal branch star, residing in the lower stellar halo at a distance of 4.4 kpc from the Sun. We measure the H & K lines of singly-ionized calcium and find two isolated velocity components, one originating in the disk, and one associated with high-velocity cloud complex WD. This detection demonstrated that complex WD is closer than ~4.4 kpc and is the first distance constraint on the +100 km/s Galactic complex of clouds. We find that Complex WD is not in corotation with the Galactic disk as has been assumed for decades. We examine a number of scenarios, and find that the most likely is that Complex WD was ejected from the solar neighborhood and is only a few kpc from the Sun. △ Less

Submitted 21 July, 2016; originally announced July 2016.

Comments: 6 pages, 4 figures, submitted to ApJL