Search

The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index \gamma. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, \gamma_run(f) = \gamma + \beta \ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter \beta consistent with no running, \beta \in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero \beta. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA., Comment: 17 pages, 4 figures, 2 tables

Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $\sigma$ and 3\ $\sigma$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav., Comment: 20 pages, 14 Figures

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., Comment: Accepted for publication in ApJ, comments welcome

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\sigma$ 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.

Dark matter may exist today in the form of ultraheavy composite bound states. Collisions between such dark matter states can release intense bursts of radiation that includes gamma-rays among the final products. Thus, indirect-detection signals of dark matter may include unconventional gamma-ray bursts. Such bursts may have been missed not necessarily because of their low arriving gamma-ray fluxes, but rather their briefness and rareness. We point out that intense bursts whose non-detection thus far are due to the latter can be detected in the near future with existing and planned facilities. In particular, we propose that, with slight experimental adjustments and suitable data analyses, imaging atmospheric Cherenkov telescopes (IACTs) and Pulsed All-sky Near-infrared and Optical Search for Extra-Terrestrial Intelligence (PANOSETI) are promising tools for detecting such rare, brief, but intense bursts. We also show that if we assume these bursts originate from collisions of dark matter states, IACTs and PANOSETI can probe a large dark matter parameter space beyond existing limits. Additionally, we present a concrete model of dark matter that produces bursts potentially detectable in these instruments.

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., Comment: 28 pages, 10 figures, 6 tables. ApJ, submitted

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.

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$\sigma$ 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., Comment: 37 pages, 37 figures, accepted to ApJ

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 $\gamma_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $\gamma_{\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.

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., Comment: submitted

PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorrelated and correlated noise processes as well as the model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated datasets as well as a real dataset of PSR B1855+09. We also describe the new features implemented in PINT since it was first described in the literature., Comment: Accepted for publication in ApJ

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 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). 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., Comment: 10 pages, 8 figures, 1 appendix, submitted to ApJ

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., Comment: 13 pages, 8 figures, 2 tables. Accepted for publication in The Astrophysical Journal

The tentative identification of approximately ten relativistic anti-helium cosmic-ray events at AMS-02 would, if confirmed, challenge our understanding of the astrophysical synthesis of heavy anti-nuclei. We propose a novel scenario for the enhanced production of such anti-nuclei that is triggered by isolated, catastrophic injections of large quantities of energetic Standard Model (SM) anti-quarks in our galaxy by physics beyond the Standard Model (BSM). We demonstrate that SM anti-nucleosynthetic processes that occur in the resulting rapidly expanding, thermalized fireballs of SM plasma can, for a reasonable range of parameters, produce the reported tentative $\sim 2:1$ ratio of anti-helium-3 to anti-helium-4 events at AMS-02, as well as their relativistic boosts. Moreover, we show that this can be achieved without violating anti-deuterium or anti-proton flux constraints for the appropriate anti-helium fluxes. A plausible BSM paradigm for the catastrophic injections is the collision of macroscopic composite dark-matter objects carrying large net anti-baryon number. Such a scenario would require these objects to be cosmologically stable, but to destabilize upon collision, promptly releasing a fraction of their mass energy into SM anti-particles within a tiny volume. We show that, in principle, the injection rate needed to attain the necessary anti-helium fluxes and the energetic conditions required to seed the fireballs appear possible to obtain in such a paradigm. We leave open the question of constructing a BSM particle physics model to realize this, but we suggest two concrete scenarios as promising targets for further investigation., Comment: 42 pages, 10 figures, journal version

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., Comment: Accepted for publication in Science Advances. We report a new supermassive black hole phenomenon that we call quasi-periodic outflows (QPOuts)

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 $\tau \propto \lambda^{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., Comment: 11 pages, 5 figures, 4 tables. Accepted for publication in ApJ. The full dataset will be available with the ApJ article

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., Comment: 14 pages, 9 figures, submitted to MNRAS

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 $\gamma$-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 $\gamma$-ray pulsation searches, which may confirm the source nature as a pulsar., Comment: Submitted to MNRAS. 13 Pages, 10 figures, 3 tables

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., Comment: 11 pages, 5 figures

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 $\eta>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 $\eta>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., Comment: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

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., Comment: 14 pages, 8 figures, version matching published paper

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., Comment: Accepted for publication in ApJ (18 pages, 16 figures), comments are welcome

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., Comment: 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

Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems., Comment: 23 pages, 13 figures, 2 tables. Accepted for publication 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 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., Comment: 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 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., Comment: 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 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., Comment: 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

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., Comment: 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

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., Comment: 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

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\sigma$) 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\sigma$). 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., Comment: 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

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 $\mu$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., Comment: 20 pages, 11 figures; accepted for publication in MNRAS

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 \theta_{\textrm{j}} \rangle = 3^{+4}_{-1}\,$deg., Comment: 20 pages, 6 figures; accepted for publication in MNRAS

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., Comment: 32 pages, 17 figures, 9 tables. Submitted to ApJ

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., Comment: 9 pages, 5 figures. Accepted to ApJ

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$\sigma$ 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$., Comment: This paper has been accepted for publication in PASA. It consists of 13 pages, 5 figures and 4 tables

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., Comment: 37 pages, 8 figures, 4 tables

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., Comment: Accepted for publication in the Astrophysical Journal (ApJ)

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., Comment: 22 pages, 12 figures. Accepted by ApJL

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 $\gamma$-ray counterparts and also exhibit significant $\gamma$-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., Comment: 21 pages, 5 figures, 7 tables

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., Comment: Accepted for publication in Nature

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, $\nu$) nor a change in dispersion measure (DM) alone (which would produce a delay proportional to $\nu^{-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)., Comment: 18 pages, 8 figures. Submitted to ApJ. Data available at https://doi.org/10.5281/zenodo.7236459

We present the results of a radio transient and polarisation survey towards the Galactic Centre, conducted as part of the Australian Square Kilometre Array Pathfinder Variables and Slow Transients pilot survey. The survey region consisted of five fields covering $\sim265\,{\rm deg}^2$ ($350^\circ\lesssim l\lesssim10^\circ$, $\vert b\vert \lesssim 10^\circ$). Each field was observed for 12\,minutes, with between 7 and 9 repeats on cadences of between one day and four months. We detected eight highly variable sources and seven highly circularly-polarised sources (14 unique sources in total). Seven of these sources are known pulsars including the rotating radio transient PSR~J1739--2521 and the eclipsing pulsar PSR~J1723--2837. One of them is a low mass X-ray binary, 4U 1758--25. Three of them are coincident with optical or infrared sources and are likely to be stars. The remaining three may be related to the class of Galactic Centre Radio Transients (including a highly likely one, VAST~J173608.2--321634, that has been reported previously), although this class is not yet understood. In the coming years, we expect to detect $\sim$40 bursts from this kind of source with the proposed four-year VAST survey if the distribution of the source is isotropic over the Galactic fields., Comment: 17 pages, 13 figures, Accepted for publication in MNRAS

We report the independent discovery of PSR J0027-1956 with the Murchison Widefield Array (MWA) in the ongoing Southern-sky MWA Rapid Two-meter (SMART) pulsar survey. J0027-1956 has a period of ~1.306 s, a dispersion measure (DM) of ~20.869 pc cm^-3 , and a nulling fraction of ~77%. This pulsar highlights the advantages of the survey's long dwell times (~80 min), which, when fully searched, will be sensitive to the expected population of similarly bright, intermittent pulsars with long nulls. A single-pulse analysis in the MWA's 140-170 MHz band also reveals a complex sub-pulse drifting behavior, including both rapid changes of the drift rate characteristic of mode switching pulsars, as well as a slow, consistent evolution of the drift rate within modes. In some longer drift sequences, interruptions in the otherwise smooth drift rate evolution occur preferentially at a particular phase, typically lasting a few pulses. These properties make this pulsar an ideal test bed for prevailing models of drifting behavior such as the carousel model., Comment: Accepted for publication in ApJ

GW170817 is the first binary neutron star (NS) merger detected in gravitational waves (GWs) and photons, and so far remains the only GW event of its class with a definitive electromagnetic (EM) counterpart. Radio emission from the structured jet associated with GW170817 has faded below the sensitivity achievable via deep radio observations with the most sensitive radio arrays currently in operation. Hence, we now have the opportunity to probe the radio re-brightening that some models predict, should emerge at late times from the interaction of the dynamically-stripped merger ejecta with the interstellar medium. Here we present the latest results from our deep radio observations of the GW170817 field with the Karl G. Jansky Very Large Array (VLA), 4.5 years after the merger. Our new data at $3\,$GHz do not show any compelling evidence for emission in excess to the tail of the jet afterglow ($<3.3\,\mu$Jy), confirming our previous results. We thus set new constraints on the dynamical ejecta afterglow models. These constraints favor single-speed ejecta with energy $\lesssim 10^{50}\,$erg (for an ejecta speed of $\beta_0=0.5$), or steeper energy-speed distributions of the kilonova ejecta. Our results also suggest larger values of the cold, non-rotating maximum NS mass in equal mass scenarios. However, without a detection of the dynamical ejecta afterglow, obtaining precise constraints on the NS equation of state remains challenging., Comment: 6 pages, 5 figures. Published in ApJ

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution., Comment: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022

We report the discovery of a highly circularly polarized, variable, steep-spectrum pulsar in the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. The pulsar is located about $1^\circ$ from the center of the Large Magellanic Cloud, and has a significant fractional circular polarization of $\sim$20%. We discovered pulsations with a period of 322.5 ms, dispersion measure (DM) of 157.5 pc cm$^{-3}$, and rotation measure (RM) of $+456$ rad m$^{-2}$ using observations from the MeerKAT and the Parkes telescopes. This DM firmly places the source, PSR J0523$-$7125, in the Large Magellanic Cloud (LMC). This RM is extreme compared to other pulsars in the LMC (more than twice that of the largest previously reported one). The average flux density of $\sim$1 mJy at 1400 MHz and $\sim$25 mJy at 400 MHz places it among the most luminous radio pulsars known. It likely evaded previous discovery because of its very steep radio spectrum (spectral index $\alpha \approx -3$, where $S_\nu \propto \nu^\alpha$) and broad pulse profile (duty cycle $\gtrsim35$%). We discuss implications for searches for unusual radio sources in continuum images, as well as extragalactic pulsars in the Magellanic Clouds and beyond. Our result highlighted the possibility of identifying pulsars, especially extreme pulsars, from radio continuum images. Future large-scale radio surveys will give us an unprecedented opportunity to discover more pulsars and potentially the most distant pulsars beyond the Magellanic Clouds., Comment: 17 pages, 7 figures; accepted by ApJ

The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that $-21$ can contribute no more than 87% of the total, while transients brighter than $-22$ can contribute no more than 58% of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26% of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory.

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

Observations of the extragalactic ($z=0.0141$) transient AT2018cow established a new class of energetic explosions shocking a dense medium, which produce luminous emission at millimeter and sub-millimeter wavelengths. Here we present detailed millimeter- through centimeter-wave observations of a similar transient, ZTF20acigmel (AT2020xnd) at $z=0.2433$. Using observations from the NOrthern Extended Millimeter Array and the Very Large Array, we model the unusual millimeter and radio emission from AT2020xnd under several different assumptions, and ultimately favor synchrotron radiation from a thermal electron population (relativistic Maxwellian). The thermal-electron model implies a fast but sub-relativistic ($v\approx0.3c$) shock and a high ambient density ($n_e\approx4\times10^{3}$cm$^{-3}$ at $\Delta t\approx40$ days). The X-ray luminosity of $L_X\approx10^{43}$ erg sec$^{-1}$ exceeds simple predictions from the radio and UVOIR luminosity and likely has a separate physical origin, such as a central engine. Using the fact that month-long luminous ($L_\nu\approx 2\times10^{30}$ erg sec$^{-1}$ Hz$^{-1}$ at 100 GHz) millimeter emission appears to be a generic feature of transients with fast ($t_{1/2}\approx3$ days) and luminous ($M_\mathrm{peak}\approx -21 $mag) optical light curves, we estimate the rate at which transients like AT2018cow and AT2020xnd will be detected by future wide-field millimeter transient surveys such as CMB-S4, and conclude that energetic explosions in dense environments may represent a significant population of extragalactic transients in the 100 GHz sky., Comment: 27 pages, 13 figures. Accepted for publication in ApJ on 23 Jan 2022

We search NANOGrav's 12.5-year data set for evidence of a gravitational wave background (GWB) with all the spatial correlations allowed by general metric theories of gravity. We find no substantial evidence in favor of the existence of such correlations in our data. We find that scalar-transverse (ST) correlations yield signal-to-noise ratios and Bayes factors that are higher than quadrupolar (tensor transverse, TT) correlations. Specifically, we find ST correlations with a signal-to-noise ratio of 2.8 that are preferred over TT correlations (Hellings and Downs correlations) with Bayesian odds of about 20:1. However, the significance of ST correlations is reduced dramatically when we include modeling of the Solar System ephemeris systematics and/or remove pulsar J0030$+$0451 entirely from consideration. Even taking the nominal signal-to-noise ratios at face value, analyses of simulated data sets show that such values are not extremely unlikely to be observed in cases where only the usual TT modes are present in the GWB. In the absence of a detection of any polarization mode of gravity, we place upper limits on their amplitudes for a spectral index of $\gamma = 5$ and a reference frequency of $f_\text{yr} = 1 \text{yr}^{-1}$. Among the upper limits for eight general families of metric theories of gravity, we find the values of $A^{95\%}_{TT} = (9.7 \pm 0.4)\times 10^{-16}$ and $A^{95\%}_{ST} = (1.4 \pm 0.03)\times 10^{-15}$ for the family of metric spacetime theories that contain both TT and ST modes., Comment: 24 pages, 18 figures, 3 appendices. Please send any comments/questions to Nima Laal (laaln@oregonstate.edu)