Showing 1–50 of 202 results for author: Murray, N

Playing with FIRE: A Galactic Feedback-Halting Experiment Challenges Star Formation Rate Theories

Authors: Shivan Khullar, Christopher D. Matzner, Norman Murray, Michael Y. Grudić, Dávid Guszejnov, Andrew Wetzel, Philip F. Hopkins

Abstract: Stellar feedback influences the star formation rate (SFR) and the interstellar medium of galaxies in ways that are difficult to quantify numerically, because feedback is an essential ingredient of realistic simulations. To overcome this, we conduct a feedback-halting experiment starting with a Milky Way-mass galaxy in the FIRE-2 simulation framework. Terminating feedback, and comparing to a simula… ▽ More Stellar feedback influences the star formation rate (SFR) and the interstellar medium of galaxies in ways that are difficult to quantify numerically, because feedback is an essential ingredient of realistic simulations. To overcome this, we conduct a feedback-halting experiment starting with a Milky Way-mass galaxy in the FIRE-2 simulation framework. Terminating feedback, and comparing to a simulation in which feedback is maintained, we monitor how the runs diverge. We find that without feedback, interstellar turbulent velocities decay. There is a marked increase of dense material, while the SFR increases by over an order of magnitude. Importantly, this SFR boost is a factor of $\sim$15-20 larger than is accounted for by the increased free fall rate caused by higher densities. This implies that feedback moderates the star formation efficiency per free-fall time more directly than simply through the density distribution. To probe changes at the scale of giant molecular clouds (GMCs), we identify GMCs using density and virial parameter thresholds, tracking clouds as the galaxy evolves. Halting feedback stimulates rapid changes, including a proliferation of new bound clouds, a decrease of turbulent support in loosely-bound clouds, an overall increase in cloud densities, and a surge of internal star formation. Computing the cloud-integrated SFR using several theories of turbulence regulation, we show that these theories underpredict the surge in SFR by at least a factor of three. We conclude that galactic star formation is essentially feedback-regulated on scales that include GMCs, and that stellar feedback affects GMCs in multiple ways. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, under ongoing review at ApJ

COMAP Pathfinder -- Season 2 results III. Implications for cosmic molecular gas content at "Cosmic Half-past Eleven"

Authors: D. T. Chung, P. C. Breysse, K. A. Cleary, D. A. Dunne, J. G. S. Lunde, H. Padmanabhan, N. -O. Stutzer, D. Tolgay, J. R. Bond, S. E. Church, H. K. Eriksen, T. Gaier, J. O. Gundersen, S. E. Harper, A. I. Harris, R. Hobbs, H. T. Ihle, J. Kim, J. W. Lamb, C. R. Lawrence, N. Murray, T. J. Pearson, L. Philip, A. C. S. Readhead, T. J. Rennie , et al. (2 additional authors not shown)

Abstract: The Carbon monOxide Mapping Array Project (COMAP) Pathfinder survey continues to demonstrate the feasibility of line-intensity mapping using high-redshift carbon monoxide (CO) line emission traced at cosmological scales. The latest COMAP Pathfinder power spectrum analysis is based on observations through the end of Season 2, covering the first three years of Pathfinder operations. We use our lates… ▽ More The Carbon monOxide Mapping Array Project (COMAP) Pathfinder survey continues to demonstrate the feasibility of line-intensity mapping using high-redshift carbon monoxide (CO) line emission traced at cosmological scales. The latest COMAP Pathfinder power spectrum analysis is based on observations through the end of Season 2, covering the first three years of Pathfinder operations. We use our latest constraints on the CO(1-0) line-intensity power spectrum at $z\sim3$ to update corresponding constraints on the cosmological clustering of CO line emission and thus the cosmic molecular gas content at a key epoch of galaxy assembly. We first mirror the COMAP Early Science interpretation, considering how Season 2 results translate to limits on the shot noise power of CO fluctuations and the bias of CO emission as a tracer of the underlying dark matter distribution. The COMAP Season 2 results place the most stringent limits on the CO tracer bias to date, at $\langle{Tb}\rangle<4.8$ $μ$K. These limits narrow the model space significantly compared to previous CO line-intensity mapping results while maintaining consistency with small-volume interferometric surveys of resolved line candidates. The results also express a weak preference for CO emission models used to guide fiducial forecasts from COMAP Early Science, including our data-driven priors. We also consider directly constraining a model of the halo-CO connection, and show qualitative hints of capturing the total contribution of faint CO emitters through the improved sensitivity of COMAP data. With continued observations and matching improvements in analysis, the COMAP Pathfinder remains on track for a detection of cosmological clustering of CO emission. △ Less

Submitted 14 June, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 9 pages + bibliography and appendices (13 pages total); 9 figures, 1 table; v2 reflects minor changes made for version submitted to A&A, with no changes to top-line results

COMAP Pathfinder -- Season 2 results II. Updated constraints on the CO(1-0) power spectrum

Authors: N. -O. Stutzer, J. G. S. Lunde, P. C. Breysse, D. T. Chung, K. A. Cleary, D. A. Dunne, H. K. Eriksen, H. T. Ihle, H. Padmanabhan, D. Tolgay, I. K. Wehus, J. R. Bond, S. E. Church, T. Gaier, J. O. Gundersen, A. I. Harris, S. E. Harper, R. Hobbs, J. Kim, J. W. Lamb, C. R. Lawrence, N. Murray, T. J. Pearson, L. Philip, A. C. S. Readhead , et al. (2 additional authors not shown)

Abstract: We present updated constraints on the cosmological 3D power spectrum of carbon monoxide CO(1-0) emission in the redshift range $2.4$-$3.4$. The constraints are derived from the two first seasons of Carbon monOxide Mapping Array Project (COMAP) Pathfinder line-intensity mapping observations aiming to trace star-formation during the Epoch of Galaxy Assembly. These results improve on the previous Ear… ▽ More We present updated constraints on the cosmological 3D power spectrum of carbon monoxide CO(1-0) emission in the redshift range $2.4$-$3.4$. The constraints are derived from the two first seasons of Carbon monOxide Mapping Array Project (COMAP) Pathfinder line-intensity mapping observations aiming to trace star-formation during the Epoch of Galaxy Assembly. These results improve on the previous Early Science (ES) results through both increased data volume and improved data processing methodology. On the methodological side, we now perform cross-correlations between groups of detectors (''feed-groups''), as opposed to cross-correlations between single feeds, and this new feed-group pseudo power spectrum (FGPXS) is constructed to be more robust against systematic effects. In terms of data volume, the effective mapping speed is significantly increased due to an improved observational strategy as well as better data selection methodology. The updated spherically- and field-averaged FGPXS, $\tilde{C}(k)$, is consistent with zero, at a probability-to-exceed of around $34\,\%$, with an excess of $2.7\,σ$ in the most sensitive bin. Our power spectrum estimate is about an order of magnitude more sensitive in our six deepest bins across ${0.09\,\mathrm{Mpc}^{-1} < k < 0.73\,\mathrm{Mpc}^{-1}}$, as compared to the feed-feed pseudo power spectrum (FPXS) of COMAP ES. Each of these bins individually constrains the CO power spectrum to ${kP_\mathrm{CO}(k)< 2400-4900\,\mathrm{μK^2 Mpc^{2}}}$ at $95\,\%$ confidence. To monitor potential contamination from residual systematic effects, we analyze a set of 312 difference-map null tests and find that these are consistent with the instrumental noise prediction. In sum, these results provide the strongest direct constraints on the cosmological 3D CO(1-0) power spectrum published to date. △ Less

Submitted 14 June, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 17 pages, 10 figures, v2 reflects changes made for version submitted to Astronomy and Astrophysics, addition of additional figure clarifying methodology, no change to final results

COMAP Pathfinder -- Season 2 results I. Improved data selection and processing

Authors: J. G. S. Lunde, N. -O. Stutzer, P. C. Breysse, D. T. Chung, K. A. Cleary, D. A. Dunne, H. K. Eriksen, S. E. Harper, H. T. Ihle, J. W. Lamb, T. J. Pearson, L. Philip, I. K. Wehus, D. P. Woody, J. R. Bond, S. E. Church, T. Gaier, J. O. Gundersen, A. I. Harris, R. Hobbs, J. Kim, C. R. Lawrence, N. Murray, H. Padmanabhan, A. C. S. Readhead , et al. (2 additional authors not shown)

Abstract: The CO Mapping Array Project (COMAP) Pathfinder is performing line intensity mapping of CO emission to trace the distribution of unresolved galaxies at redshift $z \sim 3$. We present an improved version of the COMAP data processing pipeline and apply this to the first two seasons of observations. This analysis improves on the COMAP Early Science (ES) results in several key aspects. On the observa… ▽ More The CO Mapping Array Project (COMAP) Pathfinder is performing line intensity mapping of CO emission to trace the distribution of unresolved galaxies at redshift $z \sim 3$. We present an improved version of the COMAP data processing pipeline and apply this to the first two seasons of observations. This analysis improves on the COMAP Early Science (ES) results in several key aspects. On the observational side, all second season scans were made in constant-elevation mode, after noting that the previous Lissajous scans were associated with increased systematic errors; those scans accounted for 50% of the total Season 1 data volume. Secondly, all new observations were restricted to an elevation range of 35-65 degrees, to minimize sidelobe ground pickup. On the data processing side, more effective data cleaning in both the time- and map-domain has allowed us to eliminate all data-driven power spectrum-based cuts. This increases the overall data retention and reduces the risk of signal subtraction bias. On the other hand, due to the increased sensitivity, two new pointing-correlated systematic errors have emerged, and we introduce a new map-domain PCA filter to suppress these. Subtracting only 5 out of 256 PCA modes, we find that the standard deviation of the cleaned maps decreases by 67% on large angular scales, and after applying this filter, the maps appear consistent with instrumental noise. Combining all these improvements, we find that each hour of raw Season 2 observations yields on average 3.2 times more cleaned data compared to ES analysis. Combining this with the increase in raw observational hours, the effective amount of data available for high-level analysis is a factor of 8 higher than in ES. The resulting maps have reached an uncertainty of $25$-$50\,μK$ per voxel, providing by far the strongest constraints on cosmological CO line emission published to date. △ Less

Submitted 14 June, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 23 pages, 22 figures, for submission to Astronomy and Astrophysics

Euclid. II. The VIS Instrument

Authors: Euclid Collaboration, M. Cropper, A. Al-Bahlawan, J. Amiaux, S. Awan, R. Azzollini, K. Benson, M. Berthe, J. Boucher, E. Bozzo, C. Brockley-Blatt, G. P. Candini, C. Cara, R. A. Chaudery, R. E. Cole, P. Danto, J. Denniston, A. M. Di Giorgio, B. Dryer, J. Endicott, J. -P. Dubois, M. Farina, E. Galli, L. Genolet, J. P. D. Gow , et al. (403 additional authors not shown)

Abstract: This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift ran… ▽ More This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift range z=0.1-1.5 resulting from weak gravitational lensing, one of the two principal cosmology probes of Euclid. With photometric redshifts, the distribution of dark matter can be mapped in three dimensions, and, from how this has changed with look-back time, the nature of dark energy and theories of gravity can be constrained. The entire VIS focal plane will be transmitted to provide the largest images of the Universe from space to date, reaching m_AB>24.5 with S/N >10 in a single broad I_E~(r+i+z) band over a six year survey. The particularly challenging aspects of the instrument are the control and calibration of observational biases, which lead to stringent performance requirements and calibration regimes. With its combination of spatial resolution, calibration knowledge, depth, and area covering most of the extra-Galactic sky, VIS will also provide a legacy data set for many other fields. This paper discusses the rationale behind the VIS concept and describes the instrument design and development before reporting the pre-launch performance derived from ground calibrations and brief results from the in-orbit commissioning. VIS should reach fainter than m_AB=25 with S/N>10 for galaxies of full-width half-maximum of 0.3" in a 1.3" diameter aperture over the Wide Survey, and m_AB>26.4 for a Deep Survey that will cover more than 50 deg^2. The paper also describes how VIS works with the other Euclid components of survey, telescope, and science data processing to extract the cosmological information. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

The instability mechanism of compact multiplanet systems

Authors: Caleb Lammers, Sam Hadden, Norman Murray

Abstract: To improve our understanding of orbital instabilities in compact planetary systems, we compare suites of $N$-body simulations against numerical integrations of simplified dynamical models. We show that, surprisingly, dynamical models that account for small sets of resonant interactions between the planets can accurately recover $N$-body instability times. This points towards a simple physical pict… ▽ More To improve our understanding of orbital instabilities in compact planetary systems, we compare suites of $N$-body simulations against numerical integrations of simplified dynamical models. We show that, surprisingly, dynamical models that account for small sets of resonant interactions between the planets can accurately recover $N$-body instability times. This points towards a simple physical picture in which a handful of three-body resonances (3BRs), generated by interactions between nearby two-body mean motion resonances (MMRs), overlap and drive chaotic diffusion, leading to instability. Motivated by this, we show that instability times are well-described by a power law relating instability time to planet separations, measured in units of fractional semi-major axis difference divided by the planet-to-star mass ratio to the $1/4$ power, rather than the frequently-adopted $1/3$ power implied by measuring separations in units of mutual Hill radii. For idealized systems, the parameters of this power-law relationship depend only on the ratio of the planets' orbital eccentricities to the orbit-crossing value, and we report an empirical fit to enable quick instability time predictions. This relationship predicts that observed systems comprised of three or more sub-Neptune mass planets must be spaced with period ratios $\mathcal{P} \gtrsim 1.35$ and that tightly-spaced systems ($\mathcal{P} \lesssim 1.5$) must possess very low eccentricities ($e \lesssim 0.05$) to be stable for more than $10^9$ orbits. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 20 pages, 13 figures. Submitted to AAS Journals

Probing $H_0$ and resolving AGN disks with ultrafast photon counters

Authors: Neal Dalal, Marios Galanis, Charles Gammie, Samuel E. Gralla, Norman Murray

Abstract: Intensity interferometry is a technique developed many decades ago, that has recently enjoyed a renaissance thanks in part to advances in photodetector technology. We investigate the potential for long-baseline optical intensity interferometry to observe bright, active galactic nuclei (AGN) associated with rapidly accreting supermassive black holes. We argue that realistic telescope arrays similar… ▽ More Intensity interferometry is a technique developed many decades ago, that has recently enjoyed a renaissance thanks in part to advances in photodetector technology. We investigate the potential for long-baseline optical intensity interferometry to observe bright, active galactic nuclei (AGN) associated with rapidly accreting supermassive black holes. We argue that realistic telescope arrays similar in area to existing Cherenkov arrays, if equipped with modern high-precision single photon detectors, can achieve a sufficiently high signal to noise ratio not only to detect distant AGN, but also to study them in great detail. We explore the science potential of such observations by considering two examples. First, we find that intensity interferometric observations of bright nearby AGN can allow detailed studies of the central accretion disks powering the AGN, allowing reconstruction of many disk properties like the radial profile. Next, we argue that intensity interferometers can spatially resolve the broad-line regions of AGN at cosmological distances, and thereby provide a geometric determination of the angular diameter distances to those AGN when combined with reverberation mapping. Since this measurement can be performed for AGN at distances of hundreds of megaparsecs, this directly measures the Hubble expansion rate $H_0$, with a precision adequate to resolve the recent Hubble tension. Finally, we speculate on future applications that would be enabled by even larger intensity interferometer arrays. △ Less

Submitted 23 March, 2024; originally announced March 2024.

Circumgalactic Environments around Distant Quasars 3C 9 and 4C 05.84

Authors: Sanchit Sabhlok, Shelley A. Wright, Andrey Vayner, Sonata Simonaitis-Boyd, Norman Murray, Lee Armus, Maren Cosens, James Wiley, Mariska Kriek

Abstract: We present results from the ``Quasar hosts Unveiled by high Angular Resolution Techniques" (QUART) survey studying the Circumgalactic Medium (CGM) by observing rest-frame UV emission lines Ly$α$, C IV and He II around two radio-loud quasars, 3C 9 (z=2.02) and 4C 05.84 (z=2.32), using Keck Cosmic Web Imager (KCWI). We detect large-scale Ly$α$ nebulae around both quasars with projected diameters… ▽ More We present results from the ``Quasar hosts Unveiled by high Angular Resolution Techniques" (QUART) survey studying the Circumgalactic Medium (CGM) by observing rest-frame UV emission lines Ly$α$, C IV and He II around two radio-loud quasars, 3C 9 (z=2.02) and 4C 05.84 (z=2.32), using Keck Cosmic Web Imager (KCWI). We detect large-scale Ly$α$ nebulae around both quasars with projected diameters $\sim$ 100 kpc, with spatially resolved, embedded 15-30 kpc He II and C IV nebulae around both quasars as well as kinematically distinct He II and C IV nebulae at a physical separation of $\sim$ 15 kpc from both quasars. Observations of H$α$, H$β$, and [O III] emission using Keck MOSFIRE spectroscopically confirm that the Ly$α$ nebulae extend to companion galaxies and that these quasars are in a protogroup/protocluster environment. We confirm that the He II and C IV emission is kinematically and spatially coincident with the companion galaxies. We estimate the virial masses of the companion galaxies, their metallicities, and star formation rates, and investigate the sources of ionization. We measure the dynamical mass of the host dark matter halos and estimate that the dark matter halos of these systems will grow to a mass of 2 $\times 10^{14}$ M$_{\odot}$ (3C 9) and 2 $\times 10^{13}$ M$_{\odot}$ (4C 05.84) by z=0. The combined CGM and companion galaxies observations indicate Ly$α$ substructure can indicate the presence of companion galaxies in the CGM. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Dissipative Dark Substructure: The Consequences of Atomic Dark Matter on Milky Way Analog Subhalos

Authors: Caleb Gemmell, Sandip Roy, Xuejian Shen, David Curtin, Mariangela Lisanti, Norman Murray, Philip F. Hopkins

Abstract: Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a sub-component of Atomic Dark Matter (ADM). ADM differs from Cold Dark Matter (CDM) due to the presence of self interactions that lead to energy dissipation and bound-state formation, analogous to Standard Model baryons. This model can arise in complex dark sectors that a… ▽ More Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a sub-component of Atomic Dark Matter (ADM). ADM differs from Cold Dark Matter (CDM) due to the presence of self interactions that lead to energy dissipation and bound-state formation, analogous to Standard Model baryons. This model can arise in complex dark sectors that are natural and theoretically-motivated extensions to the Standard Model. The simulations used in this work were carried out using GIZMO and utilize the FIRE-2 galaxy formation physics in the Standard Model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo's central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host visible satellite galaxies have cuspier density profiles and smaller stellar-half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is much more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 18 pages, 10 figures, 8 pages of appendices

An Analytic Model For Magnetically-Dominated Accretion Disks

Authors: Philip F. Hopkins, Jonathan Squire, Eliot Quataert, Norman Murray, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Claude-Andre Faucher-Giguere, Sarah Wellons

Abstract: Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and m… ▽ More Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and magnetically-arrested disks: the midplane pressure is dominated by toroidal magnetic fields with plasma $β\ll 1$ powered by advection of magnetic flux from the interstellar medium (ISM), and they are super-sonically and trans-Alfvenically turbulent with cooling times short compared to dynamical times yet remain gravitationally stable owing to magnetic support. In this paper, we present a simple analytic similarity model for such disks. For reasonable assumptions, the model is entirely specified by the boundary conditions (inflow rate at the BH radius of influence [BHROI]). We show that the scalings from this model are robust to various detailed assumptions, agree remarkably well with the simulations (given their simplicity), and demonstrate the self-consistency and gravitational stability of such disks even in the outer accretion disk (approaching the BHROI) at hyper-Eddington accretion rates. △ Less

Submitted 12 March, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 12 pages, 1 figure. Replaced with accepted version. Companion paper to 'FORGE'D IN FIRE II' (arXiv:2310.04506, part of a series with arXiv:2309.13115 -- animations of the simulations referred to here can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html)

FORGE'd in FIRE II: The Formation of Magnetically-Dominated Quasar Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Jonathan Squire, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Yanlong Shi, Michael Y. Grudic, Sarah Wellons, Claude-Andre Faucher-Giguere, Daniel Angles-Alcazar, Norman Murray, Eliot Quataert

Abstract: In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc s… ▽ More In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc scales. The disks are optically thick, with radiative cooling balancing accretion, but with properties that are distinct from those assumed in most previous accretion disk models. The pressure is strongly dominated by (primarily toroidal) magnetic fields, with a plasma $β\sim 10^{-4}$ even in the disk midplane. They are qualitatively distinct from magnetically elevated or arrested disks. The disks are strongly turbulent, with trans-Alfvenic and highly super-sonic turbulence, and balance this via a cooling time that is short compared to the disk dynamical time, and can sustain highly super-Eddington accretion rates. Their surface and 3D densities at $\sim 10^{3}-10^{5}$ gravitational radii are much lower than in a Shakura-Sunyaev disk, with important implications for their thermo-chemistry and stability. We show how the magnetic field strengths and geometries arise from rapid advection of flux with the inflow from much weaker galaxy-scale fields in these 'flux-frozen' disks, and how this stabilizes the disk and gives rise to efficient torques. Re-simulating without magnetic fields produces catastrophic fragmentation with a vastly smaller, lower-$\dot{M}$ Shakura-Sunyaev-like disk. △ Less

Submitted 18 January, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 46 pages, 30 figures. Updated to match accepted version from The Open Journal of Astrophysics. Part of a series with arXiv:2309.13115 . Animations of the simulations can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

FORGE'd in FIRE: Resolving the End of Star Formation and Structure of AGN Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Michael Y. Grudic, Kung-Yi Su, Sarah Wellons, Daniel Angles-Alcazar, Ulrich P. Steinwandel, David Guszejnov, Norman Murray, Claude-Andre Faucher-Giguere, Eliot Quataert, Dusan Keres

Abstract: It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydro… ▽ More It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydrodynamic (RMHD) simulation which self-consistently combines the FIRE physics (relevant on galactic/ISM scales where SF/FB are ensemble-averaged) and STARFORGE physics (relevant on small scales where we track individual (proto)stellar formation and evolution), together with explicit RMHD (including non-ideal MHD and multi-band M1-RHD) which self-consistently treats both optically-thick and thin regimes. This allows us to span scales from ~100 Mpc down to <100 au (~300 Schwarzschild radii) around a SMBH at a time where it accretes as a bright quasar, in a single simulation. We show that accretion rates up to $\sim 10-100\,{\rm M_{\odot}\,yr^{-1}}$ can be sustained into the accretion disk at $\ll 10^{3}\,R_{\rm schw}$, with gravitational torques between stars and gas dominating on sub-kpc scales until star formation is shut down on sub-pc scales by a combination of optical depth to cooling and strong magnetic fields. There is an intermediate-scale, flux-frozen disk which is gravitoturbulent and stabilized by magnetic pressure sustaining strong turbulence and inflow with persistent spiral modes. In this paper we focus on how gas gets into the small-scale disk, and how star formation is efficiently suppressed. △ Less

Submitted 12 March, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 38 pages, 18 figures. Replaced with accepted version

On the Degree of Dynamical Packing in the Kepler Multi-planet Systems

Authors: Alysa Obertas, Daniel Tamayo, Norm Murray

Abstract: Current planet formation theories rely on initially compact orbital configurations undergoing a (possibly extended) phase of giant impacts following the dispersal of the dissipative protoplanetary disk. The orbital architectures of observed mature exoplanet systems have likely been strongly sculpted by chaotic dynamics, instabilities, and giant impacts. One possible signature of systems continuall… ▽ More Current planet formation theories rely on initially compact orbital configurations undergoing a (possibly extended) phase of giant impacts following the dispersal of the dissipative protoplanetary disk. The orbital architectures of observed mature exoplanet systems have likely been strongly sculpted by chaotic dynamics, instabilities, and giant impacts. One possible signature of systems continually reshaped by instabilities and mergers is their dynamical packing. Early Kepler data showed that many multi-planet systems are maximally packed - placing an additional planet between an observed pair would make the system unstable. However, this result relied on placing the inserted planet in the most optimistic configuration for stability (e.g., circular orbits). While this would be appropriate in an ordered and dissipative picture of planet formation (i.e. planets dampen into their most stable configurations), we argue that this best-case scenario for stability is rarely realized due to the strongly chaotic nature of planet formation. Consequently, the degree of dynamical packing in multi-planet systems under a realistic formation model is likely significantly higher than previously realized. We examine the full Kepler multi planet sample through this new lens, showing that ~60-95% of Kepler multi-planet systems are strongly packed and that dynamical packing increases with multiplicity. This may be a signature of dynamical sculpting or of undetected planets, showing that dynamical packing is an important metric that can be incorporated into planet formation modelling or when searching for unseen planets. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 15 pages, 4 figures. Accepted for publication in MNRAS

Simulating Atomic Dark Matter in Milky Way Analogues

Authors: Sandip Roy, Xuejian Shen, Mariangela Lisanti, David Curtin, Norman Murray, Philip F. Hopkins

Abstract: Dark sector theories naturally lead to multi-component scenarios for dark matter where a sub-component can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogues where the majority of dark matter is collisionless Cold Dark Matter (CDM), but a sub-component (6%) is strongly diss… ▽ More Dark sector theories naturally lead to multi-component scenarios for dark matter where a sub-component can dissipate energy through self-interactions, allowing it to efficiently cool inside galaxies. We present the first cosmological hydrodynamical simulations of Milky Way analogues where the majority of dark matter is collisionless Cold Dark Matter (CDM), but a sub-component (6%) is strongly dissipative minimal Atomic Dark Matter (ADM). The simulations, implemented in GIZMO and utilizing FIRE-2 galaxy formation physics to model the standard baryonic sector, demonstrate that the addition of even a small fraction of dissipative dark matter can significantly impact galactic evolution despite being consistent with current cosmological constraints. We show that ADM gas with roughly Standard-Model-like masses and couplings can cool to form a rotating "dark disk" with angular momentum closely aligned with the visible stellar disk. The morphology of the disk depends sensitively on the parameters of the ADM model, which affect the cooling rates in the dark sector. The majority of the ADM gas gravitationally collapses into dark "clumps" (regions of black hole or mirror star formation), which form a prominent bulge and a rotating thick disk in the central galaxy. These clumps form early and quickly sink to the inner ~kpc of the galaxy, affecting the galaxy's star-formation history and present-day baryonic and CDM distributions. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 9 pages, 4 figures, 12 pages of appendices and supplementary figures

COMAP Early Science: VIII. A Joint Stacking Analysis with eBOSS Quasars

Authors: Delaney A. Dunne, Kieran A. Cleary, Patrick C. Breysse, Dongwoo T. Chung, Havard T. Ihle, J. Richard Bond, Hans Kristian Eriksen, Joshua Ott Gundersen, Laura C. Keating, Junhan Kim, Jonas Gahr Sturtzel Lunde, Norman Murray, Hamsa Padmanabhan, Liju Philip, Nils-Ole Stutzer, Doga Tolgay, Ingunn Katherine Wehus, Sarah E. Church, Todd Gaier, Andrew I. Harris, Richard Hobbs, James W. Lamb, Charles R. Lawrence, Anthony C. S. Readhead, David P. Woody

Abstract: We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 243 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.129 Jy km/s. De… ▽ More We present a new upper limit on the cosmic molecular gas density at $z=2.4-3.4$ obtained using the first year of observations from the CO Mapping Array Project (COMAP). COMAP data cubes are stacked on the 3D positions of 243 quasars selected from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalog, yielding a 95% upper limit for flux from CO(1-0) line emission of 0.129 Jy km/s. Depending on the balance of the emission between the quasar host and its environment, this value can be interpreted as an average CO line luminosity $L'_\mathrm{CO}$ of eBOSS quasars of $\leq 1.26\times10^{11}$ K km pc$^2$ s$^{-1}$, or an average molecular gas density $ρ_\mathrm{H_2}$ in regions of the universe containing a quasar of $\leq 1.52\times10^8$ M$_\odot$ cMpc$^{-3}$. The $L'_\mathrm{CO}$ upper limit falls among CO line luminosities obtained from individually-targeted quasars in the COMAP redshift range, and the $ρ_\mathrm{H_2}$ value is comparable to upper limits obtained from other Line Intensity Mapping (LIM) surveys and their joint analyses. Further, we forecast the values obtainable with the COMAP/eBOSS stack after the full 5-year COMAP Pathfinder survey. We predict that a detection is probable with this method, depending on the CO properties of the quasar sample. Based on the achieved sensitivity, we believe that this technique of stacking LIM data on the positions of traditional galaxy or quasar catalogs is extremely promising, both as a technique for investigating large galaxy catalogs efficiently at high redshift and as a technique for bolstering the sensitivity of LIM experiments, even with a fraction of their total expected survey data. △ Less

Submitted 26 February, 2024; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: 26 pages, 22 figures; Version 2 is as accepted by ApJ

Intra-system uniformity: a natural outcome of dynamical sculpting

Authors: Caleb Lammers, Sam Hadden, Norman Murray

Abstract: There is evidence that exoplanet systems display intra-system uniformity in mass, radius, and orbital spacing (like ``peas in a pod'') when compared with the system-to-system variations of planetary systems. This has been interpreted as the outcome of the early stages of planet formation, indicative of a picture in which planets form at characteristic mass scales with uniform separations. In this… ▽ More There is evidence that exoplanet systems display intra-system uniformity in mass, radius, and orbital spacing (like ``peas in a pod'') when compared with the system-to-system variations of planetary systems. This has been interpreted as the outcome of the early stages of planet formation, indicative of a picture in which planets form at characteristic mass scales with uniform separations. In this paper, we argue instead that intra-system uniformity in planet sizes and orbital spacings likely arose from the dynamical sculpting of initially-overly-packed planetary systems (in other words, the giant impact phase). With a suite of $N$-body simulations, we demonstrate that systems with random initial masses and compact planet spacings naturally develop intra-system uniformity, in quantitative agreement with observations, due to collisions between planets. Our results suggest that the pre-giant impact planet mass distribution is fairly wide and provide evidence for the prevalence of dynamical sculpting in shaping the observed population of exoplanets. △ Less

Submitted 17 July, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: 6 pages, 5 figures. Accepted for publication in MNRAS Letters

$\rm [C_{II}]$ 158 $\rm μm$ emission as an indicator of galaxy star formation rate

Authors: Lichen Liang, Robert Feldmann, Norman Murray, Desika Narayanan, Christopher C. Hayward, Daniel Anglés-Alcázar, Luigi Bassini, Alexander J. Richings, Claude-André Faucher-Giguère, Dongwoo T. Chung, Jennifer Y. H. Chan, Doǧa Tolgay, Onur Çatmabacak, Dušan Kereš, Philip F. Hopkins

Abstract: Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-lumin… ▽ More Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this `$\rm [C_{II}]$ deficit' is unclear. In this work, we study the $L_{\rm [C_{II}]}$-SFR relation of galaxies using a sample of $z=0-8$ galaxies with $M_*\approx10^7-5\times10^{11}\,M_\odot$ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (FIRE) project. We find a simple analytic expression for $L_{\rm [C_{II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C_{II}]$-emitting gas ($f_{\rm [C_{II}]}$), gas metallicity ($Z_{\rm gas}$), gas density ($n_{\rm gas}$) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes, where $t_{\rm dep}$ ($Z_{\rm gas}$) is the main driver of the $\rm [C_{II}]$ deficit in $\rm H_2$-rich ($\rm H_2$-poor) galaxies. The observed $\rm [C_{II}]$ deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that $\rm [C_{II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C_{II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C_{II}]$ line intensity mapping experiments. △ Less

Submitted 6 December, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: 44 page, 24 figures, 8 tables and 8 appendices, accepted for publication in MNRAS. Key figures: fig. 16 & 17. Comments are welcome

Cold Mode Gas Accretion on Two Galaxy Groups at z$\sim$2

Authors: Andrey Vayner, Nadia L. Zakamska, Sanchit Sabhlok, Shelley A. Wright, Lee Armus, Norman Murray, Gregory Walth, Yuzo Ishikawa

Abstract: We present Keck Cosmic Web Imager (KCWI) integral field spectroscopy (IFS) observations of rest-frame UV emission lines $\rm Lyα$, C IV $λλ$ 1548 Å, 1550Åand He II 1640 Åobserved in the circumgalactic medium (CGM) of two $z=2$ radio-loud quasar host galaxies. We detect extended emission on 80-90 kpc scale in $\rm Lyα$ in both systems with C IV, and He II emission also detected out to 30-50 kpc. Al… ▽ More We present Keck Cosmic Web Imager (KCWI) integral field spectroscopy (IFS) observations of rest-frame UV emission lines $\rm Lyα$, C IV $λλ$ 1548 Å, 1550Åand He II 1640 Åobserved in the circumgalactic medium (CGM) of two $z=2$ radio-loud quasar host galaxies. We detect extended emission on 80-90 kpc scale in $\rm Lyα$ in both systems with C IV, and He II emission also detected out to 30-50 kpc. All emission lines show kinematics with a blue and redshifted gradient pattern consistent with velocities seen in massive dark matter halos and similar to kinematic patterns of inflowing gas seen in hydrodynamical simulations. Using the kinematics of both resolved $\rm Lyα$ emission and absorption, we can confirm that both kinematic structures are associated with accretion. Combining the KCWI data with molecular gas observations with Atacama Large Millimeter/submillimeter Array (ALMA) and high spatial resolution of ionized gas with Keck OSIRIS, we find that both quasar host galaxies reside in proto-group environments at $z=2$. We estimate $1-6\times10^{10}$M$_\odot$ of warm-ionized gas within 30-50 kpc from the quasar that is likely accreting onto the galaxy group. We estimate inflow rates of 60-200 M$_\odot$yr$^{-1}$, within an order of magnitude of the outflow rates in these systems. In the 4C 09.17 system, we detect narrow gas streams associated with satellite galaxies, potentially reminiscent of ram-pressure stripping seen in local galaxy groups and clusters. We find that the quasar host galaxies reside in dynamically complex environments, with ongoing mergers, gas accretion, ISM stripping, and outflows likely playing an important role in shaping the assembly and evolution of massive galaxies at cosmic noon. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: 24 pages, 11 figures, 6 tabes. Accepted for publication in MNRAS

Boundless baryons: how diffuse gas contributes to anisotropic tSZ signal around simulated Three Hundred clusters

Authors: Martine Lokken, Weiguang Cui, J. Richard Bond, Renée Hložek, Norman Murray, Romeel Davé, Alexander van Engelen

Abstract: Upcoming advances in galaxy surveys and cosmic microwave background data will enable measurements of the anisotropic distribution of diffuse gas in filaments and superclusters at redshift $z=1$ and beyond, observed through the thermal Sunyaev-Zel'dovich (tSZ) effect. These measurements will help distinguish between different astrophysical feedback models, account for baryons that appear to be 'mis… ▽ More Upcoming advances in galaxy surveys and cosmic microwave background data will enable measurements of the anisotropic distribution of diffuse gas in filaments and superclusters at redshift $z=1$ and beyond, observed through the thermal Sunyaev-Zel'dovich (tSZ) effect. These measurements will help distinguish between different astrophysical feedback models, account for baryons that appear to be 'missing' from the cosmic census, and present opportunities for using locally-anisotropic tSZ statistics as cosmological probes. This study seeks to guide such future measurements by analysing whether diffuse intergalactic gas is a major contributor to anisotropic tSZ signal in The Three Hundred Gizmo-Simba hydrodynamic simulations. We apply multiple different halo boundary and temperature criteria to divide concentrated from diffuse gas at $z=1$, then create mock Compton-$y$ maps for the separated components. The maps from 98 simulation snapshots are centred on massive galaxy clusters, oriented by the most prominent galaxy filament axis, and stacked. Results vary significantly depending on the definition used for diffuse gas, indicating that assumptions should be clearly defined when claiming observations of the warm-hot intergalactic medium. In all cases, the diffuse gas is important, contributing 25-60% of the tSZ signal in the far field ($>4 h^{-1}$ comoving Mpc) from the stacked clusters. The gas 1-2 virial radii from halo centres is especially key. Oriented stacking and environmental selections help to amplify the signal from the warm-hot intergalactic medium, which is aligned but less concentrated along the filament axis than the hot halo gas. △ Less

Submitted 12 June, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Published in MNRAS. 18 pages, 13 figures, 2 tables. Minor revisions for clarity; results unchanged

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 523, Issue 1, pp.1346-1363 (2023)

Kinematics and Feedback in H II regions in the Dwarf Starburst Galaxy IC 10

Authors: Maren Cosens, Shelley A. Wright, Norman Murray, Lee Armus, Karin Sandstrom, Tuan Do, Kirsten Larson, Gregory Martinez, Sanchit Sabhlok, Andrey Vayner, James Wiley

Abstract: We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with th… ▽ More We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI we map the oxygen abundance and with the high spectral resolution mode we identify 46 individual H II regions. These H II regions have an average radius of 4.0 pc, star formation rate $\sim1.3\times10^{-4}$ M$_\odot$ yr$^{-1}$, and velocity dispersion $\sim$16 km s$^{-1}$. None of the H II regions appear to be virialized ($\rm α_{vir}>>1$), and, on average, they show evidence of ongoing expansion. IC 10's H II regions are offset from the star forming region size-luminosity scaling relationships, as well as Larson's Law that relates size and velocity dispersion. We investigate the balance of inward and outward pressure, $\rm P_{in}$ and $\rm P_{out}$, finding $\rm P_{out}>P_{in}$ in 89% of H II regions, indicating feedback driven expansion even in these low mass H II regions. We find warm gas pressure ($\rm P_{gas}$) provides the dominant contribution to the outward pressure ($\rm P_{out}$). This counteracts the inward pressure which is dominated by turbulence in the surrounding gas rather than self-gravity. Five H II regions show evidence of outflows which are most likely supported by either stellar winds (2 regions) or champagne flows (3 regions). These observations provide new insights into the state of the star-forming regions in IC 10 and negative feedback from low mass clusters. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 57 pages, 35 figures, 9 tables, accepted for publication in The Astrophysical Journal

Simulated Bars May Be Shorter But Are Not Slower Than Observed: TNG50 vs. MaNGA

Authors: Neige Frankel, Annalisa Pillepich, Hans-Walter Rix, Vicente Rodriguez-Gomez, Jason Sanders, Jo Bovy, Juna Kollmeier, Norm Murray, Ted Mackereth

Abstract: Galactic bars are prominent dynamical structures within disk galaxies whose size, formation time, strength, and pattern speed influence the dynamical evolution of their hosts galaxies. Yet, their formation and evolution in a cosmological context is not well understood, as cosmological simulation studies have been limited by the classic trade off between simulation volume and resolution. Here we an… ▽ More Galactic bars are prominent dynamical structures within disk galaxies whose size, formation time, strength, and pattern speed influence the dynamical evolution of their hosts galaxies. Yet, their formation and evolution in a cosmological context is not well understood, as cosmological simulation studies have been limited by the classic trade off between simulation volume and resolution. Here we analyze barred disk galaxies in the cosmological magneto-hydrodynamical simulation TNG50 and quantitatively compare the distributions of bar size and pattern speed to those from MaNGA observations at $z=0$. TNG50 galaxies are selected to match the stellar mass and size distributions of observed galaxies, to account for observational selection effects. We find that the high-resolution of TNG50 yields bars with a wide range of pattern speeds (including those with $\geq 40~\mathrm{km}\,\mathrm{s}^{-1}$\,$\mathrm{kpc}^{-1}$) and a mean value of $\sim36~\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$ consistent with observations within $6\,\mathrm{km}\,\mathrm{s}^{-1}$\,$\mathrm{kpc}^{-1}$, in contrast with previous lower-resolution cosmological simulations that produced bars that were too slow. We find, however, that bars in TNG50 are on average $\sim 35\%$ shorter than observed, although this discrepancy may partly reflect remaining inconsistencies in the simulation-data comparison. This leads to higher values of $\mathcal{R} = R_\mathrm{corot}/R_\mathrm{bar}$ in TNG50, but points to simulated bars being `too short' rather than `too slow'. After repeating the analysis on the lower-resolution run of the same simulation (with the same physical model), we qualitatively reproduce the results obtained in previous studies: this implies that, along with physical model variations, numerical resolution effects may explain the previously found `slowness' of simulated bars. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 14 pages, key figures are Fig. 5 and 6. Submitted to ApJ, comments are welcome

The Mass of the Milky Way from the H3 Survey

Authors: Jeff Shen, Gwendolyn M. Eadie, Norman Murray, Dennis Zaritsky, Joshua S. Speagle, Yuan-Sen Ting, Charlie Conroy, Phillip A. Cargile, Benjamin D. Johnson, Rohan P. Naidu, Jiwon Jesse Han

Abstract: The mass of the Milky Way is a critical quantity which, despite decades of research, remains uncertain within a factor of two. Until recently, most studies have used dynamical tracers in the inner regions of the halo, relying on extrapolations to estimate the mass of the Milky Way. In this paper, we extend the hierarchical Bayesian model applied in Eadie & Jurić (2019) to study the mass distributi… ▽ More The mass of the Milky Way is a critical quantity which, despite decades of research, remains uncertain within a factor of two. Until recently, most studies have used dynamical tracers in the inner regions of the halo, relying on extrapolations to estimate the mass of the Milky Way. In this paper, we extend the hierarchical Bayesian model applied in Eadie & Jurić (2019) to study the mass distribution of the Milky Way halo; the new model allows for the use of all available 6D phase-space measurements. We use kinematic data of halo stars out to $142~{\rm kpc}$, obtained from the H3 Survey and $\textit{Gaia}$ EDR3, to infer the mass of the Galaxy. Inference is carried out with the No-U-Turn sampler, a fast and scalable extension of Hamiltonian Monte Carlo. We report a median mass enclosed within $100~{\rm kpc}$ of $\rm M(<100 \; kpc) = 0.69_{-0.04}^{+0.05} \times 10^{12} \; M_\odot$ (68% Bayesian credible interval), or a virial mass of $\rm M_{200} = M(<216.2_{-7.5}^{+7.5} \; kpc) = 1.08_{-0.11}^{+0.12} \times 10^{12} \; M_\odot$, in good agreement with other recent estimates. We analyze our results using posterior predictive checks and find limitations in the model's ability to describe the data. In particular, we find sensitivity with respect to substructure in the halo, which limits the precision of our mass estimates to $\sim 15\%$. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: 30 pages, 11 figures, 6 tables. Accepted for publication in ApJ

Faint objects in motion: the new frontier of high precision astrometry

Authors: Fabien Malbet, Céline Boehm, Alberto Krone-Martins, Antonio Amorim, Guillem Anglada-Escudé, Alexis Brandeker, Frédéric Courbin, Torsten Enßlin, Antonio Falcão, Katherine Freese, Berry Holl, Lucas Labadie, Alain Léger, Gary Mamon, Barbara Mcarthur, Alcione Mora, Mike Shao, Alessandro Sozzetti, Douglas Spolyar, Eva Villaver, Ummi Abbas, Conrado Albertus, João Alves, Rory Barnes, Aldo Stefano Bonomo , et al. (61 additional authors not shown)

Abstract: Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the front… ▽ More Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: arXiv admin note: substantial text overlap with arXiv:1910.08028, arXiv:1707.01348

Journal ref: Experimental Astronomy, Springer Link, 2021, 51 (3), pp.845-886

COMAP Early Science: VII. Prospects for CO Intensity Mapping at Reionization

Authors: Patrick C. Breysse, Dongwoo T. Chung, Kieran A. Cleary, Håvard T. Ihle, Hamsa Padmanabhan, Marta B. Silva, J. Richard Bond, Jowita Borowska, Morgan Catha, Sarah E. Church, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Joshua Ott Gundersen, Andrew I. Harris, Richard Hobbs, Laura Keating, James W. Lamb, Charles R. Lawrence, Jonas G. S. Lunde, Norman Murray, Timothy J. Pearson, Liju Philip, Maren Rasmussen , et al. (7 additional authors not shown)

Abstract: We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts (… ▽ More We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail. △ Less

Submitted 12 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 7 of 7 in series. 19 pages, 10 figures, to be submitted to ApJ

COMAP Early Science: V. Constraints and Forecasts at $z \sim 3$

Authors: Dongwoo T. Chung, Patrick C. Breysse, Kieran A. Cleary, Håvard T. Ihle, Hamsa Padmanabhan, Marta B. Silva, J. Richard Bond, Jowita Borowska, Morgan Catha, Sarah E. Church, Delaney A. Dunne, Hans Kristian Eriksen, Marie Kristine Foss, Todd Gaier, Joshua Ott Gundersen, Stuart E. Harper, Andrew I. Harris, Brandon Hensley, Richard Hobbs, Laura C. Keating, Junhan Kim, James W. Lamb, Charles R. Lawrence, Jonas Gahr Sturtzel Lunde, Norman Murray , et al. (12 additional authors not shown)

Abstract: We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pat… ▽ More We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pathfinder early science data spanning wavenumbers $k=0.051$-$0.62\,$Mpc$^{-1}$ represent the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude $A_{\rm clust}\lesssim70\,μ$K$^2$ greatly improves on the indirect upper limit of $420\,μ$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$Mpc$^{-1}$. The COMAP limit excludes a subset of models from previous literature, and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity-bias product, $\langle{Tb}\rangle^2\lesssim50\,μ$K$^2$, and the cosmic molecular gas density, $ρ_\text{H2}<2.5\times10^8\,M_\odot\,$Mpc$^{-3}$ (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the five-year Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise of 9-17. Between then and now, we also expect to detect the CO-galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star-formation and galaxy-evolution history. △ Less

Submitted 4 March, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 5 of 7 in series. 17 pages + appendix and bibliography (30 pages total); 15 figures, 6 tables; accepted for publication in ApJ; v3 reflects the accepted version with minor changes and additions to text

Journal ref: ApJ, 933, 186 (2022)

COMAP Early Science: IV. Power Spectrum Methodology and Results

Authors: Håvard T. Ihle, Jowita Borowska, Kieran A. Cleary, Hans Kristian Eriksen, Marie K. Foss, Stuart E. Harper, Junhan Kim, Jonas G. S. Lunde, Liju Philip, Maren Rasmussen, Nils-Ole Stutzer, Bade D. Uzgil, Duncan J. Watts, Ingunn Kathrine Wehus, J. Richard Bond, Patrick C. Breysse, Morgan Catha, Sarah E. Church, Dongwoo T. Chung, Clive Dickinson, Delaney A. Dunne, Todd Gaier, Joshua Ott Gundersen, Andrew I. Harris, Richard Hobbs , et al. (8 additional authors not shown)

Abstract: We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrum… ▽ More We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales $k=0.051-0.62 \,\mathrm{Mpc}^{-1}$ we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4μ\textrm{K}^2\mathrm{Mpc}^3$, the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum in the literature. △ Less

Submitted 6 April, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 4 of 7 in series. 18 pages, 11 figures, as accepted in ApJ

COMAP Early Science: III. CO Data Processing

Authors: Marie K. Foss, Håvard T. Ihle, Jowita Borowska, Kieran A. Cleary, Hans Kristian Eriksen, Stuart E. Harper, Junhan Kim, James W. Lamb, Jonas G. S. Lunde, Liju Philip, Maren Rasmussen, Nils-Ole Stutzer, Bade D. Uzgil, Duncan J. Watts, Ingunn K. Wehus, David P. Woody, J. Richard Bond, Patrick C. Breysse, Morgan Catha, Sarah E. Church, Dongwoo T. Chung, Clive Dickinson, Delaney A. Dunne, Todd Gaier, Joshua Ott Gundersen , et al. (8 additional authors not shown)

Abstract: We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency… ▽ More We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including $χ^2$ and multi-scale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a dataset with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers. △ Less

Submitted 30 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 3 of 7 in series. 26 pages, 23 figures, submitted to ApJ

COMAP Early Science: I. Overview

Authors: Kieran A. Cleary, Jowita Borowska, Patrick C. Breysse, Morgan Catha, Dongwoo T. Chung, Sarah E. Church, Clive Dickinson, Hans Kristian Eriksen, Marie Kristine Foss, Joshua Ott Gundersen, Stuart E. Harper, Andrew I. Harris, Richard Hobbs, Håvard, T. Ihle, Junhan Kim, Jonathon Kocz, James W. Lamb, Jonas G. S. Lunde, Hamsa Padmanabhan, Timothy J. Pearson, Liju Philip, Travis W. Powell, Maren Rasmussen, Anthony C. S. Readhead , et al. (18 additional authors not shown)

Abstract: The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainte… ▽ More The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainter contribution from CO(2-1) at $z=6$-8, the Pathfinder is surveying $12$ deg$^2$ in a 5-year observing campaign to detect the CO signal from $z\sim3$. Using data from the first 13 months of observing, we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4μ\mathrm{K}^2 \mathrm{Mpc}^3$ on scales $k=0.051-0.62 \mathrm{Mpc}^{-1}$ - the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature-bias product) of $\langle Tb\rangle^2<49$ $μ$K$^2$ - nearly an order-of-magnitude improvement on the previous best measurement. These constraints allow us to rule out two models from the literature. We forecast a detection of the power spectrum after 5 years with signal-to-noise ratio (S/N) 9-17. Cross-correlation with an overlapping galaxy survey will yield a detection of the CO-galaxy power spectrum with S/N of 19. We are also conducting a 30 GHz survey of the Galactic plane and present a preliminary map. Looking to the future of COMAP, we examine the prospects for future phases of the experiment to detect and characterize the CO signal from the EoR. △ Less

Submitted 29 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: Paper 1 of 7 in series. 18 pages, 16 figures, submitted to ApJ

Multi-phase outflows in high redshift quasar host galaxies

Authors: Andrey Vayner, Nadia Zakamska, Shelley A. Wright, Lee Armus, Norman Murray, Gregory Walth

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of six radio-loud quasar host galaxies at $z=1.4-2.3$. We combine the kpc-scale resolution ALMA observations with high spatial-resolution adaptive optics integral field spectrograph data of the ionized gas. We detect molecular gas emission in five quasar host galaxies and resolve the molecular interstellar medium using the… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of six radio-loud quasar host galaxies at $z=1.4-2.3$. We combine the kpc-scale resolution ALMA observations with high spatial-resolution adaptive optics integral field spectrograph data of the ionized gas. We detect molecular gas emission in five quasar host galaxies and resolve the molecular interstellar medium using the CO (3-2) or CO (4-3) rotational transitions. Clumpy molecular outflows are detected in four quasar host galaxies and in a merger system 21 kpc away from one quasar. Between the ionized and cold-molecular gas phases, the majority of the outflowing mass is in a molecular phase, while for three out of four detected multi-phase gas outflows, the majority of the kinetic luminosity and momentum flux is in the ionized phase. Combining the energetics of the multi-phase outflows, we find that their driving mechanism is consistent with energy-conserving shocks produced by the impact of the quasar jets with the gas in the galaxy. By assessing the molecular gas mass to the dynamics of the outflows, we estimate a molecular gas depletion time scale of a few Myr. The gas outflow rates exceed the star formation rates, suggesting that quasar feedback is a major mechanism of gas depletion at the present time. The coupling efficiency between the kinetic luminosity of the outflows and the bolometric luminosity of the quasar of 0.1-1% is consistent with theoretical predictions. Studying multi-phase gas outflows at high redshift is important for quantifying the impact of negative feedback in shaping the evolution of massive galaxies. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 23 pages, 8 figures, 6 tabes. Accepted for publication in ApJ

CCAT-prime Collaboration: Science Goals and Forecasts with Prime-Cam on the Fred Young Submillimeter Telescope

Authors: CCAT-Prime collaboration, M. Aravena, J. E. Austermann, K. Basu, N. Battaglia, B. Beringue, F. Bertoldi, F. Bigiel, J. R. Bond, P. C. Breysse, C. Broughton, R. Bustos, S. C. Chapman, M. Charmetant, S. K. Choi, D. T. Chung, S. E. Clark, N. F. Cothard, A. T. Crites, A. Dev, K. Douglas, C. J. Duell, R. Dunner, H. Ebina, J. Erler , et al. (62 additional authors not shown)

Abstract: We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Corn… ▽ More We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Cornell University and sited at more than 5600 meters on Cerro Chajnantor in northern Chile. Prime-Cam is one of two instruments planned for FYST and will provide unprecedented spectroscopic and broadband measurement capabilities to address important astrophysical questions ranging from Big Bang cosmology through reionization and the formation of the first galaxies to star formation within our own Milky Way galaxy. Prime-Cam on the FYST will have a mapping speed that is over ten times greater than existing and near-term facilities for high-redshift science and broadband polarimetric imaging at frequencies above 300 GHz. We describe details of the science program enabled by this system and our preliminary survey strategies. △ Less

Submitted 8 August, 2022; v1 submitted 21 July, 2021; originally announced July 2021.

Comments: 61 pages, 16 figures. Resubmitted to ApJSS July 11, 2022

A Criterion for the Onset of Chaos in Compact, Eccentric Multiplanet Systems

Authors: Daniel Tamayo, Norman Murray, Scott Tremaine, Joshua Winn

Abstract: We derive a semi-analytic criterion for the presence of chaos in compact, eccentric multiplanet systems. Beyond a minimum semimajor-axis separation, below which the dynamics are chaotic at all eccentricities, we show that (i) the onset of chaos is determined by the overlap of two-body mean motion resonances (MMRs), like it is in two-planet systems; (ii) secular evolution causes the MMR widths to e… ▽ More We derive a semi-analytic criterion for the presence of chaos in compact, eccentric multiplanet systems. Beyond a minimum semimajor-axis separation, below which the dynamics are chaotic at all eccentricities, we show that (i) the onset of chaos is determined by the overlap of two-body mean motion resonances (MMRs), like it is in two-planet systems; (ii) secular evolution causes the MMR widths to expand and contract adiabatically, so that the chaotic boundary is established where MMRs overlap at their greatest width. For closely spaced two-planet systems, a near-symmetry strongly suppresses this secular modulation, explaining why the chaotic boundaries for two-planet systems are qualitatively different from cases with more than two planets. We use these results to derive an improved angular-momentum-deficit (AMD) stability criterion, i.e., the critical system AMD below which stability should be guaranteed. This introduces an additional factor to the expression from Laskar and Petit (2017) that is exponential in the interplanetary separations, which corrects the AMD threshold toward lower eccentricities by a factor of several for tightly packed configurations. We make routines for evaluating the chaotic boundary available to the community through the open-source SPOCK package. △ Less

Submitted 23 August, 2021; v1 submitted 28 June, 2021; originally announced June 2021.

Comments: Accepted in AJ. Routines for estimating whether compact planetary configurations are chaotic, together with example notebooks are available as part of the SPOCK package: https://github.com/dtamayo/spock

A Spatially-Resolved Survey of Distant Quasar Host Galaxies: I. Dynamics of galactic outflows

Authors: Andrey Vayner, Shelley A. Wright, Norman Murray, Lee Armus, Anna Boehle, Maren Cosens, James E. Larkin, Etsuko Mieda, Gregory Walth

Abstract: We present observations of ionized gas outflows in eleven z$ =1.39-2.59$ radio-loud quasar host galaxies. Data was taken with the integral field spectrograph (IFS) OSIRIS and the adaptive optics system at the W.M. Keck Observatory targeting nebular emission lines (H$β$, [OIII], H$α$, [NII] and [SII]) redshifted into the near-infrared (1-2.4 \micron). Outflows with velocities of 500 - 1700 km\,s… ▽ More We present observations of ionized gas outflows in eleven z$ =1.39-2.59$ radio-loud quasar host galaxies. Data was taken with the integral field spectrograph (IFS) OSIRIS and the adaptive optics system at the W.M. Keck Observatory targeting nebular emission lines (H$β$, [OIII], H$α$, [NII] and [SII]) redshifted into the near-infrared (1-2.4 \micron). Outflows with velocities of 500 - 1700 km\,s$^{-1}$ are detected in 10 systems on scales ranging from $<1$ kpc to 10 kpc with outflow rates from 8-2400 M$_\odot$yr$^{-1}$. For five sources, the outflow momentum rates are 4-80 times $L_{AGN}$/c, consistent with outflows being driven by an energy conserving shock. The five other outflows are either driven by radiation pressure or an isothermal shock. The outflows are the dominant source of gas depletion, and we find no evidence for star formation along the outflow paths. For eight objects, the outflow paths are consistent with the orientation of the jets. Yet, given the calculated pressures, we find no evidence of the jets currently doing work on these galactic-scale ionized outflows. We find that galactic-scale feedback occurs well before galaxies establish a substantial fraction of their stellar mass, as expected from local scaling relationships. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: 39 pages, 18 figures, 7 tables, accepted for publication in ApJ

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Authors: Cameron Trapp, Dusan Keres, T. K. Chan, Ivanna Escala, Cameron Hummels, Philip F. Hopkins, Claude-Andre Faucher-Giguere, Norman Murray, Eliot Quataert, Andrew Wetzel

Abstract: Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investig… ▽ More Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies, we analyze the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (M_halo ~ 10E12 solar masses), focusing on simulations with cosmic ray physics. We find that at z~0, gas approaches the disk with angular momentum similar to the gaseous disk edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominantly parallel to the disk with small but nonzero vertical velocity components, and joins the disk largely in the outskirts as opposed to "raining" down onto the disk. Once in the disk, gas trajectories are complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show clear but slow net radial infall with transport speeds of 1-3 km/s and net mass fluxes through the disk on the order of one solar mass per year, comparable to the star formation rates of the galaxies and decreasing towards galactic center as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1-7 km/s, ~4-5 solar masses per year). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies. △ Less

Submitted 17 August, 2022; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: 22 pages, 15 figures, published by MNRAS

Journal ref: 2022MNRAS.509.4149T

Mapping the Universe in HD

Authors: Patrick C. Breysse, Simon Foreman, Laura C. Keating, Joel Meyers, Norman Murray

Abstract: Hydrogen deuteride (HD) is prevalent in a wide variety of astrophysical environments, and measuring its large-scale distribution at different epochs can in principle provide information about the properties of these environments. In this paper, we explore the prospects for accessing this distribution using line intensity mapping of emission from the lowest rotational transition in HD, focusing on… ▽ More Hydrogen deuteride (HD) is prevalent in a wide variety of astrophysical environments, and measuring its large-scale distribution at different epochs can in principle provide information about the properties of these environments. In this paper, we explore the prospects for accessing this distribution using line intensity mapping of emission from the lowest rotational transition in HD, focusing on observations of the epoch of reionization ($z\sim6-10$) and earlier. We find the signal from the epoch of reionization to be strongest most promising, through cross-correlations within existing [CII] intensity mapping surveys. While the signal we predict is out of reach for current-generation projects, planned future improvements should be able to detect reionization-era HD without any additional observations, and would help to constrain the properties of the star-forming galaxies thought to play a key role in reionization. We also investigate several avenues for measuring HD during "cosmic dawn" ($z\sim10-30$), a period in which HD could provide one of the only complementary observables to 21$\,$cm intensity maps. We conclude that existing and planned facilities are poorly matched to the specifications desirable for a significant detection, though such a measurement may be achievable with sustained future effort. Finally, we explain why HD intensity mapping of the intergalactic medium during the cosmic dark ages ($z\gtrsim 30$) appears to be out of reach of any conceivable experiment. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: 19 pages, 9 figures

Characterizing mass, momentum, energy and metal outflow rates of multi-phase galactic winds in the FIRE-2 cosmological simulations

Authors: Viraj Pandya, Drummond B. Fielding, Daniel Anglés-Alcázar, Rachel S. Somerville, Greg L. Bryan, Christopher C. Hayward, Jonathan Stern, Chang-Goo Kim, Eliot Quataert, John C. Forbes, Claude-André Faucher-Giguère, Robert Feldmann, Zachary Hafen, Philip F. Hopkins, Dušan Kereš, Norman Murray, Andrew Wetzel

Abstract: We characterize mass, momentum, energy and metal outflow rates of multi-phase galactic winds in a suite of FIRE-2 cosmological "zoom-in" simulations from the Feedback in Realistic Environments (FIRE) project. We analyze simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass halos, and high-redshift massive halos. Consistent with previous work, we find that dwarfs eject about 100… ▽ More We characterize mass, momentum, energy and metal outflow rates of multi-phase galactic winds in a suite of FIRE-2 cosmological "zoom-in" simulations from the Feedback in Realistic Environments (FIRE) project. We analyze simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass halos, and high-redshift massive halos. Consistent with previous work, we find that dwarfs eject about 100 times more gas from their interstellar medium (ISM) than they form in stars, while this mass "loading factor" drops below one in massive galaxies. Most of the mass is carried by the hot phase ($>10^5$ K) in massive halos and the warm phase ($10^3-10^5$ K) in dwarfs; cold outflows ($<10^3$ K) are negligible except in high-redshift dwarfs. Energy, momentum and metal loading factors from the ISM are of order unity in dwarfs and significantly lower in more massive halos. Hot outflows have $2-5\times$ higher specific energy than needed to escape from the gravitational potential of dwarf halos; indeed, in dwarfs, the mass, momentum, and metal outflow rates increase with radius whereas energy is roughly conserved, indicating swept up halo gas. Burst-averaged mass loading factors tend to be larger during more powerful star formation episodes and when the inner halo is not virialized, but we see effectively no trend with the dense ISM gas fraction. We discuss how our results can guide future controlled numerical experiments that aim to elucidate the key parameters governing galactic winds and the resulting associated preventative feedback. △ Less

Submitted 17 September, 2021; v1 submitted 11 March, 2021; originally announced March 2021.

Comments: Accepted to MNRAS with minor revisions, main body is 25 pages with 14 figures

A Spatially-Resolved Survey of Distant Quasar Host Galaxies: II. Photoionization and Kinematics of the ISM

Authors: Andrey Vayner, Shelley A. Wright, Norman Murray, Lee Armus, Anna Boehle, Maren Cosens, James E. Larkin, Etsuko Mieda, Gregory Walth

Abstract: We present detailed observations of photoionization conditions and galaxy kinematics in eleven z$=1.39-2.59$ radio-loud quasar host galaxies. Data was taken with OSIRIS integral field spectrograph (IFS) and the adaptive optics system at the W.M. Keck Observatory that targeted nebular emission lines (H$β$,[OIII],H$α$,[NII]) redshifted into the near-infrared (1-2.4 \micron). We detect extended ioniz… ▽ More We present detailed observations of photoionization conditions and galaxy kinematics in eleven z$=1.39-2.59$ radio-loud quasar host galaxies. Data was taken with OSIRIS integral field spectrograph (IFS) and the adaptive optics system at the W.M. Keck Observatory that targeted nebular emission lines (H$β$,[OIII],H$α$,[NII]) redshifted into the near-infrared (1-2.4 \micron). We detect extended ionized emission on scales ranging from 1-30 kpc photoionized by stars, shocks, and active galactic nuclei (AGN). Spatially resolved emission-line ratios indicate that our systems reside off the star formation and AGN-mixing sequence on the Baldwin, Phillips $\&$ Terlevich (BPT) diagram at low redshift. The dominant cause of the difference between line ratios of low redshift galaxies and our sample is due to lower gas-phase metallicities, which are 2-5$\times$ less compared to galaxies with AGN in the nearby Universe. Using gas velocity dispersion as a proxy to stellar velocity dispersion and dynamical mass measurement through inclined disk modeling we find that the quasar host galaxies are under-massive relative to their central supermassive black hole (SMBH) mass, with all systems residing off the local scaling ($M_{\bullet}-σ~$,$M_{\bullet}-M_{*}~$) relationship. These quasar host galaxies require substantial growth, up to an order of magnitude in stellar mass, to grow into present-day massive elliptical galaxies. Combining these results with part I of our sample paper (Vayner et al. 2021) we find evidence for winds capable of causing feedback before the AGN host galaxies land on the local scaling relation between black hole and galaxy stellar mass, and before the enrichment of the ISM to a level observed in local galaxies with AGN. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 31 pages, 18 figures, 7 tables. Accepted for publication in ApJ

Probing the CGM of Low-redshift Dwarf Galaxies Using FIRE Simulations

Authors: Fei Li, Mubdi Rahman, Norman Murray, Zachary Hafen, Claude-André Faucher-Giguère, Jonathan Stern, Cameron B. Hummels, Philip F. Hopkins, Kareem El-Badry, Dušan Kereš

Abstract: Observations of UV metal absorption lines have provided insight into the structure and composition of the circumgalactic medium (CGM) around galaxies. We compare these observations with the low-redshift ($z \leq 0.3$) CGM around dwarf galaxies in high-resolution cosmological zoom-in runs in the FIRE-2 simulation suite. We select simulated galaxies that match the halo mass, stellar mass, and redshi… ▽ More Observations of UV metal absorption lines have provided insight into the structure and composition of the circumgalactic medium (CGM) around galaxies. We compare these observations with the low-redshift ($z \leq 0.3$) CGM around dwarf galaxies in high-resolution cosmological zoom-in runs in the FIRE-2 simulation suite. We select simulated galaxies that match the halo mass, stellar mass, and redshift of the observed samples. We produce absorption measurements using Trident for UV transitions of C IV, O VI, Mg II and Si III. The FIRE equivalent width (EW) distributions and covering fractions for the C IV ion are broadly consistent with observations inside $0.5 R_{vir}$, but are under-predicted for O VI, Mg II, and Si III. The absorption strengths of the ions in the CGM are moderately correlated with the masses and star formation activity of the galaxies. The correlation strengths increase with the ionization potential of the ions. The structure and composition of the gas from the simulations exhibit three zones around dwarf galaxies characterized by distinct ion column densities: the disky ISM, the inner CGM (the wind-dominated regime), and the outer CGM (the IGM accretion-dominated regime). We find that the outer CGM in the simulations is nearly but not quite supported by thermal pressure, so it is not in hydrostatic equilibrium (HSE), resulting in halo-scale bulk inflow and outflow motions. The net gas inflow rates are comparable to the SFR of the galaxy, but the bulk inflow and outflow rates are greater by an order of magnitude, with velocities comparable to the virial velocity of the halo. These roughly virial velocities (${\sim} 100 km s^{-1}$) produce large EWs in the simulations. This supports a picture for dwarf galaxies in which the dynamics of the CGM at large scales are coupled to the small-scale star formation activity near the centre of their halos. △ Less

Submitted 3 November, 2020; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: To appear in MNRAS

Virialization of the inner CGM in the FIRE simulations and implications for galaxy discs, star formation and feedback

Authors: Jonathan Stern, Claude-André Faucher-Giguère, Drummond Fielding, Eliot Quataert, Zachary Hafen, Alexander B. Gurvich, Xiangcheng Ma, Lindsey Byrne, Kareem El-Badry, Daniel Anglés-Alcázar, T. K. Chan, Robert Feldmann, Dušan Kereš, Andrew Wetzel, Norman Murray, Philip F. Hopkins

Abstract: We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the v… ▽ More We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0<z<5, and how the formation of this virialized phase affects the evolution of galactic discs. We demonstrate that when the halo mass crosses ~10^12 M_sun, the cooling time of shocked gas in the inner CGM (~0.1 R_vir, where R_vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average sub-virial temperatures (T<<T_vir), large pressure fluctuations and supersonic inflow/outflow velocities, to virial temperatures (T~T_vir), uniform pressures and subsonic velocities. This transition occurs when the outer CGM (~0.5 R_vir) is already subsonic and has a temperature ~T_vir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that inner CGM virialization coincides with abrupt changes in the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disc, star formation transitions from `bursty' to `steady,' and stellar-driven galaxy-scale outflows are suppressed. Our results thus suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic discs, rather than with the quenching of star formation as often assumed. △ Less

Submitted 29 December, 2020; v1 submitted 24 June, 2020; originally announced June 2020.

Comments: Accepted for publication in ApJ. Online versions of figures 13 and 15 are available at https://sites.northwestern.edu/jonathanstern/online-figures/

Survey of Extremely High-velocity Outflows in Sloan Digital Sky Survey Quasars

Authors: Paola Rodríguez Hidalgo, Abdul Moiz Khatri, Patrick B. Hall, Sean Haas, Carla Quintero, Viraja Khatu, Griffin Kowash, Norm Murray

Abstract: We present a survey of extremely high-velocity outflows (EHVOs) in quasars, defined by speeds between 0.1c and 0.2c. This region of the parameter space has not been included in previous surveys, but it might present the biggest challenge for theoretical models and it might be a large contributor to feedback due to the outflows' potentially large kinetic power. Using the Sloan Digital Sky Survey, w… ▽ More We present a survey of extremely high-velocity outflows (EHVOs) in quasars, defined by speeds between 0.1c and 0.2c. This region of the parameter space has not been included in previous surveys, but it might present the biggest challenge for theoretical models and it might be a large contributor to feedback due to the outflows' potentially large kinetic power. Using the Sloan Digital Sky Survey, we find 40 quasar spectra with broad EHVO CIV absorption, 10 times more than the number of previously known cases. We characterize the EHVO absorption and find that in 26 cases, CIV is accompanied by NV and/or OVI absorption. We find that EHVO quasars lack HeII emission and have overall larger bolometric luminosities and black hole masses than those of their parent sample and BALQSOs, while we do not find significant differences in their Eddington ratios. We also report a trend toward larger black hole masses as the velocity of the outflowing gas increases in the BALQSOs in our sample. The overall larger Lbol and lack of HeII emission of EHVO quasars suggest that radiation is likely driving these outflows. We find a potential evolutionary effect as EHVO quasars seem to be more predominant at large redshifts. We estimate that the kinetic power of these outflows may be similar to or even larger than that of the outflows from BALQSOs as the velocity factor increases this parameter by 1-2.5 orders of magnitude. Further study of EHVO quasars will help improve our understanding of quasar physics. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: 26 pages, accepted for publication in the Astrophysical Journal

Pressure balance in the multiphase ISM of cosmologically simulated disk galaxies

Authors: Alexander B. Gurvich, Claude-André Faucher-Giguère, Alexander J. Richings, Philip F. Hopkins, Michael Y. Grudić, Zachary Hafen, Sarah Wellons, Jonathan Stern, Eliot Quataert, T. K. Chan, Matthew E. Orr, Dušan Kereš, Andrew Wetzel, Christopher C. Hayward, Sarah R. Loebman, Norman Murray

Abstract: Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. W… ▽ More Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure. △ Less

Submitted 16 November, 2020; v1 submitted 26 May, 2020; originally announced May 2020.

Comments: 22 pages, 16 figures, accepted in MNRAS

Reproducing the CO-to-H$_2$ conversion factor in cosmological simulations of Milky Way-mass galaxies

Authors: Laura C. Keating, Alexander J. Richings, Norman Murray, Claude-Andre Faucher-Giguere, Philip F. Hopkins, Andrew Wetzel, Dusan Keres, Samantha Benincasa, Robert Feldmann, Sarah Loebman, Matthew E. Orr

Abstract: We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results… ▽ More We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds. △ Less

Submitted 28 September, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: Accepted to MNRAS. 16 pages, 10 figures

Variations in the slope of the resolved star-forming main sequence: a tool for constraining the mass of star-forming regions

Authors: Maan H. Hani, Christopher C. Hayward, Matthew E. Orr, Sara L. Ellison, Paul Torrey, Norm Murray, Andrew Wetzel, Claude-André Faucher-Giguère

Abstract: The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS'; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in… ▽ More The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS'; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies' rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps. △ Less

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

Comments: Accepted for publication in MNRAS Letters

Direct Detection of Black Hole-Driven Turbulence in the Centers of Galaxy Clusters

Authors: Yuan Li, Marie-Lou Gendron-Marsolais, Irina Zhuravleva, Siyao Xu, Aurora Simionescu, Grant R. Tremblay, Cassandra Lochhaas, Greg L. Bryan, Eliot Quataert, Norman W. Murray, Alessandro Boselli, Julie Hlavacek-Larrondo, Yong Zheng, Matteo Fossati, Miao Li, Eric Emsellem, Marc Sarzi, Lev Arzamasskiy, Ethan T. Vishniac

Abstract: Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide ra… ▽ More Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma. △ Less

Submitted 14 November, 2019; originally announced November 2019.

Comments: 10 pages, 5 figures. Submitted to ApJL

Live Fast, Die Young: GMC lifetimes in the FIRE cosmological simulations of Milky Way-mass galaxies

Authors: Samantha M. Benincasa, Sarah R. Loebman, Andrew Wetzel, Philip F. Hopkins, Norman Murray, Matthew A. Bellardini, Claude-André Faucher-Giguère, Dávid Guszejnov, Matthew Orr

Abstract: We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution… ▽ More We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 2$\%$ of clouds living longer than 20 Myr. We find decreasing GMC lifetimes with increasing virial parameter, and weakly increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting. △ Less

Submitted 2 December, 2020; v1 submitted 12 November, 2019; originally announced November 2019.

Comments: accepted to MNRAS

Swirls of FIRE: Spatially Resolved Gas Velocity Dispersions and Star Formation Rates in FIRE-2 Disk Environments

Authors: Matthew E. Orr, Christopher C. Hayward, Anne M. Medling, Philip F. Hopkins, Norman Murray, Jorge L. Pineda, Claude-André Faucher-Giguère, Dušan Kereš, Kung-Yi Su

Abstract: We study the spatially resolved (sub-kpc) gas velocity dispersion ($σ$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $σ\approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show th… ▽ More We study the spatially resolved (sub-kpc) gas velocity dispersion ($σ$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $σ\approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant $σ$. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher $σ$ at constant SFR. The limits of the $σ$-$Σ_{\rm SFR}$ relation correspond to the onset of strong outflows. We see evidence of "on-off" cycles of star formation in the simulations, corresponding to feedback injection timescales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disk (Toomre's $Q =1$) model predictions, and the data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e., ${\rm 1\%}$ per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of $σ$, it appears to be strongly subdominant to stellar feedback in the simulated galaxy disks. △ Less

Submitted 31 October, 2019; originally announced November 2019.

Comments: 18 pages, 12 figures, submitted to MNRAS

The Fates of the Circumgalactic Medium in the FIRE Simulations

Authors: Z. Hafen, C. -A. Faucher-Giguere, D. Angles-Alcazar, J. Stern, D. Keres, C. Esmerian, A. Wetzel, K. El-Badry, T. K. Chan, N. Murray

Abstract: We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as its location at z=… ▽ More We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 - 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its "proximate" fate), as well as its location at z=0 (its "ultimate" fate). Of the CGM at z=2, about half is found in the ISM or stars of the central galaxy by z=0 in Mh(z=2) ~ 5e11 Msun halos, but most of the CGM in lower-mass halos ends up in the IGM. This is so even though most of the CGM in M_h(z=2) ~ 5e10 Msun halos first accretes onto the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z=0.25 remains in the CGM by z=0 at all halo masses analyzed. Of the CGM gas that subsequently accretes onto the central galaxy in the progenitors of Mh(z=0) ~10^12 Msun halos, most of it is cool (T~10^4 K) at z=2 but hot (~Tvir) at z=0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z=0.25 and z=2 >~80% of the cool gas in Mh >~ 10^11 Msun halos will accrete onto a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: 16 pages, 10 figures

The CCAT-Prime Submillimeter Observatory

Authors: Manuel Aravena, Jason Austermann, Kaustuv Basu, Nicholas Battaglia, Benjamin Beringue, Frank Bertoldi, J. Richard Bond, Patrick Breysse, Ricardo Bustos, Scott Chapman, Steve Choi, Dongwoo Chung, Nicholas Cothard, Bradley Dober, Cody Duell, Shannon Duff, Rolando Dunner, Jens Erler, Michel Fich, Laura Fissel, Simon Foreman, Patricio Gallardo, Jiansong Gao, Riccardo Giovanelli, Urs Graf , et al. (31 additional authors not shown)

Abstract: The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capa… ▽ More The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys. △ Less

Submitted 5 September, 2019; originally announced September 2019.

Comments: Astro2020 APC White Paper

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

Authors: Suoqing Ji, T. K. Chan, Cameron B. Hummels, Philip F. Hopkins, Jonathan Stern, Dušan Kereš, Eliot Quataert, Claude-André Faucher-Giguère, Norman Murray

Abstract: We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity… ▽ More We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $κ\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $γ$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these "CR-dominated" halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the "low" and "mid" ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase. △ Less

Submitted 6 August, 2020; v1 submitted 30 August, 2019; originally announced September 2019.

Comments: 17 pages, 14 figures, published version on MNRAS

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

Authors: Xiangcheng Ma, Michael Y. Grudić, Eliot Quataert, Philip F. Hopkins, Claude-André Faucher-Giguère, Michael Boylan-Kolchin, Andrew Wetzel, Ji-hoon Kim, Norman Murray, Dušan Kereš

Abstract: We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at… ▽ More We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies. △ Less

Submitted 19 February, 2020; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 18 pages, 17 figures, accepted to MNRAS, high-quality images and/or animations for Figs 3,4,5,8,9 are available at http://www.tapir.caltech.edu/~xchma/HiZFIRE/globular/