Search for: All records

ABSTRACT We present the deepest and highest-resolution radio continuum imaging of the Galactic globular cluster Terzan 5, one of the most crowded locations in the radio sky. In these new 2–4 GHz Karl G. Jansky Very Large Array images, we detect 38 of the 49 confirmed pulsars, including extensive multi-frequency eclipse mapping of the luminous redback Ter5A. Nonetheless, there is still a large amount of diffuse residual flux from pulsars that are fainter than our 2.5 GHz continuum detection limit of $$\sim 11\, \mu$$Jy. Using a range of approaches including image-based simulations, we model the fluxes of the detected pulsars together with the residual flux. We find a minimum total population of $$N\sim 250$$ detectable pulsars in Terzan 5 and perhaps substantially more, though the luminosity function remains very uncertain. Consideration of the $$\gamma$$-ray properties of the cluster, though also not unambiguous to interpret, leads to consistent conclusions. These pulsar population estimates are larger than inferred from previous work and highlight Terzan 5 as a keystone target for next-generation radio facilities. 

ABSTRACT Nulling is a phenomenon where the emission from a pulsar becomes undetectable (or significantly weaker) for a relatively short period of time, followed by a return to a normal emission state. The time-scale of nulling ranges from a few pulse periods to many hours or even days. The fraction of time a nulling pulsar spends in a null state varies across the population of canonical pulsars, from 0 to 95 per cent. The long-term behaviour of a pulsar’s nulling fraction, however, is currently unknown, as published values have typically been obtained through single observations. Here, we present the first long-term analysis of nulling behaviour in eight pulsars observed in the Parkes Multibeam Pulsar Survey over the course of eight to ten years. We also apply a new Bayesian method for pulse-energy analysis, yielding posterior estimates of the nulling fraction per observation. In several cases, the nulling affects only specific components of the pulse profile, rather than the entirety of the emission. Our analysis reveals that, while most pulsars show no significant trend in their nulling fraction over time, a subset exhibit some evidence for non-zero gradients in nulling fraction. In particular, PSRs J1048–3832, J1745–3040, and J1825–0935 show statistically significant trends over the span of the data. Studying the behaviour of nulling over years and decades is valuable as it can provide insights into the physical emission processes within pulsars. Studying how nulling evolves also provides valuable insights into pulsar evolution and the characterization of the broader pulsar population. 

Abstract Recent findings from several Pulsar Timing Array (PTA) collaborations point to the existence of a Gravitational Wave Background (GWB) at nanohertz frequencies. A key next step towards characterizing this signal and identifying its origin is to map the sky distribution of its power. Several strategies have been proposed to reconstruct this distribution using PTA data. In this work, we compare these different strategies to determine which one is best suited to detect GWB anisotropies of different topologies. We find that, for both localized and large-scale anisotropies, reconstruction methods based on pixel and radiometer maps are the most promising. However, in both scenarios, even the optimistically large anisotropic signals discussed in this work remain challenging to detect with near-future PTA sensitivities. For example, we find that for a GWB hotspot contributing to 80% of the GWB power in the second frequency bin, detection probabilities reach at most 𝒪(10%) for a PTA with noise properties comparable with the ones of the upcoming IPTA third data release. Finally, we consider the fundamental limitations that cosmic variance poses to these kinds of searches by deriving the smallest deviations from isotropy that could be detected by an idealized PTA with no experimental or pulsar noise. 

Abstract Cosmic strings represent an attractive source of gravitational waves (GWs) from the early Universe. However, numerical computation of the GW signal from cosmic strings requires the evaluation of complicated integral and sum expressions, which can become computationally costly in large parameter scans. This motivates us to rederive the GW signal from a network of local stable cosmic strings in the Nambu-Goto approximation and based on the velocity-dependent one-scale model from a “pedestrian” perspective. That is, we derive purely analytical expressions for the total GW spectrum, which remain exact wherever possible and whose error can be tracked and reduced in a controlled way in crucial situations in which we are forced to introduce approximations. In this way, we obtain powerful formulas that, unlike existing results in the literature, are valid across the entire frequency spectrum and across the entire conceivable range of cosmic-string tensions. We provide an in-depth discussion of the GW spectra thus obtained, including their characteristic break frequencies and approximate power-law behaviors, comment on the effect of changes in the effective number of degrees of freedom during radiation domination, and conclude with a concise summary of our main formulas that can readily be used in future studies. 

Abstract Spectroscopic observations of the quasar J0950+5128 spanning 22 yr reveal monotonic radial velocity variations in its broad Hβemission line. Moreover, the line profile becomes broader over time, necessitating careful measurements. We present robust Hβvelocity shift measurements obtained via cross correlation, applied to both the full spectra and to isolated broad Hβcomponents derived from spectral decomposition. We also examine the light curves for variability consistent with the spectroscopic trends. Using Lomb–Scargle periodogram analysis, we find no significant periodic signal. We consider several interpretations for the observed changes, including a binary supermassive black hole, dust-cloud obscuration, outflows, a recoiling black hole, and a single perturbed, disk-like broad-line region (BLR). We deem the binary and perturbed BLR scenarios to be physically plausible. The binary interpretation is the only one for which we can immediately compare a physical model to the available data. Thus, we incorporate radial velocity “jitter” to emulate typical quasar variability and fit the radial velocity curve with a Keplerian model to examine whether it can reproduce the observations. In this context, the available observations trace only a segment of the putative orbit. The fit yields a period of 33 yr (observed frame) and an eccentricity of 0.65, with lower limits on the semimajor axis and black hole mass of 10⁻² pc and 10⁷ M_⊙, respectively. Thus, J0950+5128 is a binary candidate deserving further study. The single, perturbed BLR interpretation remains viable but requires additional observations and modeling for further evaluation. Continued monitoring is, therefore, essential. 

Abstract Supermassive black hole binary (SMBHB) systems are expected to form as a consequence of galaxy mergers. At subparsec separations, SMBHBs can be identified as quasars with periodic variability, with previous periodicity searches uncovering significant candidates. However, these searches focused primarily on sinusoidal signals, while theoretical models and hydrodynamical simulations predict that binaries produce more complex non-sinusoidal pulse shapes. Here we examine the efficacy of the Lomb–Scargle periodogram (LSP; one of the most popular tools for periodicity searches in unevenly sampled lightcurves) to detect periodicities with a sawtooth shape mimicking results of hydrodynamical simulations. We simulate idealized well-sampled lightcurves, lightcurves that mimic the data in the Palomar Transient Factory (PTF) analyzed in M. Charisi et al. (2016), and lightcurves that resemble our expectations for single-band data in the upcoming Legacy Survey of Space and Time (LSST) of the Rubin Observatory. We approximate quasar variability with a damped random walk (DRW) model, inject sinusoidal and sawtooth pulse shapes, and assess their statistical significance. We find that in the presence of red noise, the LSP detects a relatively low fraction of the sinusoidal signals (∼45%, ∼24%, and ∼23%, in the PTF-like, idealized, and LSST-like lightcurves, respectively). The fraction is significantly reduced for sawtooth periodicity (with only ∼9% in PTF-like and ∼1% in idealized and LSST-like lightcurves). These low recovery rates imply that previous searches have missed the large majority of binaries. They also have significant implications for the detection of SMBHBs in upcoming LSST necessitating the development of advanced tools that go beyond the simple LSP. 

Abstract Binary supermassive black holes (SMBHs) are consequences of galaxy mergers and dominate the low-frequency gravitational-wave background. Finding binary SMBHs in existing time-domain observations has proven difficult, as their periodic, electromagnetic signals can be confused with the natural variability of single quasars. In this work, we investigate the effects of host-galaxy contamination and survey design (cadence and duration) on the detectability of binary SMBHs with the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We simulate millions of LSST light curves of single and binary quasars, with a distribution of quasar and host-galaxy properties motivated by empirical observations and the anticipated LSST detection space. We then apply simple sinusoidal curve fits as a potential computationally inexpensive detection method. We find that host-galaxy contamination will increase false-positive rates and decrease binary parameter recovery rates. Lower-mass, lower-luminosity binary systems are most likely to be negatively affected by host-galaxy contamination. We also find that monitoring duration affects binary detection more than survey effective cadence for this detection method. As the light-curve duration increases, false-positive rates are suppressed and binary parameter recovery rates, especially for binary periods, are improved. Increasing the light-curve duration from 5 to 10 yr shows the most dramatic improvement for successful binary detection and false-positive rejection, with additional improvement from extending the light-curve duration to 20 yr. The observation duration increase is especially critical for recovering binary periods that are longer than a decade. 

Abstract This review is focused on tests of Einstein’s theory of general relativity with gravitational waves that are detectable by ground-based interferometers and pulsar-timing experiments. Einstein’s theory has been greatly constrained in the quasi-linear, quasi-stationary regime, where gravity is weak and velocities are small. Gravitational waves are allowing us to probe a complimentary, yet previously unexplored regime: the non-linear and dynamicalextreme gravity regime. Such a regime is, for example, applicable to compact binaries coalescing, where characteristic velocities can reach fifty percent the speed of light and gravitational fields are large and dynamical. This review begins with the theoretical basis and the predicted gravitational-wave observables of modified gravity theories. The review continues with a brief description of the detectors, including both gravitational-wave interferometers and pulsar-timing arrays, leading to a discussion of the data analysis formalism that is applicable for such tests. The review then discusses gravitational-wave tests using compact binary systems, and ends with a description of the first gravitational wave observations by advanced LIGO, the stochastic gravitational wave background observations by pulsar timing arrays, and the tests that can be performed with them. 

ABSTRACT The ionized interstellar medium disperses pulsar radio signals, resulting in a stochastic time-variable delay known as the dispersion measure (DM) noise. In the wideband paradigm of pulsar timing, we measure a DM together with a time of arrival from a pulsar observation to handle frequency-dependent profile evolution, interstellar scintillation, and radio frequency interference more robustly, and to reduce data volumes. In this paper, we derive a method to incorporate arbitrary models of DM variation, including Gaussian process models, in pulsar timing and noise analysis and pulsar timing array analysis. This generalizes the existing method for handling DM noise in wideband data sets. 

Abstract We present timing solutions spanning nearly two decades for five redback (RB) systems found in globular clusters (GCs), created using a novel technique that effectively “isolates” the pulsar. By accurately measuring the time of passage through periastron (T₀) at points over the timing baseline, we use a piecewise-continuous binary model to get local solutions of the orbital variations that we pair with long-term orbital information to remove the orbital timing delays. The isolated pulse times of arrival can then be fit to describe the spin behavior of the millisecond pulsar (MSP). The results of our timing analyses via this method are consistent with those of conventional timing methods for binaries in GCs as demonstrated by analyses of NGC 6440D. We also investigate the observed orbital phase variations for these systems. Quasiperiodic oscillations in Terzan 5P’s orbit may be the result of changes to the gravitational quadruple moment of the companion as prescribed by the Applegate model. We find a striking correlation between the standard deviation of the phase variations as a fraction of a system’s orbit ( ${\sigma }_{\mathrm{\Delta }{T}_0}$) and the MSP’s spin frequency, as well as a potential correlation between ${\sigma }_{\mathrm{\Delta }{T}_0}$and the binary’s projected semimajor axis. While long-term RB timing is fraught with large systematics, our work provides a needed alternative for studying systems with significant orbital variations, especially when high-cadence monitoring observations are unavailable.