Showing 1–50 of 459 results for author: Bower, G

AtLAST Science Overview Report

Authors: Mark Booth, Pamela Klaassen, Claudia Cicone, Tony Mroczkowski, Martin A. Cordiner, Luca Di Mascolo, Doug Johnstone, Eelco van Kampen, Minju M. Lee, Daizhong Liu, John Orlowski-Scherer, Amélie Saintonge, Matthew W. L. Smith, Alexander Thelen, Sven Wedemeyer, Kazunori Akiyama, Stefano Andreon, Doris Arzoumanian, Tom J. L. C. Bakx, Caroline Bot, Geoffrey Bower, Roman Brajša, Chian-Chou Chen, Elisabete da Cunha, David Eden , et al. (59 additional authors not shown)

Abstract: Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still… ▽ More Submillimeter and millimeter wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open questions that cannot be answered with current facilities. In this report we summarise the science that is guiding the design of the Atacama Large Aperture Submillimeter Telescope (AtLAST). We demonstrate how tranformational advances in topics including star formation in high redshift galaxies, the diffuse circumgalactic medium, Galactic ecology, cometary compositions and solar flares motivate the need for a 50m, single-dish telescope with a 1-2 degree field of view and a new generation of highly multiplexed continuum and spectral cameras. AtLAST will have the resolution to drastically lower the confusion limit compared to current single-dish facilities, whilst also being able to rapidly map large areas of the sky and detect extended, diffuse structures. Its high sensitivity and large field of view will open up the field of submillimeter transient science by increasing the probability of serendipitous detections. Finally, the science cases listed here motivate the need for a highly flexible operations model capable of short observations of individual targets, large surveys, monitoring programmes, target of opportunity observations and coordinated observations with other observatories. AtLAST aims to be a sustainable, upgradeable, multipurpose facility that will deliver orders of magnitude increases in sensitivity and mapping speeds over current and planned submillimeter observatories. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 47 pages, 12 figures. For further details on AtLAST see https://atlast.uio.no

The key science drivers for the Atacama Large Aperture Submillimeter Telescope (AtLAST)

Authors: Mark Booth, Pamela Klaassen, Claudia Cicone, Tony Mroczkowski, Sven Wedemeyer, Kazunori Akiyama, Geoffrey Bower, Martin A. Cordiner, Luca Di Mascolo, Doug Johnstone, Eelco van Kampen, Minju M. Lee, Daizhong Liu, John Orlowski-Scherer, Amélie Saintonge, Matthew Smith, Alexander E. Thelen

Abstract: Sub-mm and mm wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open ques… ▽ More Sub-mm and mm wavelengths provide a unique view of the Universe, from the gas and dust that fills and surrounds galaxies to the chromosphere of our own Sun. Current single-dish facilities have presented a tantalising view of the brightest (sub-)mm sources, and interferometers have provided the exquisite resolution necessary to analyse the details in small fields, but there are still many open questions that cannot be answered with current facilities: Where are all the baryons? How do structures interact with their environments? What does the time-varying (sub-)mm sky look like? In order to make major advances on these questions and others, what is needed now is a facility capable of rapidly mapping the sky spatially, spectrally, and temporally, which can only be done by a high throughput, single-dish observatory. An extensive design study for this new facility is currently being undertaken. In this paper, we focus on the key science drivers and the requirements they place on the observatory. As a 50m single dish telescope with a 1-2° field of view, the strength of the Atacama Large Aperture Submillimeter Telescope (AtLAST) is in science where a large field of view, highly multiplexed instrumentation and sensitivity to faint large-scale structure is important. AtLAST aims to be a sustainable, upgradeable, multipurpose facility that will deliver orders of magnitude increases in sensitivity and mapping speeds over current and planned telescopes. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 12 pages, Conference proceedings paper for the 2024 SPIE Astronomical Telescopes + Instrumentation meeting

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

An unidentified Fermi source emitting radio bursts in the Galactic bulge

Authors: Reshma Anna-Thomas, Sarah Burke-Spolaor, Casey J. Law, F. K. Schinzel, Kshitij Aggarwal, Geoffrey C. Bower, Liam Connor, Paul B. Demorest

Abstract: We report on the detection of radio bursts from the Galactic bulge using the real-time transient detection and localization system, realfast. The pulses were detected commensally on the Karl G. Jansky Very Large Array during a survey of unidentified Fermi $γ$-ray sources. The bursts were localized to subarcsecond precision using realfast fast-sampled imaging. Follow-up observations with the Green… ▽ More We report on the detection of radio bursts from the Galactic bulge using the real-time transient detection and localization system, realfast. The pulses were detected commensally on the Karl G. Jansky Very Large Array during a survey of unidentified Fermi $γ$-ray sources. The bursts were localized to subarcsecond precision using realfast fast-sampled imaging. Follow-up observations with the Green Bank Telescope detected additional bursts from the same source. The bursts do not exhibit periodicity in a search up to periods of 480s, assuming a duty cycle of < 20%. The pulses are nearly 100% linearly polarized, show circular polarization up to 12%, have a steep radio spectral index of -2.7, and exhibit variable scattering on timescales of months. The arcsecond-level realfast localization links the source confidently with the Fermi $γ$-ray source and places it nearby (though not coincident with) an XMM-Newton X-ray source. Based on the source's overall properties, we discuss various options for the nature of this object and propose that it could be a young pulsar, magnetar, or a binary pulsar system. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

Absolute Flux Density Calibration of the Greenland Telescope Data for Event Horizon Telescope Observations

Authors: J. Y. Koay, K. Asada, S. Matsushita, C. -Y. Kuo, C. -W. L. Huang, C. Romero-Cañizales, S. Koyama, J. Park, W. -P. Lo, G. Bower, M. -T. Chen, S. -H. Chang, C. -C. Chen, R. Chilson, C. C. Han, P. T. P. Ho, Y. -D. Huang, M. Inoue, B. Jeter, H. Jiang, P. M. Koch, D. Kubo, C. -T. Li, C. -T. Liu, K. -Y. Liu , et al. (13 additional authors not shown)

Abstract: Starting from the observing campaign in April 2018, the Greenland Telescope (GLT) has been added as a new station of the Event Horizon Telescope (EHT) array. Visibilities on baselines to the GLT, particularly in the North-South direction, potentially provide valuable new constraints for the modeling and imaging of sources such as M87*. The GLT's location at high Northern latitudes adds unique chal… ▽ More Starting from the observing campaign in April 2018, the Greenland Telescope (GLT) has been added as a new station of the Event Horizon Telescope (EHT) array. Visibilities on baselines to the GLT, particularly in the North-South direction, potentially provide valuable new constraints for the modeling and imaging of sources such as M87*. The GLT's location at high Northern latitudes adds unique challenges to its calibration strategies. Additionally, the performance of the GLT was not optimal during the 2018 observations due to it being only partially commissioned at the time. This document describes the steps taken to estimate the various parameters (and their uncertainties) required for the absolute flux calibration of the GLT data as part of the EHT. In particular, we consider the non-optimized status of the GLT in 2018, as well as its improved performance during the 2021 EHT campaign. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 17 pages, 4 figures, EHT Memo Series 2023-L1-02

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

The Greenland Telescope: Construction, Commissioning, and Operations in Pituffik

Authors: Ming-Tang Chen, Keiichi Asada, Satoki Matsushita, Philippe Raffin, Makoto Inoue, Paul T. P. Ho, Chih-Chiang Han, Derek Kubo, Timothy Norton, Nimesh A. Patel, George Nystrom, Chih-Wei L. Huang, Pierre Martin-Cocher, Jun Yi Koay, Cristina Romero-Cañizales, Ching-Tang Liu, Teddy Huang, Kuan-Yu Liu, Tashun Wei, Shu-Hao Chang, Ryan Chilson, Peter Oshiro, Homin Jiang, Chao-Te Li, Geoffrey Bower , et al. (29 additional authors not shown)

Abstract: In 2018, the Greenland Telescope (GLT) started scientific observation in Greenland. Since then, we have completed several significant improvements and added new capabilities to the telescope system. This paper presents a full review of the GLT system, a summary of our observation activities since 2018, the lessons learned from the operations in the Arctic regions, and the prospect of the telescope… ▽ More In 2018, the Greenland Telescope (GLT) started scientific observation in Greenland. Since then, we have completed several significant improvements and added new capabilities to the telescope system. This paper presents a full review of the GLT system, a summary of our observation activities since 2018, the lessons learned from the operations in the Arctic regions, and the prospect of the telescope. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 26 pages, 11 figures, and 8 tables. This is the version of the article before publication editing, as submitted by an author to Publications of the Astronomical Society of the Pacific. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record will be added when it becomes available

Testing the linear relationship between black hole mass and variability timescale in low-luminosity AGN at submillimeter wavelengths

Authors: Bo-Yan Chen, Geoffrey C. Bower, Jason Dexter, Sera Markoff, Anthony Ridenour, Mark A. Gurwell, Ramprasad Rao, Sofia H. J. Wallström

Abstract: The variability of submillimeter emission provides a useful tool to probe the accretion physics in low-luminosity active galactic nuclei. We accumulate four years of observations using Submillimeter Array for Centaurus A, NGC 4374, NGC 4278, and NGC 5077 and one year of observations for NGC 4552 and NGC 4579. All sources are variable. We measure the characteristic timescale at which the variabilit… ▽ More The variability of submillimeter emission provides a useful tool to probe the accretion physics in low-luminosity active galactic nuclei. We accumulate four years of observations using Submillimeter Array for Centaurus A, NGC 4374, NGC 4278, and NGC 5077 and one year of observations for NGC 4552 and NGC 4579. All sources are variable. We measure the characteristic timescale at which the variability is saturated by modeling these sources' light curve as a damped random walk. We detect a timescale for all the sources except NGC 4552. The detected timescales are comparable to the orbital timescale at the event horizon scale for most sources. Combined with previous studies, we show a correlation between the the timescale and the black hole mass over three orders of magnitude. This discovery suggests the sub-mm emission is optically thin with the emission originating from the event horizon. The mass scaling relationship further suggests that a group of radio sources with a broadband spectrum that peaks at submillimeter wavelengths have similar inner accretion physics. Sources that follow this relationship may be good targets for high-resolution imaging with the Event Horizon Telescope. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 18 pages, 7 figures, 4 tables, accepted for publication by ApJ

KURVS: The outer rotation curve shapes and dark matter fractions of $z \sim 1.5 $ star-forming galaxies

Authors: Annagrazia Puglisi, Ugnė Dudzevičiūtė, Mark Swinbank, Steven Gillman, Alfred L. Tiley, Richard G. Bower, Michele Cirasuolo, Luca Cortese, Karl Glazebrook, Chris Harrison, Edo Ibar, Juan Molina, Danail Obreschkow, Kyle A. Oman, Matthieu Schaller, Francesco Shankar, Ray M. Sharples

Abstract: We present first results from the KMOS Ultra-deep Rotation Velocity Survey (KURVS), aimed at studying the outer rotation curves shape and dark matter content of 22 star-forming galaxies at $z\sim1.5$. These galaxies represent `typical' star-forming discs at $z \sim 1.5$, being located within the star-forming main sequence and stellar mass-size relation with stellar masses $9.5\leqslant$log… ▽ More We present first results from the KMOS Ultra-deep Rotation Velocity Survey (KURVS), aimed at studying the outer rotation curves shape and dark matter content of 22 star-forming galaxies at $z\sim1.5$. These galaxies represent `typical' star-forming discs at $z \sim 1.5$, being located within the star-forming main sequence and stellar mass-size relation with stellar masses $9.5\leqslant$log$(M_{\star}/\mathrm{M_{\odot}})\leqslant11.5$. We extract individual rotation curves out to 4 times the effective radius, on average, or $\sim 10-15$ kpc. Most rotation curves are flat or rising between three- and six-disc scale radii. Only three objects with dispersion-dominated dynamics ($v_{\rm rot}/σ_0\sim0.2$) have declining outer rotation curves at more than 5$σ$ significance. After accounting for seeing and pressure support, the nine rotation-dominated discs with $v_{\rm rot}/σ_0\geqslant1.5$ have average dark matter fractions of $50 \pm 20\%$ at the effective radius, similar to local discs. Together with previous observations of star-forming galaxies at cosmic noon, our measurements suggest a trend of declining dark matter fraction with increasing stellar mass and stellar mass surface density at the effective radius. Simulated EAGLE galaxies are in quantitative agreement with observations up to log$(M_{\star}R_{\rm eff}^{-2}/\mathrm{M_{\odot}kpc^{-2}}) \sim 9.2$, and over-predict the dark matter fraction of galaxies with higher mass surface densities by a factor of $\sim 3$. We conclude that the dynamics of typical rotationally-supported discs at $z \sim 1.5$ is dominated by dark matter from effective radius scales, in broad agreement with cosmological models. The tension with observations at high stellar mass surface density suggests that the prescriptions for baryonic processes occurring in the most massive galaxies (such as bulge growth and quenching) need to be reassessed. △ Less

Submitted 7 May, 2023; originally announced May 2023.

Comments: 23 pages, 9 figures. Resubmitted to MNRAS after addressing the referee's comments. Abstract slightly modified to compile with the arXiv formatting

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

EAGLE-like simulation models do not solve the entropy core problem in groups and clusters of galaxies

Authors: Edoardo Altamura, Scott T. Kay, Richard G. Bower, Matthieu Schaller, Yannick M. Bahé, Joop Schaye, Josh Borrow, Imogen Towler

Abstract: Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to inve… ▽ More Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to investigate the sensitivity of the central entropy and the shape of the profiles to changes in the sub-grid model applied to a suite of zoom-in cosmological simulations of a group of mass $M_{500} = 8.8 \times 10^{12}~{\rm M}_\odot$ and a cluster of mass $2.9 \times 10^{14}~{\rm M}_\odot$. Using our reference model, calibrated to match the stellar mass function of field galaxies, we confirm that our simulated groups and clusters contain hot gas with too high entropy in their cores. Additional simulations run without artificial conduction, metal cooling or AGN feedback produce lower entropy levels but still fail to reproduce observed profiles. Conversely, the two objects run without supernova feedback show a significant entropy increase which can be attributed to excessive cooling and star formation. Varying the AGN heating temperature does not greatly affect the profile shape, but only the overall normalisation. Finally, we compared runs with four AGN heating schemes and obtained similar profiles, with the exception of bipolar AGN heating, which produces a higher and more uniform entropy distribution. Our study leaves open the question of whether the entropy core problem in simulations, and particularly the lack of power-law cool-core profiles, arise from incorrect physical assumptions, missing physical processes, or insufficient numerical resolution. △ Less

Submitted 9 February, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS on 2023 February 06. The code used in the analysis is publicly available on the corresponding author's GitHub repository: https://github.com/edoaltamura/entropy-core-problem

The JCMT Transient Survey: Single Epoch Transients and Variability of Faint Sources

Authors: Doug Johnstone, Bhavana Lalchand, Steve Mairs, Hsien Shang, Wen Ping Chen, Geoffrey C. Bower, Gregory J. Herczeg, Jeong-Eun Lee, Jan Forbrich, Bo-Yan Chen, Carlos Contreras-Pena, Yong-Hee Lee, Wooseok Park, Colton Broughton, Spencer Plovie, The JCMT Transient Team

Abstract: Short-duration flares at millimeter wavelengths provide unique insights into the strongest magnetic reconnection events in stellar coronae, and combine with longer-term variability to introduce complications to next-generation cosmology surveys. We analyze 5.5 years of JCMT Transient Survey 850 micron submillimeter monitoring observations toward eight Gould Belt star-forming regions to search for… ▽ More Short-duration flares at millimeter wavelengths provide unique insights into the strongest magnetic reconnection events in stellar coronae, and combine with longer-term variability to introduce complications to next-generation cosmology surveys. We analyze 5.5 years of JCMT Transient Survey 850 micron submillimeter monitoring observations toward eight Gould Belt star-forming regions to search for evidence of transient events or long-duration variability from faint sources. The eight regions (30 arcmin diameter fields), including ~1200 infrared-selected YSOs, have been observed on average 47 times with integrations of approximately half an hour, or one day total spread over 5.5 years. Within this large data set, only two robust faint source detections are recovered: JW 566 in OMC 2/3 and MGM12 2864 in NGC 2023. JW 566, a Class II TTauri binary system previously identified as an extraordinary submillimeter flare, remains unique, the only clear single-epoch transient detection in this sample with a flare eight times bright than our ~4.5 sigma detection threshold of 55 mJy/beam. The lack of additional recovered flares intermediate between JW 566 and our detection limit is puzzling, if smaller events are more common than larger events. In contrast, the other submillimeter variable identified in our analysis, Source 2864, is highly variable on all observed timescales. Although Source 2864 is occasionally classified as a YSO, the source is most likely a blazar. The degree of variability across the electromagnetic spectrum may be used to aid source classification. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: Accepted for publication in the Astrophysicall Journal. 17 pages, 6 figures, 2 tables

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

Authors: Maciek Wielgus, Nicola Marchili, Ivan Marti-Vidal, Garrett K. Keating, Venkatessh Ramakrishnan, Paul Tiede, Ed Fomalont, Sara Issaoun, Joey Neilsen, Michael A. Nowak, Lindy Blackburn, Charles F. Gammie, Ciriaco Goddi, Daryl Haggard, Daeyoung Lee, Monika Moscibrodzka, Alexandra J. Tetarenko, Geoffrey C. Bower, Chi-Kwan Chan, Koushik Chatterjee, Paul M. Chesler, Jason Dexter, Sheperd S. Doeleman, Boris Georgiev, Mark Gurwell , et al. (6 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data s… ▽ More The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Journal ref: The Astrophysical Journal Letters, 930:L19 (2022)

The Science Performance of JWST as Characterized in Commissioning

Authors: Jane Rigby, Marshall Perrin, Michael McElwain, Randy Kimble, Scott Friedman, Matt Lallo, René Doyon, Lee Feinberg, Pierre Ferruit, Alistair Glasse, Marcia Rieke, George Rieke, Gillian Wright, Chris Willott, Knicole Colon, Stefanie Milam, Susan Neff, Christopher Stark, Jeff Valenti, Jim Abell, Faith Abney, Yasin Abul-Huda, D. Scott Acton, Evan Adams, David Adler , et al. (601 additional authors not shown)

Abstract: This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries f… ▽ More This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies. △ Less

Submitted 10 April, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 5th version as accepted to PASP; 31 pages, 18 figures; https://iopscience.iop.org/article/10.1088/1538-3873/acb293

Journal ref: PASP 135 048001 (2023)

BURSTT: Bustling Universe Radio Survey Telescope in Taiwan

Authors: Hsiu-Hsien Lin, Kai-yang Lin, Chao-Te Li, Yao-Huan Tseng, Homin Jiang, Jen-Hung Wang, Jen-Chieh Cheng, Ue-Li Pen, Ming-Tang Chen, Pisin Chen, Yaocheng Chen, Tomotsugu Goto, Tetsuya Hashimoto, Yuh-Jing Hwang, Sun-Kun King, Derek Kubo, Chung-Yun Kuo, Adam Mills, Jiwoo Nam, Peter Oshiro, Chang-Shao Shen, Hsien-Chun Tseng, Shih-Hao Wang, Vigo Feng-Shun Wu, Geoffrey Bower , et al. (22 additional authors not shown)

Abstract: Fast Radio Bursts (FRBs) are bright millisecond-duration radio transients that appear about 1,000 times per day, all-sky, for a fluence threshold 5 Jy ms at 600 MHz. The FRB radio-emission physics and the compact objects involved in these events are subjects of intense active debate. To better constrain source models, the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is optimized to… ▽ More Fast Radio Bursts (FRBs) are bright millisecond-duration radio transients that appear about 1,000 times per day, all-sky, for a fluence threshold 5 Jy ms at 600 MHz. The FRB radio-emission physics and the compact objects involved in these events are subjects of intense active debate. To better constrain source models, the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is optimized to discover and localize a large sample of rare, high-fluence, nearby FRBs. This is the population most amenable to multi-messenger, multi-wavelength follow-up, allowing deeper understanding of source mechanisms. BURSTT will provide horizon-to-horizon sky coverage with a half power field-of-view (FoV) of $\sim$10$^{4}$ deg$^{2}$, a 400 MHz effective bandwidth between 300-800 MHz, and sub-arcsecond localization, made possible using outrigger stations hundreds to thousands of km from the main array. Initially, BURSTT will employ 256 antennas. After tests of various antenna designs and optimization of system performance we plan to expand to 2048 antennas. We estimate that BURSTT-256 will detect and localize $\sim$100 bright ($\geq$100 Jy ms) FRBs per year. Another advantage of BURSTT's large FoV and continuous operation will be greatly enhanced monitoring of FRBs for repetition. The current lack of sensitive all-sky observations likely means that many repeating FRBs are currently cataloged as single-event FRBs. △ Less

Submitted 26 September, 2022; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: Published version

Journal ref: PASP 134 094106 (2022)

A sparse regression approach for populating dark matter halos and subhalos with galaxies

Authors: M. Icaza-Lizaola, Richard G. Bower, Peder Norberg, Shaun Cole, Matthieu Schaller

Abstract: We use sparse regression methods (SRM) to build accurate and explainable models that predict the stellar mass of central and satellite galaxies as a function of properties of their host dark matter halos. SRM are machine learning algorithms that provide a framework for modelling the governing equations of a system from data. In contrast with other machine learning algorithms, the solutions of SRM… ▽ More We use sparse regression methods (SRM) to build accurate and explainable models that predict the stellar mass of central and satellite galaxies as a function of properties of their host dark matter halos. SRM are machine learning algorithms that provide a framework for modelling the governing equations of a system from data. In contrast with other machine learning algorithms, the solutions of SRM methods are simple and depend on a relatively small set of adjustable parameters. We collect data from 35,459 galaxies from the EAGLE simulation using 19 redshift slices between $z=0$ and $z=4$ to parameterize the mass evolution of the host halos. Using an appropriate formulation of input parameters, our methodology can model satellite and central halos using a single predictive model that achieves the same accuracy as when predicted separately. This allows us to remove the somewhat arbitrary distinction between those two galaxy types and model them based only on their halo growth history. Our models can accurately reproduce the total galaxy stellar mass function and the stellar mass-dependent galaxy correlation functions ($ξ(r)$) of EAGLE. We show that our SRM model predictions of $ξ(r)$ is competitive with those from sub-halo abundance matching and might be comparable to results from extremely randomized trees. We suggest SRM as an encouraging approach for populating the halos of dark matter only simulations with galaxies and for generating mock catalogues that can be used to explore galaxy evolution or analyse forthcoming large-scale structure surveys. △ Less

Submitted 8 November, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 19 pages, 12 figures

Non-thermal radio continuum emission from young nearby stars

Authors: R. Launhardt, L. Loinard, S. A. Dzib, J. Forbrich, G. C. Bower, Th. K. Henning, A. J. Mioduszewski, S. Reffert

Abstract: Young and magnetically active low-mass stars often exhibit non-thermal coronal radio emission due to the gyration of electrons in their magnetized chromospheres. This emission is easily detectable at centimeter wavelengths with the current sensitivity of large radio interferometers like the VLA. With the aim of identifying nearby stars adequate for future accurate radio astrometric monitoring usin… ▽ More Young and magnetically active low-mass stars often exhibit non-thermal coronal radio emission due to the gyration of electrons in their magnetized chromospheres. This emission is easily detectable at centimeter wavelengths with the current sensitivity of large radio interferometers like the VLA. With the aim of identifying nearby stars adequate for future accurate radio astrometric monitoring using Very Long Baseline Interferometry (VLBI), we have used the VLA in its B\, configuration to search for radio emission at \mbox{$ν\simeq 6$\,GHz} \mbox{($λ\simeq 5$\,cm)} toward a sample of 170 nearby ($<$130\,pc) mostly young (5\,--\,500\,Myr) stars of spectral types between F4 and M2. At our mean $3σ$ detection limit of $\simeq$\,50\,$μ$Jy, we identify 31 young stars with coronal radio emission (an 18\% system detection rate) and more than 600 background (most likely extra-galactic) sources. Among the targeted stars, we find a significant decline of the detection rate with age from 56$\pm$20\% for stars with ages $\le10$\,Myr to 10$\pm$3\% for stars with ages 100\,--\,200\,Myr. No star older than 200\,Myr was detected. The detection rate also declines with $T_{\rm eff}$ from 36$\pm$10\% for stars with $T_{\rm eff}<4000$\,K to 13$\pm$3\% for earlier spectral types with $T_{\rm eff}>5000$\,K. The binarity fraction among the radio-bright stars is at least twice as high as among the radio-quiet stars. The radio-bright nearby young stars identified here provide an interesting sample for future astrometric studies using VLBI arrays aimed at searching for hitherto unknown tight binary components or even exoplanets. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 24 pages, 15 figures plus electronic figures with 31 figures, 3 tables, accepted by ApJ

The resolved chemical abundance properties within the interstellar medium of star-forming galaxies at $\mathbf{ \textit{z} \approx 1.5}$

Authors: S. Gillman, A. Puglisi, U. Dudzevičiūtė, A. M. Swinbank, A. L. Tiley, C. M. Harrison, J. Molina, R. M. Sharples, R. G. Bower, M. Cirasuolo, Edo Ibar, D. Obreschkow

Abstract: We exploit the unprecedented depth of integral field data from the KMOS Ultra-deep Rotational Velocity Survey (KURVS) to analyse the strong (H$α$) and forbidden ([NII], [SII]) emission line ratios in 22 main-sequence galaxies at $z\approx1.5$. Using the [NII]/H$α$ emission-line ratio we confirm the presence of the stellar mass $-$ gas-phase metallicity relation at this epoch, with galaxies exhibit… ▽ More We exploit the unprecedented depth of integral field data from the KMOS Ultra-deep Rotational Velocity Survey (KURVS) to analyse the strong (H$α$) and forbidden ([NII], [SII]) emission line ratios in 22 main-sequence galaxies at $z\approx1.5$. Using the [NII]/H$α$ emission-line ratio we confirm the presence of the stellar mass $-$ gas-phase metallicity relation at this epoch, with galaxies exhibiting on average 0.13$\pm$0.04 dex lower gas-phase metallicity (12+log(O/H)$_{\rm M13}$=8.40$\pm$0.03) for a given stellar mass ($\log_{10}$($M_{\rm *}$[$M_{\odot}$]=10.1$\pm$0.1) than local main-sequence galaxies. We determine the galaxy-integrated [SII] doublet ratio, with a median value of [SII]$λ$6716/$λ$6731=1.26$\pm$0.14 equivalent to an electron density of log$_{10}$($n_{\rm e}$[cm$^{-3}$])=1.95$\pm$0.12. Utilising CANDELS $HST$ multi-band imaging we define the pixel surface-mass and star-formation rate density in each galaxy and spatially resolve the fundamental metallicity relation at $z\approx1.5$, finding an evolution of 0.05$\pm$0.01 dex compared to the local relation. We quantify the intrinsic gas-phase metallicity gradient within the galaxies using the [NII]/H$α$ calibration, finding a median annuli-based gradient of $Δ$Z/$Δ$R=$-$0.015$\pm$0.005 dex kpc$^{-1}$. Finally we examine the azimuthal variations in gas-phase metallicity, which show a negative correlation with the galaxy integrated star-formation rate surface density ($r_{\rm s}$=$-$0.40, $p_{\rm s}$=0.07) but no connection to the galaxies kinematic or morphological properties nor radial variations in stellar mass surface density or star formation rate surface density. This suggests both the radial and azimuthal variations in interstellar medium properties are connected to the galaxy integrated density of recent star formation. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: Submitted for publication in MNRAS. Comments welcome. 18 pages

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

The importance of black hole repositioning for galaxy formation simulations

Authors: Yannick M. Bahé, Joop Schaye, Matthieu Schaller, Richard G. Bower, Josh Borrow, Evgenii Chaikin, Roi Kugel, Folkert Nobels, Sylvia Ploeckinger

Abstract: Active galactic nucleus (AGN) feedback from accreting supermassive black holes (SMBHs) is an essential ingredient of galaxy formation simulations. The orbital evolution of SMBHs is affected by dynamical friction that cannot be predicted self-consistently by contemporary simulations of galaxy formation in representative volumes. Instead, such simulations typically use a simple "repositioning" of SM… ▽ More Active galactic nucleus (AGN) feedback from accreting supermassive black holes (SMBHs) is an essential ingredient of galaxy formation simulations. The orbital evolution of SMBHs is affected by dynamical friction that cannot be predicted self-consistently by contemporary simulations of galaxy formation in representative volumes. Instead, such simulations typically use a simple "repositioning" of SMBHs, but the effects of this approach on SMBH and galaxy properties have not yet been investigated systematically. Based on a suite of smoothed particle hydrodynamics simulations with the SWIFT code and a Bondi-Hoyle-Lyttleton subgrid gas accretion model, we investigate the impact of repositioning on SMBH growth and on other baryonic components through AGN feedback. Across at least a factor ~1000 in mass resolution, SMBH repositioning (or an equivalent approach) is a necessary prerequisite for AGN feedback; without it, black hole growth is negligible. Limiting the effective repositioning speed to $\lesssim$ 10 km/s delays the onset of AGN feedback and severely limits its impact on stellar mass growth in the centre of massive galaxies. Repositioning has three direct physical consequences. It promotes SMBH mergers and thus accelerates their initial growth. In addition, it raises the peak density of the ambient gas and reduces the SMBH velocity relative to it, giving a combined boost to the accretion rate that can reach many orders of magnitude. Our results suggest that a more sophisticated and/or better calibrated treatment of SMBH repositioning is a critical step towards more predictive galaxy formation simulations. △ Less

Submitted 16 June, 2022; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 17 pages, 11 figures (including appendix); accepted by MNRAS. Minor updates to match accepted version

The JCMT Transient Survey: Four Year Summary of Monitoring the Submillimeter Variability of Protostars

Authors: Yong-Hee Lee, Doug Johnstone, Jeong-Eun Lee, Gregory Herczeg, Steve Mairs, Carlos Contreras-Peña, Jennifer Hatchell, Tim Naylor, Graham S. Bell, Tyler L. Bourke, Colton Broughton, Logan Francis, Aashish Gupta, Daniel Harsono, Sheng-Yuan Liu, Geumsook Park, Spencer Plovie, Gerald H. Moriarty-Schieven, Aleks Scholz, Tanvi Sharma, Paula Stella Teixeira, Yao-Te Wang, Yuri Aikawa, Geoffrey C. Bower, Huei-Ru Vivien Chen , et al. (27 additional authors not shown)

Abstract: We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources.… ▽ More We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted for publication in the Astrophysical Journal

A repeating fast radio burst source in a globular cluster

Authors: F. Kirsten, B. Marcote, K. Nimmo, J. W. T. Hessels, M. Bhardwaj, S. P. Tendulkar, A. Keimpema, J. Yang, M. P. Snelders, P. Scholz, A. B. Pearlman, C. J. Law, W. M. Peters, M. Giroletti, Z. Paragi, C. Bassa, D. M. Hewitt, U. Bach, V. Bezrukovs, M. Burgay, S. T. Buttaccio, J. E. Conway, A. Corongiu, R. Feiler, O. Forssén , et al. (41 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (… ▽ More Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (Margalit & Metzger 2018, CHIME/FRB Collaboration et al. 2020). Recently, CHIME/FRB announced the discovery (Bhardwaj et al. 2021) of the repeating FRB 20200120E, coming from the direction of the nearby grand design spiral galaxy M81. Four potential counterparts at other observing wavelengths were identified (Bhardwaj et al. 2021) but no definitive association with these sources, or M81, could be made. Here we report an extremely precise localisation of FRB 20200120E, which allows us to associate it with a globular cluster (GC) in the M81 galactic system and to place it ~2pc offset from the optical center of light of the GC. This confirms (Bhardwaj et al. 2021) that FRB 20200120E is 40 times closer than any other known extragalactic FRB. Because such GCs host old stellar populations, this association strongly challenges FRB models that invoke young magnetars formed in a core-collapse supernova as powering FRB emission. We propose, instead, that FRB 20200120E is a highly magnetised neutron star formed via either accretion-induced collapse of a white dwarf or via merger of compact stars in a binary system (Margalit et al. 2019). Alternative scenarios involving compact binary systems, efficiently formed inside globular clusters, could also be responsible for the observed bursts. △ Less

Submitted 29 September, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: Submitted. Comments welcome

Constraints on black-hole charges with the 2017 EHT observations of M87*

Authors: Prashant Kocherlakota, Luciano Rezzolla, Heino Falcke, Christian M. Fromm, Michael Kramer, Yosuke Mizuno, Antonios Nathanail, Hector Olivares, Ziri Younsi, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell, Wilfred Boland , et al. (212 additional authors not shown)

Abstract: Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87*… ▽ More Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 15 pages, 4 figures, published in PRD on May 19

ALMA and NOEMA constraints on synchrotron nebular emission from embryonic superluminous supernova remnants and radio-gamma-ray connection

Authors: Kohta Murase, Conor M. B. Omand, Deanne L. Coppejans, Hiroshi Nagai, Geoffrey C. Bower, Ryan Chornock, Derek B. Fox, Kazumi Kashiyama, Casey Law, Raffaella Margutti, Peter Meszaros

Abstract: Fast-rotating pulsars and magnetars have been suggested as the central engines of super-luminous supernovae (SLSNe) and fast radio bursts, and this scenario naturally predicts non-thermal synchrotron emission from their nascent pulsar wind nebulae (PWNe). We report results of high-frequency radio observations with ALMA and NOEMA for three SLSNe (SN 2015bn, SN 2016ard, and SN 2017egm), and present… ▽ More Fast-rotating pulsars and magnetars have been suggested as the central engines of super-luminous supernovae (SLSNe) and fast radio bursts, and this scenario naturally predicts non-thermal synchrotron emission from their nascent pulsar wind nebulae (PWNe). We report results of high-frequency radio observations with ALMA and NOEMA for three SLSNe (SN 2015bn, SN 2016ard, and SN 2017egm), and present a detailed theoretical model to calculate non-thermal emission from PWNe with an age of about 1-3 yr. We find that the ALMA data disfavors a PWN model motivated by the Crab nebula for SN 2015bn and SN 2017egm, and argue that this tension can be resolved if the nebular magnetization is very high or very low. Such models can be tested by future MeV-GeV gamma-ray telescopes such as AMEGO. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: 7 pages, 4 figures, 2 tables

Journal ref: Mon.Not.Roy.Astron.Soc. 508 (2021) 44

Polarimetric properties of Event Horizon Telescope targets from ALMA

Authors: Ciriaco Goddi, Ivan Marti-Vidal, Hugo Messias, Geoffrey C. Bower, Avery E. Broderick, Jason Dexter, Daniel P. Marrone, Monika Moscibrodzka, Hiroshi Nagai, Juan Carlos Algaba, Keiichi Asada, Geoffrey B. Crew, Jose L. Gomez, C. M. Violette Impellizzeri, Michael Janssen, Matthias Kadler, Thomas P. Krichbaum, Rocco Lico, Lynn D. Matthews, Antonios Nathanail, Angelo Ricarte, Eduardo Ros, Ziri Younsi, The Event Horizon Telescope Collaboration, Gabriele Bruni , et al. (9 additional authors not shown)

Abstract: We present the results from a full polarization study carried out with ALMA during the first VLBI campaign, which was conducted in Apr 2017 in the $λ$3mm and $λ$1.3mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud AGN.… ▽ More We present the results from a full polarization study carried out with ALMA during the first VLBI campaign, which was conducted in Apr 2017 in the $λ$3mm and $λ$1.3mm bands, in concert with the Global mm-VLBI Array (GMVA) and the Event Horizon Telescope (EHT), respectively. We determine the polarization and Faraday properties of all VLBI targets, including Sgr A*, M87, and a dozen radio-loud AGN. We detect high linear polarization fractions (2-15%) and large rotation measures (RM $>10^{3.3}-10^{5.5}$ rad m$^{-2}$). For Sgr A* we report a mean RM of $(-4.2\pm0.3) \times10^5$ rad m$^{-2}$ at 1.3 mm, consistent with measurements over the past decade, and, for the first time, an RM of $(-2.1\pm0.1) \times10^5$ rad m$^{-2}$ at 3 mm, suggesting that about half of the Faraday rotation at 1.3 mm may occur between the 3 mm photosphere and the 1.3 mm source. We also report the first unambiguous measurement of RM toward the M87 nucleus at mm wavelengths, which undergoes significant changes in magnitude and sign reversals on a one year time-scale, spanning the range from -1.2 to 0.3 $\times\,10^5$ rad m$^{-2}$ at 3 mm and -4.1 to 1.5 $\times\,10^5$ rad m$^{-2}$ at 1.3 mm. Given this time variability, we argue that, unlike the case of Sgr A*, the RM in M87 does not provide an accurate estimate of the mass accretion rate onto the black hole. We put forward a two-component model, comprised of a variable compact region and a static extended region, that can simultaneously explain the polarimetric properties observed by both the EHT and ALMA. These measurements provide critical constraints for the calibration, analysis, and interpretation of simultaneously obtained VLBI data with the EHT and GMVA. △ Less

Submitted 5 May, 2021; originally announced May 2021.

Comments: 62 pages, 26 figures, published in The Astrophysical Journal Letters on March 24, 2021

Journal ref: The Astrophysical Journal Letters, Volume 910, Issue 1, id.L14, 54 pp., 2021 March 24

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

An 86-GHz search for Pulsars in the Galactic Center with the Atacama Large Millimeter/submillimeter Array

Authors: Kuo Liu, Gregory Desvignes, Ralph P. Eatough, Ramesh Karuppusamy, Michael Kramer, Pablo Torne, Robert Wharton, Shami Chatterjee, James M. Cordes, Geoffrey B. Crew, Ciriaco Goddi, Scott M. Ransom, Helge Rottmann, Federico Abbate, Geoffrey C. Bower, Christiaan D. Brinkerink, Heino Falcke, Aristeidis Noutsos, Antonio Hernandez-Gomez, Wu Jiang, Michael D. Johnson, Ru-Sen Lu, Yurii Pidopryhora, Luciano Rezzolla, Lijing Shao , et al. (2 additional authors not shown)

Abstract: We report on the first pulsar and transient survey of the Galactic Center (GC) with the Atacama Large Millimeter/submillimeter Array (ALMA). The observations were conducted during the Global Millimeter VLBI Array campaign in 2017 and 2018. We carry out searches using timeseries of both total intensity and other polarization components in the form of Stokes parameters. We incorporate acceleration a… ▽ More We report on the first pulsar and transient survey of the Galactic Center (GC) with the Atacama Large Millimeter/submillimeter Array (ALMA). The observations were conducted during the Global Millimeter VLBI Array campaign in 2017 and 2018. We carry out searches using timeseries of both total intensity and other polarization components in the form of Stokes parameters. We incorporate acceleration and its derivative in the pulsar search, and also search in segments of the entire observation to compensate for potential orbital motion of the pulsar. While no new pulsar is found, our observations yield the polarization profile of the GC magnetar PSR J1745-2900 at mm-wavelength for the first time, which turns out to be nearly 100 % linearly polarized. Additionally, we estimate the survey sensitivity placed by both system and red noise, and evaluate its capability of finding pulsars in orbital motion with either Sgr A* or a binary companion. We show that the survey is sensitive to only the most luminous pulsars in the known population, and future observations with ALMA in Band-1 will deliver significantly deeper survey sensitivity on the GC pulsar population. △ Less

Submitted 18 April, 2021; originally announced April 2021.

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

Persistent Non-Gaussian Structure in the Image of Sagittarius A* at 86 GHz

Authors: S. Issaoun, M. D. Johnson, L. Blackburn, A. Broderick, P. Tiede, M. Wielgus, S. S. Doeleman, H. Falcke, K. Akiyama, G. C. Bower, C. D. Brinkerink, A. Chael, I. Cho, J. L. Gómez, A. Hernández-Gómez, D. Hughes, M. Kino, T. P. Krichbaum, E. Liuzzo, L. Loinard, S. Markoff, D. P. Marrone, Y. Mizuno, J. M. Moran, Y. Pidopryhora , et al. (4 additional authors not shown)

Abstract: Observations of the Galactic Center supermassive black hole Sagittarius A* (Sgr A*) with very long baseline interferometry (VLBI) are affected by interstellar scattering along our line of sight. At long radio observing wavelengths ($\gtrsim1\,$cm), the scattering heavily dominates image morphology. At 3.5 mm (86 GHz), the intrinsic source structure is no longer sub-dominant to scattering, and thus… ▽ More Observations of the Galactic Center supermassive black hole Sagittarius A* (Sgr A*) with very long baseline interferometry (VLBI) are affected by interstellar scattering along our line of sight. At long radio observing wavelengths ($\gtrsim1\,$cm), the scattering heavily dominates image morphology. At 3.5 mm (86 GHz), the intrinsic source structure is no longer sub-dominant to scattering, and thus the intrinsic emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). Long-baseline detections to the phased Atacama Large Millimeter/submillimeter Array (ALMA) in 2017 provided new constraints on the intrinsic and scattering properties of Sgr A*, but the stochastic nature of the scattering requires multiple observing epochs to reliably estimate its statistical properties. We present new observations with the GMVA+ALMA, taken in 2018, which confirm non-Gaussian structure in the scattered image seen in 2017. In particular, the ALMA-GBT baseline shows more flux density than expected for an anistropic Gaussian model, providing a tight constraint on the source size and an upper limit on the dissipation scale of interstellar turbulence. We find an intrinsic source extent along the minor axis of $\sim100\,μ$as both via extrapolation of longer wavelength scattering constraints and direct modeling of the 3.5 mm observations. Simultaneously fitting for the scattering parameters, we find an at-most modestly asymmetrical (major-to-minor axis ratio of $1.5\pm 0.2$) intrinsic source morphology for Sgr A*. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 18 pages, 10 figures, submitted to ApJ

Robust Assessment of Clustering Methods for Fast Radio Transient Candidates

Authors: Kshitij Aggarwal, Sarah Burke-Spolaor, Casey J. Law, Geoffrey C. Bower, Bryan J. Butler, Paul B. Demorest, T. Joseph W. Lazio, Justin Linford, Jessica Sydnor, Reshma Anna-Thomas

Abstract: Fast radio transient search algorithms identify signals of interest by iterating and applying a threshold on a set of matched filters. These filters are defined by properties of the transient such as time and dispersion. A real transient can trigger hundreds of search trials, each of which has to be post-processed for visualization and classification tasks. In this paper, we have explored a range… ▽ More Fast radio transient search algorithms identify signals of interest by iterating and applying a threshold on a set of matched filters. These filters are defined by properties of the transient such as time and dispersion. A real transient can trigger hundreds of search trials, each of which has to be post-processed for visualization and classification tasks. In this paper, we have explored a range of unsupervised clustering algorithms to cluster these redundant candidate detections. We demonstrate this for Realfast, the commensal fast transient search system at the Very Large Array. We use four features for clustering: sky position (l, m), time and dispersion measure (DM). We develop a custom performance metric that makes sure that the candidates are clustered into a small number of pure clusters, i.e, clusters with either astrophysical or noise candidates. We then use this performance metric to compare eight different clustering algorithms. We show that using sky location along with DM/time improves clustering performance by $\sim$10% as compared to the traditional DM/time-based clustering. Therefore, positional information should be used during clustering if it can be made available. We conduct several tests to compare the performance and generalisability of clustering algorithms to other transient datasets and propose a strategy that can be used to choose an algorithm. Our performance metric and clustering strategy can be easily extended to different single-pulse search pipelines and other astronomy and non-astronomy-based applications. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 17 pages, 7 figures, Submitted to AAS Journals (ApJ)

The Galactic Faraday rotation sky 2020

Authors: Sebastian Hutschenreuter, Craig S. Anderson, Sarah Betti, Geoffrey C. Bower, Jo-Anne Brown, Marcus Brüggen, Ettore Carretti, Tracy Clarke, Andrew Clegg, Allison Costa, Steve Croft, Cameron Van Eck, B. M. Gaensler, Francesco de Gasperin, Marijke Haverkorn, George Heald, Charles L. H. Hull, Makoto Inoue, Melanie Johnston-Hollitt, Jane Kaczmarek, Casey Law, Yik Ki Ma, David MacMahon, Sui Ann Mao, Christopher Riseley , et al. (12 additional authors not shown)

Abstract: This work gives an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. Observations of extra-Galactic sources in recent years have, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This h… ▽ More This work gives an update to existing reconstructions of the Galactic Faraday rotation sky by processing almost all Faraday rotation data sets available at the end of the year 2020. Observations of extra-Galactic sources in recent years have, among other regions, further illuminated the previously under-constrained southern celestial sky, as well as parts of the inner disc of the Milky Way. This has culminated in an all-sky data set of 55,190 data points, which is a significant expansion on the 41,330 used in previous works, hence making an updated separation of the Galactic component a promising venture. The increased source density allows us to present our results in a resolution of about $1.3\cdot 10^{-2}\, \mathrm{deg}^2$ ($46.8\,\mathrm{arcmin}^2$), which is a twofold increase compared to previous works. As for previous Faraday rotation sky reconstructions, this work is based on information field theory, a Bayesian inference scheme for field-like quantities which handles noisy and incomplete data. In contrast to previous reconstructions, we find a significantly thinner and pronounced Galactic disc with small-scale structures exceeding values of several thousand $\mathrm{rad}\,\mathrm{m}^{-2}$. The improvements can mainly be attributed to the new catalog of Faraday data, but are also supported by advances in correlation structure modeling within numerical information field theory. We furthermore give a detailed discussion on statistical properties of the Faraday rotation sky and investigate correlations to other data sets. △ Less

Submitted 8 October, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Comments: accepted in A&A; 15 pages, 12 Figures; results at https://wwwmpa.mpa-garching.mpg.de/~ensslin/research/data/faraday2020.html and http://cutouts.cirada.ca/rmcutout

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

A sparse regression approach to modeling the relation between galaxy stellar masses and their host halos

Authors: M. Icaza-Lizaola, Richard G. Bower, Peder Norberg, Shaun Cole, Matthieu Schaller, Stefan Egan

Abstract: Sparse regression algorithms have been proposed as the appropriate framework to model the governing equations of a system from data, without needing prior knowledge of the underlying physics. In this work, we use sparse regression to build an accurate and explainable model of the stellar mass of central galaxies given properties of their host dark matter (DM) halo. Our data set comprises 9,521 cen… ▽ More Sparse regression algorithms have been proposed as the appropriate framework to model the governing equations of a system from data, without needing prior knowledge of the underlying physics. In this work, we use sparse regression to build an accurate and explainable model of the stellar mass of central galaxies given properties of their host dark matter (DM) halo. Our data set comprises 9,521 central galaxies from the EAGLE hydrodynamic simulation. By matching the host halos to a DM-only simulation, we collect the halo mass and specific angular momentum at present time and for their main progenitors in 10 redshift bins from $z=0$ to $z=4$. The principal component of our governing equation is a third-order polynomial of the host halo mass, which models the stellar-mass halo-mass relation. The scatter about this relation is driven by the halo mass evolution and is captured by second and third-order correlations of the halo mass evolution with the present halo mass. An advantage of sparse regression approaches is that unnecessary terms are removed. Although we include information on halo specific angular momentum, these parameters are discarded by our methodology. This suggests that halo angular momentum has little connection to galaxy formation efficiency. Our model has a root mean square error (RMSE) of $0.167 \log_{10}(M^*/M_\odot)$, and accurately reproduces both the stellar mass function and central galaxy correlation function of EAGLE. The methodology appears to be an encouraging approach for populating the halos of DM-only simulations with galaxies, and we discuss the next steps that are required. △ Less

Submitted 13 August, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: 20 pages, 15 figures

Sphenix: Smoothed Particle Hydrodynamics for the next generation of galaxy formation simulations

Authors: Josh Borrow, Matthieu Schaller, Richard G. Bower, Joop Schaye

Abstract: Smoothed Particle Hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We introduce the Sphenix SPH scheme, which was designed with three key goals in mind: to work well with sub-grid physics modules that inject energy, be highly computationally efficient (both in terms of compute an… ▽ More Smoothed Particle Hydrodynamics (SPH) is a ubiquitous numerical method for solving the fluid equations, and is prized for its conservation properties, natural adaptivity, and simplicity. We introduce the Sphenix SPH scheme, which was designed with three key goals in mind: to work well with sub-grid physics modules that inject energy, be highly computationally efficient (both in terms of compute and memory), and to be Lagrangian. Sphenix uses a Density-Energy equation of motion, along with variable artificial viscosity and conduction, including limiters designed to work with common sub-grid models of galaxy formation. In particular, we present and test a novel limiter that prevents conduction across shocks, preventing spurious radiative losses in feedback events. Sphenix is shown to solve many difficult test problems for traditional SPH, including fluid mixing and vorticity conservation, and it is shown to produce convergent behaviour in all tests where this is appropriate. Crucially, we use the same parameters within Sphenix for the various switches throughout, to demonstrate the performance of the scheme as it would be used in production simulations. △ Less

Submitted 28 October, 2021; v1 submitted 7 December, 2020; originally announced December 2020.

Comments: Accepted for publication in MNRAS

Inconsistencies arising from the coupling of galaxy formation sub-grid models to Pressure-Smoothed Particle Hydrodynamics

Authors: Josh Borrow, Matthieu Schaller, Richard G. Bower

Abstract: Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method for solving the fluid equations that is commonplace in astrophysics, prized for its natural adaptivity and stability. The choice of variable to smooth in SPH has been the topic of contention, with smoothed pressure (P-SPH) being introduced to reduce errors at contact discontinuities relative to smoothed density schemes. Smoothed pressure… ▽ More Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method for solving the fluid equations that is commonplace in astrophysics, prized for its natural adaptivity and stability. The choice of variable to smooth in SPH has been the topic of contention, with smoothed pressure (P-SPH) being introduced to reduce errors at contact discontinuities relative to smoothed density schemes. Smoothed pressure schemes produce excellent results in isolated hydrodynamics tests; in more complex situations however, especially when coupling to the `sub-grid' physics and multiple time-stepping used in many state-of-the-art astrophysics simulations, these schemes produce large force errors that can easily evade detection as they do not manifest as energy non-conservation. Here two scenarios are evaluated: the injection of energy into the fluid (common for stellar feedback) and radiative cooling. In the former scenario, force and energy conservation errors manifest (of the same order as the injected energy), and in the latter large force errors that change rapidly over a few timesteps lead to instability in the fluid (of the same order as the energy lost to cooling). Potential ways to remedy these issues are explored with solutions generally leading to large increases in computational cost. Schemes using a Density-based formulation do not create these instabilities and as such it is recommended that they are preferred over Pressure-based solutions when combined with an energy diffusion term to reduce errors at contact discontinuities. △ Less

Submitted 23 November, 2020; originally announced November 2020.

Comments: Submitted to MNRAS; comments to be sent to joshua.borrow@durham.ac.uk

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

Monitoring the Morphology of M87* in 2009-2017 with the Event Horizon Telescope

Authors: Maciek Wielgus, Kazunori Akiyama, Lindy Blackburn, Chi-kwan Chan, Jason Dexter, Sheperd S. Doeleman, Vincent L. Fish, Sara Issaoun, Michael D. Johnson, Thomas P. Krichbaum, Ru-Sen Lu, Dominic W. Pesce, George N. Wong, Geoffrey C. Bower, Avery E. Broderick, Andrew Chael, Koushik Chatterjee, Charles F. Gammie, Boris Georgiev, Kazuhiro Hada, Laurent Loinard, Sera Markoff, Daniel P. Marrone, Richard Plambeck, Jonathan Weintroub , et al. (4 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature - a ring with azimuthal brightness asymmetry - and to quantify the imag… ▽ More The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature - a ring with azimuthal brightness asymmetry - and to quantify the image variability on event horizon scales. To address this need, we analyze M87* data collected with prototype EHT arrays in 2009, 2011, 2012, and 2013. While these observations do not contain enough information to produce images, they are sufficient to constrain simple geometric models. We develop a modeling approach based on the framework utilized for the 2017 EHT data analysis and validate our procedures using synthetic data. Applying the same approach to the observational data sets, we find the M87* morphology in 2009-2017 to be consistent with a persistent asymmetric ring of ~40 uas diameter. The position angle of the peak intensity varies in time. In particular, we find a significant difference between the position angle measured in 2013 and 2017. These variations are in broad agreement with predictions of a subset of general relativistic magnetohydrodynamic simulations. We show that quantifying the variability across multiple observational epochs has the potential to constrain the physical properties of the source, such as the accretion state or the black hole spin. △ Less

Submitted 24 September, 2020; originally announced September 2020.

Journal ref: The Astrophysical Journal, 901:67 (2020)

An Intensity Mapping Detection of Aggregate CO Line Emission at 3 mm

Authors: Garrett K. Keating, Daniel P. Marrone, Geoffrey C. Bower, Ryan P. Keenan

Abstract: We present a detection of molecular gas emission at $z\sim1-5$ using the technique of line intensity mapping. We make use of a pair of 3 mm interferometric data sets, the first from the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS), and the second from a series of Atacama Compact Array (ACA) observations conducted between 2016 and 2018, targeting the COSMOS field. At 100 GHz, w… ▽ More We present a detection of molecular gas emission at $z\sim1-5$ using the technique of line intensity mapping. We make use of a pair of 3 mm interferometric data sets, the first from the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS), and the second from a series of Atacama Compact Array (ACA) observations conducted between 2016 and 2018, targeting the COSMOS field. At 100 GHz, we measure non-zero power at 97.8% and 99.9% confidence in the ACA and ALMA data sets, respectively. In the joint result, we reject the zero-power hypothesis at 99.99% confidence, finding $\tilde{I}^{2}_{s}(ν)=770\pm210\ μ\textrm{K}^2\ \textrm{Hz}\ \textrm{sr}$. After accounting for sample variance effects, the estimated spectral shot power is $\tilde{I}^{2}_{s}(ν)=1010_{-390}^{+550}\ μ\textrm{K}^2\ \textrm{Hz}\ \textrm{sr}$. We derive a model for the various line species our measurement is expected to be sensitive to, and estimate the shot power to be $120_{-40}^{+80}\ μ\textrm{K}^2\ h^{-3}\,\textrm{Mpc}^{3}$, $200^{+120}_{-70}\ μ\textrm{K}^2\ h^{-3}\,\textrm{Mpc}^{3}$, and $90^{+70}_{-40}\ μ\textrm{K}^2\ h^{-3}\,\textrm{Mpc}^{3}$ for CO(2-1) at $z=1.3$, CO(3-2) at $z=2.5$, and CO(4-3) at $z=3.6$, respectively. Using line ratios appropriate for high-redshift galaxies, we find these results to be in good agreement with those from the CO Power Spectrum Survey (COPSS). Adopting $α_{\rm CO}=3.6\ M_{\odot}\ (\textrm{K}\ \textrm{km}\ \textrm{s}^{-1}\ \textrm{pc}^{2})^{-1}$, we estimate a cosmic molecular gas density of $ρ_{\textrm{H}_2}(z)\sim 10^{8}\ M_{\odot}\ \textrm{Mpc}^{-3}$ between $z=1-3$. △ Less

Submitted 18 August, 2020; originally announced August 2020.

Comments: 25 pages, 12 figures, 6 tables, 2 appendices. Accepted for publication in ApJ

The Star-Formation Properties of the Observed and Simulated AGN Universe: BAT vs EAGLE

Authors: Thomas M. Jackson, D. J. Rosario, D. M. Alexander, J. Scholtz, Stuart McAlpine, R. G. Bower

Abstract: In this paper we present data from 72 low redshift, hard X-ray selected AGN taken from the {\it Swift}-BAT 58 month catalogue. We utilise spectral energy distribution fitting to the optical to IR photometry in order to estimate host galaxy properties. We compare this observational sample to a volume and flux matched sample of AGN from the EAGLE hydrodynamical simulations in order to verify how acc… ▽ More In this paper we present data from 72 low redshift, hard X-ray selected AGN taken from the {\it Swift}-BAT 58 month catalogue. We utilise spectral energy distribution fitting to the optical to IR photometry in order to estimate host galaxy properties. We compare this observational sample to a volume and flux matched sample of AGN from the EAGLE hydrodynamical simulations in order to verify how accurately the simulations can reproduce observed AGN host galaxy properties. After correcting for the known +0.2 dex offset in the SFRs between EAGLE and previous observations, we find agreement in the SFR and X-ray luminosity distributions; however we find that the stellar masses in EAGLE are $0.2 - 0.4$ dex greater than the observational sample, which consequently leads to lower sSFRs. We compare these results to our previous study at high redshift, finding agreement in both the observations and simulations, whereby the widths of sSFR distributions are similar ($\sim0.4-0.6$ dex) and the median of the SFR distributions lie below the star forming main sequence by $\sim0.3-0.5$ dex across all samples. We also use EAGLE to select a sample of AGN host galaxies at high and low redshift and follow their characteristic evolution from $z=8$ to $z=0$. We find similar behaviour between these two samples, whereby star formation is quenched when the black hole goes through its phase of most rapid growth. Utilising EAGLE we find that 23\% of AGN selected at $z\sim0$ are also AGN at high redshift, and that their host galaxies are among the most massive objects in the simulation. Overall we find EAGLE reproduces the observations well, with some minor inconsistencies ($\sim$ 0.2 dex in stellar masses and $\sim$ 0.4 dex in sSFRs). △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 14 pages, 6 figures, 3 tables, accepted for publication in MNRAS

A Distant Fast Radio Burst Associated to its Host Galaxy with the Very Large Array

Authors: C. J. Law, B. J. Butler, J. X. Prochaska, B. Zackay, S. Burke-Spolaor, A. Mannings, N. Tejos, A. Josephy, B. Andersen, P. Chawla, K. E. Heintz, K. Aggarwal, G. C. Bower, P. B. Demorest, C. D. Kilpatrick, T. J. W. Lazio, J. Linford, R. Mckinven, S. Tendulkar, S. Simha

Abstract: We present the discovery and subarcsecond localization of a new Fast Radio Burst with the Karl G. Jansky Very Large Array and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc/cm3. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hou… ▽ More We present the discovery and subarcsecond localization of a new Fast Radio Burst with the Karl G. Jansky Very Large Array and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc/cm3. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hours of VLA observing and 153 hours of CHIME/FRB observing. We describe a suite of statistical and data quality tests we used to verify the significance of the event and its localization precision. Follow-up optical/infrared photometry with Keck and Gemini associate the FRB to a pair of galaxies with $\rm{r}\sim23$ mag. The false-alarm rate for radio transients of this significance that are associated with a host galaxy is roughly $3\times10^{-4}\ \rm{hr}^{-1}$. The two putative host galaxies have similar photometric redshifts of $z_{\rm{phot}}\sim0.6$, but different colors and stellar masses. Comparing the host distance to that implied by the dispersion measure suggests a modest (~ 50 pc/cm3) electron column density associated with the FRB environment or host galaxy/galaxies. △ Less

Submitted 4 July, 2020; originally announced July 2020.

Comments: Submitted to AAS Journals (ApJ) and revised for referee comments

VLA/Realfast Detection of a Burst from FRB180916.J0158+65 and Tests for Periodic Activity

Authors: Kshitij Aggarwal, Casey J. Law, Sarah Burke-Spolaor, Geoffrey Bower, Bryan J. Butler, Paul Demorest, Justin Linford, T. J. W. Lazio

Abstract: We report on the detection of a burst from FRB180916 by realfast/VLA and present software for interpreting fast radio bursts (FRB) periodicity. We demonstrate a range of periodicity analyses with bursts from FRB180916, FRB121102 and FRB180814. Our results for FRB180916 and FRB121102 are consistent with published results. For FRB180814, we did not detect any significant periodic episodes. The realf… ▽ More We report on the detection of a burst from FRB180916 by realfast/VLA and present software for interpreting fast radio bursts (FRB) periodicity. We demonstrate a range of periodicity analyses with bursts from FRB180916, FRB121102 and FRB180814. Our results for FRB180916 and FRB121102 are consistent with published results. For FRB180814, we did not detect any significant periodic episodes. The realfast-detected and other high-frequency bursts for FRB180916 tend to lie at the beginning of the activity window, indicating a possible phase-frequency relation. The python package $\texttt{frbpa}$ can be used to reproduce and expand on this analysis to test models for repeating FRBs. △ Less

Submitted 26 June, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: Published in Research Notes of the AAS

Journal ref: Research Notes of the AAS, 2020, Volume 4, Issue 6, id.94

Quasar Sightline and Galaxy Evolution (QSAGE) survey -- II. Galaxy overdensities around UV luminous quasars at z=1-2

Authors: J. P. Stott, R. M. Bielby, F. Cullen, J. N. Burchett, N. Tejos, M. Fumagalli, R. A. Crain, S. L. Morris, N. Amos, R. G. Bower, J. X. Prochaska

Abstract: We demonstrate that the UV brightest quasars at z=1-2 live in overdense environments. This is based on an analysis of deep Hubble Space Telescope WFC3 G141 grism spectroscopy of the galaxies along the lines-of-sight to UV luminous quasars in the redshift range z=1-2. This constitutes some of the deepest grism spectroscopy performed by WFC3, with 4 roll angles spread over a year of observations to… ▽ More We demonstrate that the UV brightest quasars at z=1-2 live in overdense environments. This is based on an analysis of deep Hubble Space Telescope WFC3 G141 grism spectroscopy of the galaxies along the lines-of-sight to UV luminous quasars in the redshift range z=1-2. This constitutes some of the deepest grism spectroscopy performed by WFC3, with 4 roll angles spread over a year of observations to mitigate the effect of overlapping spectra. Of the 12 quasar fields studied, 8 display evidence for a galaxy overdensity at the redshift of the quasar. One of the overdensities, PG0117+213 at z=1.50, has potentially 36 spectroscopically confirmed members, consisting of 19 with secure redshifts and 17 with single-line redshifts, within a cylinder of radius ~700 kpc. Its halo mass is estimated to be log (M/Msol)=14.7. This demonstrates that spectroscopic and narrow-band observations around distant UV bright quasars may be an excellent route for discovering protoclusters. Our findings agree with previous hints from statistical observations of the quasar population and theoretical works, as feedback regulated black hole growth predicts a correlation between quasar luminosity and halo mass. We also present the high signal-to-noise rest-frame optical spectral and photometric properties of the quasars themselves. △ Less

Submitted 19 August, 2020; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: Published in MNRAS. This is the author accepted version. A small number of typos were corrected after journal proof stage

Journal ref: 2020, MNRAS, 497, 3083

A Demonstration of Extremely Low Latency $γ$-ray, X-Ray & UV Follow-Up of a Millisecond Radio Transient

Authors: Aaron Tohuvavohu, Casey J. Law, Jamie A. Kennea, Elizabeth A. K. Adams, Kshitij Aggarwal, Geoffrey Bower, Sarah Burke-Spolaor, Bryan J. Butler, John M. Cannon, S. Bradley Cenko, James DeLaunay, Paul Demorest, Maria R. Drout, Philip A. Evans, Alec S. Hirschauer, T. J. W. Lazio, Justin Linford, Francis E. Marshall, K. McQuinn, Emily Petroff, Evan D. Skillman

Abstract: We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved aro… ▽ More We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5σ$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3σ$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: Technical note and capability update for the community. We encourage low latency FRB alerts from relevant facilities to enable this science

Redshift evolution of the hot intracluster gas metallicity in the C-EAGLE cluster simulations

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Yannick M. Bahe, Richard G. Bower

Abstract: The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investi… ▽ More The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investigate the redshift evolution of metallicity, with particular focus on the cluster outskirts. The early enrichment model, in which the majority of metals are produced in the core of cluster progenitors at high redshift, suggests that metals in cluster outskirts have not significantly evolved since $z=2$. With the C-EAGLE sample, we find reasonable agreement with the early enrichment model as there is very little scatter in the metallicity abundance at large radius across the whole sample, out to at least $z=2$. The exception is Fe for which the radial dependence of metallicity was found to evolve at low redshift as a result of being mainly produced by Type Ia supernovae, which are more likely to be formed at later times than core-collapse supernovae. We also found considerable redshift evolution of metal abundances in the cores of the C-EAGLE clusters which has not been seen in other simulations or observation based metallicity studies. Since we find this evolution to be driven by accretion of low metallicity gas, it suggests that the interaction between outflowing, AGN heated material and the surrounding gas is important for determining the core abundances in clusters. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 18 pages, 16 figures, submitted to MNRAS

Galactic inflow and wind recycling rates in the EAGLE simulations

Authors: Peter D. Mitchell, Joop Schaye, Richard G. Bower

Abstract: The role of galactic wind recycling represents one of the largest unknowns in galaxy evolution, as any contribution of recycling to galaxy growth is largely degenerate with the inflow rates of first-time infalling material, and the rates with which outflowing gas and metals are driven from galaxies. We present measurements of the efficiency of wind recycling from the EAGLE cosmological simulation… ▽ More The role of galactic wind recycling represents one of the largest unknowns in galaxy evolution, as any contribution of recycling to galaxy growth is largely degenerate with the inflow rates of first-time infalling material, and the rates with which outflowing gas and metals are driven from galaxies. We present measurements of the efficiency of wind recycling from the EAGLE cosmological simulation project, leveraging the statistical power of large-volume simulations that reproduce a realistic galaxy population. We study wind recycling at the halo scale, i.e. gas that has been ejected beyond the halo virial radius, and at the galaxy scale, i.e. gas that has been ejected from the ISM to at least $\approx 10 \, \%$ of the virial radius (thus excluding smaller-scale galactic fountains). Galaxy-scale wind recycling is generally inefficient, with a characteristic return timescale that is comparable or longer than a Hubble time, and with an efficiency that clearly peaks at the characteristic halo mass of $M_{200} = 10^{12} \, \mathrm{M_\odot}$. Correspondingly, the majority of gas being accreted onto galaxies in EAGLE is infalling for the first time. At the halo scale, the efficiency of recycling onto haloes differs by orders of magnitude from values assumed by semi-analytic galaxy formation models. Differences in the efficiency of wind recycling with other hydrodynamical simulations are currently difficult to assess, but are likely smaller. We are able to show that the fractional contribution of wind recycling to galaxy growth is smaller in EAGLE than in some other simulations. We find that cumulative first-time gas accretion rates at the virial radius are reduced relative to the expectation from dark matter accretion for haloes with mass, $M_{200} < 10^{12} \, \mathrm{M_\odot}$, indicating efficient preventative feedback on halo scales. △ Less

Submitted 7 August, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: MNRAS accepted, 22 pages

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabili… ▽ More Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 15 figures, accepted for publication in A&A

A repeating fast radio burst source localised to a nearby spiral galaxy

Authors: B. Marcote, K. Nimmo, J. W. T. Hessels, S. P. Tendulkar, C. G. Bassa, Z. Paragi, A. Keimpema, M. Bhardwaj, R. Karuppusamy, V. M. Kaspi, C. J. Law, D. Michilli, K. Aggarwal, B. Andersen, A. M. Archibald, K. Bandura, G. C. Bower, P. J. Boyle, C. Brar, S. Burke-Spolaor, B. J. Butler, T. Cassanelli, P. Chawla, P. Demorest, M. Dobbs , et al. (29 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the ho… ▽ More Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 61 pages, 11 figures, 3 tables. Published in Nature

From Peculiar Morphologies to Hubble-type Spirals: The relation between galaxy dynamics and morphology in star-forming galaxies at z~1.5

Authors: S. Gillman, A. L. Tiley, A. M. Swinbank, C. M. Harrison, Ian Smail, U. Dudzevičiūtė, R. M. Sharples, L. Cortese, D. Obreschkow, R. G. Bower, T. Theuns, M. Cirasuolo, D. B. Fisher, K. Glazebrook, Edo Ibar, J. Trevor Mendel, Sarah M. Sweet

Abstract: We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$α$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular m… ▽ More We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$α$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$α$ star formation rate surface density $Σ_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. Sérsic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium. △ Less

Submitted 5 February, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Comments: Accepted for publication in MNRAS. Comments welcome. 25 pages, 14 figures. The full pdf, which includes an additional 17 pages of supplementary figures is available on this http://astro.dur.ac.uk/KROSS/publications.html

Hydrostatic mass estimates of massive galaxy clusters: a study with varying hydrodynamics flavours and non-thermal pressure support

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Richard G. Bower, Matthieu Schaller

Abstract: We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but s… ▽ More We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but still producing clusters which broadly match observations. At the same mass resolution as BAHAMAS, for the most massive clusters ($M_{500} > 10^{15}$ M$_{\odot}$), we find changes in the SPH method produce the greatest differences in the final halo, while the subgrid models dominate at lower mass. By calculating the mass of all of the clusters using different permutations of the pressure, temperature and density profiles, created with either the true simulated data or mock spectroscopic data, we find that the spectroscopic temperature causes a bias in the hydrostatic mass estimates which increases with the mass of the cluster, regardless of the SPH flavour used. For the most massive clusters, the estimated mass of the cluster using spectroscopic density and temperature profiles is found to be as low as 50 per cent of the true mass compared to $\sim$ 90 per cent for low mass clusters. When including a correction for non-thermal pressure, the spectroscopic hydrostatic mass estimates are less biased on average and the mass dependence of the bias is reduced, although the scatter in the measurements does increase. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 21 pages, 14 figures. Accepted for publication by MNRAS

Galactic outflow rates in the EAGLE simulations

Authors: Peter D. Mitchell, Joop Schaye, Richard G. Bower, Robert A. Crain

Abstract: We present measurements of galactic outflow rates from the EAGLE suite of cosmological simulations. We find that gas is removed from the interstellar medium (ISM) of central galaxies with a dimensionless mass loading factor that scales approximately with circular velocity as $V_{\mathrm{c}}^{-3/2}$ in the low-mass regime where stellar feedback dominates. Feedback from active galactic nuclei (AGN)… ▽ More We present measurements of galactic outflow rates from the EAGLE suite of cosmological simulations. We find that gas is removed from the interstellar medium (ISM) of central galaxies with a dimensionless mass loading factor that scales approximately with circular velocity as $V_{\mathrm{c}}^{-3/2}$ in the low-mass regime where stellar feedback dominates. Feedback from active galactic nuclei (AGN) causes an upturn in the mass loading for halo masses $> 10^{12} \, \mathrm{M_\odot}$. We find that more gas outflows through the halo virial radius than is removed from the ISM of galaxies, particularly at low redshift, implying substantial mass loading within the circum-galactic medium (CGM). Outflow velocities span a wide range at a given halo mass/redshift, and on average increase positively with redshift and halo mass up to $M_{200} \sim 10^{12} \, \mathrm{M_\odot}$. Outflows exhibit a bimodal flow pattern on circum-galactic scales, aligned with the galactic minor axis. We present a number of like-for-like comparisons to outflow rates from other recent cosmological hydrodynamical simulations, and show that comparing the propagation of galactic winds as a function of radius reveals substantial discrepancies between different models. Relative to some other simulations, EAGLE favours a scenario for stellar feedback where agreement with the galaxy stellar mass function is achieved by removing smaller amounts of gas from the ISM, but with galactic winds that then propagate and entrain ambient gas out to larger radii. △ Less

Submitted 1 April, 2020; v1 submitted 21 October, 2019; originally announced October 2019.

Comments: MNRAS accepted