Showing 1–50 of 164 results for author: Cooke, J

A two-minute burst of highly polarised radio emission originating from low Galactic latitude

Authors: Dougal Dobie, Andrew Zic, Lucy S. Oswald, Joshua Pritchard, Marcus E. Lower, Ziteng Wang, Hao Qiu, Natasha Hurley-Walker, Yuanming Wang, Emil Lenc, David L. Kaplan, Akash Anumarlapudi, Katie Auchettl, Matthew Bailes, Andrew D. Cameron, Jeffrey Cooke, Adam Deller, Laura N. Driessen, James Freeburn, Tara Murphy, Ryan M. Shannon, Adam J. Stewart

Abstract: Several sources of repeating coherent bursts of radio emission with periods of many minutes have now been reported in the literature. These ``ultra-long period'' (ULP) sources have no clear multi-wavelength counterparts and challenge canonical pulsar emission models, leading to debate regarding their nature. In this work we report the discovery of a bright, highly-polarised burst of radio emission… ▽ More Several sources of repeating coherent bursts of radio emission with periods of many minutes have now been reported in the literature. These ``ultra-long period'' (ULP) sources have no clear multi-wavelength counterparts and challenge canonical pulsar emission models, leading to debate regarding their nature. In this work we report the discovery of a bright, highly-polarised burst of radio emission at low Galactic latitude as part of a wide-field survey for transient and variable radio sources. ASKAP\,J175534.9$-$252749.1 does not appear to repeat, with only a single intense two-minute $\sim 200$\,mJy burst detected from 60~hours of observations. The burst morphology and polarisation properties are comparable to those of classical pulsars but the duration is more than one hundred times longer, analogous to ULPs. No comparable bursts are detected in the rest of our widefield survey to date. Combined with the existing ULP population, this suggests that these sources have a strong Galactic latitude dependence and hints at an unexplored population of transient and variable radio sources in the thin disk of the Milky Way. The resemblance of this burst with both ULPs and pulsars calls for a unified coherent emission model for objects with spin periods from milliseconds to tens of minutes. However, whether or not these are all neutron stars or have the same underlying power source remains open for debate. △ Less

Submitted 18 June, 2024; originally announced June 2024.

The Evolution of Turbulence Producing Motions in the ABL Across a Natural Roughness Transition

Authors: Justin P. Cooke, Douglas J. Jerolmack, George I. Park

Abstract: Landforms such as sand dunes act as roughness elements to Atmospheric Boundary Layer (ABL) flows, triggering the development of new scales of turbulent motions. These turbulent motions, in turn, energize and kick-up sand particles, influencing sediment transport and ultimately the formation and migration of dunes -- with knock on consequences for dust emission. While feedbacks between flow and for… ▽ More Landforms such as sand dunes act as roughness elements to Atmospheric Boundary Layer (ABL) flows, triggering the development of new scales of turbulent motions. These turbulent motions, in turn, energize and kick-up sand particles, influencing sediment transport and ultimately the formation and migration of dunes -- with knock on consequences for dust emission. While feedbacks between flow and form have been studied at the scale of dunes, research has not examined how the development of an Internal Boundary Layer (IBL) over the entire dune field influences sediment-transporting turbulence. Here, we deploy large-eddy simulation of an ABL encountering a natural roughness transition: the sand dunes at White Sands National Park, New Mexico. We analyze turbulence producing motions and how they change as the IBL grows over the dune field. Frequency spectrum and Reynolds shear stress profiles show that IBL thickness determines the largest scales of turbulence. More, the developing IBL enhances the frequency, magnitude and duration of sweep and ejection events -- turbulence producing motions whose peaks systematically migrate away from the wall as the IBL thickens. Because sweep and ejection events are known to drive sediment transport, our findings provide a mechanism for coupling the co-evolution of the landscape and the ABL flow over it. More broadly, our results have implications for how roughness transitions influence the transport of pollutants, particulates, heat, and moisture. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 31 pages, 15 figures, 2 tables

Lethal surface ozone concentrations are possible on habitable zone exoplanets

Authors: G. J. Cooke, D. R. Marsh, C. Walsh, F. Sainsbury-Martinez

Abstract: Ozone ($\textrm{O}_3$) is important for the survival of life on Earth because it shields the surface from ionising ultraviolet (UV) radiation. However, the existence of $\textrm{O}_3$ in Earth's atmosphere is not always beneficial. Resulting from anthropogenic activity, $\textrm{O}_3$ exists as a biologically harmful pollutant at the surface when it forms in the presence of sunlight and other poll… ▽ More Ozone ($\textrm{O}_3$) is important for the survival of life on Earth because it shields the surface from ionising ultraviolet (UV) radiation. However, the existence of $\textrm{O}_3$ in Earth's atmosphere is not always beneficial. Resulting from anthropogenic activity, $\textrm{O}_3$ exists as a biologically harmful pollutant at the surface when it forms in the presence of sunlight and other pollutants. As a strong oxidiser, $\textrm{O}_3$ can be lethal to several different organisms; thus, when assessing the potential habitability of an exoplanet, a key part is determining whether toxic gases could be present at its surface. Using the Whole Atmosphere Community Climate Model version 6 (WACCM6; a three-dimensional chemistry-climate model), twelve atmospheric simulations of the terrestrial exoplanet TRAPPIST-1 e are performed with a variety of $\textrm{O}_2$ concentrations and assuming two different stellar spectra proposed in the literature. Four atmospheric simulations of the exoplanet Proxima Centauri b are also included. Some scenarios for both exoplanets exhibit time-averaged surface $\textrm{O}_3$ mixing ratios exceeding harmful levels of 40 ppbv, with 2200 ppbv the maximum concentration found in the cases simulated. These concentrations are toxic and can be fatal to most life on Earth. In other scenarios $\textrm{O}_3$ remains under harmful limits over a significant fraction of the surface, despite there being present regions which may prove inhospitable. In the case that $\textrm{O}_3$ is detected in a terrestrial exoplanet's atmosphere, determining the surface concentration is an important step when evaluating a planet's habitability. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 24 pages, 8 figures. Accepted in The Planetary Science Journal

A Fast-cadenced Search for Gamma-Ray Burst Orphan Afterglows with the Deeper, Wider, Faster Programme

Authors: James Freeburn, Jeff Cooke, Anais Möller, Dougal Dobie, Jielai Zhang, Om Sharan Salafia, Karelle Siellez, Katie Auchettl, Simon Goode, Timothy M. C. Abbott, Igor Andreoni, Rebecca Allen, Natasha Van Bemmel, Sara Webb

Abstract: The relativistic outflows that produce Long GRBs (LGRBs) can be described by a structured jet model where prompt $γ$-ray emission is restricted to a narrow region in the jet's core. Viewing the jet off-axis from the core, a population of afterglows without an associated GRB detection can be predicted. In this work, we conduct an archival search for these `orphan' afterglows (OAs) with minute-caden… ▽ More The relativistic outflows that produce Long GRBs (LGRBs) can be described by a structured jet model where prompt $γ$-ray emission is restricted to a narrow region in the jet's core. Viewing the jet off-axis from the core, a population of afterglows without an associated GRB detection can be predicted. In this work, we conduct an archival search for these `orphan' afterglows (OAs) with minute-cadence, deep ($g\sim23$) data from the Dark Energy Camera (DECam) taken as part of the Deeper, Wider, Faster programme (DWF). We introduce a method to select fast-evolving OA candidates within DWF data that comprises a machine learning model, based on a realistic synthetic population of OAs. Using this classifier, we recover 51 OA candidates. Of these candidates, 42 are likely flare events from M-class stars. The remaining nine possess quiescent, coincident sources in archival data with angular profiles consistent with a star and are inconsistent with the expected population of LGRB host galaxies. We therefore conclude that these are likely Galactic events. We calculate an upper limit on the rate of OAs down to $g<22$ AB mag of 7.46\,deg$^{-2}$yr$^{-1}$ using our criteria and constrain possible jet structures. We also place an upper limit of the characteristic angle between the $γ$-ray emitting region and the jet's half opening angle. For a smooth power-law and a power-law with core jet model respectively, these values are $58.3^{\circ}$ and $56.6^{\circ}$, for a power-law index of 0.8 and $75.3^{\circ}$ and $76.8^{\circ}$ for a power-law index of 1.2. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS

Latest Developments and Opportunities in Sky Survey

Authors: Anthony Brown, Federica Bianco, Varun Bhalerao, Shri Kulkarni, Jeffery Cooke, David H. Reitze, Pranav Sharma, Ashish Mahabal

Abstract: Policy Brief on "Latest Developments and Opportunities in Sky Survey", distilled from the corresponding panel that was part of the discussions during S20 Policy Webinar on Astroinformatics for Sustainable Development held on 6-7 July 2023. Sky surveys have been a crucial tool in advancing our understanding of the Universe. The last few decades have seen an explosion in the number and scope of sk… ▽ More Policy Brief on "Latest Developments and Opportunities in Sky Survey", distilled from the corresponding panel that was part of the discussions during S20 Policy Webinar on Astroinformatics for Sustainable Development held on 6-7 July 2023. Sky surveys have been a crucial tool in advancing our understanding of the Universe. The last few decades have seen an explosion in the number and scope of sky surveys, both ground-based and space-based. This growth has led to a wealth of data that has enabled us to make significant advances in many areas of astronomy, and help understand the physics of the universe. They have helped us discover new astronomical objects, the origin of the elements, dark matter and dark energy, the accelerated expansion of the universe, and gravitational waves. They have helped us study the distribution of neutral and ionized matter in the Universe and test our theories about the origin and evolution of galaxies, stars, and planets. We explore recent advances and potential avenues in sky surveys, and examine how these developments may impact the field of international astronomical research. The policy webinar took place during the G20 presidency in India (2023). A summary based on the seven panels can be found here: arxiv:2401.04623. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 5 pages. The panel videos including keynotes and the white papers are available on the S20 site at: https://s20india.org/science-policy-webinar-astroinformatics-for-sustainable-development/

FLIMFLAM DR1: The First Constraints on the Cosmic Baryon Distribution from 8 FRB sightlines

Authors: Ilya S. Khrykin, Metin Ata, Khee-Gan Lee, Sunil Simha, Yuxin Huang, J. Xavier Prochaska, Nicolas Tejos, Keith W. Bannister, Jeff Cooke, Cherie K. Day, Adam Deller, Marcin Glowacki, Alexa C. Gordon, Clancy W. James, Lachlan Marnoch, Ryan. M. Shannon, Jielai Zhang, Lucas Bernales-Cortes

Abstract: The dispersion measure of fast radio bursts (FRBs), arising from the interactions of the pulses with free electrons along the propagation path, constitutes a unique probe of the cosmic baryon distribution. Their constraining power is further enhanced in combination with observations of the foreground large-scale structure and intervening galaxies. In this work, we present the first constraints on… ▽ More The dispersion measure of fast radio bursts (FRBs), arising from the interactions of the pulses with free electrons along the propagation path, constitutes a unique probe of the cosmic baryon distribution. Their constraining power is further enhanced in combination with observations of the foreground large-scale structure and intervening galaxies. In this work, we present the first constraints on the partition of the cosmic baryons between the intergalactic medium (IGM) and circumgalactic medium (CGM), inferred from the FLIMFLAM spectroscopic survey. In its first data release, the FLIMFLAM survey targeted galaxies in the foreground of 8 localized FRBs. Using Bayesian techniques, we reconstruct the underlying ~Mpc-scale matter density field that is traced by the IGM gas. Simultaneously, deeper spectroscopy of intervening foreground galaxies (at impact parameters $b_\perp \lesssim r_{200}$) and the FRB host galaxies constrains the contribution from the CGM. Applying Bayesian parameter inference to our data and assuming a fiducial set of priors, we infer the IGM cosmic baryon fraction to be $f_{\rm igm}=0.59^{+0.11}_{-0.10}$, and a CGM gas fraction of $f_{\rm gas} = 0.55^{+0.26}_{-0.29}$ for $10^{10}\,M_\odot \lesssim M_{\rm halo}\lesssim 10^{13}\,M_\odot$ halos. The mean FRB host dispersion measure (rest-frame) in our sample is $\langle \rm{DM_{host}}\rangle = 90^{+29}_{-19}\rm{pc~cm^{-3}}$, of which $\langle{\rm DM_{host}^{unk}}\rangle =69^{+28}_{-19}~\rm{pc~cm^{-3}}$ arises from the host galaxy ISM and/or the FRB progenitor environment. While our current $f_{\rm igm}$ and $f_{\rm gas}$ uncertainties are too broad to constrain most galactic feedback models, this result marks the first measurement of the IGM and CGM baryon fractions, as well as the first systematic separation of the FRB host dispersion measure into two components: arising from the halo and from the inner ISM/FRB engine. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 21 pages, 9 figures, 5 tables; Submitted to ApJ

AstroInformatics: Recommendations for Global Cooperation

Authors: Ashish Mahabal, Pranav Sharma, Rana Adhikari, Mark Allen, Stefano Andreon, Varun Bhalerao, Federica Bianco, Anthony Brown, S. Bradley Cenko, Paula Coehlo, Jeffery Cooke, Daniel Crichton, Chenzhou Cui, Reinaldo de Carvalho, Richard Doyle, Laurent Eyer, Bernard Fanaroff, Christopher Fluke, Francisco Forster, Kevin Govender, Matthew J. Graham, Renée Hložek, Puji Irawati, Ajit Kembhavi, Juna Kollmeier , et al. (23 additional authors not shown)

Abstract: Policy Brief on "AstroInformatics, Recommendations for Global Collaboration", distilled from panel discussions during S20 Policy Webinar on Astroinformatics for Sustainable Development held on 6-7 July 2023. The deliberations encompassed a wide array of topics, including broad astroinformatics, sky surveys, large-scale international initiatives, global data repositories, space-related data, regi… ▽ More Policy Brief on "AstroInformatics, Recommendations for Global Collaboration", distilled from panel discussions during S20 Policy Webinar on Astroinformatics for Sustainable Development held on 6-7 July 2023. The deliberations encompassed a wide array of topics, including broad astroinformatics, sky surveys, large-scale international initiatives, global data repositories, space-related data, regional and international collaborative efforts, as well as workforce development within the field. These discussions comprehensively addressed the current status, notable achievements, and the manifold challenges that the field of astroinformatics currently confronts. The G20 nations present a unique opportunity due to their abundant human and technological capabilities, coupled with their widespread geographical representation. Leveraging these strengths, significant strides can be made in various domains. These include, but are not limited to, the advancement of STEM education and workforce development, the promotion of equitable resource utilization, and contributions to fields such as Earth Science and Climate Science. We present a concise overview, followed by specific recommendations that pertain to both ground-based and space data initiatives. Our team remains readily available to furnish further elaboration on any of these proposals as required. Furthermore, we anticipate further engagement during the upcoming G20 presidencies in Brazil (2024) and South Africa (2025) to ensure the continued discussion and realization of these objectives. The policy webinar took place during the G20 presidency in India (2023). Notes based on the seven panels will be separately published. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 7 pages

Lyman-alpha at Cosmic Noon II: The relationship between kinematics and Lyman-alpha in z~2-3 Lyman Break Galaxies

Authors: Garry Foran, Jeff Cooke, Emily Wisnioski, Naveen Reddy, Charles Steidel

Abstract: We report for the first time a relationship between galaxy kinematics and net Lyman-alpha equivalent width (net Lya EW) in star forming galaxies during the epoch of peak cosmic star formation. Building on the previously reported broadband imaging segregation of Lya-emitting and Lya-absorbing Lyman break galaxies (LBGs) at z~2 (Paper I in this series) and previously at z~3, we use the Lya spectral… ▽ More We report for the first time a relationship between galaxy kinematics and net Lyman-alpha equivalent width (net Lya EW) in star forming galaxies during the epoch of peak cosmic star formation. Building on the previously reported broadband imaging segregation of Lya-emitting and Lya-absorbing Lyman break galaxies (LBGs) at z~2 (Paper I in this series) and previously at z~3, we use the Lya spectral type classification method to study the relationship between net Lya EW and nebular emission-line kinematics determined from IFU spectroscopy. We show that z~2 and z~3 LBGs segregate in colour-magnitude space according to their kinematic properties and Lyman-alpha spectral type, and conclude that LBGs with Lya dominant in absorption are almost exclusively rotation-dominated (presumably disc-like) systems, and LBGs with Lya dominant in emission characteristically have dispersion-dominated kinematics. We quantify the relationship between the strength of rotational dynamic support and net Lya EW, and demonstrate the consistency of our result with other properties that scale with net Lya EW and kinematics. Based on these findings, we suggest a method by which large samples of rotation- and dispersion-dominated galaxies might be selected using broadband imaging in as few as three filters and/or net Lya EW alone. Application of this method will enable an understanding of galaxy kinematic behaviour over large scales in datasets from current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs in redshift ranges from z~2-6 across many hundreds to thousands of Mpc. Finally, we speculate that the combination of our result linking net Lya EW and nebular emission-line kinematics with the known large-scale clustering behaviour of Lya-absorbing and Lya-emitting LBGs is evocative of a nascent morphology-density relation at z~2-3. △ Less

Submitted 7 January, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: Published in PASA (Version 2: Typos corrected and minor clarifying edits in S3.1 & S3.2.1)

Journal ref: Publications of the Astronomical Society of Australia (2023), 1-20

Lyman-alpha at Cosmic Noon I: Ly-alpha Spectral Type Selection of z ~ 2-3 Lyman Break Galaxies with Broadband Imaging

Authors: Garry Foran, Jeff Cooke, Naveen Reddy, Charles Steidle, Alice Shapley

Abstract: High-redshift Lyman break galaxies (LBGs) are efficiently selected in deep images using as few as three broadband filters, and have been shown to have multiple intrinsic and small- to large-scale environmental properties related to Lyman-alpha. In this paper we demonstrate a statistical relationship between net Lyman-alpha equivalent width (net Lya EW) and the optical broadband photometric propert… ▽ More High-redshift Lyman break galaxies (LBGs) are efficiently selected in deep images using as few as three broadband filters, and have been shown to have multiple intrinsic and small- to large-scale environmental properties related to Lyman-alpha. In this paper we demonstrate a statistical relationship between net Lyman-alpha equivalent width (net Lya EW) and the optical broadband photometric properties of LBGs at z~2. We show that LBGs with the strongest net Lya EW in absorption (aLBGs) and strongest net Lya EW in emission (eLBGs) separate into overlapping but discrete distributions in $(U_n-R)$ colour and $R$-band magnitude space, and use this segregation behaviour to determine photometric criteria by which sub-samples with a desired Lya spectral type can be selected using data from as few as three broadband optical filters. We propose application of our result to current and future large-area and all-sky photometric surveys that will select hundreds of millions of LBGs across many hundreds to thousands of Mpc, and for which spectroscopic follow-up to obtain Lya spectral information is prohibitive. To this end, we use spectrophotometry of composite spectra derived from a sample of 798 LBGs divided into quartiles on the basis of net Lya EW to calculate criteria for the selection of Lya absorbing and Lya emitting populations of z~3 LBGs using $ugri$ broadband photometric data from the Vera Rubin Observatory Legacy Survey of Space and Time (LSST). △ Less

Submitted 18 November, 2023; originally announced November 2023.

Comments: Published in Publications of the Astronomical Society of Australia 16 pages, 6 figures, 8 tables

Journal ref: Publications of the Astronomical Society of Australia 2023, 40, e052

Minutes-duration Optical Flares with Supernova Luminosities

Authors: Anna Y. Q. Ho, Daniel A. Perley, Ping Chen, Steve Schulze, Vik Dhillon, Harsh Kumar, Aswin Suresh, Vishwajeet Swain, Michael Bremer, Stephen J. Smartt, Joseph P. Anderson, G. C. Anupama, Supachai Awiphan, Sudhanshu Barway, Eric C. Bellm, Sagi Ben-Ami, Varun Bhalerao, Thomas de Boer, Thomas G. Brink, Rick Burruss, Poonam Chandra, Ting-Wan Chen, Wen-Ping Chen, Jeff Cooke, Michael W. Coughlin , et al. (52 additional authors not shown)

Abstract: In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Seve… ▽ More In recent years, certain luminous extragalactic optical transients have been observed to last only a few days. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae) whose timescale is weeks. Some short-duration transients, most notably AT2018cow, display blue optical colours and bright radio and X-ray emission. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source, such as X-ray variability, prolonged ultraviolet emission, a tentative X-ray quasiperiodic oscillation, and large energies coupled to fast (but subrelativistic) radio-emitting ejecta. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the "Tasmanian Devil"). The flares occur over a period of months, are highly energetic, and are likely nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that in some AT2018cow-like transients the embedded energy source is a compact object, either a magnetar or an accreting black hole. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 79 pages, 3 figures (main text) + 7 figures (extended data) + 2 figures (supplementary information). Published online in Nature on 15 November 2023

Mesoscale structure of the atmospheric boundary layer across a natural roughness transition

Authors: Justin P. Cooke, Douglas J. Jerolmack, George I. Park

Abstract: The structure and intensity of turbulence in the atmospheric boundary layer (ABL) drives fluxes of sediment, contaminants, heat, moisture and CO$_2$ at the Earth's surface. Where ABL flows encounter changes in roughness -- such as cities, wind farms, forest canopies and landforms -- a new mesoscopic flow scale is introduced: the internal boundary layer (IBL), which represents a near-bed region of… ▽ More The structure and intensity of turbulence in the atmospheric boundary layer (ABL) drives fluxes of sediment, contaminants, heat, moisture and CO$_2$ at the Earth's surface. Where ABL flows encounter changes in roughness -- such as cities, wind farms, forest canopies and landforms -- a new mesoscopic flow scale is introduced: the internal boundary layer (IBL), which represents a near-bed region of transient flow adjustment that develops over kilometers. This important scale lies within a gap in present observational capabilities of ABL flows, and simplified models fail to capture the sensitive dependence of turbulence on roughness geometry. Here we use large-eddy simulations, run over high-resolution topographic data and validated against field observations, to examine the structure of the ABL across a natural roughness transition: the emergent sand dunes at White Sands National Park. We observe that development of the IBL is triggered by the abrupt transition from smooth playa surface to dunes; however, continuous changes in the size and spacing of dunes over several kilometers influence the downwind patterns of boundary stress and near-bed turbulence. Coherent flow structures grow and merge over the entire $\sim$10-km distance of the dune field, and modulate the influence of large-scale atmospheric turbulence on the bed. Simulated boundary stresses in the developing IBL explain the observed downwind decrease in dune migration, demonstrating a mesoscale coupling between flow and form that governs landscape dynamics. More broadly, our findings demonstrate the importance of resolving both turbulence and realistic roughness for understanding fluid-boundary interactions in environmental flows. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 11 Pages, 4 Figures, 76 References

The Golden Meteorite Fall: Fireball Trajectory, Orbit and Meteorite Characterization

Authors: P. G. Brown, P. J. A. McCausland, A. R Hildebrand, L. T. J. Hanton, L. M. Eckart, H. Busemann, D. Krietsch, C. Maden, K. Welten, M. W. Caffee, M. Laubenstein, D. Vida, F. Ciceri, E. Silber, C. D. K. Herd, P. Hill, H. Devillepoix, Eleanor K. Sansom, Martin Cupák, Seamus Anderson, R. L. Flemming, A. J. Nelson, M. Mazur, D. E. Moser, W. J. Cooke , et al. (4 additional authors not shown)

Abstract: The Golden (British Columbia, Canada) meteorite fall occurred on Oct 4, 2021 at 0534 UT with the first recovered fragment (1.3 kg) landing on an occupied bed. The meteorite is an unbrecciated, low-shock (S2) ordinary chondrite of intermediate composition, typed as an L/LL5. From noble gas measurements the cosmic ray exposure age is 25 Ma while gas retention ages are all >2 Ga. Short-lived radionuc… ▽ More The Golden (British Columbia, Canada) meteorite fall occurred on Oct 4, 2021 at 0534 UT with the first recovered fragment (1.3 kg) landing on an occupied bed. The meteorite is an unbrecciated, low-shock (S2) ordinary chondrite of intermediate composition, typed as an L/LL5. From noble gas measurements the cosmic ray exposure age is 25 Ma while gas retention ages are all >2 Ga. Short-lived radionuclides and noble gas measurements of the pre-atmospheric size overlap with estimates from infrasound and lightcurve modelling producing a preferred pre-atmospheric mass of 70-200 kg. The orbit of Golden has a high inclination (23.5 degs) and is consistent with delivery from the inner main belt. The highest probability (60%) of an origin is from the Hungaria group. We propose that Golden may originate among the background S-type asteroids found interspersed in the Hungaria region. The current collection of 18 L and LL chondrite orbits shows a strong preference for origins in the inner main belt, suggesting multiple parent bodies may be required to explain the diversity in CRE ages and shock states. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 92 Pages, 20 Tables, 21 Figures, plus 3 appendices, accepted in Meteoritics and Planetary Science Oct 26 2023

Ready for O4 II: GRANDMA Observations of Swift GRBs during eight-weeks of Spring 2022

Authors: I. Tosta e Melo, J. -G. Ducoin, Z. Vidadi, C. Andrade, V. Rupchandani, S. Agayeva, J. Abdelhadi, L. Abe, O. Aguerre-Chariol, V. Aivazyan, S. Alishov, S. Antier, J. -M. Bai, A. Baransky, S. Bednarz, Ph. Bendjoya, Z. Benkhaldoun, S. Beradze, M. A. Bizouard, U. Bhardwaj, M. Blazek, M. Boër, E. Broens, O. Burkhonov, N. Christensen , et al. (84 additional authors not shown)

Abstract: We present a campaign designed to train the GRANDMA network and its infrastructure to follow up on transient alerts and detect their early afterglows. In preparation for O4 II campaign, we focused on GRB alerts as they are expected to be an electromagnetic counterpart of gravitational-wave events. Our goal was to improve our response to the alerts and start prompt observations as soon as possible… ▽ More We present a campaign designed to train the GRANDMA network and its infrastructure to follow up on transient alerts and detect their early afterglows. In preparation for O4 II campaign, we focused on GRB alerts as they are expected to be an electromagnetic counterpart of gravitational-wave events. Our goal was to improve our response to the alerts and start prompt observations as soon as possible to better prepare the GRANDMA network for the fourth observational run of LIGO-Virgo-Kagra (which started at the end of May 2023), and future missions such as SM. To receive, manage and send out observational plans to our partner telescopes we set up dedicated infrastructure and a rota of follow-up adcates were organized to guarantee round-the-clock assistance to our telescope teams. To ensure a great number of observations, we focused on Swift GRBs whose localization errors were generally smaller than the GRANDMA telescopes' field of view. This allowed us to bypass the transient identification process and focus on the reaction time and efficiency of the network. During 'Ready for O4 II', 11 Swift/INTEGRAL GRB triggers were selected, nine fields had been observed, and three afterglows were detected (GRB 220403B, GRB 220427A, GRB 220514A), with 17 GRANDMA telescopes and 17 amateur astronomers from the citizen science project Kilonova-Catcher. Here we highlight the GRB 220427A analysis where our long-term follow-up of the host galaxy allowed us to obtain a photometric redshift of $z=0.82\pm0.09$, its lightcurve elution, fit the decay slope of the afterglows, and study the properties of the host galaxy. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Air Traffic Controller Workload Level Prediction using Conformalized Dynamical Graph Learning

Authors: Yutian Pang, Jueming Hu, Christopher S. Lieber, Nancy J. Cooke, Yongming Liu

Abstract: Air traffic control (ATC) is a safety-critical service system that demands constant attention from ground air traffic controllers (ATCos) to maintain daily aviation operations. The workload of the ATCos can have negative effects on operational safety and airspace usage. To avoid overloading and ensure an acceptable workload level for the ATCos, it is important to predict the ATCos' workload accura… ▽ More Air traffic control (ATC) is a safety-critical service system that demands constant attention from ground air traffic controllers (ATCos) to maintain daily aviation operations. The workload of the ATCos can have negative effects on operational safety and airspace usage. To avoid overloading and ensure an acceptable workload level for the ATCos, it is important to predict the ATCos' workload accurately for mitigation actions. In this paper, we first perform a review of research on ATCo workload, mostly from the air traffic perspective. Then, we briefly introduce the setup of the human-in-the-loop (HITL) simulations with retired ATCos, where the air traffic data and workload labels are obtained. The simulations are conducted under three Phoenix approach scenarios while the human ATCos are requested to self-evaluate their workload ratings (i.e., low-1 to high-7). Preliminary data analysis is conducted. Next, we propose a graph-based deep-learning framework with conformal prediction to identify the ATCo workload levels. The number of aircraft under the controller's control varies both spatially and temporally, resulting in dynamically evolving graphs. The experiment results suggest that (a) besides the traffic density feature, the traffic conflict feature contributes to the workload prediction capabilities (i.e., minimum horizontal/vertical separation distance); (b) directly learning from the spatiotemporal graph layout of airspace with graph neural network can achieve higher prediction accuracy, compare to hand-crafted traffic complexity features; (c) conformal prediction is a valuable tool to further boost model prediction accuracy, resulting a range of predicted workload labels. The code used is available at \href{https://github.com/ymlasu/para-atm-collection/blob/master/air-traffic-prediction/ATC-Workload-Prediction/}{$\mathsf{Link}$}. △ Less

Submitted 22 July, 2023; v1 submitted 19 July, 2023; originally announced July 2023.

Numerical and Experimental Study on the Addition of Surface Roughness to Micro-Propellers

Authors: Justin P Cooke, Matthew F Campbell, Edward B Steager, Igor Bargatin, Mark H Yim, George I Park

Abstract: Micro aerial vehicles are making a large impact in applications such as search-and-rescue, package delivery, and recreation. Unfortunately, these diminutive drones are currently constrained to carrying small payloads, in large part because they use propellers optimized for larger aircraft and inviscid flow regimes. Fully realizing the potential of emerging microflyers requires next-generation prop… ▽ More Micro aerial vehicles are making a large impact in applications such as search-and-rescue, package delivery, and recreation. Unfortunately, these diminutive drones are currently constrained to carrying small payloads, in large part because they use propellers optimized for larger aircraft and inviscid flow regimes. Fully realizing the potential of emerging microflyers requires next-generation propellers that are specifically designed for low-Reynolds number conditions and that include new features advantageous in highly viscous flows. One aspect that has received limited attention in the literature is the addition of roughness to propeller blades as a method of reducing drag and increasing thrust. To investigate this possibility, we used large eddy simulation to conduct a numerical investigation of smooth and rough propellers. Our results indicate that roughness produces a 2% increase in thrust and a 5% decrease in power relative to a baseline smooth propeller operating at the same Reynolds number of Rec = 6500, held constant by rotational speed. We corroborated our numerical findings using thrust-stand-based experiments of 3D-printed propellers identical to those of the numerical simulations. Our study confirms that surface roughness is an additional parameter within the design space for micro-propellers that will lead to unprecedented drone efficiencies and payloads. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 23 Pages, 9 Figures

The FRB20190520B Sightline Intersects Foreground Galaxy Clusters

Authors: Khee-Gan Lee, Ilya S. Khrykin, Sunil Simha, Metin Ata, Yuxin Huang, J. Xavier Prochaska, Nicolas Tejos, Jeff Cooke, Kentaro Nagamine, Jielai Zhang

Abstract: The repeating fast radio burst FRB20190520B is an anomaly of the FRB population thanks to its high dispersion measure (DM$=1205\,$pc/cc) despite its low redshift of $z_\mathrm{frb}=0.241$. This excess has been attributed to a large host contribution of $DM_{host}\approx 900\,$pc/cc, far larger than any other known FRB. In this paper, we describe spectroscopic observations of the FRB20190520B field… ▽ More The repeating fast radio burst FRB20190520B is an anomaly of the FRB population thanks to its high dispersion measure (DM$=1205\,$pc/cc) despite its low redshift of $z_\mathrm{frb}=0.241$. This excess has been attributed to a large host contribution of $DM_{host}\approx 900\,$pc/cc, far larger than any other known FRB. In this paper, we describe spectroscopic observations of the FRB20190520B field obtained as part of the FLIMFLAM survey, which yielded 701 galaxy redshifts in the field. We find multiple foreground galaxy groups and clusters, for which we then estimated halo masses by comparing their richness with numerical simulations. We discover two separate $M_{halo} >10^{14}\,M_\odot$ galaxy clusters, at $z=0.1867$ and $z=0.2170$, respectively, that are directly intersected by the FRB sightline within their characteristic halo radius \rvir{}. Subtracting off their estimated DM contributions as well that of the diffuse intergalactic medium, we estimate a host contribution of $DM_{host}=430^{+140}_{-220}\,$pc/cc or $DM_{host}=280^{+140}_{-170}\,$pc/cc (observed frame) depending on whether we assume the halo gas extends to $r_{200}$ or $2\times r_{200}$. This significantly smaller $DM_{host}$ -- no longer the largest known value -- is now consistent with H$α$ emission measures of the host galaxy without invoking unusually high gas temperatures. Combined with the observed FRB scattering timescale, we estimate the turbulent fluctuation and geometric amplification factor of the scattering layer to be $\tilde{F} G\approx4.5 - 11\,(\mathrm{pc^2\;km})^{-1/3}$, suggesting most of the gas is close to the FRB host. This result illustrates the importance of incorporating foreground data for FRB analyses, both for understanding the nature of FRBs and to realize their potential as a cosmological probe. △ Less

Submitted 14 August, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Accepted by ApJL. Interactive figure included (link in text). Note numerical values have changed from v1

Searching for the sources of excess extragalactic dispersion of FRBs

Authors: Sunil Simha, Khee-Gan Lee, J. Xavier Prochaska, Ilya S. Khrykin, Yuxin Huang, Nicolas Tejos, Lachlan Marnoch, Metin Ata, Lucas Bernales, Shivani Bhandari, Jeff Cooke, Adam T. Deller, Suart Ryder, Jielai Zhang

Abstract: The FLIMFLAM survey is collecting spectroscopic data of field galaxies near fast radio burst (FRB) sightlines to constrain key parameters describing the distribution of matter in the Universe. In this work, we leverage the survey data to determine the source of the excess extragalactic dispersion measure (DM), compared to the Macquart relation estimate of four FRBs: FRB20190714A, FRB20200430A, FRB… ▽ More The FLIMFLAM survey is collecting spectroscopic data of field galaxies near fast radio burst (FRB) sightlines to constrain key parameters describing the distribution of matter in the Universe. In this work, we leverage the survey data to determine the source of the excess extragalactic dispersion measure (DM), compared to the Macquart relation estimate of four FRBs: FRB20190714A, FRB20200430A, FRB20200906A, and FRB20210117A. By modeling the gas distribution around the foreground galaxy halos and galaxy groups of the sightlines, we estimate $\rm DM_{halos}$, their contribution to the FRB dispersion measures. The FRB20190714A sightline shows a clear excess of foreground halos which contribute roughly 2/3$^{rd}$ of the observed excess DM, thus implying a baryon-dense sightline. FRB20200906A shows a smaller but non-negligible foreground halo contribution, and further analysis of the IGM is necessary to ascertain the true cosmic contribution to its DM. RB20200430A and FRB20210117A show negligible foreground contributions, implying a large host galaxy excess and/or progenitor environment excess. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 14 pages, 6 figures, 3 tables. Data access available here: https://drive.google.com/drive/folders/1gVVwTAmiHJtFB_GD56hgmDylXHYbOTFS

The restframe ultraviolet of superluminous supernovae -- I. Potential as cosmological probes

Authors: Nandita Khetan, Jeff Cooke, Marica Branchesi

Abstract: Superluminous supernovae (SLSNe) have been detected to $z\sim4$ and can be detected to $z\gtrsim15$ using current and upcoming facilities. SLSNe are extremely UV luminous, and hence objects at $z\gtrsim7$ are detected exclusively via their rest-frame UV using optical and infrared facilities. SLSNe have great utility in multiple areas of stellar and galactic evolution. Here, we explore the potentia… ▽ More Superluminous supernovae (SLSNe) have been detected to $z\sim4$ and can be detected to $z\gtrsim15$ using current and upcoming facilities. SLSNe are extremely UV luminous, and hence objects at $z\gtrsim7$ are detected exclusively via their rest-frame UV using optical and infrared facilities. SLSNe have great utility in multiple areas of stellar and galactic evolution. Here, we explore the potential use of SLSNe type-I as high-redshift cosmological distance indicators in their rest-frame UV. Using a SLSNe-I sample in the redshift range $1\lesssim z\lesssim 3$, we investigate correlations between the peak absolute magnitude in a synthetic UV filter centered at 250 nm and rise time, colour and decline rate of SLSNe-I light curves. We observe a linear correlation between $M_0(250)$ and the rise time with an intrinsic scatter of 0.29. Interestingly, this correlation is further tightened ($σ_{int} \approx 0.2$) by eliminating those SLSNe which show a pre-peak bump in their light curve. This result hints at the possibility that the "bumpy" SLSNe could belong to a different population. Weak correlations are observed between the peak luminosity and colour indices. No relationship is found between UV peak magnitude and the decline rate in contrast to what is typically found in optical band. The correlations found here are promising, and give encouraging insights for the use of SLSNe as cosmological probes at high redshifts using standardising relations in the UV. We also highlight the importance of early, and consistent, photometric data for constraining the light curve properties. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS

Direct measurement of decimeter-sized rocky material in the Oort cloud

Authors: Denis Vida, Peter G. Brown, Hadrien A. R. Devillepoix, Paul Wiegert, Danielle E. Moser, Pavol Matlovič, Christopher D. K. Herd, Patrick J. A. Hill, Eleanor K. Sansom, Martin C. Towner, Juraj Tóth, William J. Cooke, Donald W. Hladiuk

Abstract: The Oort cloud is thought to be a reservoir of icy planetesimals and the source of long-period comets (LPCs) implanted from the outer Solar System during the time of giant planet formation. The abundance of rocky ice-free bodies is a key diagnostic of Solar System formation models as it can distinguish between ``massive" and ``depleted" proto-asteroid belt scenarios and thus disentangle competing… ▽ More The Oort cloud is thought to be a reservoir of icy planetesimals and the source of long-period comets (LPCs) implanted from the outer Solar System during the time of giant planet formation. The abundance of rocky ice-free bodies is a key diagnostic of Solar System formation models as it can distinguish between ``massive" and ``depleted" proto-asteroid belt scenarios and thus disentangle competing planet formation models. Here we report a direct observation of a decimeter-sized ($\sim2$ kg) rocky meteoroid on a retrograde LPC orbit ($e \approx 1.0$, i = $121^{\circ}$). During its flight, it fragmented at dynamic pressures similar to fireballs dropping ordinary chondrite meteorites. A numerical ablation model fit produces bulk density and ablation properties also consistent with asteroidal meteoroids. We estimate the flux of rocky objects impacting Earth from the Oort cloud to be $1.08^{+2.81}_{-0.95} \mathrm{meteoroids/10^6 km^2/yr}$ to a mass limit of 10 g. This corresponds to an abundance of rocky meteoroids of $\sim6^{+13}_{-5}$\% of all objects originating in the Oort cloud and impacting Earth to these masses. Our result gives support to migration-based dynamical models of the formation of the Solar System which predict that significant rocky material is implanted in the Oort cloud, a result not explained by traditional Solar System formation models. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Nature Astronomy

Enhancing the Electron Mobility in Si-doped (010) $β$-Ga$_2$O$_3$ films with Low-Temperature Buffer Layers

Authors: Arkka Bhattacharyya, Carl Peterson, Takeki Itoh, Saurav Roy, Jacqueline Cooke, Steve Rebollo, Praneeth Ranga, Berardi Sensale-Rodriguez, Sriram Krishnamoorthy

Abstract: We demonstrate a new substrate cleaning and buffer growth scheme in $β$-Ga$_2$O$_3$ epitaxial thin films using metalorganic vapor phase epitaxy (MOVPE). For the channel structure, a low-temperature (LT, 600 $^\circ$C) undoped Ga$_2$O$_3$ buffer is grown followed by transition layers to a high-temperature (HT, 810 $^\circ$C) Si-doped Ga$_2$O$_3$ channel layers without growth interruption. The (010)… ▽ More We demonstrate a new substrate cleaning and buffer growth scheme in $β$-Ga$_2$O$_3$ epitaxial thin films using metalorganic vapor phase epitaxy (MOVPE). For the channel structure, a low-temperature (LT, 600 $^\circ$C) undoped Ga$_2$O$_3$ buffer is grown followed by transition layers to a high-temperature (HT, 810 $^\circ$C) Si-doped Ga$_2$O$_3$ channel layers without growth interruption. The (010) Ga$_2$O$_3$ Fe-doped substrate cleaning uses solvent cleaning followed by an additional HF (49% in water) treatment for 30 mins before the epilayer growth. This step is shown to compensate the parasitic Si channel at the epilayer-substrate interface that originates from the substrate polishing process or contamination from the ambient. SIMS analysis shows the Si peak atomic density at the substrate interface is several times lower than the Fe atomic density in the substrate - indicating full compensation. The elimination of the parasitic electron channel at the epi-substrate interface was also verified by electrical (capacitance-voltage profiling) measurements. In the LT-grown buffer layers, it is seen that the Fe forward decay tail from the substrate is very sharp with a decay rate of $\sim$ 9 nm$/$dec. These channels show record high electron mobility in the range of 196 - 85 cm$^2$/Vs in unintentionally doped and Si-doped films in the doping range of 2$\times$10$^{16}$ to 1$\times$10$^{20}$ cm$^{-3}$. Si delta-doped channels were also grown utilizing this substrate cleaning and the hybrid LT-buffers. Record high electron Hall mobility of 110 cm$^2$/Vs was measured for sheet charge density of 9.2$\times$10$^{12}$ cm$^{-2}$. This substrate cleaning combined with the LT-buffer scheme shows the potential of designing Si-doped $β$-Ga$_2$O$_3$ channels with exceptional transport properties for high performance gallium oxide-based electron devices. △ Less

Submitted 2 February, 2023; v1 submitted 5 December, 2022; originally announced December 2022.

Comments: 20 pages, 10 figures,

A very luminous jet from the disruption of a star by a massive black hole

Authors: Igor Andreoni, Michael W. Coughlin, Daniel A. Perley, Yuhan Yao, Wenbin Lu, S. Bradley Cenko, Harsh Kumar, Shreya Anand, Anna Y. Q. Ho, Mansi M. Kasliwal, Antonio de Ugarte Postigo, Ana Sagues-Carracedo, Steve Schulze, D. Alexander Kann, S. R. Kulkarni, Jesper Sollerman, Nial Tanvir, Armin Rest, Luca Izzo, Jean J. Somalwar, David L. Kaplan, Tomas Ahumada, G. C. Anupama, Katie Auchettl, Sudhanshu Barway , et al. (56 additional authors not shown)

Abstract: Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jett… ▽ More Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Nature

Radio Transients and Variables in the Tenth Deeper, Wider, Faster Observing Run

Authors: D. Dobie, J. Pritchard, Y. Wang, L. W. Graham, J. Freeburn, H. Qiu, T. R. White, A. O'Brien, E. Lenc, J. K. Leung, C. Lynch, Tara Murphy, A. J. Stewart, Z. Wang, A. Zic, T. M. C. Abbott, C. Cai, J. Cooke, M. Dobiecki, S. Goode, S. Jia, C. Li, A. Möller, S. Webb, J. Zhang , et al. (1 additional authors not shown)

Abstract: The Deeper, Wider, Faster (DWF) program coordinates observations with telescopes across the electromagnetic spectrum, searching for transients on timescales of milliseconds to days. The tenth DWF observing run was carried out in near real-time during September 2021 and consisted of six consecutive days of observations of the NGC 6744 galaxy group and a field containing the repeating fast radio bur… ▽ More The Deeper, Wider, Faster (DWF) program coordinates observations with telescopes across the electromagnetic spectrum, searching for transients on timescales of milliseconds to days. The tenth DWF observing run was carried out in near real-time during September 2021 and consisted of six consecutive days of observations of the NGC 6744 galaxy group and a field containing the repeating fast radio burst FRB190711 with the Australian Square Kilometre Array Pathfinder, the Dark Energy Camera, the Hard X-ray Modulation Telescope and the Parkes 64m "Murriyang" radio telescope. In this work we present the results of an image-domain search for transient, variable and circularly polarised sources carried out with ASKAP using data from the observing run, along with test observations prior to the run and follow-up observations carried out during and after the run. We identified eight variable radio sources, consisting of one pulsar, six stellar systems (five of which exhibit circularly polarised emission) and one previously uncatalogued source. Of particular interest is the detection of pulses from the ultra-cool dwarf SCR J1845-6357 with a period of $14.2\pm 0.3$ h, in good agreement with the known optical rotation period, making this the slowest rotating radio-loud ultra-cool dwarf discovered. △ Less

Submitted 13 November, 2022; originally announced November 2022.

MUSE Analysis of Gas around Galaxies (MAGG) -- IV: The gaseous environment of $z\sim$ 3-4 Lyman-alpha emitting galaxies

Authors: Emma K. Lofthouse, Michele Fumagalli, Matteo Fossati, Rajeshwari Dutta, Marta Galbiati, Fabrizio Arrigoni Battaia, Sebastiano Cantalupo, Lise Christensen, Ryan J. Cooke, Alessia Longobardi, Michael T. Murphy, J. Xavier. Prochaska

Abstract: We study the link between galaxies and HI-selected absorption systems at z~3-4 in the MUSE Analysis of Gas around Galaxies (MAGG) survey, an ESO large programme consisting of integral field pectroscopic observations of 28 quasar fields hosting 61 strong absorbers with $\rm N_{\rm HI}\gtrsim 10^{16.5}~\rm cm^{-2}$. We identify 127 Ly$α$ emitting galaxies (LAEs) around the absorbers, corresponding t… ▽ More We study the link between galaxies and HI-selected absorption systems at z~3-4 in the MUSE Analysis of Gas around Galaxies (MAGG) survey, an ESO large programme consisting of integral field pectroscopic observations of 28 quasar fields hosting 61 strong absorbers with $\rm N_{\rm HI}\gtrsim 10^{16.5}~\rm cm^{-2}$. We identify 127 Ly$α$ emitting galaxies (LAEs) around the absorbers, corresponding to a detection rate of 82$\pm$16 per cent. The luminosity function of these LAEs is approximately 5 times higher in normalization than the field population and we detect a significant clustering of galaxies with respect to the gas, confirming that high column density absorbers and LAEs trace each other. Between 30 and 40 per cent of the absorbers are associated with multiple LAEs, which lie preferentially along filaments. Galaxies in groups also exhibit a three times higher covering factor of optically-thick gas compared to isolated systems. No significant correlations are identified between the emission properties of LAEs and the absorption properties of optically-thick gas clouds, except for a weak preference of brighter and multiple galaxies to reside near broad absorbers. Based on the measured impact parameters and the covering factor, we conclude that the near totality of optically-thick gas in the Universe can be found in the outer circumgalactic medium (CGM) of LAEs or in the intergalactic medium (IGM) in proximity to these galaxies. Thus, LAEs act as tracers of larger scale structures within which both galaxies and optically-thick clouds are embedded. The patchy and inhomogeneous nature of the CGM and IGM explains the lack of correlations between absorption and emission properties. This implies that very large samples are needed to unveil the trends that encode the properties of the baryon cycle. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 30 pages, 19 figures, submitted to MNRAS

Science Cases for the Keck Wide-Field Imager

Authors: J. Cooke, C. Angus, K. Auchettl, J. Bally, B. Bolin, S. Brough, J. N. Burchett, R. Foley, G. Foran, D. Forbes, J. Gannon, R. Hirai, G. G. Kacprzak, R. Margutti, C. Martinez-Lombilla, U. Mestric, A. Moller, A. Rest, J. Rhodes, R. M. Rich, F. Schussler, R. Wainscoat, J. Walawender, I. Wold, J. Zhang

Abstract: The Keck Wide-Field Imager (KWFI) is a proposed 1-degree diameter field of view UV-sensitive optical camera for Keck prime focus. KWFI will be the most powerful optical wide-field camera in the world and the only such 8m-class camera sensitive down to ~3000 A for the foreseeable future. Twenty science cases are described for KWFI compiled largely during 2019-2021, preceded by a brief discussion of… ▽ More The Keck Wide-Field Imager (KWFI) is a proposed 1-degree diameter field of view UV-sensitive optical camera for Keck prime focus. KWFI will be the most powerful optical wide-field camera in the world and the only such 8m-class camera sensitive down to ~3000 A for the foreseeable future. Twenty science cases are described for KWFI compiled largely during 2019-2021, preceded by a brief discussion of the instrument, components, and capabilities for context. △ Less

Submitted 24 July, 2022; originally announced July 2022.

Comments: 56 pages, 31 figures

Sonification and Sound Design for Astronomy Research, Education and Public Engagement

Authors: A. Zanella, C. M. Harrison, S. Lenzi, J. Cooke, P. Damsma, S. W. Fleming

Abstract: Over the last ten years there has been a large increase in the number of projects using sound to represent astronomical data and concepts. Motivation for these projects includes the potential to enhance scientific discovery within complex datasets, by utilising the inherent multi-dimensionality of sound and the ability of our hearing to filter signals from noise. Other motivations include creating… ▽ More Over the last ten years there has been a large increase in the number of projects using sound to represent astronomical data and concepts. Motivation for these projects includes the potential to enhance scientific discovery within complex datasets, by utilising the inherent multi-dimensionality of sound and the ability of our hearing to filter signals from noise. Other motivations include creating engaging multi-sensory resources, for education and public engagement, and making astronomy more accessible to people who are blind or have low vision, promoting their participation in science and related careers. We describe potential benefits of sound within these contexts and provide an overview of the nearly 100 sound-based astronomy projects that we identified. We discuss current limitations and challenges of the approaches taken. Finally, we suggest future directions to help realise the full potential of sound-based techniques in general and to widen their application within the astronomy community. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Nature Astronomy. A Word document (more accessible with screen readers) is available under 'ancillary files'. This is the author's own version (it is not the Version of Record and does not reflect post-acceptance improvements, or any corrections). The Version of Record will be available with doi: 10.1038/s41550-022-01721-z

Resolving the HI in Damped Lyman-α systems that power star-formation

Authors: Rongmon Bordoloi, John M. O'Meara, Keren Sharon, Jane R. Rigby, Jeff Cooke, Ahmed Shaban, Mateusz Matuszewski, Luca Rizzi, Greg Doppmann, D. Christopher Martin, Anna M. Moore, Patrick Morrissey, James D. Neill

Abstract: Reservoirs of dense atomic gas (primarily hydrogen), contain approximately 90 percent of the neutral gas at a redshift of 3, and contribute to 2-3 percent of the total baryons in the Universe. These damped Lyman-$α$ systems (so called because they absorb Lyman-$α$ photons from within and from background sources) have been studied for decades, but only through absorption lines present in the spectr… ▽ More Reservoirs of dense atomic gas (primarily hydrogen), contain approximately 90 percent of the neutral gas at a redshift of 3, and contribute to 2-3 percent of the total baryons in the Universe. These damped Lyman-$α$ systems (so called because they absorb Lyman-$α$ photons from within and from background sources) have been studied for decades, but only through absorption lines present in the spectra of background quasars and gamma-ray bursts. Such pencil beams do not constrain the physical extent of the systems. Here, we report integral-field spectroscopy of a bright, gravitationally lensed galaxy at a redshift of 2.7 with two foreground damped Lyman-$α$ systems. These systems are $>$ 238 $kpc^2$ in extent, with column densities of neutral hydrogen varying by more than an order of magnitude on $<$ 3 kpc-scales. The mean column densities are $10^{20.46}$ - $10^{20.84} cm^{-2}$ and the total masses are $> 5.5 \times 10^{8}$ - $1.4 \times 10^{9} M_{\odot}$, showing that they contain the necessary fuel for the next generation of star formation, consistent with relatively massive, low-luminosity primeval galaxies at redshifts $>$ 2. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Comments: 31 pages, 4 Figures, 4 Extended Data figures, 2 Extended Data tables, Author's version, Accepted: 4 March 2022, Published online by Nature on May 18, 2022

Journal ref: Nature 606, 59-63 (2022)

Higher Rank Askey-Wilson Algebras as Skein Algebras

Authors: Juliet Cooke, Abel Lacabanne

Abstract: In this paper we give a topological interpretation and diagrammatic calculus for the rank $(n-2)$ Askey-Wilson algebra by proving there is an explicit isomorphism with the Kauffman bracket skein algebra of the $(n+1)$-punctured sphere. To do this we consider the Askey-Wilson algebra in the braided tensor product of $n$ copies of either the quantum group $\mathcal{U}_q{(\mathfrak{sl}_2)}$ or the re… ▽ More In this paper we give a topological interpretation and diagrammatic calculus for the rank $(n-2)$ Askey-Wilson algebra by proving there is an explicit isomorphism with the Kauffman bracket skein algebra of the $(n+1)$-punctured sphere. To do this we consider the Askey-Wilson algebra in the braided tensor product of $n$ copies of either the quantum group $\mathcal{U}_q{(\mathfrak{sl}_2)}$ or the reflection equation algebra. We then use the isomorpism of the Kauffman bracket skein algebra of the $(n+1)$-punctured sphere with the $\mathcal{U}_q{(\mathfrak{sl}_2})$ invariants of the Aleeksev moduli algebra to complete the correspondence. We also find the graded vector space dimension of the $\mathcal{U}_q{(\mathfrak{sl}_2})$ invariants of the Aleeksev moduli algebra and apply this to finding a presentation of the skein algebra of the five-punctured sphere and hence also find a presentation for the rank $2$ Askey-Wilson algebra. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 40 pages

Primordial helium-3 redux: The helium isotope ratio of the Orion nebula

Authors: Ryan J. Cooke, Pasquier Noterdaeme, James W. Johnson, Max Pettini, Louise Welsh, Celine Peroux, Michael T. Murphy, David H. Weinberg

Abstract: We report the first direct measurement of the helium isotope ratio, 3He/4He, outside of the Local Interstellar Cloud, as part of science verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph (CRIRES). Our determination of 3He/4He is based on metastable HeI* absorption along the line-of-sight towards Tet02 Ori A in the Orion Nebula. We measure a value 3He/4He=(1.77+/-0… ▽ More We report the first direct measurement of the helium isotope ratio, 3He/4He, outside of the Local Interstellar Cloud, as part of science verification observations with the upgraded CRyogenic InfraRed Echelle Spectrograph (CRIRES). Our determination of 3He/4He is based on metastable HeI* absorption along the line-of-sight towards Tet02 Ori A in the Orion Nebula. We measure a value 3He/4He=(1.77+/-0.13)x10^{-4}, which is just ~40 per cent above the primordial relative abundance of these isotopes, assuming the Standard Model of particle physics and cosmology, (3He/4He)_p = (1.257+/-0.017)x10^-4. We calculate a suite of galactic chemical evolution simulations to study the Galactic build up of these isotopes, using the yields from Limongi & Chieffi (2018) for stars in the mass range M=8-100 M_sun and Lagarde (2011,2012) for M=0.8-8 M_sun. We find that these simulations simultaneously reproduce the Orion and protosolar 3He/4He values if the calculations are initialized with a primordial ratio (3He/4He)_p=(1.043+/-0.089)x10^-4. Even though the quoted error does not include the model uncertainty, this determination agrees with the Standard Model value to within ~2sigma. We also use the present-day Galactic abundance of deuterium (D/H), helium (He/H), and 3He/4He to infer an empirical limit on the primordial 3He abundance, (3He/H)_p < (1.09+/-0.18)x10^-5, which also agrees with the Standard Model value. We point out that it is becoming increasingly difficult to explain the discrepant primordial 7Li/H abundance with non-standard physics, without breaking the remarkable simultaneous agreement of three primordial element ratios (D/H, 4He/H, and 3He/4He) with the Standard Model values. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 23 pages, 9 figures, Resubmitted to the Astrophysical Journal after addressing referee comments

A Cautionary Tale of LyC Escape Fraction Estimates from High Redshift Galaxies

Authors: R. Bassett, E. V. Ryan-Weber, J. Cooke, U. Mestric, L. J. Prichard, M. Rafelski, I. Iwata, M. Sawicki, S. Gwyn, S. Arnouts

Abstract: Measuring the escape fraction, $f_{\rm esc}$, of ionizing, Lyman Continuum (LyC) radiation is key to our understanding of the process of cosmic reionization. In this paper we provide a methodology for recovering the posterior probability distribution of the LyC escape fraction, $f_{\rm esc}^{\rm PDF}$, considering both the observational uncertainties and ensembles of simulated transmission functio… ▽ More Measuring the escape fraction, $f_{\rm esc}$, of ionizing, Lyman Continuum (LyC) radiation is key to our understanding of the process of cosmic reionization. In this paper we provide a methodology for recovering the posterior probability distribution of the LyC escape fraction, $f_{\rm esc}^{\rm PDF}$, considering both the observational uncertainties and ensembles of simulated transmission functions through the intergalactic medium (IGM). We present an example of this method applied to a VUDS galaxy at $z=3.64$ and find $f_{\rm esc}^{\rm PDF}$ = 0.51$^{+0.33}_{-0.34}$ and compare this to the values computed assuming averaged IGM transmission with and without consideration of detection bias along average sightlines yielding $f_{\rm esc}^{\langle T \rangle}$ = 1.40$^{+0.80}_{-0.42}$, and $f_{\rm esc}^{\rm bias}$ = 0.82$^{+0.33}_{-0.16}$. Our results highlight the limitations of methods assuming average, smooth transmission functions. We also present MOSFIRE data for a sample of seven LyC candidates selected based on photometric redshifts at $z > 3.4$, but find that all seven have overestimated photometric redshifts by $Δz \sim 0.2$ making them unsuitable for LyC measurements. This results likely due to a bias induced by our selection criteria. △ Less

Submitted 6 February, 2022; originally announced February 2022.

Comments: 12 pages, 9 figues, 2 tables, accepted for publication in MNRAS

Lyman Continuum Galaxy Candidates in COSMOS

Authors: Laura J. Prichard, Marc Rafelski, Jeff Cooke, Uros Mestric, Robert Bassett, Emma V. Ryan-Weber, Ben Sunnquist, Anahita Alavi, Nimish Hathi, Xin Wang, Mitchell Revalski, Varun Bajaj, John M. O'Meara, Lee Spitler

Abstract: Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$,… ▽ More Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$, $\sim29\text{--}30$~mag depth at 5$σ$), we attempt to identify a sample of Lyman continuum galaxies (LCGs). These are individual sources that emit ionizing flux below the Lyman break ($<912~\textÅ$). This population would allow us to constrain cosmic reionization parameters such as the number density and escape fraction ($f_{\rm esc}$) of ionizing sources. We compile a comprehensive parent sample that does not rely on the Lyman-break technique for redshifts. We present three new spectroscopic candidates at $z\sim3.7\text{--}4.4$, and 32 new photometric candidates. The high-resolution multi-band HST imaging and new Keck/Low Resolution Imaging Spectrometer (LRIS) redshifts make these promising spectroscopic LCG candidates. Using both a traditional and probabilistic approach, we find the most likely $f_{\rm esc}$ values for the three spectroscopic LCG candidates are $>100\%$, and therefore not physical. We are unable to confirm the true nature of these sources with the best available imaging and direct blue Keck/LRIS spectroscopy. More spectra, especially from the new class of 30 m telescopes, will be required to build a statistical sample of LCGs to place firm observational constraints on cosmic reionization. △ Less

Submitted 14 February, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ. 38 pages (28 in main body), 14 figures

SOAR/Goodman Spectroscopic Assessment of Candidate Counterparts of the LIGO-Virgo Event GW190814

Authors: Douglas Tucker, Matthew Wiesner, Sahar Allam, Marcelle Soares-Santos, Clecio de Bom, Melissa Butner, Alyssa Garcia, Robert Morgan, Felipe Olivares, Antonella Palmese, Luidhy Santana-Silva, Anushka Shrivastava, James Annis, Juan Garcia-Bellido, Mandeep Gill, Kenneth Herner, Charles Kilpatrick, Martin Makler, Nora Sherman, Adam Amara, Huan Lin, Mathew Smith, Elizabeth Swann, Iair Arcavi, Tristan Bachmann , et al. (118 additional authors not shown)

Abstract: On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity in… ▽ More On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond. △ Less

Submitted 2 June, 2022; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: 32 pages, 18 figures, accepted for publication by ApJ

Report number: DES-2020-601, FERMILAB-PUB-21-454-AE-E-SCD

Journal ref: ApJ, 929, 115 (2022)

Finding Fast Transients in Real Time Using Novel Light Curve Analysis Algorithm

Authors: Robert Strausbaugh, Antonino Cucchiara, Michael Dow Jr., Sara Webb, Jielai Zhang, Simon Goode, Jeff Cooke

Abstract: The current data acquisition rate of astronomical transient surveys and the promise for significantly higher rates during in the next decade necessitate the development of novel approaches to analyze astronomical data sets and promptly detect objects of interest. The Deeper, Wider, Faster (DWF) program is a survey focused on the identification of fast evolving transients, such as fast radio bursts… ▽ More The current data acquisition rate of astronomical transient surveys and the promise for significantly higher rates during in the next decade necessitate the development of novel approaches to analyze astronomical data sets and promptly detect objects of interest. The Deeper, Wider, Faster (DWF) program is a survey focused on the identification of fast evolving transients, such as fast radio bursts, gamma-ray bursts, and supernova shock breakouts. It employs a multi-frequency simultaneous coverage of the same part of the sky over several orders of magnitude. Using the Dark Energy Camera mounted on the 4-meter Blanco telescope, DWF captures a 20 second g-band exposure every minute, at a typical seeing of ~ 1" and an airmass of ~ 1.5. These optical data are collected simultaneously with observations conducted over the entire electromagnetic spectrum - from radio to gamma-rays - as well as cosmic ray observations. In this paper, we present a novel real-time light curve analysis algorithm, designed to detect transients in the DWF optical data; this algorithm functions independently from, or in conjunction with, image subtraction. We present a sample of fast transients detected by our algorithm, as well as a false-positive analysis. Our algorithm is customizable and can be tuned to be sensitive to transients evolving over different timescales and flux ranges. △ Less

Submitted 6 December, 2021; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: AJ accepted, 19 pages, 11 figures

Upper limits on the escape fraction of ionizing radiation from galaxies at $2\lesssim z < 6$

Authors: U. Meštrić, E. V. Ryan-Weber, J. Cooke, R. Bassett, L. J. Prichard, M. Rafelski

Abstract: In this work, we investigate upper limits on the global escape fraction of ionizing photons ($f_{\rm esc/global}^{\rm abs}$) from a sample of galaxies probed for Lyman-continuum (LyC) emission characterized as non-LyC and LyC leakers. We present a sample of 9 clean non-contaminated (by low redshift interlopers, CCD problems and internal reflections of the instrument) galaxies which do not show sig… ▽ More In this work, we investigate upper limits on the global escape fraction of ionizing photons ($f_{\rm esc/global}^{\rm abs}$) from a sample of galaxies probed for Lyman-continuum (LyC) emission characterized as non-LyC and LyC leakers. We present a sample of 9 clean non-contaminated (by low redshift interlopers, CCD problems and internal reflections of the instrument) galaxies which do not show significant ($>$ $3σ$) LyC flux between 880Å $<λ_{rest}<$ 910Å. The 9 galaxy stacked spectrum reveals no significant LyC flux with an upper limit of $f_{\rm esc}^{\rm abs} \leq 0.06$. In the next step of our analysis, we join all estimates of $f_{\rm esc}^{\rm abs}$ upper limits derived from different samples of $2\lesssim z < 6$ galaxies from the literature reported in last $\sim$20 years and include the sample presented in this work. We find the $f_{\rm esc}^{\rm abs}$ upper limit $\leq$ 0.084 for the galaxies recognized as non-LyC leakers. After including all known detections from literature $f_{\rm esc/global}^{\rm abs}$ upper limit $\leq$ 0.088 for all galaxies examined for LyC flux. Furthermore, $f_{\rm esc}^{\rm abs}$ upper limits for different groups of galaxies indicate that the strongest LyC emitters could be galaxies classified as Lyman alpha emitters. We also discuss the possible existence of a correlation among the observed flux density ratio $(F_ν^{LyC}/F_ν^{UV})_{\rm obs}$ and Lyman alpha equivalent width EW(Ly$α)$, where we confirm the existence of moderately significant correlation among galaxies classified as non-LyC leakers. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: 17 pages, 5 Figures

A comprehensive search for the radio counterpart of GW190814 with the Australian Square Kilometre Array Pathfinder

Authors: D. Dobie, A. Stewart, K. Hotokezaka, Tara Murphy, D. L. Kaplan, D. A. H. Buckley, J. Cooke, A. Y. Q. Ho, E. Lenc, J. K. Leung, M. Gromadzki, A. O'Brien, S. Pintaldi, J. Pritchard, Y. Wang, Z. Wang

Abstract: We present results from a search for the radio counterpart to the possible neutron star-black hole merger GW190814 with the Australian Square Kilometre Array Pathfinder. We have carried out 10 epochs of observation spanning 2-655 days post-merger at a frequency of 944 MHz. Each observation covered 30 deg$^2$, equivalent to 87% of the event localisation. We conducted an untargeted search for radio… ▽ More We present results from a search for the radio counterpart to the possible neutron star-black hole merger GW190814 with the Australian Square Kilometre Array Pathfinder. We have carried out 10 epochs of observation spanning 2-655 days post-merger at a frequency of 944 MHz. Each observation covered 30 deg$^2$, equivalent to 87% of the event localisation. We conducted an untargeted search for radio transients in the field, as well as a targeted search for transients associated with known galaxies. We find one radio transient, ASKAP J005022.3-230349, but conclude that it is unlikely to be associated with the merger. We use our observations to place constraints on the inclination angle of the merger and the density of the surrounding environment by comparing our non-detection to model predictions for radio emission from compact binary coalescences. This survey is also the most comprehensive widefield search (in terms of sensitivity and both areal and temporal coverage) for radio transients to-date and we calculate the radio transient surface density at 944 MHz. △ Less

Submitted 26 September, 2021; v1 submitted 17 September, 2021; originally announced September 2021.

Constraining the Cosmic Baryon Distribution with Fast Radio Burst Foreground Mapping

Authors: Khee-Gan Lee, Metin Ata, Ilya S. Khrykin, Yuxin Huang, J. Xavier Prochaska, Jeff Cooke, Jielai Zhang, Adam Batten

Abstract: The dispersion measures (DM) of fast radio bursts (FRBs) encode the integrated electron density along the line-of-sight, which is typically dominated by the intergalactic medium (IGM) contribution in the case of extragalactic FRBs. In this paper, we show that incorporating wide-field spectroscopic galaxy survey data in the foreground of localized FRBs can significantly improve constraints on the p… ▽ More The dispersion measures (DM) of fast radio bursts (FRBs) encode the integrated electron density along the line-of-sight, which is typically dominated by the intergalactic medium (IGM) contribution in the case of extragalactic FRBs. In this paper, we show that incorporating wide-field spectroscopic galaxy survey data in the foreground of localized FRBs can significantly improve constraints on the partition of diffuse cosmic baryons. Using mock DMs and realistic lightcone galaxy catalogs derived from the Millennium simulation, we define spectroscopic surveys that can be carried out with 4m and 8m-class wide field spectroscopic facilities. On these simulated surveys, we carry out Bayesian density reconstructions in order to estimate the foreground matter density field. In comparison with the `true' matter density field, we show that these can help reduce the uncertainties in the foreground structures by $\sim 2-3\times$ compared to cosmic variance. We calculate the Fisher matrix to forecast that $N=30\: (96)$ localized FRBs should be able to constrain the diffuse cosmic baryon fraction to $\sim 10\%\: (\sim 5\%) $, and parameters governing the size and baryon fraction of galaxy circumgalactic halos to within $\sim 20-25\%\: (\sim 8-12\%)$. From the Fisher analysis, we show that the foreground data increases the sensitivity of localized FRBs toward our parameters of interest by $\sim 25\times$. We briefly introduce FLIMFLAM, an ongoing galaxy redshift survey that aims to obtain foreground data on $\sim 30$ localized FRB fields. △ Less

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

Comments: 27 pages, 12 Figures. Accepted by ApJ

The Deeper, Wider, Faster Program: Exploring stellar flare activity with deep, fast cadenced DECam imaging via machine learning

Authors: Sara Webb, Chris Flynn, Jeff Cooke, Jielai Zhang, Ashish Mahabal, Tim Abbott, Rebecca Allen, Igor Andreoni, Sarah Bird, Simon Goode, Michelle Lochner, Tyler Pritchard

Abstract: We present our 500 pc distance-limited study of stellar fares using the Dark Energy Camera as part of the Deeper, Wider, Faster Program. The data was collected via continuous 20-second cadence g band imaging and we identify 19,914 sources with precise distances from Gaia DR2 within twelve, ~3 square-degree, fields over a range of Galactic latitudes. An average of ~74 minutes is spent on each field… ▽ More We present our 500 pc distance-limited study of stellar fares using the Dark Energy Camera as part of the Deeper, Wider, Faster Program. The data was collected via continuous 20-second cadence g band imaging and we identify 19,914 sources with precise distances from Gaia DR2 within twelve, ~3 square-degree, fields over a range of Galactic latitudes. An average of ~74 minutes is spent on each field per visit. All light curves were accessed through a novel unsupervised machine learning technique designed for anomaly detection. We identify 96 flare events occurring across 80 stars, the majority of which are M dwarfs. Integrated are energies range from $\sim 10^{31}-10^{37}$ erg, with a proportional relationship existing between increased are energy with increased distance from the Galactic plane, representative of stellar age leading to declining yet more energetic are events. In agreement with previous studies we observe an increase in flaring fraction from M0 -> M6 spectral types. Furthermore, we find a decrease in the flaring fraction of stars as vertical distance from the galactic plane is increased, with a steep decline present around ~100 pc. We find that ~70% of identified flares occur on short timescales of ~8 minutes. Finally we present our associated are rates, finding a volumetric rate of $2.9 \pm 0.3 \times 10^{-6}$ flares pc$^{-3}$ hr$^{-1}$. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 20 pages, 14 figures, 6 tables

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

Authors: Charles D. Kilpatrick, David A. Coulter, Iair Arcavi, Thomas G. Brink, Georgios Dimitriadis, Alexei V. Filippenko, Ryan J. Foley, D. Andrew Howell, David O. Jones, Martin Makler, Anthony L. Piro, César Rojas-Bravo, David J. Sand, Jonathan J. Swift, Douglas Tucker, WeiKang Zheng, Sahar S. Allam, James T. Annis, Juanita Antilen, Tristan G. Bachmann, Joshua S. Bloom, Clecio R. Bom, K. Azalee Bostroem, Dillon Brout, Jamison Burke , et al. (57 additional authors not shown)

Abstract: We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localiz… ▽ More We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources. △ Less

Submitted 12 June, 2021; originally announced June 2021.

Comments: 86 pages, 9 figures

MUSE Analysis of Gas around Galaxies (MAGG) -- III: The gas and galaxy environment of z = 3-4.5 quasars

Authors: Matteo Fossati, Michele Fumagalli, Emma K. Lofthouse, Rajeshwari Dutta, Sebastiano Cantalupo, Fabrizio Arrigoni Battaia, Johan P. U. Fynbo, Elisabeta Lusso, Michael T. Murphy, J. Xavier Prochaska, Tom Theuns, Ryan J. Cooke

Abstract: We present a study of the environment of 27 z=3-4.5 bright quasars from the MUSE Analysis of Gas around Galaxies (MAGG) survey. With medium-depth MUSE observations (4 hours on target per field), we characterise the effects of quasars on their surroundings by studying simultaneously the properties of extended gas nebulae and Lyalpha emitters (LAEs) in the quasar host haloes. We detect extended (up… ▽ More We present a study of the environment of 27 z=3-4.5 bright quasars from the MUSE Analysis of Gas around Galaxies (MAGG) survey. With medium-depth MUSE observations (4 hours on target per field), we characterise the effects of quasars on their surroundings by studying simultaneously the properties of extended gas nebulae and Lyalpha emitters (LAEs) in the quasar host haloes. We detect extended (up to ~ 100 kpc) Lyalpha emission around all MAGG quasars, finding a very weak redshift evolution between z=3 and z=6. By stacking the MUSE datacubes, we confidently detect extended emission of CIV and only marginally detect extended HeII up to ~40 kpc, implying that the gas is metal enriched. Moreover, our observations show a significant overdensity of LAEs within 300 km/s from the quasar systemic redshifts estimated from the nebular emission. The luminosity functions and equivalent width distributions of these LAEs show similar shapes with respect to LAEs away from quasars suggesting that the Lyalpha emission of the majority of these sources is not significantly boosted by the quasar radiation or other processes related to the quasar environment. Within this framework, the observed LAE overdensities and our kinematic measurements imply that bright quasars at z=3-4.5 are hosted by haloes in the mass range ~ 10^{12.0}-10^{12.5} Msun. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: 23 pages, 20 figures, accepted for publication in MNRAS

IGM Transmission Bias for $z$ $\geq$ 2.9 Lyman Continuum Detected Galaxies

Authors: R. Bassett, E. V. Ryan-Weber, J. Cooke, U. Meštrić, K. Kakiichi, L. Prichard, M. Rafelski

Abstract: Understanding the relationship between the underlying escape fraction of Lyman continuum (LyC) photons ($f_{\rm esc}$) emitted by galaxies and measuring the distribution of observed $f_{\rm esc}$ values at high redshift is fundamental to the interpretation of the reionization process. In this paper we perform a statistical exploration of the attenuation of LyC photons by neutral hydrogen in the in… ▽ More Understanding the relationship between the underlying escape fraction of Lyman continuum (LyC) photons ($f_{\rm esc}$) emitted by galaxies and measuring the distribution of observed $f_{\rm esc}$ values at high redshift is fundamental to the interpretation of the reionization process. In this paper we perform a statistical exploration of the attenuation of LyC photons by neutral hydrogen in the intergalactic medium using ensembles of simulated transmission functions. We show that LyC detected galaxies are more likely to be found in sightlines with higher-than-average transmission of LyC photons. This means that adopting a mean transmission at a given redshift leads to an overestimate of the true $f_{\rm esc}$ for LyC detected galaxies. We note, however, that mean values are appropriate for $f_{\rm esc}$ estimates of larger parent samples that include LyC non-detected galaxies. We quantify this IGM transmission bias for LyC detections in photometric and spectroscopic surveys in the literature and show that the bias is stronger for both shallower observations and for fainter parent samples (i.e. Lyman $α$ emitters versus Lyman break galaxies). We also explore the effects of varying the underlying probability distribution function (PDF) of $f_{\rm esc}$ on recovered values, showing that the underlying $f_{\rm esc}$ PDF may depend on sample selection by comparing with observational surveys. This work represents a first step in improved interpretation of LyC detections in the context of understanding $f_{\rm esc}$ from high redshift galaxies. △ Less

Submitted 3 January, 2021; originally announced January 2021.

Comments: 20 pages, 17 figures, resubmitted to MNRAS after minor revision

Constraints on the rate of supernovae lasting for more than a year from Subaru/Hyper Suprime-Cam

Authors: Takashi J. Moriya, Ji-an Jiang, Naoki Yasuda, Mitsuru Kokubo, Kojiro Kawana, Keiichi Maeda, Yen-Chen Pan, Robert M. Quimby, Nao Suzuki, Ichiro Takahashi, Masaomi Tanaka, Nozomu Tominaga, Ken'ichi Nomoto, Jeff Cooke, Lluis Galbany, Santiago Gonzalez-Gaitan, Chien-Hsiu Lee, Giuliano Pignata

Abstract: Some supernovae such as pair-instability supernovae are predicted to have the duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2m Subaru telescope. HSC is a wide-field (a 1.75 deg2 field-of-view) camera and it can efficiently conduct transient… ▽ More Some supernovae such as pair-instability supernovae are predicted to have the duration of more than a year in the observer frame. To constrain the rates of supernovae lasting for more than a year, we conducted a long-term deep transient survey using Hyper Suprime-Cam (HSC) on the 8.2m Subaru telescope. HSC is a wide-field (a 1.75 deg2 field-of-view) camera and it can efficiently conduct transient surveys. We observed the same 1.75 deg2 field repeatedly using the g, r, i, and z band filters with the typical depth of 26 mag for 4 seasons (from late 2016 to early 2020). Using these data, we searched for transients lasting for more than a year. Two supernovae were detected in 2 continuous seasons, one supernova was detected in 3 continuous seasons, but no transients lasted for all 4 seasons searched. The discovery rate of supernovae lasting for more than a year with the typical limiting magnitudes of 26 mag is constrained to be 1.4^{+1.3}_{-0.7}(stat.)^{+0.2}_{-0.3}(sys.) events deg-2 yr-1. All the long-lasting supernovae we found are likely Type IIn supernovae and our results indicate that about 40% of Type IIn supernovae have long-lasting light curves. No plausible pair-instability supernova candidates lasting for more than a year are discovered. By comparing the survey results and survey simulations, we constrain the luminous pair-instability supernova rate up to z ~ 3 should be of the order of 100 Gpc-3 yr-1 at most, which is 0.01 - 0.1 per cent of the core-collapse supernova rate. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: 18 pages, 9 figures, 6 tables, accepted by The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 908, Issue 2, id.249, 13 pp. (2021)

Inverse Designed THz Spectral Splitters

Authors: Sourangsu Banerji, Yu Shi, Vivian Song-En Su, Udayan Ghosh, Jacqueline Cooke, Yong Lin Kong, Lei Liu, Berardi Sensale-Rodriguez

Abstract: This letter reports proof-of-principle demonstration of 3D printable, low-cost, and compact THz spectral splitters based on diffractive optical elements (DOEs) designed to disperse the incident collimated broadband THz radiation (0.5 THz - 0.7 THz) at a pre-specified distance. Via inverse design, we show that it is possible to design such a diffractive optic, which can split the broadband incident… ▽ More This letter reports proof-of-principle demonstration of 3D printable, low-cost, and compact THz spectral splitters based on diffractive optical elements (DOEs) designed to disperse the incident collimated broadband THz radiation (0.5 THz - 0.7 THz) at a pre-specified distance. Via inverse design, we show that it is possible to design such a diffractive optic, which can split the broadband incident spectrum in any desired fashion, as is evidenced from both FDTD simulations and measured intensity profiles using a 500-750 GHz VNA. Due to its straightforward construction without the usage of any movable parts, our approach, in principle, can have various applications such as in portable, low-cost spectroscopy as well as in wireless THz communication systems as a THz demultiplexer. △ Less

Submitted 28 October, 2020; originally announced November 2020.

Comments: This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

MUSE Analysis of Gas around Galaxies (MAGG) -- II: Metal-enriched halo gas around z~1 galaxies

Authors: Rajeshwari Dutta, Michele Fumagalli, Matteo Fossati, Emma K. Lofthouse, J. Xavier Prochaska, Fabrizio Arrigoni Battaia, Richard M. Bielby, Sebastiano Cantalupo, Ryan J. Cooke, Michael T. Murphy, John M. O'Meara

Abstract: We present a study of the metal-enriched cool halo gas traced by MgII absorption around 228 galaxies at z~0.8-1.5 within 28 quasar fields from the MUSE Analysis of Gas around Galaxies (MAGG) survey. We observe no significant evolution in the MgII equivalent width versus impact parameter relation and in the MgII covering fraction compared to surveys at z<~0.5. The stellar mass, along with distance… ▽ More We present a study of the metal-enriched cool halo gas traced by MgII absorption around 228 galaxies at z~0.8-1.5 within 28 quasar fields from the MUSE Analysis of Gas around Galaxies (MAGG) survey. We observe no significant evolution in the MgII equivalent width versus impact parameter relation and in the MgII covering fraction compared to surveys at z<~0.5. The stellar mass, along with distance from galaxy centre, appears to be the dominant factor influencing the MgII absorption around galaxies. With a sample that is 90% complete down to a star formation rate of ~0.1 Msun/yr and up to impact parameters ~250-350 kpc from quasars, we find that the majority (67^{+12}_{-15}% or 14/21) of the MgII absorption systems are associated with more than one galaxy. The complex distribution of metals in these richer environments adds substantial scatter to previously-reported correlations. Multiple galaxy associations show on average five times stronger absorption and three times higher covering fraction within twice the virial radius than isolated galaxies. The dependence of MgII absorption on galaxy properties disfavours the scenario in which a widespread intra-group medium dominates the observed absorption. This leaves instead gravitational interactions among group members or hydrodynamic interactions of the galaxy haloes with the intra-group medium as favoured mechanisms to explain the observed enhancement in the MgII absorption strength and cross section in rich environments. △ Less

Submitted 16 November, 2020; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 499, Issue 4, pp.5022-5046, October 2020

Role of Fabrication Errors and Refractive Index on Multilevel Diffractive Lens Performance

Authors: Sourangsu Banerji, Jacqueline Cooke, Berardi Sensale-Rodriguez

Abstract: Multilevel diffractive lenses (MDLs) have emerged as an alternative to both conventional diffractive optical elements (DOEs) and metalenses for applications ranging from imaging to holographic and immersive displays. Recent work has shown that by harnessing structural parametric optimization of DOEs, one can design MDLs to enable multiple functionalities like achromaticity, depth of focus, wide-an… ▽ More Multilevel diffractive lenses (MDLs) have emerged as an alternative to both conventional diffractive optical elements (DOEs) and metalenses for applications ranging from imaging to holographic and immersive displays. Recent work has shown that by harnessing structural parametric optimization of DOEs, one can design MDLs to enable multiple functionalities like achromaticity, depth of focus, wide-angle imaging, etc. with great ease in fabrication. Therefore, it becomes critical to understand how fabrication errors still do affect the performance of MDLs and numerically evaluate the trade-off between efficiency and initial parameter selection, right at the onset of designing an MDL, i.e., even before putting it into fabrication. Here, we perform a statistical simulation-based study on MDLs (primarily operating in the THz regime) to analyze the impact of various fabrication imperfections (single and multiple) on the final structure as a function of the number of ring height levels. Furthermore, we also evaluate the performance of these same MDLs with the change in the refractive index of the constitutive material. We use focusing efficiency as the evaluation criterion in our numerical analysis; since it is the most fundamental property that can be used to compare and assess the performance of lenses (and MDLs) in general designed for any application with any specific functionality. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: Initial manuscript version of the published article. Some comments were removed in the final published version at the discretion of the reviewers and not of the authors. However the authors would still like to include the discussion for the benefit of the optics community

Journal ref: Sci Rep 10, 14608 (2020)

On the genus two skein algebra

Authors: Juliet Cooke, Peter Samuelson

Abstract: We study the skein algebra of the genus 2 surface and its action on the skein module of the genus 2 handlebody. We compute this action explicitly, and we describe how the module decomposes over certain subalgebras in terms of polynomial representations of double affine Hecke algebras. Finally, we show that this algebra is isomorphic to the $t=q$ specialisation of the genus two spherical double aff… ▽ More We study the skein algebra of the genus 2 surface and its action on the skein module of the genus 2 handlebody. We compute this action explicitly, and we describe how the module decomposes over certain subalgebras in terms of polynomial representations of double affine Hecke algebras. Finally, we show that this algebra is isomorphic to the $t=q$ specialisation of the genus two spherical double affine Hecke algebra recently defined by Arthamonov and Shakirov. △ Less

Submitted 26 October, 2021; v1 submitted 28 August, 2020; originally announced August 2020.

Comments: 35 pages; Journal of the London Mathematical Society (in press), Early view online

Unsupervised machine learning for transient discovery in Deeper, Wider, Faster light curves

Authors: Sara Webb, Michelle Lochner, Daniel Muthukrishna, Jeff Cooke, Chris Flynn, Ashish Mahabal, Simon Goode, Igor Andreoni, Tyler Pritchard, Timothy M. C. Abbott

Abstract: Identification of anomalous light curves within time-domain surveys is often challenging. In addition, with the growing number of wide-field surveys and the volume of data produced exceeding astronomers ability for manual evaluation, outlier and anomaly detection is becoming vital for transient science. We present an unsupervised method for transient discovery using a clustering technique and the… ▽ More Identification of anomalous light curves within time-domain surveys is often challenging. In addition, with the growing number of wide-field surveys and the volume of data produced exceeding astronomers ability for manual evaluation, outlier and anomaly detection is becoming vital for transient science. We present an unsupervised method for transient discovery using a clustering technique and the Astronomaly package. As proof of concept, we evaluate 85553 minute-cadenced light curves collected over two 1.5 hour periods as part of the Deeper, Wider, Faster program, using two different telescope dithering strategies. By combining the clustering technique HDBSCAN with the isolation forest anomaly detection algorithm via the visual interface of Astronomaly, we are able to rapidly isolate anomalous sources for further analysis. We successfully recover the known variable sources, across a range of catalogues from within the fields, and find a further 7 uncatalogued variables and two stellar flare events, including a rarely observed ultra fast flare (5 minute) from a likely M-dwarf. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: Accepted 7 Aug 2020, 19 pages, 8 figures,

Constraining the Kilonova Rate with Zwicky Transient Facility Searches Independent of Gravitational Wave and Short Gamma-ray Burst Triggers

Authors: Igor Andreoni, Erik C. Kool, Ana Sagues Carracedo, Mansi M. Kasliwal, Mattia Bulla, Tomas Ahumada, Michael W. Coughlin, Shreya Anand, Jesper Sollerman, Ariel Goobar, David L. Kaplan, Tegan T. Loveridge, Viraj Karambelkar, Jeff Cooke, Ashot Bagdasaryan, Eric C. Bellm, S. Bradley Cenko, David O. Cook, Kishalay De, Richard Dekany, Alexandre Delacroix, Andrew Drake, Dmitry A. Duev, Christoffer Fremling, V. Zach Golkhou , et al. (12 additional authors not shown)

Abstract: The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of r… ▽ More The first binary neutron star merger, GW170817, was accompanied by a radioactivity-powered optical/infrared transient called a kilonova. To date, no compelling kilonova has been found during optical surveys of the sky, independent of gravitational-wave triggers. In this work, we searched the first 23 months of the Zwicky Transient Facility (ZTF) data stream for candidate kilonovae in the form of rapidly evolving transients. We combined ZTF alert queries with forced point-spread-function photometry and nightly flux stacking to increase our sensitivity to faint and fast transients. Automatic queries yielded $>11,200$ candidates, 24 of which passed quality checks and strict selection criteria based on a grid of kilonova models tailored for both binary neutron star and neutron star-black hole mergers. None of the candidates in our sample was deemed a possible kilonova after thorough vetting, catalog cross-matching, and study of their color evolution. The sources that passed our selection criteria are dominated by Galactic cataclysmic variables. In addition, we identified two fast transients at high Galactic latitude, one of which is the confirmed afterglow of long-duration GRB190106A, and the other is a possible cosmological afterglow. Using a survey simulation code, we constrained the kilonova rate for a range of models including top-hat and linearly decaying light curves and synthetic light curves obtained with radiative transfer simulations. For prototypical GW170817-like kilonovae, we constrain the rate to be $R < 1775$ Gpc$^{-3}$ yr$^{-1}$ at 95% confidence level by requiring at least 2 high-significance detections. By assuming a population of kilonovae with the same geometry and composition of GW170817 observed under a uniform viewing angle distribution, we obtained a constraint on the rate of $R < 4029$ Gpc$^{-3}$ yr$^{-1}$. △ Less

Submitted 31 July, 2020; originally announced August 2020.

Comments: Submitted for publication in ApJ

Neutron Star Extreme Matter Observatory: A kilohertz-band gravitational-wave detector in the global network

Authors: K. Ackley, V. B. Adya, P. Agrawal, P. Altin, G. Ashton, M. Bailes, E. Baltinas, A. Barbuio, D. Beniwal, C. Blair, D. Blair, G. N. Bolingbroke, V. Bossilkov, S. Shachar Boublil, D. D. Brown, B. J. Burridge, J. Calderon Bustillo, J. Cameron, H. Tuong Cao, J. B. Carlin, S. Chang, P. Charlton, C. Chatterjee, D. Chattopadhyay, X. Chen , et al. (139 additional authors not shown)

Abstract: Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provid… ▽ More Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica. △ Less

Submitted 5 November, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: Accepted for publication in PASA

Journal ref: PASA (2020) 37, e047

A DESGW Search for the Electromagnetic Counterpart to the LIGO/Virgo Gravitational Wave Binary Neutron Star Merger Candidate S190510g

Authors: DES Collaboration, A. Garcia, R. Morgan, K. Herner, A. Palmese, M. Soares-Santos, J. Annis, D. Brout, A. K. Vivas, A. Drlica-Wagner, L. Santana-Silva, D. L. Tucker, S. Allam, M. Wiesner, J. García-Bellido, M. S. S. Gill, M. Sako, R. Kessler, T. M. Davis, D. Scolnic, F. Olivares E., F. Paz-Chinchón, N. Sherman, C. Conselice, H. Chen , et al. (65 additional authors not shown)

Abstract: We present the results from a search for the electromagnetic counterpart of the LIGO/Virgo event S190510g using the Dark Energy Camera (DECam). S190510g is a binary neutron star (BNS) merger candidate of moderate significance detected at a distance of 227$\pm$92 Mpc and localized within an area of 31 (1166) square degrees at 50\% (90\%) confidence. While this event was later classified as likely n… ▽ More We present the results from a search for the electromagnetic counterpart of the LIGO/Virgo event S190510g using the Dark Energy Camera (DECam). S190510g is a binary neutron star (BNS) merger candidate of moderate significance detected at a distance of 227$\pm$92 Mpc and localized within an area of 31 (1166) square degrees at 50\% (90\%) confidence. While this event was later classified as likely non-astrophysical in nature within 30 hours of the event, our short latency search and discovery pipeline identified 11 counterpart candidates, all of which appear consistent with supernovae following offline analysis and spectroscopy by other instruments. Later reprocessing of the images enabled the recovery of 6 more candidates. Additionally, we implement our candidate selection procedure on simulated kilonovae and supernovae under DECam observing conditions (e.g., seeing, exposure time) with the intent of quantifying our search efficiency and making informed decisions on observing strategy for future similar events. This is the first BNS counterpart search to employ a comprehensive simulation-based efficiency study. We find that using the current follow-up strategy, there would need to be 19 events similar to S190510g for us to have a 99\% chance of detecting an optical counterpart, assuming a GW170817-like kilonova. We further conclude that optimization of observing plans, which should include preference for deeper images over multiple color information, could result in up to a factor of 1.5 reduction in the total number of followups needed for discovery. △ Less

Submitted 30 June, 2020; originally announced July 2020.

Comments: Paper submitted to ApJL

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

Authors: Mansi M. Kasliwal, Shreya Anand, Tomas Ahumada, Robert Stein, Ana Sagues Carracedo, Igor Andreoni, Michael W. Coughlin, Leo P. Singer, Erik C. Kool, Kishalay De, Harsh Kumar, Mouza AlMualla, Yuhan Yao, Mattia Bulla, Dougal Dobie, Simeon Reusch, Daniel A. Perley, S. Bradley Cenko, Varun Bhalerao, David L. Kaplan, Jesper Sollerman, Ariel Goobar, Christopher M. Copperwheat, Eric C. Bellm, G. C. Anupama , et al. (78 additional authors not shown)

Abstract: We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GR… ▽ More We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged) △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ

The PHLEK Survey: A New Determination of the Primordial Helium Abundance

Authors: Tiffany Hsyu, Ryan J. Cooke, J. Xavier Prochaska, Michael Bolte

Abstract: We present Keck NIRSPEC and Keck NIRES spectroscopy of sixteen metal-poor galaxies that have pre-existing optical observations. The near-infrared (NIR) spectroscopy specifically targets the HeI 10830A emission line, due to its sensitivity to the physical conditions of the gas in HII regions. We use these NIR observations, combined with optical spectroscopy, to determine the helium abundance of six… ▽ More We present Keck NIRSPEC and Keck NIRES spectroscopy of sixteen metal-poor galaxies that have pre-existing optical observations. The near-infrared (NIR) spectroscopy specifically targets the HeI 10830A emission line, due to its sensitivity to the physical conditions of the gas in HII regions. We use these NIR observations, combined with optical spectroscopy, to determine the helium abundance of sixteen galaxies across a metallicity range 12 + log (O/H) = 7.13 - 8.00. This data set is combined with two other samples where metallicity and helium abundance measurements can be secured: star-forming galaxies selected from the Sloan Digital Sky Survey spectroscopic database and existing low-metallicity systems in the literature. We calculate a linear fit to these measurements, accounting for intrinsic scatter, and report a new determination of the primordial helium number abundance, y_P = 0.0805+/-0.0017, which corresponds to a primordial helium mass fraction Y_P = 0.2436+/-0.0040. Using our determination of the primordial helium abundance in combination with the latest primordial deuterium measurement, (D/H)_P x 10^5 =2.527+/-0.030, we place a bound on the baryon density Omega_bh^2 = 0.0215+/-0.0005 and the effective number of neutrino species Neff = 2.85+/-0.28. These values are in 1.3-sigma agreement with those deduced from the Planck satellite observations of the temperature fluctuations imprinted on the Cosmic Microwave Background. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 32 pages, 11 figures, Accepted for publication in the Astrophysical Journal