Showing 1–50 of 340 results for author: Kaplan, D L

An emission state switching radio transient with a 54 minute period

Authors: M. Caleb, E. Lenc, D. L. Kaplan, T. Murphy, Y. P. Men, R. M. Shannon, L. Ferrario, K. M. Rajwade, T. E. Clarke, S. Giacintucci, N. Hurley-Walker, S. D. Hyman, M. E. Lower, Sam McSweeney, V. Ravi, E. D. Barr, S. Buchner, C. M. L. Flynn, J. W. T. Hessels, M. Kramer, J. Pritchard, B. W. Stappers

Abstract: Long-period radio transients are an emerging class of extreme astrophysical events of which only three are known. These objects emit highly polarised, coherent pulses of typically a few tens of seconds duration and minutes to hour-long periods. While magnetic white dwarfs and magnetars, either isolated or in binary systems, have been invoked to explain these objects, a consensus has not emerged. H… ▽ More Long-period radio transients are an emerging class of extreme astrophysical events of which only three are known. These objects emit highly polarised, coherent pulses of typically a few tens of seconds duration and minutes to hour-long periods. While magnetic white dwarfs and magnetars, either isolated or in binary systems, have been invoked to explain these objects, a consensus has not emerged. Here we report on the discovery of ASKAP J193505.1+214841.0 (henceforth ASKAPJ1935+2148) with a period of 53.8 minutes exhibiting three distinct emission states - a bright pulse state with highly linearly polarised pulses with widths of 10-50 seconds; a weak pulse state which is about 26 times fainter than the bright state with highly circularly polarised pulses of widths of approximately 370 milliseconds; and a quiescent or quenched state with no pulses. The first two states have been observed to progressively evolve over the course of 8 months with the quenched state interspersed between them suggesting physical changes in the region producing the emission. A constraint on the radius of the source for the observed period rules out a magnetic white dwarf origin. Unlike other long-period sources, ASKAPJ1935+2148 is the first to exhibit drastic variations in emission modes reminiscent of neutron stars. However, its radio properties challenge our current understanding of neutron star emission and evolution. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: Published in Nature Astronomy

Radio afterglows from tidal disruption events: An unbiased sample from ASKAP RACS

Authors: Akash Anumarlapudi, Dougal Dobie, David L. Kaplan, Tara Murphy, Assaf Horesh, Emil Lenc, Laura N. Driessen, Stefan W. Duchesne, Ms. Hannah Dykaar, Bryan M. Gaensler, Timothy J. Galvin, J. A. Grundy, George Heald, Aidan Hotan, Minh Huynh, James Leung, David McConnell, Vanessa A. Moss, Joshua Pritchard, Wasim Raja, Kovi Rose, Gregory R. Sivakoff, Yuanming Wang, Ziteng Wang, Mark Wieringa , et al. (1 additional authors not shown)

Abstract: Late-time ($\sim$ year) radio follow-up of optically-discovered tidal disruption events (TDEs) is increasingly resulting in detections at radio wavelengths, and there is growing evidence for this late-time radio activity to be common to the broad class of sub-relativistic TDEs. Detailed studies of some of these TDEs at radio wavelengths are also challenging the existing models for radio emission.… ▽ More Late-time ($\sim$ year) radio follow-up of optically-discovered tidal disruption events (TDEs) is increasingly resulting in detections at radio wavelengths, and there is growing evidence for this late-time radio activity to be common to the broad class of sub-relativistic TDEs. Detailed studies of some of these TDEs at radio wavelengths are also challenging the existing models for radio emission. Using all-sky multi-epoch data from the Australian Square Kilometre Array Pathfinder (ASKAP), taken as a part of the Rapid ASKAP Continuum Survey (RACS), we searched for radio counterparts to a sample of optically-discovered TDEs. We detected late-time emission at RACS frequencies (742-1032\,MHz) in five TDEs, reporting the independent discovery of radio emission from TDE AT2019ahk and extending the time baseline out to almost 3000\,days for some events. Overall, we find that at least $22^{+15}_{-11}$\% of the population of optically-discovered TDEs has detectable radio emission in the RACS survey, while also noting that the true fraction can be higher given the limited cadence (2 epochs separated by $\sim 3\,$ years) of the survey. Finally, we project that the ongoing higher-cadence ($\sim 2$\,months) ASKAP Variable and Slow Transients (VAST) survey can detect $\sim 20$ TDEs in its operational span (4\,yrs), given the current rate from optical surveys. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: Accepted for publication in ApJ, comments welcome

The Anomalous Acceleration of PSR J2043+1711: Long-Period Orbital Companion or Stellar Flyby?

Authors: Thomas Donlon II, Sukanya Chakrabarti, Michael T. Lam, Daniel Huber, Daniel Hey, Enrico Ramirez-Ruiz, Benjamin Shappee, David L. Kaplan, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Paul R. Brook, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels, Peter A. Gentile , et al. (31 additional authors not shown)

Abstract: Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4σ$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr… ▽ More Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4σ$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in \textit{Gaia} DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6\%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 80 kyr, a companion mass of 0.3 $M_\odot$ (indicative of a white dwarf or low-mass star), and a semi-major axis of 2000 AU. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios. △ Less

Submitted 8 July, 2024; originally announced July 2024.

An Image-Based Search for Pulsar Candidates in the MeerKAT Bulge Survey

Authors: Dale A. Frail, Emil Polisensky, Scott D. Hyman, W. M. Cotton, Namir E. Kassim, Michele L. Silverstein, Rahul Sengar, David L. Kaplan, Francesca Calore, Joanna Berteaud, Maica Clavel, Marisa Geyer, Samuel Legodi, Vasaant Krishnan, Sarah Buchner, Fernando Camilo

Abstract: We report on the results of an image-based search for pulsar candidates toward the Galactic bulge. We used mosaic images from the MeerKAT radio telescope, that were taken as part of a 173 deg**2 survey of the bulge and Galactic center of our Galaxy at L band (856-1712 MHz) in all four Stokes I, Q, U and V. The image root-mean-square noise levels of 12-17 uJy/ba represent a significant increase in… ▽ More We report on the results of an image-based search for pulsar candidates toward the Galactic bulge. We used mosaic images from the MeerKAT radio telescope, that were taken as part of a 173 deg**2 survey of the bulge and Galactic center of our Galaxy at L band (856-1712 MHz) in all four Stokes I, Q, U and V. The image root-mean-square noise levels of 12-17 uJy/ba represent a significant increase in sensitivity over past image-based pulsar searches. Our primary search criterion was circular polarization, but we used other criteria including linear polarization, in-band spectral index, compactness, variability and multi-wavelength counterparts to select pulsar candidates. We first demonstrate the efficacy of this technique by searching for polarized emission from known pulsars, and comparing our results with measurements from the literature. Our search resulted in a sample of 75 polarized pulsar candidates. Bright stars or young stellar objects were associated with 28 of these sources, including a small sample of highly polarized dwarf stars with pulsar-like steep spectra. Comparing the properties of this sample with the known pulsars, we identified 30 compelling candidates for pulsation follow-up, including two sources with both strong circular and linear polarization. The remaining 17 sources are either pulsars or stars, but we cannot rule out an extragalactic origin or image artifacts among the brighter, flat spectrum objects. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 28 pages, 10 figures, 6 tables. ApJ, submitted

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.

An Untargeted Search for Radio-Emitting Tidal Disruption Events in the VAST Pilot Survey

Authors: Hannah Dykaar, Maria R. Drout, B. M. Gaensler, David L. Kaplan, Tara Murphy, Assaf Horesh, Akash Anumarlapudi, Dougal Dobie, Laura N. Driessen, Emil Lenc, Adam Stewart

Abstract: We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts… ▽ More We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts through follow-up observations. With systematic time-domain radio surveys becoming available, we can now identify new TDEs in the radio regime. A population of radio-discovered TDEs has the potential to provide several key insights including an independent constraint on their volumetric rate. We conducted a search to select variable radio sources with a single prominent radio flare and a position consistent within 2$σ$ of the nucleus of a known galaxy. While TDEs were the primary target of our search, sources identified in this search may also be consistent with active galactic nuclei exhibiting unusual flux density changes at the timescales probed, uncharacteristically bright supernovae, or a population of gamma-ray bursts. We identify a sample of 12 radio-bright candidate TDEs. The timescales and luminosities range from ~6 to 230 days and ~10$^{38}$ to 10$^{41}$ erg s$^{-1}$, consistent with models of radio emission from TDEs that launch relativistic jets. After calculating the detection efficiency of our search using a Monte Carlo simulation of TDEs, and assuming all 12 sources are jetted TDEs, we derive a volumetric rate for jetted TDEs of 0.80$^{+0.31}_{-0.23}$ Gpc$^{-3}$ yr$^{-1}$, consistent with previous empirically estimated rates. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 37 pages, 37 figures, accepted to ApJ

The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars

Authors: Bjorn Larsen, Chiara M. F. Mingarelli, Jeffrey S. Hazboun, Aurelien Chalumeau, Deborah C. Good, Joseph Simon, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Paul R. Brook, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels, Peter A. Gentile, Joseph Glaser, Ross J. Jennings , et al. (39 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characteriz… ▽ More Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect achromatic noise properties in each pulsar. The models we compare include existing DM models used by NANOGrav and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar dataset. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from $\log_{10}A_{\text{RN}} = -14.1^{+0.1}_{-0.1}$ to $-14.7^{+0.3}_{-0.5}$, and the spectral index broadens from $γ_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $γ_{\text{RN}} = 3.5^{+1.2}_{-0.9}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Searching for gravitational wave optical counterparts with the Zwicky Transient Facility: summary of O4a

Authors: Tomás Ahumada, Shreya Anand, Michael W. Coughlin, Vaidehi Gupta, Mansi M. Kasliwal, Viraj R. Karambelkar, Robert D. Stein, Gaurav Waratkar, Vishwajeet Swain, Theophile Jegou du Laz, Akash Anumarlapudi, Igor Andreoni, Mattia Bulla, Gokul P. Srinivasaragavan, Andrew Toivonen, Avery Wold, Eric C. Bellm, S. Bradley Cenko, David L. Kaplan, Jesper Sollerman, Varun Bhalerao, Daniel Perley, Anirudh Salgundi, Aswin Suresh, K-Ryan Hinds , et al. (27 additional authors not shown)

Abstract: During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candida… ▽ More During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations (ToO) and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg$^2$ within the skymaps of the high-significance GW events, reaching a median depth of g~20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than -16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH), along with our O4a follow-ups. Following Kasliwal et al. (2020), we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag/day can peak at a magnitude brighter than -17.5 mag. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: submitted

The Sydney Radio Star Catalogue: properties of radio stars at megahertz to gigahertz frequencies

Authors: Laura N. Driessen, Joshua Pritchard, Tara Murphy, George Heald, Jan Robrade, Barnali Das, Stefan Duchesne, David L. Kaplan, Emil Lenc, Christene R. Lynch, Benjamin J. S. Pope, Kovi Rose, Beate Stelzer, Yuanming Wang, Andrew Zic

Abstract: We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and pro… ▽ More We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and proper motion searches as well as presenting hundreds of newly detected stars from our search of Australian SKA Pathfinder observations. The focus of this first version of the catalogue is on objects detected in surveys using SKA precursor instruments; however we will expand this scope in future versions. The 839 objects in the Sydney Radio Star Catalogue are distributed across the whole sky and range from ultracool dwarfs to Wolf-Rayet stars. We find that the radio luminosities of cool dwarfs are lower than the radio luminosities of more evolved sub-giant and giant stars. We use X-ray detections of 530 radio stars by the eROSITA soft X-ray instrument onboard the SRG spacecraft to show that almost all of the radio stars in the catalogue are over-luminous in the radio, indicating that the majority of stars at these radio frequencies are coherent radio emitters. The Sydney Radio Star Catalogue can be found in Vizier or at https://radiostars.org. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 21 pages, 8 Figures, 7 Tables, Submitted to PASA

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

Authors: Gabriella Agazie, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Lucas Brown, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels , et al. (58 additional authors not shown)

Abstract: The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analyt… ▽ More The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures, 1 appendix, submitted to ApJ

The GMRT High-Resolution Southern Sky Survey for pulsars and transients -- VII: Timing of Spider MSP J1242-4712, A Bridge Between Redback and Black Widow Pulsars

Authors: Ankita Ghosh, Bhaswati Bhattacharyya, Andrew Lyne, David L. Kaplan, Jayanta Roy, Paul S. Ray, Ben Stappers, Sangita Kumari, Shubham Singh, Rahul Sharan

Abstract: We present the timing solution for the 5.31-ms spider millisecond pulsar (MSP) J1242-4712, discovered with the GMRT. PSR J1242-4712 orbits a companion of minimum mass 0.08 M$_{\odot}$ with an orbital period of 7.7 hrs and occupies a relatively unexplored region in the orbital period versus companion mass space. We did not detect gamma-ray pulsations for this MSP, and also could not identify the op… ▽ More We present the timing solution for the 5.31-ms spider millisecond pulsar (MSP) J1242-4712, discovered with the GMRT. PSR J1242-4712 orbits a companion of minimum mass 0.08 M$_{\odot}$ with an orbital period of 7.7 hrs and occupies a relatively unexplored region in the orbital period versus companion mass space. We did not detect gamma-ray pulsations for this MSP, and also could not identify the optical counterpart for PSR J1242-4712 in available optical/near-infrared data. The profile of J1242-4712 evolves with frequency showing a clear single component at lower frequencies and a three-component profile at 650 MHz. PSR J1242-4712 eclipses for a very short duration near superior conjunction (orbital phase ~ 0.23-0.25) below 360 MHz. Moreover, significant DM delays and errors in pulse times of arrivals are observed near inferior conjunction (orbital phase ~ 0.7), along with an observed eclipse in one epoch at 650 MHz. Observed eclipses and significant orbital period variability suggest that PSR J1242-4712 is possibly not a He-WD binary, but has a semi or non-degenerate companion, indicating that this is a ``spider" MSP lying in a region between typical black widows and redbacks. This system may represent a distinct category of spider MSPs, displaying characteristics that bridge the gap between known black widow and redback MSPs. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 13 pages, 8 figures, 2 tables. Accepted for publication in The Astrophysical Journal

Results of the follow-up of ANTARES neutrino alerts

Authors: A. Albert, S. Alves, M. André, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, Y. Becherini, B. Belhorma, M. Bendahman, F. Benfenati, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzas, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (166 additional authors not shown)

Abstract: High-energy neutrinos could be produced in the interaction of charged cosmic rays with matter or radiation surrounding astrophysical sources. To look for transient sources associated with neutrino emission, a follow-up program of neutrino alerts has been operating within the ANTARES Collaboration since 2009. This program, named TAToO, has triggered robotic optical telescopes (MASTER, TAROT, ROTSE… ▽ More High-energy neutrinos could be produced in the interaction of charged cosmic rays with matter or radiation surrounding astrophysical sources. To look for transient sources associated with neutrino emission, a follow-up program of neutrino alerts has been operating within the ANTARES Collaboration since 2009. This program, named TAToO, has triggered robotic optical telescopes (MASTER, TAROT, ROTSE and the SVOM ground based telescopes) immediately after the detection of any relevant neutrino candidate and scheduled several observations in the weeks following the detection. A subset of ANTARES events with highest probabilities of being of cosmic origin has also been followed by the Swift and the INTEGRAL satellites, the Murchison Widefield Array radio telescope and the H.E.S.S. high-energy gamma-ray telescope. The results of twelve years of observations are reported. No optical counterpart has been significantly associated with an ANTARES candidate neutrino signal during image analysis. Constraints on transient neutrino emission have been set. In September 2015, ANTARES issued a neutrino alert and during the follow-up, a potential transient counterpart was identified by Swift and MASTER. A multi-wavelength follow-up campaign has allowed to identify the nature of this source and has proven its fortuitous association with the neutrino. The return of experience is particularly important for the design of the alert system of KM3NeT, the next generation neutrino telescope in the Mediterranean Sea. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 27 pages, 14 figures, submitted to JCAP

A Case for a Binary Black Hole System Revealed via Quasi-Periodic Outflows

Authors: Dheeraj R. Pasham, Francesco Tombesi, Petra Sukova, Michal Zajacek, Suvendu Rakshit, Eric Coughlin, Peter Kosec, Vladimir Karas, Megan Masterson, Andrew Mummery, Thomas W. -S. Holoien, Muryel Guolo, Jason Hinkle, Bart Ripperda, Vojtech Witzany, Ben Shappee, Erin Kara, Assaf Horesh, Sjoert van Velzen, Itai Sfaradi, David L. Kaplan, Noam Burger, Tara Murphy, Ronald Remillard, James F. Steiner , et al. (11 additional authors not shown)

Abstract: Binaries containing a compact object orbiting a supermassive black hole are thought to be precursors of gravitational wave events, but their identification has been extremely challenging. Here, we report quasi-periodic variability in X-ray absorption which we interpret as quasi-periodic outflows (QPOuts) from a previously low-luminosity active galactic nucleus after an outburst, likely caused by a… ▽ More Binaries containing a compact object orbiting a supermassive black hole are thought to be precursors of gravitational wave events, but their identification has been extremely challenging. Here, we report quasi-periodic variability in X-ray absorption which we interpret as quasi-periodic outflows (QPOuts) from a previously low-luminosity active galactic nucleus after an outburst, likely caused by a stellar tidal disruption. We rule out several models based on observed properties and instead show using general relativistic magnetohydrodynamic simulations that QPOuts, separated by roughly 8.3 days, can be explained with an intermediate-mass black hole secondary on a mildly eccentric orbit at a mean distance of about 100 gravitational radii from the primary. Our work suggests that QPOuts could be a new way to identify intermediate/extreme-mass ratio binary candidates. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: Accepted for publication in Science Advances. We report a new supermassive black hole phenomenon that we call quasi-periodic outflows (QPOuts)

Intensive Swift and LCO monitoring of PG 1302$-$102: AGN disk reverberation mapping of a supermassive black hole binary candidate

Authors: Tingting Liu, Rick Edelson, Juan V. Hernández Santisteban, Erin Kara, John Montano, Jonathan Gelbord, Keith Horne, Aaron J. Barth, Edward M. Cackett, David L. Kaplan

Abstract: We present an intensive multiwavelength monitoring campaign of the quasar PG 1302$-$102 with Swift and the Las Cumbres Observatory network telescopes. At $z\sim0.3$, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time… ▽ More We present an intensive multiwavelength monitoring campaign of the quasar PG 1302$-$102 with Swift and the Las Cumbres Observatory network telescopes. At $z\sim0.3$, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302$-$102 was observed at a $\sim$daily cadence for $\sim 9$ months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the inter-band time lags which are consistent with a $τ\propto λ^{4/3}$ relation as expected from standard disk reprocessing, albeit with large errors. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped AGN in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 11 pages, 5 figures, 4 tables. Accepted for publication in ApJ. The full dataset will be available with the ApJ article

Multi-epoch sampling of the radio star population with the Australian SKA Pathfinder

Authors: Joshua Pritchard, Tara Murphy, George Heald, Michael S. Wheatland, David L. Kaplan, Emil Lenc, Andrew O'Brien, Ziteng Wang

Abstract: The population of radio-loud stars has to date been studied primarily through either targeted observations of a small number of highly active stars or widefield, single-epoch surveys that cannot easily distinguish stellar emission from background extra-Galactic sources. As a result it has been difficult to constrain population statistics such as the surface density and fraction of the population p… ▽ More The population of radio-loud stars has to date been studied primarily through either targeted observations of a small number of highly active stars or widefield, single-epoch surveys that cannot easily distinguish stellar emission from background extra-Galactic sources. As a result it has been difficult to constrain population statistics such as the surface density and fraction of the population producing radio emission in a particular variable or spectral class. In this paper we present a sample of 36 radio stars detected in a circular polarisation search of the multi-epoch Variables and Slow Transients (VAST) pilot survey with ASKAP at 887.5~MHz. Through repeat sampling of the VAST pilot survey footprint we find an upper limit to the duty cycle of M-dwarf radio bursts of 8.5 per cent, and that at least $10 \pm 3$ per cent of the population should produce radio bursts more luminous than $10^{15}$ erg s$^{-1}$ Hz$^{-1}$. We infer a lower limit on the long-term surface density of such bursts in a shallow 1.25 mJy PSF$^{-1}$ sensitivity survey of $9^{+11}_{-7} \times 10^{-3}$ deg$^{-2}$ and an instantaneous radio star surface density of $1.7 \pm 0.2 \times 10^{-3}$ deg$^{-2}$ on 12 min timescales. Based on these rates we anticipate ${\sim}200 \pm 50$ new radio star detections per year over the full VAST survey and ${\sim}41\,000^{+10\,000}_{-9\,000}$ in next-generation all-sky surveys with the Square Kilometre Array. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 14 pages, 9 figures, submitted to MNRAS

The Green Bank North Celestial Cap Survey IX: Timing Follow-up for 128 Pulsars

Authors: A. E. McEwen, J. K. Swiggum, D. L. Kaplan, C. M. Tan, B. W. Meyers, E. Fonseca, G. Y. Agazie, P. Chawla, K. Crowter, M. E. DeCesar, T. Dolch, F. A. Dong, W. Fiore, E. Fonseca, D. C. Good, A. G. Istrate, V. M. Kaspi, V. I. Kondratiev, J. van Leeuwen, L. Levin, E. F. Lewis, R. S. Lynch, K. W. Masui, J. W. McKee, M. A. McLaughlin , et al. (6 additional authors not shown)

Abstract: The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations take… ▽ More The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including 9 previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6)\Msun and 1.412(6)\Msun, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves: $\dot{P}_{\rm B} = -1.37(7)\times10^{-12}$, which is in agreement with the rate predicted by general relativity for these masses. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: accepted for publication in The Astrophysical Journal

Discovery of radio eclipses from 4FGL J1646.5$-$4406: a new candidate redback pulsar binary

Authors: Andrew Zic, Ziteng Wang, Emil Lenc, David L. Kaplan, Tara Murphy, Alessandro Ridolfi, Rahul Sengar, Natasha Hurley-Walker, Dougal Dobie, James K. Leung, Joshua Pritchard, Yuanming Wang

Abstract: Large widefield surveys make possible the serendipitous discovery of rare sub-classes of pulsars. One such class are "spider"-type pulsar binaries, comprised of a pulsar in a compact orbit with a low-mass (sub)stellar companion. In a search for circularly-polarized radio sources in ASKAP Pilot Survey observations, we discovered highly variable and circularly polarized emission from a radio source… ▽ More Large widefield surveys make possible the serendipitous discovery of rare sub-classes of pulsars. One such class are "spider"-type pulsar binaries, comprised of a pulsar in a compact orbit with a low-mass (sub)stellar companion. In a search for circularly-polarized radio sources in ASKAP Pilot Survey observations, we discovered highly variable and circularly polarized emission from a radio source within the error region of the $γ$-ray source {4FGL}~J1646.5$-$4406. The variability is consistent with the eclipse of a compact, steep-spectrum source behind ablated material from a companion in a $\sim 5.3\,$h binary orbit. Based on the eclipse properties and spatial coincidence with {4FGL} J1646.5$-$4406, we argue that the source is likely a recycled pulsar in a "redback" binary system. Using properties of the eclipses from ASKAP and Murchison Widefield Array observations, we provide broad constraints on the properties of the eclipse medium. We identified a potential optical/infra-red counterpart in archival data consistent with a variable low-mass star. Using the Parkes Radio Telescope "Murriyang" and MeerKAT, we searched extensively for radio pulsations but yielded no viable detections of pulsed emission. We suggest that the non-detection of pulses is due to scattering in the intra-binary material, but scattering from the ISM can also plausibly explain the pulse non-detections if the interstellar dispersion measure exceeds $\sim$600$\,$pc$\,$cm$^{-3}$. Orbital constraints derived from optical observations of the counterpart would be highly valuable for future $γ$-ray pulsation searches, which may confirm the source nature as a pulsar. △ Less

Submitted 30 November, 2023; originally announced December 2023.

Comments: Submitted to MNRAS. 13 Pages, 10 figures, 3 tables

Discovery of a young, highly scattered pulsar PSR J1032-5804 with the Australian SKA Pathfinder

Authors: Ziteng Wang, David L. Kaplan, Rahul Sengar, Emil Lenc, Andrew Zic, Akash Anumarlapudi, B. M. Gaensler, Natasha Hurley-Walker, Tara Murphy, Yuanming Wang

Abstract: We report the discovery of a young, highly scattered pulsar in a search for highly circularly polarized radio sources as part of the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. In follow-up observations with Murriyang/Parkes, we identified PSR J1032-5804 and measured a period of 78.7 ms, dispersion measure (DM) of 819$\pm$4 pc cm$^{-3}$, rotati… ▽ More We report the discovery of a young, highly scattered pulsar in a search for highly circularly polarized radio sources as part of the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. In follow-up observations with Murriyang/Parkes, we identified PSR J1032-5804 and measured a period of 78.7 ms, dispersion measure (DM) of 819$\pm$4 pc cm$^{-3}$, rotation measure of -2000$\pm$1 rad m$^{-2}$, and a characteristic age of 34.6 kyr. We found a pulse scattering timescale at 3 GHz of ~22 ms, implying a timescale at 1 GHz of ~3845 ms, which is the third most scattered pulsar known and explains its non-detection in previous pulsar surveys. We discuss the identification of a possible pulsar wind nebula and supernova remnant in the pulsar's local environment by analyzing the pulsar spectral energy distribution and the surrounding extended emission from multiwavelength images. Our result highlights the possibility of identifying extremely scattered pulsars from radio continuum images. Ongoing and future large-scale radio continuum surveys will offer us an unprecedented opportunity to find more extreme pulsars (e.g., highly scattered, highly intermittent, highly accelerated), which will enhance our understanding of the characteristics of pulsars and the interstellar medium. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 15 pages, 10 figures, 2 tables. Accepted for publication in ApJ

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

ForceGen: End-to-end de novo protein generation based on nonlinear mechanical unfolding responses using a protein language diffusion model

Authors: Bo Ni, David L. Kaplan, Markus J. Buehler

Abstract: Through evolution, nature has presented a set of remarkable protein materials, including elastins, silks, keratins and collagens with superior mechanical performances that play crucial roles in mechanobiology. However, going beyond natural designs to discover proteins that meet specified mechanical properties remains challenging. Here we report a generative model that predicts protein designs to m… ▽ More Through evolution, nature has presented a set of remarkable protein materials, including elastins, silks, keratins and collagens with superior mechanical performances that play crucial roles in mechanobiology. However, going beyond natural designs to discover proteins that meet specified mechanical properties remains challenging. Here we report a generative model that predicts protein designs to meet complex nonlinear mechanical property-design objectives. Our model leverages deep knowledge on protein sequences from a pre-trained protein language model and maps mechanical unfolding responses to create novel proteins. Via full-atom molecular simulations for direct validation, we demonstrate that the designed proteins are novel, and fulfill the targeted mechanical properties, including unfolding energy and mechanical strength, as well as the detailed unfolding force-separation curves. Our model offers rapid pathways to explore the enormous mechanobiological protein sequence space unconstrained by biological synthesis, using mechanical features as target to enable the discovery of protein materials with superior mechanical properties. △ Less

Submitted 15 December, 2023; v1 submitted 16 October, 2023; originally announced October 2023.

The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Belinda D. Cheeseboro, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Lankeswar Dey, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara , et al. (63 additional authors not shown)

Abstract: The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for contin… ▽ More The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $η>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. △ Less

Submitted 15 January, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

MWA rapid follow-up of gravitational wave transients: prospects for detecting prompt radio counterparts

Authors: J. Tian, G. E. Anderson, A. J. Cooper, K. Gourdji, M. Sokolowski, A. Rowlinson, A. Williams, G. Sleap, D. Dobie, D. L. Kaplan, Tara Murphy, S. J. Tingay, F. H. Panther, P. D. Lasky, A. Bahramian, J. C. A. Miller-Jones, C. W. James, B. W. Meyers, S. J. McSweeney, P. J. Hancock

Abstract: We present and evaluate the prospects for detecting coherent radio counterparts to gravitational wave (GW) events using Murchison Widefield Array (MWA) triggered observations. The MWA rapid-response system, combined with its buffering mode ($\sim4$ minutes negative latency), enables us to catch any radio signals produced from seconds prior to hours after a binary neutron star (BNS) merger. The lar… ▽ More We present and evaluate the prospects for detecting coherent radio counterparts to gravitational wave (GW) events using Murchison Widefield Array (MWA) triggered observations. The MWA rapid-response system, combined with its buffering mode ($\sim4$ minutes negative latency), enables us to catch any radio signals produced from seconds prior to hours after a binary neutron star (BNS) merger. The large field of view of the MWA ($\sim1000\,\text{deg}^2$ at 120\,MHz) and its location under the high sensitivity sky region of the LIGO-Virgo-KAGRA (LVK) detector network, forecast a high chance of being on-target for a GW event. We consider three observing configurations for the MWA to follow up GW BNS merger events, including a single dipole per tile, the full array, and four sub-arrays. We then perform a population synthesis of BNS systems to predict the radio detectable fraction of GW events using these configurations. We find that the configuration with four sub-arrays is the best compromise between sky coverage and sensitivity as it is capable of placing meaningful constraints on the radio emission from 12.6\% of GW BNS detections. Based on the timescales of four BNS merger coherent radio emission models, we propose an observing strategy that involves triggering the buffering mode to target coherent signals emitted prior to, during or shortly following the merger, which is then followed by continued recording for up to three hours to target later time post-merger emission. We expect MWA to trigger on $\sim5\text{--}22$ BNS merger events during the LVK O4 observing run, which could potentially result in two detections of predicted coherent emission. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in PASA

Generative modeling, design and analysis of spider silk protein sequences for enhanced mechanical properties

Authors: Wei Lu, David L. Kaplan, Markus J. Buehler

Abstract: Spider silks are remarkable materials characterized by superb mechanical properties such as strength, extensibility and lightweightedness. Yet, to date, limited models are available to fully explore sequence-property relationships for analysis and design. Here we propose a custom generative large-language model to enable design of novel spider silk protein sequences to meet complex combinations of… ▽ More Spider silks are remarkable materials characterized by superb mechanical properties such as strength, extensibility and lightweightedness. Yet, to date, limited models are available to fully explore sequence-property relationships for analysis and design. Here we propose a custom generative large-language model to enable design of novel spider silk protein sequences to meet complex combinations of target mechanical properties. The model, pretrained on a large set of protein sequences, is fine-tuned on ~1,000 major ampullate spidroin (MaSp) sequences for which associated fiber-level mechanical properties exist, to yield an end-to-end forward and inverse generative strategy. Performance is assessed through: (1), a novelty analysis and protein type classification for generated spidroin sequences through BLAST searches, (2) property evaluation and comparison with similar sequences, (3) comparison of molecular structures, as well as, and (4) a detailed sequence motif analyses. We generate silk sequences with property combinations that do not exist in nature, and develop a deep understanding the mechanistic roles of sequence patterns in achieving overarching key mechanical properties (elastic modulus, strength, toughness, failure strain). The model provides an efficient approach to expand the silkome dataset, facilitating further sequence-structure analyses of silks, and establishes a foundation for synthetic silk design and optimization. △ Less

Submitted 18 September, 2023; originally announced September 2023.

How to Detect an Astrophysical Nanohertz Gravitational-Wave Background

Authors: Bence Bécsy, Neil J. Cornish, Patrick M. Meyers, Luke Zoltan Kelley, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Tyler Cohen, James M. Cordes, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch , et al. (71 additional authors not shown)

Abstract: Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam… ▽ More Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them. △ Less

Submitted 1 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 14 pages, 8 figures, version matching published paper

Journal ref: ApJ 959 9 (2023)

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

Characterizing Pulsars Detected in the Rapid ASKAP Continuum Survey

Authors: Akash Anumarlapudi, Anna Ehlke, Megan L. Jones, David L. Kaplan, Dougal Dobie, Emil Lenc, James K. Leung, Tara Murphy, Joshua Pritchard, Adam J. Stewart, Rahul Sengar, Craig Anderson, Julie Banfield, George Heald, Aidan W. Hotan, David McConnell, Vanessa A. Moss, Wasim Raja, Matthew T. Whiting

Abstract: We present the detection of 661 known pulsars observed with the Australian SKA Pathfinder (ASKAP) telescope at 888 MHz as a part of the Rapid ASKAP Continuum Survey (RACS). Detections were made through astrometric coincidence and we estimate the false alarm rate of our sample to be ~0.5%. Using archival data at 400 and 1400 MHz, we estimate the power law spectral indices for the pulsars in our sam… ▽ More We present the detection of 661 known pulsars observed with the Australian SKA Pathfinder (ASKAP) telescope at 888 MHz as a part of the Rapid ASKAP Continuum Survey (RACS). Detections were made through astrometric coincidence and we estimate the false alarm rate of our sample to be ~0.5%. Using archival data at 400 and 1400 MHz, we estimate the power law spectral indices for the pulsars in our sample and find that the mean spectral index is -1.78 +/- 0.6. However, we also find that a single power law is inadequate to model all the observed spectra. With the addition of the flux densities between 150 MHz and 3 GHz from various imaging surveys, we find that up to 40% of our sample shows deviations from a simple power law model. Using Stokes V measurements from the RACS data, we measured the circular polarization fraction for 9% of our sample and find that the mean polarization fraction is ~10% (consistent between detections and upper limits). Using the dispersion measure (DM) derived distance we estimate the pseudo luminosity of the pulsars and do not find any strong evidence for a correlation with the pulsars' intrinsic properties. △ Less

Submitted 31 July, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ (18 pages, 16 figures), comments are welcome

The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis , et al. (65 additional authors not shown)

Abstract: We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a sign… ▽ More We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately $3.3 \times 10^{-14}$. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613$-$0200. This suggests that this pulsar has some excess noise which can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 29 pages, 5 figures

The NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline

Authors: Aaron D. Johnson, Patrick M. Meyers, Paul T. Baker, Neil J. Cornish, Jeffrey S. Hazboun, Tyson B. Littenberg, Joseph D. Romano, Stephen R. Taylor, Michele Vallisneri, Sarah J. Vigeland, Ken D. Olum, Xavier Siemens, Justin A. Ellis, Rutger van Haasteren, Sophie Hourihane, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde , et al. (71 additional authors not shown)

Abstract: This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P… ▽ More This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data. △ Less

Submitted 7 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore , et al. (68 additional authors not shown)

Abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these… ▽ More The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Paul T. Baker, Bence Bécsy, Laura Blecha, Alexander Bonilla, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler , et al. (89 additional authors not shown)

Abstract: The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul… ▽ More The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations. △ Less

Submitted 18 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Authors: Adeela Afzal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Jose Juan Blanco-Pillado, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie , et al. (98 additional authors not shown)

Abstract: The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string… ▽ More The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. Decesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca , et al. (66 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated… ▽ More Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars

Authors: Gabriella Agazie, Md Faisal Alam, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Victoria Bonidie, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, Christopher Chapman, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler , et al. (75 additional authors not shown)

Abstract: We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival… ▽ More We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler, Megan E. DeCesar , et al. (89 additional authors not shown)

Abstract: We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr… ▽ More We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Periodic Radio Emission from the T8 Dwarf WISE J062309.94-045624.6

Authors: Kovi Rose, Joshua Pritchard, Tara Murphy, Manisha Caleb, Dougal Dobie, Laura Driessen, Stefan W. Duchesne, David L. Kaplan, Emil Lenc, Ziteng Wang

Abstract: We present the detection of rotationally modulated, circularly polarized radio emission from the T8 brown dwarf WISE J062309.94-045624.6 between 0.9 and 2.0 GHz. We detected this high proper motion ultracool dwarf with the Australian SKA Pathfinder in $1.36$ GHz imaging data from the Rapid ASKAP Continuum Survey. We observed WISE J062309.94-045624.6 to have a time and frequency averaged Stokes I f… ▽ More We present the detection of rotationally modulated, circularly polarized radio emission from the T8 brown dwarf WISE J062309.94-045624.6 between 0.9 and 2.0 GHz. We detected this high proper motion ultracool dwarf with the Australian SKA Pathfinder in $1.36$ GHz imaging data from the Rapid ASKAP Continuum Survey. We observed WISE J062309.94-045624.6 to have a time and frequency averaged Stokes I flux density of $4.17\pm0.41$ mJy beam$^{-1}$, with an absolute circular polarization fraction of $66.3\pm9.0\%$, and calculated a specific radio luminosity of $L_ν\sim10^{14.8}$ erg s$^{-1}$ Hz$^{-1}$. In follow-up observations with the Australian Telescope Compact Array and MeerKAT we identified a multi-peaked pulse structure, used dynamic spectra to place a lower limit of $B>0.71$ kG on the dwarf's magnetic field, and measured a $P=1.912\pm0.005$ h periodicity which we concluded to be due to rotational modulation. The luminosity and period we measured are comparable to those of other ultracool dwarfs observed at radio wavelengths. This implies that future megahertz to gigahertz surveys, with increased cadence and improved sensitivity, are likely to detect similar or later-type dwarfs. Our detection of WISE J062309.94-045624.6 makes this dwarf the coolest and latest-type star observed to produce radio emission. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ Letters; 11 pages, 3 figures and 2 tables

Radio Variable and Transient Sources on Minute Timescales in the ASKAP Pilot Surveys

Authors: Yuanming Wang, Tara Murphy, Emil Lenc, Louis Mercorelli, Laura Driessen, Joshua Pritchard, Baoqiang Lao, David L. Kaplan, Tao An, Keith W. Bannister, George Heald, 5 Shuoying Lu, Artem Tuntsov, Mark Walker, Andrew Zic

Abstract: We present results from a radio survey for variable and transient sources on 15-min timescales, using the Australian SKA Pathfinder (ASKAP) pilot surveys. The pilot surveys consist of 505 h of observations conducted at around 1 GHz observing frequency, with a total sky coverage of 1476 deg$^2$. Each observation was tracked for approximately 8-10h, with a typical rms sensitivity of $\sim$30 $μ$jy/b… ▽ More We present results from a radio survey for variable and transient sources on 15-min timescales, using the Australian SKA Pathfinder (ASKAP) pilot surveys. The pilot surveys consist of 505 h of observations conducted at around 1 GHz observing frequency, with a total sky coverage of 1476 deg$^2$. Each observation was tracked for approximately 8-10h, with a typical rms sensitivity of $\sim$30 $μ$jy/beam and an angular resolution of $\sim$12 arcsec. The variability search was conducted within each 8-10h observation on a 15-min timescale. We detected 38 variable and transient sources. Seven of them are known pulsars, including an eclipsing millisecond pulsar, PSR J2039$-$5617. Another eight sources are stars, only one of which has been previously identified as a radio star. For the remaining 23 objects, 22 are associated with active galactic nuclei or galaxies (including the five intra-hour variables that have been reported previously), and their variations are caused by discrete, local plasma screens. The remaining source has no multi-wavelength counterparts and is therefore yet to be identified. This is the first large-scale radio survey for variables and transient sources on minute timescales at a sub-mJy sensitivity level. We expect to discover $\sim$1 highly variable source per day using the same technique on the full ASKAP surveys. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 20 pages, 11 figures; accepted for publication in MNRAS

A matched-filter approach to radio variability and transients: searching for orphan afterglows in the VAST Pilot Survey

Authors: James K. Leung, Tara Murphy, Emil Lenc, Philip G. Edwards, Giancarlo Ghirlanda, David L. Kaplan, Andrew O'Brien, Ziteng Wang

Abstract: Radio transient searches using traditional variability metrics struggle to recover sources whose evolution timescale is significantly longer than the survey cadence. Motivated by the recent observations of slowly evolving radio afterglows at gigahertz frequency, we present the results of a search for radio variables and transients using an alternative matched-filter approach. We designed our match… ▽ More Radio transient searches using traditional variability metrics struggle to recover sources whose evolution timescale is significantly longer than the survey cadence. Motivated by the recent observations of slowly evolving radio afterglows at gigahertz frequency, we present the results of a search for radio variables and transients using an alternative matched-filter approach. We designed our matched-filter to recover sources with radio light curves that have a high-significance fit to power-law and smoothly broken power-law functions; light curves following these functions are characteristic of synchrotron transients, including "orphan" gamma-ray burst afterglows, which were the primary targets of our search. Applying this matched-filter approach to data from Variables and Slow Transients Pilot Survey conducted using the Australian SKA Pathfinder, we produced five candidates in our search. Subsequent Australia Telescope Compact Array observations and analysis revealed that: one is likely a synchrotron transient; one is likely a flaring active galactic nucleus, exhibiting a flat-to-steep spectral transition over $4\,$months; one is associated with a starburst galaxy, with the radio emission originating from either star formation or an underlying slowly evolving transient; and the remaining two are likely extrinsic variables caused by interstellar scintillation. The synchrotron transient, VAST J175036.1$-$181454, has a multi-frequency light curve, peak spectral luminosity, and volumetric rate that is consistent with both an off-axis afterglow and an off-axis tidal disruption event; interpreted as an off-axis afterglow would imply an average inverse beaming factor $\langle f^{-1}_{\text{b}} \rangle = 860^{+1980}_{-710}$, or equivalently, an average jet opening angle of $\langle θ_{\textrm{j}} \rangle = 3^{+4}_{-1}\,$deg. △ Less

Submitted 15 June, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 20 pages, 6 figures; accepted for publication in MNRAS

The Green Bank North Celestial Cap Survey. VIII. 21 New Pulsar Timing Solutions

Authors: William Fiore, Lina Levin, Maura A. McLaughlin, Akash Anumarlapudi, David L. Kaplan, Joseph K. Swiggum, Gabriella Y. Agazie, Robert Bavisotto, Pragya Chawla, Megan E. DeCesar, Timothy Dolch, Emmanuel Fonseca, Victoria M. Kaspi, Zachary Komassa, Vlad I. Kondratiev, Joeri van Leeuwen, Evan F. Lewis, Ryan S. Lynch, Alexander E. McEwen, Rusty Mundorf, Hind Al Noori, Emilie Parent, Ziggy Pleunis, Scott M. Ransom, Xavier Siemens , et al. (4 additional authors not shown)

Abstract: We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957-0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528-day… ▽ More We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957-0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528-day orbit, PSR J0032+6946 joins a small group of five other mildly recycled wide binary pulsars, for which the duration of recycling through accretion is limited by the length of the companion's giant phase. PSRs J0141+6303 and J1327+3423 are new disrupted recycled pulsars. We incorporate Arecibo observations from the NANOGrav pulsar timing array into our analysis of the latter. We also observed PSR J1327+3423 with the Long Wavelength Array, and our data suggest a frequency-dependent dispersion measure. PSR J0957-0619 was discovered as a rotating radio transient, but is a nulling pulsar at 820 MHz. PSR J1239+3239 is a new millisecond pulsar (MSP) in a 4-day orbit with a low-mass companion. Four of our pulsars already have published timing solutions, which we update in this work: the recycled wide binary PSR J0214+5222, the non-eclipsing black widow PSR J0636+5128, the disrupted recycled pulsar J1434+7257, and the eclipsing binary MSP J1816+4510, which is in an 8.7 hr orbit with a redback-mass companion. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 32 pages, 17 figures, 9 tables. Submitted to ApJ

Constraints on Undetected Long-Period Binaries in the Known Pulsar Population

Authors: Megan L. Jones, David L. Kaplan, Maura A. McLaughlin, Duncan R. Lorimer

Abstract: Although neutron star-black hole binaries have been identified through mergers detected in gravitational waves, a pulsar-black hole binary has yet to be detected. While short-period binaries are detectable due to a clear signal in the pulsar's timing residuals, effects from a long-period binary could be masked by other timing effects, allowing them to go undetected. In particular, a long-period bi… ▽ More Although neutron star-black hole binaries have been identified through mergers detected in gravitational waves, a pulsar-black hole binary has yet to be detected. While short-period binaries are detectable due to a clear signal in the pulsar's timing residuals, effects from a long-period binary could be masked by other timing effects, allowing them to go undetected. In particular, a long-period binary measured over a small subset of its orbital period could manifest via time derivatives of the spin-frequency incompatible with isolated pulsar properties. We assess the possibility of pulsars having unknown companions in long-period binaries and put constraints on the range of binary properties that may remain undetected in current data, but that may be detectable with further observations. We find that for 35% of canonical pulsars with published higher order derivatives, the precision of measurements is not enough to confidently reject binarity (period greater than ~2 kyr), and that a black-hole binary companion could not be ruled out for a sample of pulsars without published constraints if the period is greater than 1 kyr. While we find no convincing cases in the literature, we put more stringent limits on orbital period and longitude of periastron for the few pulsars with published higher-order frequency derivatives (n greater than 3). We discuss the detectability of candidates and find that a sample pulsar in a 100 yr orbit could be detectable within 5-10 yr. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 9 pages, 5 figures. Accepted to ApJ

Classical Novae in the ASKAP Pilot Surveys

Authors: Ashna Gulati, Tara Murphy, David L. Kaplan, Roberto Soria, James K. Leung, Yuanming Wang, Joshua Pritchard, Emil Lenc, Stefan W. Duchesne, Andrew O'Brien

Abstract: We present a systematic search for radio counterparts of novae using the Australian Square Kilometer Array Pathfinder (ASKAP). Our search used the Rapid ASKAP Continuum Survey, which covered the entire sky south of declination $+41^{\circ}$ ($\sim34,000$ square degrees) at a central frequency of 887.5 MHz, the Variables and Slow Transients Pilot Survey, which covered $\sim5,000$ square degrees per… ▽ More We present a systematic search for radio counterparts of novae using the Australian Square Kilometer Array Pathfinder (ASKAP). Our search used the Rapid ASKAP Continuum Survey, which covered the entire sky south of declination $+41^{\circ}$ ($\sim34,000$ square degrees) at a central frequency of 887.5 MHz, the Variables and Slow Transients Pilot Survey, which covered $\sim5,000$ square degrees per epoch (887.5 MHz), and other ASKAP pilot surveys, which covered $\sim200-2000$ square degrees with 2-12 hour integration times. We crossmatched radio sources found in these surveys over a two-year period, from April 2019 to August 2021, with 440 previously identified optical novae, and found radio counterparts for four novae: V5668 Sgr, V1369 Cen, YZ Ret, and RR Tel. Follow-up observations with the Australian Telescope Compact Array confirm the ejecta thinning across all observed bands with spectral analysis indicative of synchrotron emission in V1369 Cen and YZ Ret. Our light-curve fit with the Hubble Flow model yields a value of $1.65\pm 0.17 \times 10^{-4} \rm \:M_\odot$ for the mass ejected in V1369 Cen. We also derive a peak surface brightness temperature of $250\pm80$ K for YZ Ret. Using Hubble Flow model simulated radio lightcurves for novae, we demonstrate that with a 5$σ$ sensitivity limit of 1.5 mJy in 15-min survey observations, we can detect radio emission up to a distance of 4 kpc if ejecta mass is in the range $10^{-3}\rm \:M_\odot$, and upto 1 kpc if ejecta mass is in the range $10^{-5}-10^{-3}\rm \:M_\odot$. Our study highlights ASKAP's ability to contribute to future radio observations for novae within a distance of 1 kpc hosted on white dwarfs with masses $0.4-1.25\:\rm M_\odot$ , and within a distance of 4 kpc hosted on white dwarfs with masses $0.4-1.0\:\rm M_\odot$. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: This paper has been accepted for publication in PASA. It consists of 13 pages, 5 figures and 4 tables

A LOFAR sample of luminous compact sources coincident with nearby dwarf galaxies

Authors: D. Vohl, H. K. Vedantham, J. W. T. Hessels, C. G. Bassa, D. O. Cook, D. L. Kaplan, T. W. Shimwell, C. Zhang

Abstract: The vast majority of extragalactic compact continuum radio sources are associated with star formation or jets from (super)massive black holes and, as such, are more likely to be found in association with starburst galaxies or early-type galaxies. Two new populations of radio sources were recently identified: (a) compact and persistent sources (PRSs) associated with fast radio bursts (FRBs) in dwar… ▽ More The vast majority of extragalactic compact continuum radio sources are associated with star formation or jets from (super)massive black holes and, as such, are more likely to be found in association with starburst galaxies or early-type galaxies. Two new populations of radio sources were recently identified: (a) compact and persistent sources (PRSs) associated with fast radio bursts (FRBs) in dwarf galaxies and (b) compact sources in dwarf galaxies that could belong to the long-sought population of intermediate-mass black holes. Despite the interesting aspects of these newly found sources, the current sample size is small, limiting scrutiny of the underlying population. Here, we present a search for compact radio sources coincident with dwarf galaxies. We search the LOFAR Two-meter Sky Survey (LoTSS), the most sensitive low-frequency (144 MHz central frequency) large-area survey for optically thin synchrotron emission to date. Exploiting the high spatial resolution (6 arcsec) and low astrometric uncertainty (~0.2 arcsec) of the LoTSS, we match its compact sources to the compiled sample of dwarf galaxies in the Census of the Local Universe, an H alpha survey with the Palomar Observatory 48 inch Samuel Oschin Telescope. We identify 29 over-luminous compact radio sources, evaluate the probability of chance alignment within the sample, investigate the potential nature of these sources, and evaluate their volumetric density and volumetric rate. While optical line-ratio diagnostics on the nebular lines from the host galaxies support a star-formation origin rather than an AGN origin, future high-angular-resolution radio data are necessary to ascertain the origin of the radio sources. We discuss planned strategies to differentiate between candidate FRB hosts and intermediate-mass black holes. △ Less

Submitted 24 October, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: Accepted to A&A

Journal ref: A&A 680, A98 (2023)

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

The Southern-sky MWA Rapid Two-metre (SMART) pulsar survey -- II. Survey status, pulsar census, and first pulsar discoveries

Authors: N. D. R. Bhat, N. A. Swainston, S. J. McSweeney, M. Xue, B. W. Meyers, S. Kudale, S. Dai, S. E. Tremblay, W. van Straten, R. M. Shannon, K. R. Smith, M. Sokolowski, S. M. Ord, G. Sleap, A. Williams, P. J. Hancock, R. Lange, J. Tocknell, M. Johnston-Hollitt, D. L. Kaplan, S. J. Tingay, M. Walker

Abstract: In Paper I, we presented an overview of the Southern-sky MWA Rapid Two-metre (SMART) survey, including the survey design and search pipeline. While the combination of MWA's large field-of-view and the voltage capture system brings a survey speed of ~450 square degrees per hour, the survey progression relies on the availability of compact configuration of the Phase II array. Over the past few years… ▽ More In Paper I, we presented an overview of the Southern-sky MWA Rapid Two-metre (SMART) survey, including the survey design and search pipeline. While the combination of MWA's large field-of-view and the voltage capture system brings a survey speed of ~450 square degrees per hour, the survey progression relies on the availability of compact configuration of the Phase II array. Over the past few years, by taking advantage of multiple windows of opportunity when the compact configuration was available, we have advanced the survey to 75% completion. To date, about 10% of the data collected thus far have been processed for a first-pass search, where 10 minutes of observation is processed for dispersion measures out to 250 ${\rm pc\,cm^{-3}}$, to realise a shallow survey for long-period pulsars. The ongoing analysis has led to two new pulsar discoveries, as well as an independent discovery and a rediscovery of a previously incorrectly characterised pulsar, all from ~3% of the data for which candidate scrutiny is completed. Here we describe the strategies for further detailed follow-up including improved sky localisation and convergence to timing solution, and illustrate them using example pulsar discoveries. The processing has also led to re-detection of 120 pulsars in the SMART observing band, bringing the total number of pulsars detected to date with the MWA to 180, and these are used to assess the search sensitivity of current processing pipelines. The planned second-pass (deep survey) processing is expected to yield a three-fold increase in sensitivity for long-period pulsars, and a substantial improvement to millisecond pulsars by adopting optimal de-dispersion plans. The SMART survey will complement the highly successful Parkes High Time Resolution Universe survey at 1.2-1.5 GHz, and inform future large survey efforts such as those planned with the low-frequency Square Kilometre Array. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 22 pages, 9 figures, 7 tables, Accepted for publication in PASA

The Southern-sky MWA Rapid Two-metre (SMART) pulsar survey -- I. Survey design and processing pipeline

Authors: N. D. R. Bhat, N. A. Swainston, S. J. McSweeney, M. Xue, B. W. Meyers, S. Kudale, S. Dai, S. E. Tremblay, W. van Straten, R. M. Shannon, K. R. Smith, M. Sokolowski, S. M. Ord, G. Sleap, A. Williams, P. J. Hancock, R. Lange, J. Tocknell, M. Johnston-Hollitt, D. L. Kaplan, S. J. Tingay, M. Walker

Abstract: We present an overview of the Southern-sky MWA Rapid Two-metre (SMART) pulsar survey that exploits the MWA's large field of view and voltage capture system to survey the sky south of 30 degree in declination for pulsars and fast transients in the 140-170 MHz band. The survey is enabled by the advent of the Phase II MWA's compact configuration, which offers an enormous efficiency in beam-forming an… ▽ More We present an overview of the Southern-sky MWA Rapid Two-metre (SMART) pulsar survey that exploits the MWA's large field of view and voltage capture system to survey the sky south of 30 degree in declination for pulsars and fast transients in the 140-170 MHz band. The survey is enabled by the advent of the Phase II MWA's compact configuration, which offers an enormous efficiency in beam-forming and processing costs, thereby making an all-sky survey of this magnitude tractable with the MWA. Even with the long dwell times of the survey (4800 s), data collection can be completed in < 100 hours of telescope time, while still retaining the ability to reach a limiting sensitivity of ~2-3 mJy. Each observation is processed to generate ~5000-8000 tied-array beams that tessellate the full ~610 square degree field of view, which are then processed to search for pulsars. The voltage-capture recording allows a multitude of post hoc processing options including the reprocessing of data for higher time resolution. Due to the substantial computational cost in pulsar searches at low frequencies, processing is undertaken in multiple passes: in the first pass, a shallow survey is performed, where 10 minutes of each observation is processed, reaching about one-third of the full search sensitivity. Here we present the system overview and initial results. Further details including first pulsar discoveries and a census of low-frequency detections are presented in a companion paper. Future plans include deeper searches to reach the full sensitivity and acceleration searches to target binary and millisecond pulsars. Simulation analysis forecasts ~300 new pulsars upon the completion of full processing. The SMART survey will also generate a complete digital record of the low-frequency sky, which will serve as a valuable reference for future pulsar searches planned with the low-frequency Square Kilometre Array. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 22 pages, 12 figures, 2 tables, Accepted for publication in PASA

Collapsars as Sites of r-process Nucleosynthesis: Systematic Near-Infrared Follow-up of Type Ic-BL Supernovae

Authors: Shreya Anand, Jennifer Barnes, Sheng Yang, Mansi M. Kasliwal, Michael W. Coughlin, Jesper Sollerman, Kishalay De, Christoffer Fremling, Alessandra Corsi, Anna Y. Q. Ho, Arvind Balasubramanian, Conor Omand, Gokul P. Srinivasaragavan, S. Bradley Cenko, Tomas Ahumada, Igor Andreoni, Aishwarya Dahiwale, Kaustav Kashyap Das, Jacob Jencson, Viraj Karambelkar, Harsh Kumar, Brian D. Metzger, Daniel Perley, Nikhil Sarin, Tassilo Schweyer , et al. (19 additional authors not shown)

Abstract: One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star… ▽ More One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of $r$-process nucleosynthesis in the binary neutron star merger GW170817 was its long-lasting near-infrared emission, thus motivating a systematic photometric study of the light curves of broadlined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL -- including 18 observed with the Zwicky Transient Facility and 7 from the literature -- in the optical/near-infrared bands to determine what quantity of $r$-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for $r$-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the $r$-process mass for these SNe. We also perform independent light curve fits to models without $r$-process. We find that the $r$-process-free models are a better fit to the light curves of the objects in our sample. Thus we find no compelling evidence of $r$-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of $r$-process ejecta mass or indicate whether all collapsars are completely devoid of $r$-process nucleosynthesis. △ Less

Submitted 12 February, 2024; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 37 pages, 8 figures, 4 tables

Journal ref: ApJ, 962, 68 (2024)

A Pilot Study of Nulling in 22 Pulsars Using Mixture Modeling

Authors: Akash Anumarlapudi, Joseph K. Swiggum, David L. Kaplan, Travis D. J. Fichtenbauer

Abstract: The phenomenon of pulsar nulling, observed as the temporary inactivity of a pulsar, remains poorly understood both observationally and theoretically. Most observational studies that quantify nulling employ a variant of Ritchings (1976)'s algorithm which can suffer significant biases for pulsars where the emission is weak. Using a more robust mixture model method, we study pulsar nulling in a sampl… ▽ More The phenomenon of pulsar nulling, observed as the temporary inactivity of a pulsar, remains poorly understood both observationally and theoretically. Most observational studies that quantify nulling employ a variant of Ritchings (1976)'s algorithm which can suffer significant biases for pulsars where the emission is weak. Using a more robust mixture model method, we study pulsar nulling in a sample of 22 recently discovered pulsars, for which we publish the nulling fractions for the first time. These data clearly demonstrate biases of the former approach and show how an otherwise non-nulling pulsar can be classified as having significant nulls. We show that the population-wide studies that find a positive correlation of nulling with pulsar period/characteristic age can similarly be biased because of the bias in estimating the nulling fraction. We use our probabilistic approach to find the evidence for periodicity in the nulls in a subset of three pulsars in our sample. In addition, we also provide improved timing parameters for 17 of the 22 pulsars that had no prior follow-up. △ Less

Submitted 13 September, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: Accepted for publication in the Astrophysical Journal (ApJ)

The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, E. C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman , et al. (53 additional authors not shown)

Abstract: Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data s… ▽ More Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of $h_0 < $ $(6.82 \pm 0.35) \times 10^{-15}$, and $h_0 <$ $(2.66 \pm 0.15) \times 10^{-15}$ in our most sensitive sky location. Finally, we placed a multi-messenger limit of $\mathcal{M} <$ $(1.41 \pm 0.02) \times 10^9 M_\odot$ on the chirp mass of the supermassive black hole binary candidate 3C~66B. △ Less

Submitted 6 June, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures. Accepted by ApJL

The Green Bank North Celestial Cap Survey. VII. 12 New Pulsar Timing Solutions

Authors: Joseph K. Swiggum, Ziggy Pleunis, Emilie Parent, David L. Kaplan, Maura A. McLaughlin, Ingrid H. Stairs, Renée Spiewak, Gabriella Y. Agazie, Pragya Chawla, Megan E. DeCesar, Timothy Dolch, William Fiore, Emmanuel Fonseca, Alina G. Istrate, Victoria M. Kaspi, Vlad I. Kondratiev, Joeri van Leeuwen, Lina Levin, Evan F. Lewis, Ryan S. Lynch, Alex E. McEwen, Hind Al Noori, Scott M. Ransom, Xavier Siemens, Mayuresh Surnis

Abstract: We present timing solutions for 12 pulsars discovered in the Green Bank North Celestial Cap (GBNCC) 350 MHz pulsar survey, including six millisecond pulsars (MSPs), a double neutron star (DNS) system, and a pulsar orbiting a massive white dwarf companion. Timing solutions presented here include 350 and 820 MHz Green Bank Telescope data from initial confirmation and follow-up as well as a dedicated… ▽ More We present timing solutions for 12 pulsars discovered in the Green Bank North Celestial Cap (GBNCC) 350 MHz pulsar survey, including six millisecond pulsars (MSPs), a double neutron star (DNS) system, and a pulsar orbiting a massive white dwarf companion. Timing solutions presented here include 350 and 820 MHz Green Bank Telescope data from initial confirmation and follow-up as well as a dedicated timing campaign spanning one year. PSR J1122$-$3546 is an isolated MSP, PSRs J1221$-$0633 and J1317$-$0157 are MSPs in black widow systems and regularly exhibit eclipses, and PSRs J2022+2534 and J2039$-$3616 are MSPs that can be timed with high precision and have been included in pulsar timing array experiments seeking to detect low-frequency gravitational waves. PSRs J1221$-$0633 and J2039$-$3616 have Fermi Large Area Telescope $γ$-ray counterparts and also exhibit significant $γ$-ray pulsations. We measure proper motion for three of the MSPs in this sample and estimate their space velocities, which are typical compared to those of other MSPs. We have detected the advance of periastron for PSR J1018$-$1523 and therefore measure the total mass of the double neutron star system, $m_{\rm tot}=2.3\pm0.3$ M$_{\odot}$. Long-term pulsar timing with data spanning more than one year is critical for classifying recycled pulsars, carrying out detailed astrometry studies, and shedding light on the wealth of information in these systems post-discovery. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: 21 pages, 5 figures, 7 tables

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

An unusual pulse shape change event in PSR J1713+0747 observed with the Green Bank Telescope and CHIME

Authors: Ross J. Jennings, James M. Cordes, Shami Chatterjee, Maura A. McLaughlin, Paul B. Demorest, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Paul R. Brook, Tyler Cohen, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Timothy Dolch, Elizabeth C. Ferrara, Emmanuel Fonseca, Deborah C. Good, Jeffrey S. Hazboun, Megan L. Jones, David L. Kaplan, Michael T. Lam, T. Joseph W. Lazio, Duncan R. Lorimer, Jing Luo, Ryan S. Lynch , et al. (19 additional authors not shown)

Abstract: The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter G… ▽ More The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter Green Bank Telescope (GBT) in a three-year period encompassing the shape change event, between February 2020 and February 2023. As of February 2023, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying TOA residuals display a strong non-monotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency, $ν$) nor a change in dispersion measure (DM) alone (which would produce a delay proportional to $ν^{-2}$). However, it does bear some resemblance to the two previous "chromatic timing events" observed in J1713+0747 (Demorest et al. 2013; Lam et al. 2016), as well as to a similar event observed in PSR J1643-1224 in 2015 (Shannon et al. 2016). △ Less

Submitted 31 January, 2024; v1 submitted 21 October, 2022; originally announced October 2022.

Comments: 18 pages, 8 figures. Submitted to ApJ. Data available at https://doi.org/10.5281/zenodo.7236459