Showing 1–50 of 186 results for author: Stevens, A

Astronomy's climate emissions: Global travel to scientific meetings in 2019

Authors: Andrea Gokus, Knud Jahnke, Paul M Woods, Vanessa A Moss, Volker Ossenkopf-Okada, Elena Sacchi, Adam R H Stevens, Leonard Burtscher, Cenk Kayhan, Hannah Dalgleish, Victoria Grinberg, Travis A Rector, Jan Rybizki, Jacob White

Abstract: Travel to academic conferences -- where international flights are the norm -- is responsible for a sizeable fraction of the greenhouse gas (GHG) emissions associated with academic work. In order to provide a benchmark for comparison with other fields, as well as for future reduction strategies and assessments, we estimate the CO2-equivalent emissions for conference travel in the field of astronomy… ▽ More Travel to academic conferences -- where international flights are the norm -- is responsible for a sizeable fraction of the greenhouse gas (GHG) emissions associated with academic work. In order to provide a benchmark for comparison with other fields, as well as for future reduction strategies and assessments, we estimate the CO2-equivalent emissions for conference travel in the field of astronomy for the prepandemic year 2019. The GHG emission of the international astronomical community's 362 conferences and schools in 2019 amounted to 42,500 tCO2e, assuming a radiative-forcing index factor of 1.95 for air travel. This equates to an average of 1.0 $\pm$ 0.6 tCO2e per participant per meeting. The total travel distance adds up to roughly 1.5 Astronomical Units, that is, 1.5 times the distance between the Earth and the Sun. We present scenarios for the reduction of this value, for instance with virtual conferencing or hub models, while still prioritizing the benefits conferences bring to the scientific community. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: Supplementary material is available at PNAS Nexus online: https://academic.oup.com/pnasnexus/article/3/5/pgae143/7659884

Journal ref: PNAS Nexus, Volume 3, Issue 5, May 2024, pgae143

New constraints on ultraheavy dark matter from the LZ experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer, C. A. J. Brew , et al. (174 additional authors not shown)

Abstract: Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal f… ▽ More Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal freeze-out, but they have generally been less explored experimentally. In this work, we present a re-analysis of the first science run (SR1) of the LZ experiment, with an exposure of $0.9$ tonne$\times$year, to search for ultraheavy particle dark matter. The signal topology consists of multiple energy deposits in the active region of the detector forming a straight line, from which the velocity of the incoming particle can be reconstructed on an event-by-event basis. Zero events with this topology were observed after applying the data selection calibrated on a simulated sample of signal-like events. New experimental constraints are derived, which rule out previously unexplored regions of the dark matter parameter space of spin-independent interactions beyond a mass of 10$^{17}$ GeV/$c^2$. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 9 pages, 7 figures

Journal ref: Phys. Rev. D 109, 112010 (2024)

Dark Sage: Next-generation semi-analytic galaxy evolution with multidimensional structure and minimal free parameters

Authors: Adam R. H. Stevens, Manodeep Sinha, Alexander Rohl, Mawson W. Sammons, Boryana Hadzhiyska, César Hernández-Aguayo, Lars Hernquist

Abstract: After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a… ▽ More After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a mass-loading factor for every galaxy. Building on the model's resolved angular-momentum structure of galaxies, we now consider the heating of stellar discs, delivering predictions for disc structure both radially and vertically. We add a further dimension to stellar discs by tracking the distribution of stellar ages in each annulus. Each annulus--age bin has its own velocity dispersion and metallicity evolved in the model. This allows Dark Sage to make structural predictions for galaxies that previously only hydrodynamic simulations could. We present the model as run on the merger trees of the highest-resolution gravity-only simulation of the MillenniumTNG suite. Despite its additional complexity relative to other SAMs, Dark Sage only has three free parameters, the least of any SAM, which we calibrate exclusively against the cosmic star formation history and the $z=0$ stellar and HI mass functions using a particle-swarm optimisation method. The Dark Sage codebase, written in C and Python, is publicly available at https://github.com/arhstevens/DarkSage △ Less

Submitted 14 February, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Updates post referee's report. Get in touch for data access. 34 pages before references; 14 figures, 2 tables, 130 equations

Journal ref: PASA, 2024

First Constraints on WIMP-Nucleon Effective Field Theory Couplings in an Extended Energy Region From LUX-ZEPLIN

Authors: LZ Collaboration, J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger , et al. (175 additional authors not shown)

Abstract: Following the first science results of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time projection chamber operating from the Sanford Underground Research Facility in Lead, South Dakota, USA, we report the initial limits on a model-independent non-relativistic effective field theory describing the complete set of possible interactions of a weakly interacting massive particle (WIMP) with a n… ▽ More Following the first science results of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time projection chamber operating from the Sanford Underground Research Facility in Lead, South Dakota, USA, we report the initial limits on a model-independent non-relativistic effective field theory describing the complete set of possible interactions of a weakly interacting massive particle (WIMP) with a nucleon. These results utilize the same 5.5 t fiducial mass and 60 live days of exposure collected for the LZ spin-independent and spin-dependent analyses while extending the upper limit of the energy region of interest by a factor of 7.5 to 270 keVnr. No significant excess in this high energy region is observed. Using a profile-likelihood ratio analysis, we report 90% confidence level exclusion limits on the coupling of each individual non-relativistic WIMP-nucleon operator for both elastic and inelastic interactions in the isoscalar and isovector bases. △ Less

Submitted 26 February, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 17 pages 11 figures

Journal ref: Phys. Rev. D 109, 092003 (2024)

Why do semi-analytic models predict higher scatter in the stellar mass-halo mass relation than cosmological hydrodynamic simulations?

Authors: Antonio J. Porras-Valverde, John C. Forbes, Rachel S. Somerville, Adam R. H. Stevens, Kelly Holley-Bockelmann, Andreas A. Berlind, Shy Genel

Abstract: Semi-analytic models (SAMs) systematically predict higher stellar-mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAM Dark Sage. We design two black hole formation models that approximate results from the IllustrisTNG 300-1 hydrodynamical simu… ▽ More Semi-analytic models (SAMs) systematically predict higher stellar-mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAM Dark Sage. We design two black hole formation models that approximate results from the IllustrisTNG 300-1 hydrodynamical simulation. In the first model, we assign a fixed black hole mass of $10^{6}\, \mathrm{M}_{\odot}$ to every halo that reaches $10^{10.5}\, \mathrm{M}_{\odot}$. In the second model, we disregard any black hole growth as implemented in the standard Dark Sage model. Instead, we force all black hole masses to follow the median black hole mass-halo mass relation in IllustrisTNG 300-1 with a fixed scatter. We find that each model on its own does not significantly reduce the scatter in stellar mass. To do this, we replace the native Dark Sage AGN feedback model with a simple model where we turn off cooling for galaxies with black hole masses above $10^{8}\, \mathrm{M}_{\odot}$. With this additional modification, the SMBH seeding and fixed conditional distribution models find a significant reduction in the scatter in stellar mass at halo masses between $10^{11-14}\, \mathrm{M}_{\odot}$. These results suggest that AGN feedback in SAMs acts in a qualitatively different way than feedback implemented in cosmological simulations. Either or both may require substantial modification to match the empirically inferred scatter in the Stellar Mass Halo Mass Relation (SMHMR). △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 21 pages, 16 figures

VERTICO and IllustrisTNG: The spatially resolved effects of environment on galactic gas

Authors: Adam R. H. Stevens, Toby Brown, Benedikt Diemer, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Yannick M. Bahé, Alessandro Boselli, Timothy A. Davis, Pascal J. Elahi, Sara L. Ellison, María J. Jiménez-Donaire, Ian D. Roberts, Kristine Spekkens, Vicente Villanueva, Adam B. Watts, Christine D. Wilson, Nikki Zabel

Abstract: It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies w… ▽ More It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, accepted in ApJL

VERTICO VII: Environmental quenching caused by suppression of molecular gas content and star formation efficiency in Virgo Cluster galaxies

Authors: Toby Brown, Ian D. Roberts, Mallory Thorp, Sara L. Ellison, Nikki Zabel, Christine D. Wilson, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Alessandro Boselli, Aeree Chung, Luca Cortese, Barbara Catinella, Timothy A. Davis, María J. Jiménez-Donaire, Claudia D. P. Lagos, Bumhyun Lee, Laura C. Parker, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Vicente Villanueva, Adam B. Watts

Abstract: We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with resp… ▽ More We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with respect to both HI-normal cluster and field galaxies (0.5 dex), suggesting that mechanisms regulating the global HI content are responsible for quenching local star formation. We demonstrate that the observed quenching in HI-poor galaxies is caused by environmental processes such as ram pressure stripping (RPS) simultaneously reducing molecular gas surface density and star formation efficiency (SFE), compared to regions in HI-normal systems (by 0.38 and 0.22 dex, respectively). We observe systematically elevated SFRs that are driven by increased molecular gas surface densities at fixed stellar mass surface density in the outskirts of early-stage RPS galaxies, while SFE remains unchanged with respect to the field sample. We quantify how RPS and starvation affect the star formation cycle of inner and outer galaxy discs as they are processed by the cluster. We show both are effective quenching mechanisms with the key difference being that RPS acts upon the galaxy outskirts while starvation regulates the star formation cycle throughout disc, including within the truncation radius. For both processes, the quenching is caused by a simultaneous reduction in molecular gas surface densities and SFE at fixed stellar mass surface density. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

The Galaxy Number Density Profile of Haloes

Authors: Fei Qin, David Parkinson, Adam R. H. Stevens, Cullan Howlett

Abstract: More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new galaxy number density profile of haloes, which makes predictions for t… ▽ More More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new galaxy number density profile of haloes, which makes predictions for the positions of galaxies in the host halo, different to the widely adopted Navarro-Frenk-White (NFW) profile, since galaxies tend to be found more in the outskirts of halos (nearer the virial radius) than an NFW profile. The parameterised galaxy number density profile model of haloes is fit and tested using the DARKSAGE semi-analytic model of galaxy formation. We find that our galaxy number density profile model of haloes can accurately reproduce the halo occupation distribution and galaxy two-point correlation function of the DARKSAGE simulation. We also derive the analytic expressions for the circular velocity and gravitational potential energy for this profile model. We use the SDSS DR10 galaxy group catalogue to validate this galaxy number density profile model of haloes. Compared to the NFW profile, we find that our model more accurately predicts the positions of galaxies in their host halo and the galaxy two-point correlation function. △ Less

Submitted 25 September, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 13 pages. 10 figures. Appear on ApJ

A search for new physics in low-energy electron recoils from the first LZ exposure

Authors: The LZ Collaboration, J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, P. Beltrame, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, G. M. Blockinger , et al. (178 additional authors not shown)

Abstract: The LUX-ZEPLIN (LZ) experiment is a dark matter detector centered on a dual-phase xenon time projection chamber. We report searches for new physics appearing through few-keV-scale electron recoils, using the experiment's first exposure of 60 live days and a fiducial mass of 5.5t. The data are found to be consistent with a background-only hypothesis, and limits are set on models for new physics inc… ▽ More The LUX-ZEPLIN (LZ) experiment is a dark matter detector centered on a dual-phase xenon time projection chamber. We report searches for new physics appearing through few-keV-scale electron recoils, using the experiment's first exposure of 60 live days and a fiducial mass of 5.5t. The data are found to be consistent with a background-only hypothesis, and limits are set on models for new physics including solar axion electron coupling, solar neutrino magnetic moment and millicharge, and electron couplings to galactic axion-like particles and hidden photons. Similar limits are set on weakly interacting massive particle (WIMP) dark matter producing signals through ionized atomic states from the Migdal effect. △ Less

Submitted 9 September, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 13 pages, 10 figures. See https://tinyurl.com/LZDataReleaseRun1ER for a data release related to this paper

Journal ref: Phys. Rev. D 108, 072006 (2023)

Exploring the Angular Momentum -- Atomic Gas Content Connection with EAGLE and IllustrisTNG

Authors: Jennifer A. Hardwick, Luca Cortese, Danail Obreschkow, Claudia Lagos, Adam R. H. Stevens, Barbara Catinella, Lilian Garratt-Smithson

Abstract: We use the EAGLE (Evolution and Assembly of GaLaxies and their Environments) and IllustrisTNG (The Next Generation) cosmological simulations to investigate the properties of the baryonic specific angular momentum (j), baryonic mass (M) and atomic gas fraction ($f_{\rm{atm}}$) plane for nearby galaxies. We find EAGLE and TNG to be in excellent agreement with each other. These simulations are also c… ▽ More We use the EAGLE (Evolution and Assembly of GaLaxies and their Environments) and IllustrisTNG (The Next Generation) cosmological simulations to investigate the properties of the baryonic specific angular momentum (j), baryonic mass (M) and atomic gas fraction ($f_{\rm{atm}}$) plane for nearby galaxies. We find EAGLE and TNG to be in excellent agreement with each other. These simulations are also consistent with the results obtained with xGASS (eXtended GALEX Arecibo SDSS Survey) for gas fractions greater than 0.01. This implies that the disagreements previously identified between xGASS and predictions from simple analytical disc stability arguments also holds true for EAGLE and TNG. For lower gas fraction (the regime currently unconstrained by observations), both simulations deviate from the plane but still maintain good agreement with each other. Despite the challenges posed by resolution limits at low gas fractions, our findings suggest a potential disconnect between angular momentum and gas fraction in the gas-poor regime, implying that not all gas-poor galaxies have low specific angular momentum. △ Less

Submitted 18 September, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 12 pages, 10 figures. Accepted for publication in MNRAS

VERTICO VI: Cold-gas asymmetries in Virgo cluster galaxies

Authors: Ian D. Roberts, Toby Brown, Nikki Zabel, Christine D. Wilson, Aeree Chung, Laura C. Parker, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Ryan Chown, Luca Cortese, Timothy A. Davis, Sara Ellison, Maria Jesus Jimenez-Donaire, Bumhyun Lee, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Mallory Thorp, Vincente Villanueva, Adam B. Watts, Charlotte Welker, Hyein Yoon

Abstract: We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to th… ▽ More We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to the relatively dense molecular gas component, and that these environmental perturbations can compress the cold interstellar medium in cluster galaxies leading to elevated star formation. In this work we observationally test these predictions for star-forming satellite galaxies within the Virgo cluster. We divide our Virgo galaxy sample into HI-normal, HI-tailed, and HI-truncated classes and show, unsurprisingly, that the HI-tailed galaxies have the largest quantitative HI asymmetries. We also compare to a control sample of non-cluster galaxies and find that Virgo galaxies, on average, have HI asymmetries that are 40 +/- 10 per cent larger than the control. There is less separation between control, HI-normal, HI-tailed, and HI-truncated galaxies in terms of H2 asymmetries, and on average, Virgo galaxies have H2 asymmetries that are only marginally (20 +/- 10 per cent) larger than the control sample. We find a weak correlation between HI and H2 asymmetries over our entire sample, but a stronger correlation for those specific galaxies being strongly impacted by environmental perturbations. Finally, we divide the discs of the HI-tailed Virgo galaxies into a leading half and trailing half according to the observed tail direction. We find evidence for excess molecular gas mass on the leading halves of the disc. This excess molecular gas on the leading half is accompanied by an excess in star formation rate such that the depletion time is, on average, unchanged. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 15 pages, 8 figures, 1 table, accepted for publication in A&A

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

Galaxy And Mass Assembly (GAMA): The group HI mass as a function of halo mass

Authors: Ajay Dev, Simon P. Driver, Martin Meyer, Sambit Roychowdhury, Jonghwan Rhee, Adam R. H. Stevens, Claudia del P. Lagos, Joss Bland-Hawthorn, Barbara Catinella, A. M. Hopkins, Jonathan Loveday, Danail Obreschkow, Steven Phillipps, Aaron S. G. Robotham

Abstract: We determine the atomic hydrogen (HI) to halo mass relation (HIHM) using Arecibo Legacy Fast ALFA survey HI data at the location of optically selected groups from the Galaxy and Mass Assembly (GAMA) survey. We make direct HI detections for 37 GAMA groups. Using HI group spectral stacking of 345 groups, we study the group HI content as function of halo mass across a halo mass range of… ▽ More We determine the atomic hydrogen (HI) to halo mass relation (HIHM) using Arecibo Legacy Fast ALFA survey HI data at the location of optically selected groups from the Galaxy and Mass Assembly (GAMA) survey. We make direct HI detections for 37 GAMA groups. Using HI group spectral stacking of 345 groups, we study the group HI content as function of halo mass across a halo mass range of $10^{11} - 10^{14.7}\text{ M}_\odot$. We also correct our results for Eddington bias. We find that the group HI mass generally rises as a function of halo mass from $1.3\%$ of the halo mass at $10^{11.6} \text{M}_\odot$ to $0.4\%$ at $10^{13.7} \text{M}_\odot$ with some indication of flattening towards the high-mass end. Despite the differences in optical survey limits, group catalogues, and halo mass estimation methods, our results are consistent with previous group HI-stacking studies. Our results are also consistent with mock observations from SHARK and IllustrisTNG. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: Accepted in MNRAS; 18 pages, 12 figures

VERTICO V: The environmentally driven evolution of the inner cold gas discs of Virgo cluster galaxies

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Toby Brown, Christine D. Wilson, Nikki Zabel, Ian D. Roberts, Timothy A. Davis, Mallory Thorp, Aeree Chung, Adam R. H. Stevens, Sara L. Ellison, Kristine Spekkens, Laura C. Parker, Yannick M. Bahé, Vicente Villanueva, María Jiménez-Donaire, Dhruv Bisaria, Alessandro Boselli, Alberto D. Bolatto, Bumhyun Lee

Abstract: The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virg… ▽ More The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ($Σ_{\star}$), HI ($Σ_\mathrm{HI}$), and molecular gas ($Σ_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average $Σ_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower $Σ_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical $Σ_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of $Σ_\mathrm{HI}$ or $Σ_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of $Σ_\mathrm{mol}$-to-$Σ_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PASA. Main text is 19 pages including 12 figures and 3 tables, plus 1 appendix. A 2.5 min, high-level summary can be found at https://youtu.be/7djMmVEpDVc

Driving action on the climate crisis through Astronomers for Planet Earth and beyond

Authors: Adam R. H. Stevens, Vanessa A. Moss

Abstract: While an astronomer's job is typically to look out from Earth, the seriousness of the climate crisis has meant a shift in many astronomers' focus. Astronomers are starting to consider how our resource requirements may contribute to this crisis and how we may better conduct our research in a more environmentally sustainable fashion. Astronomers for Planet Earth is an international organisation (mor… ▽ More While an astronomer's job is typically to look out from Earth, the seriousness of the climate crisis has meant a shift in many astronomers' focus. Astronomers are starting to consider how our resource requirements may contribute to this crisis and how we may better conduct our research in a more environmentally sustainable fashion. Astronomers for Planet Earth is an international organisation (more than 1,700 members from over 70 countries as of November 2022) that seeks to answer the call for sustainability to be at the heart of astronomers' practices. In this article, we review the organisation's history, summarising the proactive, collaborative efforts and research into astronomy sustainability conducted by its members. We update the state of affairs with respect to the carbon footprint of astronomy research, noting an improvement in renewable energy powering supercomputing facilities in Australia, reducing that component of our footprint by a factor of 2--3. We discuss how, despite accelerated changes made throughout the pandemic, we still must address the format of our meetings. Using recent annual meetings of the Australian and European astronomical societies as examples, we demonstrate that the more online-focussed a meeting is, the greater its attendance and the lower its emissions. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 11 pages, 2 figures, refereed and accepted to CAPjournal for special issue on CAP2022 conference

The relationship between cluster environment and molecular gas content of star-forming galaxies in the EAGLE simulation

Authors: Aditya Manuwal, Adam R. H. Stevens

Abstract: We employ the EAGLE hydrodynamical simulation to uncover the relationship between cluster environment and $\rm H_2$ content of star-forming galaxies at redshifts spanning $0\leq z\leq 1$. To do so, we divide the star-forming sample into those that are bound to clusters and those that are not. We find that, at any given redshift, the galaxies in clusters generally have less $\rm H_2$ than their non… ▽ More We employ the EAGLE hydrodynamical simulation to uncover the relationship between cluster environment and $\rm H_2$ content of star-forming galaxies at redshifts spanning $0\leq z\leq 1$. To do so, we divide the star-forming sample into those that are bound to clusters and those that are not. We find that, at any given redshift, the galaxies in clusters generally have less $\rm H_2$ than their non-cluster counterparts with the same stellar mass (corresponding to an offset of $\lesssim 0.5$ dex), but this offset varies with stellar mass and is virtually absent at $M_\star\lesssim10^{9.3}~{\rm M}_\odot$. The $\rm H_2$ deficit in star-forming cluster galaxies can be traced back to a decline in their $\rm H_2$ content that commenced after first infall into a cluster, which occurred later than a typical cluster galaxy. Evolution of the full cluster population after infall is generally consistent with `slow-then-rapid' quenching, but galaxies with $M_\star\lesssim 10^{9.5}~{\rm M}_\odot$ exhibit rapid quenching. Unlike most cluster galaxies, star-forming ones were not pre-processed in groups prior to being accreted by clusters. For both of these cluster samples, the star formation efficiency remained oblivious to the infall. We track the particles associated with star-forming cluster galaxies and attribute the drop in $\rm H_2$ mass after infall to poor replenishment, depletion due to star formation, and stripping of $\rm H_2$ in cluster environments. These results provide predictions for future surveys, along with support and theoretical insights for existing molecular gas observations that suggest there is less $\rm H_2$ in cluster galaxies. △ Less

Submitted 23 May, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: Accepted for publication in MNRAS

VERTICO III: The Kennicutt-Schmidt relation in Virgo cluster galaxies

Authors: M. J. Jiménez-Donaire, T. Brown, C. D. Wilson, I. D. Roberts, N. Zabel, S. L. Ellison, M. Thorp, V. Villanueva, R. Chown, D. Bisaria, A. D. Bolatto, A. Boselli, B. Catinella, A. Chung, L. Cortese, T. A. Davis, C. D. P. Lagos, B. Lee, L. C. Parker, K. Spekkens, A. R. H. Stevens, J. Sun

Abstract: In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies.… ▽ More In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies. We use new ACA and TP observations from the VERTICO-ALMA Large Program at 720pc resolution to resolve the molecular gas content, as traced by the 12CO(2-1) transition, across the disks of 37 spiral galaxies in the Virgo cluster. In combination with archival observations, we estimate the parameters of the KS relation for the entire ensemble of galaxies, and within individual galaxies. We find the KS slope for the entire population to be N=0.97+/-0.07, with a characteristic molecular gas depletion time of 1.86Gyr for our full sample, in agreement with previous work in isolated star-forming galaxies. In individual galaxies, we find KS slopes ranging between 0.69 and 1.40, and typical star formation efficiencies (SFE) that can vary from galaxy to galaxy by a factor of ~4. These galaxy-to-galaxy variations account for ~0.20dex in scatter in the ensemble KS relation, which is characterized by a 0.42dex scatter. We find that the HI-deficient galaxies in the Virgo cluster show a steeper resolved KS relation and lower molecular gas efficiencies than HI-normal cluster galaxies. While the molecular gas content in Virgo cluster galaxies appears to behave similarly to that in isolated galaxies, our VERTICO sample shows that cluster environments play a key role in regulating star formation. The environmental mechanisms affecting the HI galaxy content also have a direct impact in the SFE of molecular gas in cluster galaxies, leading to longer depletion times in HI-deficient members. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 671, A3 (2023)

VERTICO IV: Environmental Effects on the Gas Distribution and Star Formation Efficiency of Virgo Cluster Spirals

Authors: Vicente Villanueva, Alberto D. Bolatto, Stuart Vogel, Tobias Brown, Christine D. Wilson, Nikki Zabel, Sara Ellison, Adam R. H. Stevens, Maria Jesus Jimenez Donaire, Kristine Spekkens, Mallory Thorp, Timothy A. Davis, Laura C. Parker, Ian D. Roberts, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Aeree Chung, Luca Cortese, Bumhyun Lee, Adam Watts

Abstract: We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating th… ▽ More We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We compute $R_{\rm mol}$ as a function of different physical quantities. While the spatially-resolved $R_{\rm mol}$ on average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radius $R_{\rm e}$, $R_{\rm mol}(r<R_{\rm e})$, shows a systematic increase with the level of H$_{\rm I}$, truncation and/or asymmetry (H$_{\rm I}$ perturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios within $R_{\rm e}$, $R^{\rm mol}_{\star}(r<R_{\rm e})$ and $R^{\rm atom}_{\star}(r<R_{\rm e})$, shows that VERTICO galaxies have increasingly lower $R^{\rm atom}_{\star}(r<R_{\rm e})$ for larger levels of H$_{\rm I}$perturbation (compared to field galaxies matched in stellar mass), but no significant change in $R^{\rm mol}_{\star}(r<R_{\rm e})$. We also measure a clear systematic decrease of the SFE$_{\rm mol}$ within $R_{\rm e}$, SFE$_{\rm mol}(r<R_{\rm e})$, with increasingly perturbed H$_{\rm I}$. Therefore, compared to galaxies from the field, VERTICO galaxies are more compact in CO emission in relation to their stellar distribution, but increasingly perturbed atomic gas increases their $R_{\rm mol}$ and decreases the efficiency with which their molecular gas forms stars. (abridged) △ Less

Submitted 1 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

HI HOD -- I. The Halo Occupation Distribution of HI Galaxies

Authors: Fei Qin, Cullan Howlett, Adam R. H. Stevens, David Parkinson

Abstract: The next generation of galaxy surveys will provide more precise measurements of galaxy clustering than have previously been possible. The 21-cm radio signals that are emitted from neutral atomic hydrogen (HI) gas will be detected by large-area radio surveys such as WALLABY and the SKA, and deliver galaxy positions and velocities that can be used to measure galaxy clustering statistics. But, to har… ▽ More The next generation of galaxy surveys will provide more precise measurements of galaxy clustering than have previously been possible. The 21-cm radio signals that are emitted from neutral atomic hydrogen (HI) gas will be detected by large-area radio surveys such as WALLABY and the SKA, and deliver galaxy positions and velocities that can be used to measure galaxy clustering statistics. But, to harness this information to improve our cosmological understanding, and learn about the physics of dark matter and dark energy, we need to accurately model the manner in which galaxies detected in HI trace the underlying matter distribution of the Universe. For this purpose, we develop a new HI-based Halo Occupation Distribution (HOD) model, which makes predictions for the number of galaxies present in dark matter halos conditional on their HI mass. The parameterised HOD model is fit and validated using the Dark Sage semi-analytic model, where we show that the HOD parameters can be modelled by simple linear and quadratic functions of HI mass. However, we also find that the clustering predicted by the HOD depends sensitively on the radial distributions of the HI galaxies within their host dark matter halos, which does not follow the NFW profile in the Dark Sage simulation. As such, this work enables -- for the first time -- a simple prescription for placing galaxies of different HI mass within dark matter halos in a way that is able to reproduce the HI mass-dependent galaxy clustering and HI mass function simultaneously and without requiring knowledge of the optical properties of the galaxies. Further efforts are required to demonstrate that this model can be used to produce large ensembles of mock galaxy catalogues for upcoming surveys. △ Less

Submitted 1 September, 2022; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 18 pages, 13 figures, 1 table. Published in ApJ

Early results from GLASS-JWST XIV: A first morphological atlas of the 1 < z < 5 Universe in the rest-frame optical

Authors: Colin Jacobs, Karl Glazebrook, Antonello Calabrò, Tommaso Treu, Themiya Nanayakkara, Tucker Jones, Emiliano Merlin, Roberto G. Abraham, Adam R H Stevens, Benedetta Vulcani, Lilan Yang, Andrea Bonchi, Marusa Bradac, Marco Castellano, Adriano Fontana, Matthew A. Malkan, Charlotte A Mason, Takahiro Morishita, Diego Paris, Michele Trenti, Danilo Marchesini, Xin Wang, Paola Santini

Abstract: We present a rest-frame optical morphological analysis of galaxies observed with the NIRCam imager on the James Webb Space Telescope (JWST) as part of the GLASS-JWST Early Release Science program. We select 388 sources at redshifts \(0.8 < z < 5.4\) and use the seven 0.9--5\micron\ NIRCam filters to generate rest-frame $gri$ composite color images, and conduct visual morphological classification.… ▽ More We present a rest-frame optical morphological analysis of galaxies observed with the NIRCam imager on the James Webb Space Telescope (JWST) as part of the GLASS-JWST Early Release Science program. We select 388 sources at redshifts \(0.8 < z < 5.4\) and use the seven 0.9--5\micron\ NIRCam filters to generate rest-frame $gri$ composite color images, and conduct visual morphological classification. Compared to HST-based work we find a higher incidence of disks and bulges than expected at $z>1.5$, revealed by rest frame optical imaging. We detect 123 clear disks (58 at $z>1.5$) of which 76 have bulges. No evolution of bulge fraction with redshift is evident: 61\% at \(z<2\) (\(N=110\)) versus 60\% at \(z\geq2\) (\(N=13\)). A stellar mass dependence is evident, with bulges visible in 80\% of all disk galaxies with mass \(> 10^{9.5}\, {\rm M}_{\odot}\) (\(N=41\)) but only 52\% at \(M < 10^{9.5}\, {\rm M}_{\odot}\) (\(N=82\)). We supplement visual morphologies with non-parametric measurements of Gini and Asymmetry coefficients in the rest-frame $i$-band. Our sources are more asymmetric than local galaxies, with slightly higher Gini values. When compared to high-z rest-frame ultraviolet measurements with Hubble Space Telescope, JWST shows more regular morphological types such as disks, bulges and spiral arms at $z>1.5$, with smoother (i.e. lower Gini) and more symmetrical light distributions. △ Less

Submitted 17 April, 2023; v1 submitted 12 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ Letters

First Dark Matter Search Results from the LUX-ZEPLIN (LZ) Experiment

Authors: J. Aalbers, D. S. Akerib, C. W. Akerlof, A. K. Al Musalhi, F. Alder, A. Alqahtani, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, S. Azadi, A. J. Bailey, A. Baker, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, M. J. Barry, J. Barthel, D. Bauer, A. Baxter , et al. (322 additional authors not shown)

Abstract: The LUX-ZEPLIN experiment is a dark matter detector centered on a dual-phase xenon time projection chamber operating at the Sanford Underground Research Facility in Lead, South Dakota, USA. This Letter reports results from LUX-ZEPLIN's first search for weakly interacting massive particles (WIMPs) with an exposure of 60~live days using a fiducial mass of 5.5 t. A profile-likelihood ratio analysis s… ▽ More The LUX-ZEPLIN experiment is a dark matter detector centered on a dual-phase xenon time projection chamber operating at the Sanford Underground Research Facility in Lead, South Dakota, USA. This Letter reports results from LUX-ZEPLIN's first search for weakly interacting massive particles (WIMPs) with an exposure of 60~live days using a fiducial mass of 5.5 t. A profile-likelihood ratio analysis shows the data to be consistent with a background-only hypothesis, setting new limits on spin-independent WIMP-nucleon, spin-dependent WIMP-neutron, and spin-dependent WIMP-proton cross sections for WIMP masses above 9 GeV/c$^2$. The most stringent limit is set for spin-independent scattering at 36 GeV/c$^2$, rejecting cross sections above 9.2$\times 10^{-48}$ cm$^2$ at the 90% confidence level. △ Less

Submitted 2 August, 2023; v1 submitted 8 July, 2022; originally announced July 2022.

Comments: 9 pages, 8 figures. See https://doi.org/10.1103/PhysRevLett.131.041002 for a data release related to this paper

Journal ref: Phys. Rev. Lett. 131, 041002 (2023)

An orbital perspective on the starvation, stripping, and quenching of satellite galaxies in the EAGLE simulations

Authors: Ruby J. Wright, Claudia del P. Lagos, Chris Power, Adam R. H. Stevens, Luca Cortese, Rhys J. J. Poulton

Abstract: Using the EAGLE suite of simulations, we demonstrate that both cold gas stripping {\it and} starvation of gas inflow play an important role in quenching satellite galaxies across a range of stellar and halo masses, $M_{\star}$ and $M_{200}$. By quantifying the balance between gas inflows, outflows, and star formation rates, we show that even at $z=2$, only $\approx30\%$ of satellite galaxies are a… ▽ More Using the EAGLE suite of simulations, we demonstrate that both cold gas stripping {\it and} starvation of gas inflow play an important role in quenching satellite galaxies across a range of stellar and halo masses, $M_{\star}$ and $M_{200}$. By quantifying the balance between gas inflows, outflows, and star formation rates, we show that even at $z=2$, only $\approx30\%$ of satellite galaxies are able to maintain equilibrium or grow their reservoir of cool gas - compared to $\approx50\%$ of central galaxies at this redshift. We find that the number of orbits completed by a satellite is a very good predictor of its quenching, even more so than the time since infall. On average, we show that intermediate-mass satellites with $M_{\star}$ between $10^{9}{\rm M}_{\odot}-10^{10}{\rm M}_{\odot}$ will be quenched at first pericenter in massive group environments, $M_{200}>10^{13.5}{\rm M}_{\odot}$; and will be quenched at second pericenter in less massive group environments, $M_{200}<10^{13.5}{\rm M}_{\odot}$. On average, more massive satellites ($M_{\star}>10^{10}{\rm M}_{\odot}$) experience longer depletion time-scales, being quenched between first and second pericenters in massive groups; while in smaller group environments, just $\approx30\%$ will be quenched even after two orbits. Our results suggest that while starvation alone may be enough to slowly quench satellite galaxies, direct gas stripping, particularly at pericenters, is required to produce the short quenching time-scales exhibited in the simulation. △ Less

Submitted 19 July, 2022; v1 submitted 17 May, 2022; originally announced May 2022.

Comments: Accepted for publication in MNRAS

VERTICO II: effects of HI-identified environmental mechanisms on molecular gas

Authors: Nikki Zabel, Toby Brown, Christine D. Wilson, Timothy A. Davis, Luca Cortese, Laura C. Parker, Alessandro Boselli, Barbara Catinella, Ryan Chown, Aeree Chung, Tirna Deb, Sara L. Ellison, María J. Jiménez-Donaire, Bumhyun Lee, Ian D. Roberts, Kristine Spekkens, Adam R. H. Stevens, Mallory Thorp, Stephanie Tonnesen, Vicente Villanueva

Abstract: In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a s… ▽ More In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a stellar mass range of 9 < log M*/Msol < 11. We study molecular gas radial profiles, isodensity radii, and surface densities as a function of galaxy HI deficiency and morphology. There is a weak correlation between global HI and H2 deficiencies, and resolved properties of molecular gas correlate with HI deficiency: galaxies that have large HI deficiencies have relatively steep and truncated molecular gas radial profiles, which is due to the removal of low-surface density molecular gas on the outskirts. Therefore, while the environmental mechanisms observed in HI also affect molecular gas reservoirs, there is only a moderate reduction of the total amount of molecular gas. △ Less

Submitted 12 January, 2023; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Published in ApJ. 22 pages, 6 figures, 1 table, 1 appendix. Erratum accepted for publication in ApJ

A Next-Generation Liquid Xenon Observatory for Dark Matter and Neutrino Physics

Authors: J. Aalbers, K. Abe, V. Aerne, F. Agostini, S. Ahmed Maouloud, D. S. Akerib, D. Yu. Akimov, J. Akshat, A. K. Al Musalhi, F. Alder, S. K. Alsum, L. Althueser, C. S. Amarasinghe, F. D. Amaro, A. Ames, T. J. Anderson, B. Andrieu, N. Angelides, E. Angelino, J. Angevaare, V. C. Antochi, D. Antón Martin, B. Antunovic, E. Aprile, H. M. Araújo , et al. (572 additional authors not shown)

Abstract: The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neut… ▽ More The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neutrinos through neutrinoless double-beta decay and through a variety of astrophysical sources. A next-generation xenon-based detector will therefore be a true multi-purpose observatory to significantly advance particle physics, nuclear physics, astrophysics, solar physics, and cosmology. This review article presents the science cases for such a detector. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 77 pages, 40 figures, 1262 references

Report number: INT-PUB-22-003

Journal ref: J. Phys. G: Nucl. Part. Phys. 50 (2023) 013001

Cold Gas in Massive Galaxies as A Critical Test of Black Hole Feedback Models

Authors: Jingjing Shi, Yingjie Peng, Benedikt Diemer, Adam R. H. Stevens, Annalisa Pillepich, Alvio Renzini, Jing Dou, Yu Gao, Qiusheng Gu, Luis C. Ho, Xu Kong, Claudia del P. Lagos, Di Li, Jiaxuan Li, Roberto Maiolino, Filippo Mannucci, Lizhi Xie, Chengpeng Zhang

Abstract: Black hole feedback has been widely implemented as the key recipe to quench star formation in massive galaxies in modern semi-analytic models and hydrodynamical simulations. As the theoretical details surrounding the accretion and feedback of black holes continue to be refined, various feedback models have been implemented across simulations, with notable differences in their outcomes. Yet, most o… ▽ More Black hole feedback has been widely implemented as the key recipe to quench star formation in massive galaxies in modern semi-analytic models and hydrodynamical simulations. As the theoretical details surrounding the accretion and feedback of black holes continue to be refined, various feedback models have been implemented across simulations, with notable differences in their outcomes. Yet, most of these simulations have successfully reproduced some observations, such as stellar mass function and star formation rate density in the local Universe. We use the recent observation on the change of neutral hydrogen gas mass (including both ${\rm H_2}$ and ${\rm HI}$) with star formation rate of massive central disc galaxies as a critical constraint of black hole feedback models across several simulations. We find that the predictions of IllustrisTNG agree with the observations much better than the other models tested in this work. This favors IllustrisTNG's treatment of active galactic nuclei - where kinetic winds are driven by black holes at low accretion rates - as more plausible amongst those we test. In turn, this also indirectly supports the idea that the massive central disc galaxy population in the local Universe was likely quenched by AGN feedback. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 15 pages, 3+4 figures, accepted for publication in ApJ

Radius Constraints from Reflection Modeling of Cygnus X-2 with NuSTAR and NICER

Authors: R. M. Ludlam, E. M. Cackett, J. A. García, J. M. Miller, A. L. Stevens, A. C. Fabian, J. Homan, M. NG, S. Guillot, D. J. K. Buisson, D. Chakrabarty

Abstract: We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color-color diagram; namely the horizontal branch, normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in t… ▽ More We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color-color diagram; namely the horizontal branch, normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in two different ways that produced a comparable quality fit. The spectra showed clear evidence of a reflection component in the form of a broadened Fe K line, as well as a lower energy emission feature near 1 keV likely due to an ionized plasma located far from the innermost accretion disk. We account for the reflection spectrum with two independent models (relxillns and rdblur*rfxconv). The inferred inclination is in agreement with earlier estimates from optical observations of ellipsoidal light curve modeling (relxillns: $i=67^{\circ}\pm4^{\circ}$, rdblur*rfxconv: $i=60^{\circ}\pm10^{\circ}$). The inner disk radius remains close to the NS ($R_{\rm in}\leq1.15\ R_{\mathrm{ISCO}}$) regardless of the source position along the Z-track or how the 1 keV feature is modeled. Given the optically determined NS mass of $1.71\pm0.21\ M_{\odot}$, this corresponds to a conservative upper limit of $R_{\rm in}\leq19.5$ km for $M=1.92\ M_{\odot}$ or $R_{\rm in}\leq15.3$ km for $M=1.5\ M_{\odot}$. We compare these radius constraints to those obtained from NS gravitational wave merger events and recent NICER pulsar light curve modeling measurements. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ, 13 pages, 5 tables, 5 figures

Cosmogenic production of $^{37}$Ar in the context of the LUX-ZEPLIN experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, A. Baker, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, K. Beattie, E. P. Bernard, A. Bhatti, A. Biekert, T. P. Biesiadzinski , et al. (183 additional authors not shown)

Abstract: We estimate the amount of $^{37}$Ar produced in natural xenon via cosmic ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth's surface. We then calculate the resulting $^{37}$Ar concentration in a 10-tonne payload~(similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage and delivery to the… ▽ More We estimate the amount of $^{37}$Ar produced in natural xenon via cosmic ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth's surface. We then calculate the resulting $^{37}$Ar concentration in a 10-tonne payload~(similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea level production rate of $^{37}$Ar in natural xenon is estimated to be 0.024~atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1~tonne/month, the average $^{37}$Ar activity after 10~tonnes are purified and transported underground is 0.058--0.090~$μ$Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic $^{37}$Ar will appear as a noticeable background in the early science data, while decaying with a 35~day half-life. This newly-noticed production mechanism of $^{37}$Ar should be considered when planning for future liquid xenon-based experiments. △ Less

Submitted 22 March, 2022; v1 submitted 8 January, 2022; originally announced January 2022.

WALLABY Pilot Survey: HI gas disc truncation and star formation of galaxies falling into the Hydra I cluster

Authors: T. N. Reynolds, B. Catinella, L. Cortese, T. Westmeier, G. R. Meurer, L. Shao, D. Obreschkow, J. Román, L. Verdes-Montenegro, N. Deg, H. Dénes, B. -Q. For, D. Kleiner, B. S. Koribalski, K. Lee-Waddell, C. Murugeshan, S. -H. Oh, J. Rhee, K. Spekkens, L. Staveley-Smith, A. R. H. Stevens, J. M. van der Hulst, J. Wang, O. I. Wong, B. W. Holwerda , et al. (3 additional authors not shown)

Abstract: We present results from our analysis of the Hydra I cluster observed in neutral atomic hydrogen (HI) as part of the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). These WALLABY observations cover a 60-square-degree field of view with uniform sensitivity and a spatial resolution of 30 arcsec. We use these wide-field observations to investigate the effect of galaxy environment on HI g… ▽ More We present results from our analysis of the Hydra I cluster observed in neutral atomic hydrogen (HI) as part of the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). These WALLABY observations cover a 60-square-degree field of view with uniform sensitivity and a spatial resolution of 30 arcsec. We use these wide-field observations to investigate the effect of galaxy environment on HI gas removal and star formation quenching by comparing the properties of cluster, infall and field galaxies extending up to $\sim5R_{200}$ from the cluster centre. We find a sharp decrease in the HI-detected fraction of infalling galaxies at a projected distance of $\sim1.5R_{200}$ from the cluster centre from $\sim0.85\%$ to $\sim0.35\%$. We see evidence for the environment removing gas from the outskirts of HI-detected cluster and infall galaxies through the decrease in the HI to $r$-band optical disc diameter ratio. These galaxies lie on the star forming main sequence, indicating that gas removal is not yet affecting the inner star-forming discs and is limited to the galaxy outskirts. Although we do not detect galaxies undergoing galaxy-wide quenching, we do observe a reduction in recent star formation in the outer disc of cluster galaxies, which is likely due to the smaller gas reservoirs present beyond the optical radius in these galaxies. Stacking of HI non-detections with HI masses below $M_{\rm{HI}}\lesssim10^{8.4}\,\rm{M}_{\odot}$ will be required to probe the HI of galaxies undergoing quenching at distances $\gtrsim60$ Mpc with WALLABY. △ Less

Submitted 30 November, 2021; originally announced December 2021.

Comments: 18 pages, 10 figures, accepted for publication in MNRAS

VERTICO: The Virgo Environment Traced In CO Survey

Authors: Toby Brown, Christine D. Wilson, Nikki Zabel, Timothy A. Davis, Alessandro Boselli, Aeree Chung, Sara L. Ellison, Claudia D. P. Lagos, Adam R. H. Stevens, Luca Cortese, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Claire R. Cashmore, Barbara Catinella, Ryan Chown, Benedikt Diemer, Pascal J. Elahi, Maan H. Hani, María J. Jiménez-Donaire, Bumhyun Lee, Katya Leidig, Angus Mok, Karen Pardos Olsen, Laura C. Parker , et al. (11 additional authors not shown)

Abstract: We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star forma… ▽ More We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, $^{12}$CO($2-1$) observations, and data reduction procedures. We characterize global $^{12}$CO($2-1$) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two non-detections, and produce resolved $^{12}$CO($2-1$) data products (median resolution $= 8^{\prime\prime} \approx 640~{\rm pc}$). Azimuthally averaged $^{12}$CO($2-1$) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size--mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size--mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy dataset for studying galaxy evolution in our closest cluster. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 68 pages, 13 Figures, 2 Figure Sets, Accepted for publication in ApJS, Online FITS versions of Tables 1, 2, and 3 are available with the journal publication

Drivers of asymmetry in synthetic H I emission-line profiles of galaxies in the EAGLE simulation

Authors: Aditya Manuwal, Aaron D. Ludlow, Adam R. H. Stevens, Ruby J. Wright, Aaron S. G. Robotham

Abstract: We study the shapes of spatially integrated H I emission-line profiles of galaxies in the EAGLE simulation using three separate measures of the profile's asymmetry. We show that the subset of EAGLE galaxies whose gas fractions and stellar masses are consistent with those in the xGASS survey also have similar H I line asymmetries. Central galaxies with symmetric H I line profiles typically correspo… ▽ More We study the shapes of spatially integrated H I emission-line profiles of galaxies in the EAGLE simulation using three separate measures of the profile's asymmetry. We show that the subset of EAGLE galaxies whose gas fractions and stellar masses are consistent with those in the xGASS survey also have similar H I line asymmetries. Central galaxies with symmetric H I line profiles typically correspond to rotationally supported H I and stellar disks, but those with asymmetric line profiles may or may not correspond to dispersion-dominated systems. Galaxies with symmetric H I emission lines are, on average, more gas rich than those with asymmetric lines, and also exhibit systematic differences in their specific star formation rates, suggesting that turbulence generated by stellar or AGN feedback may be one factor contributing to H I line asymmetry. The line asymmetry also correlates strongly with the dynamical state of a galaxy's host dark matter halo: older, more relaxed haloes host more-symmetric galaxies than those hosted by unrelaxed ones. At fixed halo mass, asymmetric centrals tend to be surrounded by a larger number of massive subhaloes than their symmetric counterparts, and also experience higher rates of gas accretion and outflow. At fixed stellar mass, central galaxies have, on average, more symmetric H I emission lines than satellites; for the latter, ram pressure and tidal stripping are significant sources of asymmetry. △ Less

Submitted 29 December, 2021; v1 submitted 23 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS

WALLABY Pre-Pilot Survey: The effects of angular momentum and environment on the HI gas and star formation properties of galaxies in the Eridanus supergroup

Authors: C. Murugeshan, V. A. Kilborn, B. -Q. For, O. I. Wong, J. Wang, T. Westmeier, A. R. H. Stevens, K. Spekkens, P. Kamphuis, L. Staveley-Smith, K. Lee-Waddell, D. Kleiner, B. S. Koribalski, M. E. Cluver, S. -H. Oh, J. Rhee, B. Catinella, T. N. Reynolds, H. Denes, A. Elagali

Abstract: We use high-resolution ASKAP observations of galaxies in the Eridanus supergroup to study their HI, angular momentum and star formation properties, as part of the WALLABY pre-pilot survey efforts. The Eridanus supergroup is composed of three sub-groups in the process of merging to form a cluster. The main focus of this study is the Eridanus (or NGC 1395) sub-group. The baryonic specific angular mo… ▽ More We use high-resolution ASKAP observations of galaxies in the Eridanus supergroup to study their HI, angular momentum and star formation properties, as part of the WALLABY pre-pilot survey efforts. The Eridanus supergroup is composed of three sub-groups in the process of merging to form a cluster. The main focus of this study is the Eridanus (or NGC 1395) sub-group. The baryonic specific angular momentum - baryonic mass ($j_{\mathrm{b}} - M_{\mathrm{b}}$) relation for the Eridanus galaxies is observed to be an unbroken power law of the form $j_{\mathrm{b}} \propto M_{\mathrm{b}}^{0.57 \pm 0.05}$, with a scatter of $\sim 0.10 \pm 0.01$ dex, consistent with previous works. We examine the relation between the atomic gas fraction, $f_{\mathrm{atm}}$, and the integrated atomic disc stability parameter $q$ (the $f_{\mathrm{atm}} - q$ relation), and find that the Eridanus galaxies deviate significantly from the relation owing to environmental processes such as tidal interactions and ram-pressure affecting their HI gas. We find that a majority of the Eridanus galaxies are HI deficient compared to normal star-forming galaxies in the field. We also find that the star formation among the Eridanus galaxies may be suppressed owing to their environment, thus hinting at significant levels of pre-processing within the Eridanus sub-group, even before the galaxies have entered a cluster-like environment. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: 19 pages, 19 figures, accepted for publication in MNRAS

WALLABY pre-pilot survey: Two dark clouds in the vicinity of NGC 1395

Authors: O. Ivy Wong, A. R. H. Stevens, B. -Q. For, T. Westmeier, M. Dixon, S. -H. Oh, G. I. G. Józsa, T. N. Reynolds, K. Lee-Waddell, J. Román, L. Verdes-Montenegro, H. M. Courtois, D. Pomarède, C. Murugeshan, M. T. Whiting, K. Bekki, F. Bigiel, A. Bosma, B. Catinella, H. Dénes, A. Elagali, B. W. Holwerda, P. Kamphuis, V. A. Kilborn, D. Kleiner , et al. (12 additional authors not shown)

Abstract: We present the Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pre-pilot observations of two `dark' HI sources (with HI masses of a few times 10^8 Msol and no known stellar counterpart) that reside within 363 kpc of NGC 1395, the most massive early-type galaxy in the Eridanus group of galaxies. We investigate whether these `dark' HI sources have resulted from past tidal interactions o… ▽ More We present the Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pre-pilot observations of two `dark' HI sources (with HI masses of a few times 10^8 Msol and no known stellar counterpart) that reside within 363 kpc of NGC 1395, the most massive early-type galaxy in the Eridanus group of galaxies. We investigate whether these `dark' HI sources have resulted from past tidal interactions or whether they are an extreme class of low surface brightness galaxies. Our results suggest that both scenarios are possible, and not mutually exclusive. The two `dark' HI sources are compact, reside in relative isolation and are more than 159 kpc away from their nearest HI-rich galaxy neighbour. Regardless of origin, the HI sizes and masses of both `dark' HI sources are consistent with the HI size-mass relationship that is found in nearby low-mass galaxies, supporting the possibility that these HI sources are an extreme class of low surface brightness galaxies. We identified three analogues of candidate primordial `dark' HI galaxies within the TNG100 cosmological, hydrodynamic simulation. All three model analogues are dark matter-dominated, have assembled most of their mass 12-13 Gyr ago, and have not experienced much evolution until cluster infall 1-2 Gyr ago. Our WALLABY pre-pilot science results suggest that the upcoming large area HI surveys will have a significant impact on our understanding of low surface brightness galaxies and the physical processes that shape them. △ Less

Submitted 9 August, 2021; originally announced August 2021.

Comments: 16 pages, 11 figures, accepted for publication in MNRAS

HI intensity mapping correlation function from UNIT simulations: BAO and observationally induced anisotropy

Authors: Santiago Avila, Bernhard Vos-Ginés, Steven Cunnington, Adam R. H. Stevens, Gustavo Yepes, Alexander Knebe, Chia-Hsun Chuang

Abstract: We study the clustering of HI intensity maps produced from simulations with a focus on baryonic acoustic oscillations (BAO) and the effects induced by telescope beam smoothing and foreground cleaning. We start by creating a HI catalogue at $z=1.321$ based on the Semi-Analytic Galaxy Evolution (SAGE) model applied to the UNIT simulations. With this catalogue we investigate the relation between mode… ▽ More We study the clustering of HI intensity maps produced from simulations with a focus on baryonic acoustic oscillations (BAO) and the effects induced by telescope beam smoothing and foreground cleaning. We start by creating a HI catalogue at $z=1.321$ based on the Semi-Analytic Galaxy Evolution (SAGE) model applied to the UNIT simulations. With this catalogue we investigate the relation between model HI and the dark matter haloes and we also study the abundance of HI, $Ω_{\rm HI}$, predicted by this model. We then create synthetic HI intensity maps with a Nearest-Grid-Point approach. In order to simulate the telescope beam effect, a Gaussian smoothing is applied on the plane perpendicular to the line of sight. The effect of foreground removal methods is simulated by exponentially damping the largest wavelength Fourier modes on the radial direction. We study the anisotropic 2-point correlation function (2PCF) $ξ(r_\perp,r_\parallel)$ and how it is affected by the aforementioned observational effects. In order to better isolate the BAO signal, we study several 2PCF $μ$-wedges (with a restricted range of orientations $μ$) tailored to address the systematics effects and we compare them with different definitions of radial 2PCFs. Finally, we discuss our findings in the context of an SKA-like survey, finding a clear BAO signal in most of the estimators here proposed. △ Less

Submitted 10 December, 2021; v1 submitted 21 May, 2021; originally announced May 2021.

Comments: published in MNRAS; 18 pages, 10 figures

An Empirical Background Model for the NICER X-ray Timing Instrument

Authors: Ronald A. Remillard, Michael Loewenstein, James F. Steiner, Gregory Y. Prigozhin, Beverly LaMarr, Teruaki Enoto, Keith C. Gendreau, Zaven Arzoumanian, Craig Markwardt, Arkadip Basak, Abigail L. Stevens, Paul S. Ray, Diego Altamirano, Douglas J. K. Buisson

Abstract: NICER has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4-12 keV;… ▽ More NICER has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4-12 keV; 50 detectors) of 0.87 c/s, but in 5 percent (1 percent) of cases, the rate exceeds 10 (300) c/s. Model residuals persist at 20-30 percent of the initial rate for the brightest GTIs, implying one or more missing model parameters. Filtering criteria are given to flag GTIs likely to have unsatisfactory background predictions. With such filtering, we estimate a detection limit, 1.20 c/s (3 sigma, single GTI) at 0.4-12 keV, equivalent to 3.6e-12 erg/cm^2/s for a Crab-like spectrum. The corresponding limit for soft X-ray sources is 0.51 c/s at 0.3-2.0 keV, or 4.3e-13 erg/cm^2/s for a 100 eV blackbody. Faint-source filtering selects 85 percent of the background GTIs, and higher rates are expected for targets scheduled more favorably. An application of the model to 1 s timescale makes it possible to distinguish source flares from possible surges in the background. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 27 pages, 14 figures; submitted to the Astronomical Journal

WALLABY Pilot Survey: First Look at the Hydra I Cluster and Ram Pressure Stripping of ESO 501-G075

Authors: T. N. Reynolds, T. Westmeier, A. Elagali, B. Catinella, L. Cortese, N. Deg, B. -Q. For, P. Kamphuis, D. Kleiner, B. S. Koribalski, K. Lee-Waddell, S. -H. Oh, J. Rhee, P. Serra, K. Spekkens, L. Staveley-Smith, A. R. H. Stevens, E. N. Taylor, J. Wang, O. I. Wong

Abstract: We present results from neutral atomic hydrogen (HI) observations of Hydra I, the first cluster observed by the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY) on the Australian Square Kilometre Array Pathfinder. For the first time we show that WALLABY can reach its final survey sensitivity. Leveraging the sensitivity, spatial resolution and wide field of view of WALLABY, we identify… ▽ More We present results from neutral atomic hydrogen (HI) observations of Hydra I, the first cluster observed by the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY) on the Australian Square Kilometre Array Pathfinder. For the first time we show that WALLABY can reach its final survey sensitivity. Leveraging the sensitivity, spatial resolution and wide field of view of WALLABY, we identify a galaxy, ESO 501-G075, that lies near the virial radius of Hydra I and displays an HI tail. ESO 501-G075 shows a similar level of morphological asymmetry as another cluster member, which lies near the cluster centre and shows signs of experiencing ram pressure. We investigate possible environmental processes that could be responsible for producing the observed disturbance in the HI morphology of ESO 501-G075. We rule out tidal interactions, as ESO 501-G075 has no nearby neighbours within $\sim0.34$Mpc. We use a simple model to determine that ram pressure can remove gas from the disc at radii $r\gtrsim25$kpc. We conclude that, as ESO 501-G075 has a typical HI mass compared to similar galaxies in the field and its morphology is compatible with a ram pressure scenario, ESO 501-G075 is likely recently infalling into the cluster and in the early stages of experiencing ram pressure. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

UNITSIM-Galaxies: data release and clustering of emission-line galaxies

Authors: Alexander Knebe, Daniel Lopez-Cano, Santiago Avila, Ginevra Favole, Adam R. H. Stevens, Violeta Gonzalez-Perez, Guillermo Reyes-Peraza, Gustavo Yepes, Chia-Hsun Chuang, Francisco-Shu Kitaura

Abstract: New surveys such as ESA's Euclid mission are planned to map with unprecedented precision the large-scale structure of the Universe by measuring the 3D positions of tens of millions of galaxies. It is necessary to develop theoretically modelled galaxy catalogues to estimate the expected performance and to optimise the analysis strategy of these surveys. We populate two pairs of (1 Gpc/h)^3 volume d… ▽ More New surveys such as ESA's Euclid mission are planned to map with unprecedented precision the large-scale structure of the Universe by measuring the 3D positions of tens of millions of galaxies. It is necessary to develop theoretically modelled galaxy catalogues to estimate the expected performance and to optimise the analysis strategy of these surveys. We populate two pairs of (1 Gpc/h)^3 volume dark-matter-only simulations from the UNIT project with galaxies using the SAGE semi-analytic model of galaxy formation, coupled to the photoionisation model get_emlines to estimate their H_alpha emission. These catalogues represent a unique suite that includes galaxy formation physics and - thanks to the fixed-pair technique used - an effective volume of ~(5 Gpc/h)^3, which is several times larger than the Euclid survey. We present the performance of these data and create five additional emission-line galaxy (ELG) catalogues by applying a dust attenuation model as well as adjusting the flux threshold as a function of redshift in order to reproduce Euclid-forecast dN/dz values. As a first application, we study the abundance and clustering of those model H_alpha ELGs: for scales greater than ~5 Mpc/h, we find a scale-independent bias with a value of $b\sim 1$ at redshift z ~ 0.5, that can increase nearly linearly to $b\sim 4$ at $z\sim 2$, depending on the ELG catalogue. Model galaxy properties, including their emission-line fluxes (with and without dust extinction) are publicly available. △ Less

Submitted 12 January, 2022; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 16 pages, revised version as accepted by MNRAS

Unveiling the atomic hydrogen-halo mass relation via spectral stacking

Authors: Garima Chauhan, Claudia del P. Lagos, Adam R. H. Stevens, Matias Bravo, Jonghwan Rhee, Chris Power, Danail Obreschkow, Martin Meyer

Abstract: Measuring the HI-halo mass scaling relation (HIHM) is fundamental to understanding the role of HI in galaxy formation and its connection to structure formation. While direct measurements of the HI mass in haloes are possible using HI-spectral stacking, the reported shape of the relation depends on the techniques used to measure it (e.g. monotonically increasing with mass versus flat, mass-independ… ▽ More Measuring the HI-halo mass scaling relation (HIHM) is fundamental to understanding the role of HI in galaxy formation and its connection to structure formation. While direct measurements of the HI mass in haloes are possible using HI-spectral stacking, the reported shape of the relation depends on the techniques used to measure it (e.g. monotonically increasing with mass versus flat, mass-independent). Using a simulated HI and optical survey produced with the SHARK semi-analytic galaxy formation model, we investigate how well different observational techniques can recover the intrinsic, theoretically predicted, HIHM relation. We run a galaxy group finder and mimic the HI stacking procedure adopted by different surveys and find we can reproduce their observationally derived HIHM relation. However, none of the adopted techniques recover the underlying HIHM relation predicted by the simulation. We find that systematic effects in halo mass estimates of galaxy groups modify the inferred shape of the HIHM relation from the intrinsic one in the simulation, while contamination by interloping galaxies, not associated with the groups, contribute to the inferred HI mass of a halo mass bin, when using large velocity windows for stacking. The effect of contamination is maximal at Mvir~10^(12-12.5)Msol. Stacking methods based on summing the HI emission spectra to infer the mean HI mass of galaxies of different properties belonging to a group suffer minimal contamination but are strongly limited by the use of optical counterparts, which miss the contribution of dwarf galaxies. Deep spectroscopic surveys will provide significant improvements by going deeper while maintaining high spectroscopic completeness; for example, the WAVES survey will recover ~52% of the total HI mass of the groups with Mvir~10^(14)Msol compared to ~21% in GAMA. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: 20 pages,16 figures. Paper submitted to MNRAS

Projected sensitivities of the LUX-ZEPLIN (LZ) experiment to new physics via low-energy electron recoils

Authors: The LZ Collaboration, D. S. Akerib, A. K. Al Musalhi, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, P. Beltrame, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch , et al. (172 additional authors not shown)

Abstract: LUX-ZEPLIN (LZ) is a dark matter detector expected to obtain world-leading sensitivity to weakly interacting massive particles (WIMPs) interacting via nuclear recoils with a ~7-tonne xenon target mass. This manuscript presents sensitivity projections to several low-energy signals of the complementary electron recoil signal type: 1) an effective neutrino magnetic moment and 2) an effective neutrino… ▽ More LUX-ZEPLIN (LZ) is a dark matter detector expected to obtain world-leading sensitivity to weakly interacting massive particles (WIMPs) interacting via nuclear recoils with a ~7-tonne xenon target mass. This manuscript presents sensitivity projections to several low-energy signals of the complementary electron recoil signal type: 1) an effective neutrino magnetic moment and 2) an effective neutrino millicharge, both for pp-chain solar neutrinos, 3) an axion flux generated by the Sun, 4) axion-like particles forming the galactic dark matter, 5) hidden photons, 6) mirror dark matter, and 7) leptophilic dark matter. World-leading sensitivities are expected in each case, a result of the large 5.6t 1000d exposure and low expected rate of electron recoil backgrounds in the $<$100keV energy regime. A consistent signal generation, background model and profile-likelihood analysis framework is used throughout. △ Less

Submitted 18 May, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: v2 updates exclusion sensitivities from single-sided to two-sided

Enhancing the sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment to low energy signals

Authors: D. S. Akerib, A. K. Al Musalhi, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, P. Beltrame, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, G. M. Blockinger , et al. (162 additional authors not shown)

Abstract: Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matt… ▽ More Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 14 pages, 6 figures

Molecular hydrogen in IllustrisTNG galaxies: carefully comparing signatures of environment with local CO & SFR data

Authors: Adam R. H. Stevens, Claudia del P. Lagos, Luca Cortese, Barbara Catinella, Benedikt Diemer, Dylan Nelson, Annalisa Pillepich, Lars Hernquist, Federico Marinacci, Mark Vogelsberger

Abstract: We examine how the post-processed content of molecular hydrogen (H$_2$) in galaxies from the TNG100 cosmological, hydrodynamic simulation changes with environment at $z\!=\!0$, assessing central/satellite status and host halo mass. We make close comparisons with the carbon monoxide (CO) emission survey xCOLD GASS where possible, having mock-observed TNG100 galaxies to match the survey's specificat… ▽ More We examine how the post-processed content of molecular hydrogen (H$_2$) in galaxies from the TNG100 cosmological, hydrodynamic simulation changes with environment at $z\!=\!0$, assessing central/satellite status and host halo mass. We make close comparisons with the carbon monoxide (CO) emission survey xCOLD GASS where possible, having mock-observed TNG100 galaxies to match the survey's specifications. For a representative sample of host haloes across $10^{11}\!\lesssim\!M_{\rm 200c}/{\rm M}_{\odot}\!<\!10^{14.6}$, TNG100 predicts that satellites with $m_*\!\geq\!10^9\,{\rm M}_{\odot}$ should have a median deficit in their H$_2$ fractions of $\sim$0.6 dex relative to centrals of the same stellar mass. Once observational and group-finding uncertainties are accounted for, the signature of this deficit decreases to $\sim$0.2 dex. Remarkably, we calculate a deficit in xCOLD GASS satellites' H$_2$ content relative to centrals of 0.2--0.3 dex, in line with our prediction. We further show that TNG100 and SDSS data exhibit continuous declines in the average star formation rates of galaxies at fixed stellar mass in denser environments, in quantitative agreement with each other. By tracking satellites from their moment of infall in TNG100, we directly show that atomic hydrogen (HI) is depleted at fractionally higher rates than H$_2$ on average. Supporting this picture, we find that the H$_2$/HI mass ratios of satellites are elevated relative to centrals in xCOLD GASS. We provide additional predictions for the effect of environment on H$_2$ -- both absolute and relative to HI -- that can be tested with spectral stacking in future CO surveys. △ Less

Submitted 17 December, 2020; v1 submitted 6 November, 2020; originally announced November 2020.

Comments: Accepted in MNRAS. Minor proofing and reference edits from v1. Main body: 19 pages, 11 figures

Global HI asymmetries in IllustrisTNG: a diversity of physical processes disturb the cold gas in galaxies

Authors: Adam B. Watts, Chris Power, Barbara Catinella, Luca Cortese, Adam R. H. Stevens

Abstract: Observations of the cold neutral atomic hydrogen (HI) in and around disc galaxies have revealed that spatial and kinematic asymmetries are commonplace, and are reflected in the global HI spectra. We use the TNG100 box from the IllustrisTNG suite of cosmological simulations to study the conditions under which these asymmetries may arise in current theoretical galaxy formation models. We find that m… ▽ More Observations of the cold neutral atomic hydrogen (HI) in and around disc galaxies have revealed that spatial and kinematic asymmetries are commonplace, and are reflected in the global HI spectra. We use the TNG100 box from the IllustrisTNG suite of cosmological simulations to study the conditions under which these asymmetries may arise in current theoretical galaxy formation models. We find that more than 50% of the sample has at least a 10% difference in integrated flux between the high- and low-velocity half of the spectrum, thus the typical TNG100 galaxy has an HI profile that is not fully symmetric. We find that satellite galaxies are a more asymmetric population than centrals, consistent with observational results. Using halo mass as a proxy for environment, this trend appears to be driven by the satellite population within the virial radius of haloes more massive than $10^{13} M_{\odot}$, typical of medium/large groups. We show that, while the excess of HI asymmetry in group satellites is likely driven by ram pressure, the bulk of the asymmetric HI profiles observed in TNG100 are driven by physical processes able to affect both the central and satellite populations. Our results highlight how asymmetries are not driven solely by environment, and multiple physical processes can produce the same asymmetric shape in global HI spectra. △ Less

Submitted 11 October, 2020; originally announced October 2020.

Comments: 16 pages, 13 figures, 1 appendix, accepted for publication in MRNRAS

The distribution and properties of DLAs at z $\leq$ 2 in the EAGLE simulations

Authors: Lilian Garratt-Smithson, Chris Power, Claudia del P. Lagos, Adam R. H. Stevens, James R. Allison, Elaine M. Sadler

Abstract: Determining the spatial distribution and intrinsic physical properties of neutral hydrogen on cosmological scales is one of the key goals of next-generation radio surveys. We use the EAGLE galaxy formation simulations to assess the properties of damped Lyman-alpha absorbers (DLAs) that are associated with galaxies and their underlying dark matter haloes between 0 $\leq$ z $\leq$ 2. We find that th… ▽ More Determining the spatial distribution and intrinsic physical properties of neutral hydrogen on cosmological scales is one of the key goals of next-generation radio surveys. We use the EAGLE galaxy formation simulations to assess the properties of damped Lyman-alpha absorbers (DLAs) that are associated with galaxies and their underlying dark matter haloes between 0 $\leq$ z $\leq$ 2. We find that the covering fraction of DLAs increases at higher redshift; a significant fraction of neutral atomic hydrogen (HI) resides in the outskirts of galaxies with stellar mass greater than or equal to 10$^{10}$ M$_\odot$; and the covering fraction of DLAs in the circumgalactic medium (CGM) is enhanced relative to that of the interstellar medium (ISM) with increasing halo mass. Moreover, we find that the mean density of the HI in galaxies increases with increasing stellar mass, while the DLAs in high- and low-halo-mass systems have higher column densities than those in galaxies with intermediate halo masses (~ 10$^{12}$ M$_\odot$ at z = 0). These high-impact CGM DLAs in high-stellar-mass systems tend to be metal-poor, likely tracing smooth accretion. Overall, our results point to the CGM playing an important role in DLA studies at high redshift (z $\geq$ 1). However, their properties are impacted both by numerical resolution and the detailed feedback prescriptions employed in cosmological simulations, particularly that of AGN. △ Less

Submitted 14 December, 2020; v1 submitted 21 August, 2020; originally announced August 2020.

Comments: 25 pages. Accepted for publication in MNRAS

The physical drivers of the atomic hydrogen-halo mass relation

Authors: Garima Chauhan, Claudia del P. Lagos, Adam R. H. Stevens, Danail Obreschkow, Chris Power, Martin Meyer

Abstract: We use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0\leq z \leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at… ▽ More We use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0\leq z \leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at $M_{\rm vir}\approx 10^{12-13} M_{\odot}$ in the latter, but otherwise great agreement below and above that mass threshold. In SHARK, we find that the HI mass increases with the halo mass up to a critical mass of $\approx 10^{11.8} M_{\odot}$; between $\sim 10^{11.8}-10^{13}M_{\odot}$, the scatter in the relation increases by 0.7 dex and the HI mass decreases with the halo mass on average; at $M_{\rm vir} \geq 10^{13} M_{\odot}$, the HI content continues to increase with halo mass. We find that the critical halo mass of $\approx 10^{12} M_{\odot}$ is largely set by feedback from Active Galactic Nuclei (AGN), and the exact shape and scatter of the HI-halo mass relation around that mass is extremely sensitive to how AGN feedback is modelled, with other physical processes playing a less significant role. We determine the main secondary parameters responsible for the scatter of the HI-halo mass relation, namely the halo spin parameter at $M_{\rm vir}\leq 10^{11.8} M_{\odot}$, and the fractional contribution from substructure to the total halo mass for $M_{\rm vir}\geq 10^{13} M_{\odot}$. The scatter at $10^{11.8}<M_{\rm vir}<10^{13} M_{\odot}$ is best described by the black-hole-to-stellar mass ratio of the central galaxy, reflecting the AGN feedback relevance. We present a numerical model to populate dark matter-only simulations with HI at $0\leq z \leq 2$ based solely on halo parameters that are measurable in such simulations. △ Less

Submitted 27 July, 2020; v1 submitted 22 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS

A NICER look at the state transitions of the black hole candidate MAXI J1535-571 during its reflares

Authors: V. A. Cúneo, K. Alabarta, L. Zhang, D. Altamirano, M. Méndez, M. Armas Padilla, R. Remillard, J. Homan, J. F. Steiner, J. A. Combi, T. Muñoz-Darias, K. C. Gendreau, Z. Arzoumanian, A. L. Stevens, M. Loewenstein, F. Tombesi, P. Bult, A. C. Fabian, D. J. K. Buisson, J. Neilsen, A. Basak

Abstract: The black hole candidate and X-ray binary MAXI J1535-571 was discovered in September 2017. During the decay of its discovery outburst, and before returning to quiescence, the source underwent at least four reflaring events, with peak luminosities of $\sim$10$^{35-36}$ erg s$^{-1}$ (d/4.1 kpc)$^2$. To investigate the nature of these flares, we analysed a sample of NICER observations taken with almo… ▽ More The black hole candidate and X-ray binary MAXI J1535-571 was discovered in September 2017. During the decay of its discovery outburst, and before returning to quiescence, the source underwent at least four reflaring events, with peak luminosities of $\sim$10$^{35-36}$ erg s$^{-1}$ (d/4.1 kpc)$^2$. To investigate the nature of these flares, we analysed a sample of NICER observations taken with almost daily cadence. In this work we present the detailed spectral and timing analysis of the evolution of the four reflares. The higher sensitivity of NICER at lower energies, in comparison with other X-ray detectors, allowed us to constrain the disc component of the spectrum at $\sim$0.5 keV. We found that during each reflare the source appears to trace out a q-shaped track in the hardness-intensity diagram similar to those observed in black hole binaries during full outbursts. MAXI J1535-571 transits between the hard state (valleys) and softer states (peaks) during these flares. Moreover, the Comptonised component is undetected at the peak of the first reflare, while the disc component is undetected during the valleys. Assuming the most likely distance of 4.1 kpc, we find that the hard-to-soft transitions take place at the lowest luminosities ever observed in a black hole transient, while the soft-to-hard transitions occur at some of the lowest luminosities ever reported for such systems. △ Less

Submitted 6 July, 2020; v1 submitted 4 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS main journal

Galaxy And Mass Assembly (GAMA): A forensic SED reconstruction of the cosmic star formation history and metallicity evolution by galaxy type

Authors: Sabine Bellstedt, Aaron S. G. Robotham, Simon P. Driver, Jessica E. Thorne, Luke J. M. Davies, Claudia del P. Lagos, Adam R. H. Stevens, Edward N. Taylor, Ivan K. Baldry, Amanda J. Moffett, Andrew M. Hopkins, Steven Phillipps

Abstract: We apply the spectral energy distribution (SED) fitting code ProSpect to multiwavelength imaging for $\sim$7,000 galaxies from the GAMA survey at $z<0.06$, in order to extract their star formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We s… ▽ More We apply the spectral energy distribution (SED) fitting code ProSpect to multiwavelength imaging for $\sim$7,000 galaxies from the GAMA survey at $z<0.06$, in order to extract their star formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We show with this approach that we are able to recover the observationally determined cosmic star formation history (CSFH), an indication that stars are being formed in the correct epoch of the Universe, on average, for the manner in which we are conducting SED fitting. We also show the contribution to the CSFH of galaxies of different present-day visual morphologies, and stellar masses. Our analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago. In other morphological types, the stellar mass formed later, up to 6 Gyr ago for present-day irregular galaxies. Similarly, the most massive galaxies in our sample were shown to have formed half their stellar mass by 11 Gyr ago, whereas the least massive galaxies reached this stage as late as 4 Gyr ago (the well-known effect of "galaxy downsizing"). Finally, our metallicity approach allows us to follow the average evolution in gas-phase metallicity for populations of galaxies, and extract the evolution of the cosmic metal mass density in stars and in gas, producing results in broad agreement with independent, higher-redshift observations of metal densities in the Universe. △ Less

Submitted 7 October, 2020; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: 24 pages, 16 figures, accepted for publication by MNRAS

Journal ref: MNRAS, 2020, 498, 5581

The 2019 super-Eddington outburst of RX J0209.6-7427: Detection of pulsations and constraints on the magnetic field strength

Authors: G. Vasilopoulos, P. S. Ray, K. C. Gendreau, P. A. Jenke, G. K. Jaisawal, C. A. Wilson-Hodge, T. E. Strohmayer, D. Altamirano, W. B. Iwakiri, M. T. Wolff, S. Guillot, C. Malacaria, A. L. Stevens

Abstract: In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provid… ▽ More In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 10 pages, 8 figures, accepted for publication in MNRAS

A rapid change in X-ray variability and a jet ejection in the black hole transient MAXI J1820+070

Authors: Jeroen Homan, Joe Bright, Sara E. Motta, Diego Altamirano, Zaven Arzoumanian, Arkadip Basak, Tomaso M. Belloni, Edward M. Cackett, Rob Fender, Keith C. Gendreau, Erin Kara, Dheeraj R. Pasham, Ronald A. Remillard, James F. Steiner, Abigail L. Stevens, Phil Uttley

Abstract: We present Neutron Star Interior Composition Explorer X-ray and Arcminute Microkelvin Imager Large Array radio observations of a rapid hard-to-soft state transition in the black hole X-ray transient MAXI J1820+070. During the transition from the hard state to the soft state a switch between two particular types of quasiperiodic oscillations (QPOs) was seen in the X-ray power density spectra, from… ▽ More We present Neutron Star Interior Composition Explorer X-ray and Arcminute Microkelvin Imager Large Array radio observations of a rapid hard-to-soft state transition in the black hole X-ray transient MAXI J1820+070. During the transition from the hard state to the soft state a switch between two particular types of quasiperiodic oscillations (QPOs) was seen in the X-ray power density spectra, from type-C to type-B, along with a drop in the strength of the broadband X-ray variability and a brief flare in the 7-12 keV band. Soon after this switch (~1.5-2.5 hr) a strong radio flare was observed that corresponded to the launch of superluminal ejecta. Although hints of a connection between QPO transitions and radio flares have been seen in other black hole X-ray transients, our observations constitute the strongest observational evidence to date for a link between the appearance of type-B QPOs and the launch of discrete jet ejections. △ Less

Submitted 2 March, 2020; originally announced March 2020.

Comments: Accepted for publication in ApJ Letters

WALLABY -- An SKA Pathfinder HI Survey

Authors: B. S. Koribalski, L. Staveley-Smith, T. Westmeier, P. Serra, K. Spekkens, O. I. Wong, C. D. P. Lagos, D. Obreschkow, E. V. Ryan-Weber, M. Zwaan, V. Kilborn, G. Bekiaris, K. Bekki, F. Bigiel, A. Boselli, A. Bosma, B. Catinella, G. Chauhan, M. E. Cluver, M. Colless, H. M. Courtois, R. A. Crain, W. J. G. de Blok, H. Dénes, A. R. Duffy , et al. (45 additional authors not shown)

Abstract: The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western A… ▽ More The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA). △ Less

Submitted 7 July, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApSS (38 pages, 14 figures), see also https://www.atnf.csiro.au/research/WALLABY/ - Contact email: Baerbel.Koribalski@csiro.au

Journal ref: Ap&SS 365, 118 (2020) - https://rdcu.be/b5Bfg

Galaxy Cold Gas Contents in Modern Cosmological Hydrodynamic Simulations

Authors: Romeel Davé, Robert A. Crain, Adam R. H. Stevens, Desika Narayanan, Amelie Saintonge, Barbara Catinella, Luca Cortese

Abstract: We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ m… ▽ More We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ mass functions (HIMF, H2MF), CO(1-0) luminosity functions (COLF), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density $μ_*$, with some quantitative differences. To compare to the COLF, we apply an H$_2$-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in $α_{\rm CO}$ and significant differences between models. Using this, predicted $z=0$ COLFs agree better with data than predicted H2MFs. Out to $z\sim 2$, EAGLE's and Simba's HIMF and COLF strongly increase, while TNG's HIMF declines and COLF evolves slowly. EAGLE and Simba reproduce high $L_{\rm CO1-0}$ galaxies at $z\sim 1-2$ as observed, owing partly to a median $α_{\rm CO}(z=2)\sim 1$ versus $α_{\rm CO}(z=0)\sim 3$. Examining \HI, H$_2$, and CO scaling relations, their trends with $M_*$ are broadly reproduced in all models, but EAGLE yields too little HI in green valley galaxies, TNG and Simba overproduce cold gas in massive galaxies, and Simba overproduces molecular gas in small systems. Using Simba variants that exclude individual AGN feedback modules, we find that Simba's AGN jet feedback is primarily responsible by lowering cold gas contents from $z\sim 1\to0$ by suppressing cold gas in $M_*> 10^{10}{\rm M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-$μ_*$ systems. △ Less

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

Comments: 22 pages, MNRAS accepted

The imperative to reduce carbon emissions in astronomy

Authors: Adam R. H. Stevens, Sabine Bellstedt, Pascal J. Elahi, Michael T. Murphy

Abstract: For astronomers to make a significant contribution to the reduction of climate change-inducing greenhouse gas emissions, we first must quantify our sources of emissions and review the most effective approaches for reducing them. Here we estimate that Australian astronomers' total greenhouse gas emissions from their regular work activities are $\gtrsim$25 ktCO$_2$-e/yr (equivalent kilotonnes of car… ▽ More For astronomers to make a significant contribution to the reduction of climate change-inducing greenhouse gas emissions, we first must quantify our sources of emissions and review the most effective approaches for reducing them. Here we estimate that Australian astronomers' total greenhouse gas emissions from their regular work activities are $\gtrsim$25 ktCO$_2$-e/yr (equivalent kilotonnes of carbon dioxide per year). This can be broken into $\sim$15 ktCO$_2$-e/yr from supercomputer usage, $\sim$4.2 ktCO$_2$-e/yr from flights (where individuals' flight emissions correlate with seniority), $>$3.3 ktCO$_2$-e/yr from the operation of observatories, and 2.6$\pm$0.4 ktCO$_2$-e/yr from powering office buildings. Split across faculty scientists, postdoctoral researchers, and PhD students, this averages to $\gtrsim$37 tCO$_2$-e/yr per astronomer, over 40% more than what the average Australian non-dependant emits in total, equivalent to $\sim$5$\times$ the global average. To combat these environmentally unsustainable practices, we suggest astronomers should strongly preference use of supercomputers, observatories, and office spaces that are predominantly powered by renewable energy sources. Where facilities that we currently use do not meet this requirement, their funders should be lobbied to invest in renewables, such as solar or wind farms. Air travel should also be reduced wherever possible, replaced primarily by video conferencing, which should also promote inclusivity. △ Less

Submitted 18 June, 2020; v1 submitted 12 December, 2019; originally announced December 2019.

Comments: Accepted Perspective in Nature Astronomy. Significant additions/updates to results from v1. Text refinement post-refereeing since v2. 9 pages (excl. refs), 2 figures

Angular momentum and morphological sequence of massive galaxies through Dark Sage

Authors: Antonio J. Porras-Valverde, Kelly Holley-Bockelmann, Andreas A. Berlind, Adam R. H. Stevens

Abstract: We study the present-day connection between galaxy morphology and angular momentum using the {\sc Dark Sage} semi-analytic model of galaxy formation. For galaxies between $ 10^{11}-10^{12} \mathrm{M}_{\odot}$ in stellar mass, the model successfully predicts the observed trend whereby galaxies with more prominent disks exhibit higher {\em stellar} disk specific angular momentum (… ▽ More We study the present-day connection between galaxy morphology and angular momentum using the {\sc Dark Sage} semi-analytic model of galaxy formation. For galaxies between $ 10^{11}-10^{12} \mathrm{M}_{\odot}$ in stellar mass, the model successfully predicts the observed trend whereby galaxies with more prominent disks exhibit higher {\em stellar} disk specific angular momentum ($j_{\rm stellar, disk}$) at fixed stellar mass. However, when we include the gas in the disk, bulge-dominated galaxies have the highest {\em total} disk specific angular momentum ($j_{\rm total, disk}$). We attribute this to a large contribution from an extended disk of cold gas in typical bulge-dominated galaxies. We find the relationship between $j_{\rm dark matter}$ and morphology to be quite complex. Surprisingly, in this stellar mass range, not only do bulge-dominated galaxies tend to live in halos with higher $j_{\rm dark matter}$ than disk-dominated galaxies, but intermediate galaxies (those with roughly equal fractions of bulge and disk mass) have the lowest $j_{\rm dark matter}$ of all. Yet, when controlling for halo mass, rather than stellar mass, the relationship between $j_{\rm dark matter}$ and morphology vanishes. Based on these results, halo mass rather than angular momentum is the main driver of the predicted morphology sequence at high masses. In fact, in our stellar mass range, disk-dominated galaxies live in dark matter halos that are roughly 1/10th the mass of their bulge-dominated counterparts. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: 11 pages, 9 figures, submitted to ApJ