Search

The stochastic gravitational wave background (SGWB) can be observed in the nanohertz band using a pulsar timing array (PTA). Here a computationally efficient state-space framework is developed for analysing SGWB data, in which the stochastic gravitational wave strain at Earth is tracked with a non-linear Kalman filter and separated simultaneously from intrinsic, achromatic pulsar spin wandering. The filter is combined with a nested sampler to estimate the parameters of the model, and to calculate a Bayes factor for selecting between models with and without a SGWB. The procedure extends previous state-space formulations of PTA data analysis applied to individually resolvable binary black hole sources. The performance of the new algorithm is tested on synthetic data from the first International PTA Mock Data Challenge. It is shown that the algorithm distinguishes a SGWB from pure noise for $A_{\rm gw} \geq 3 \times 10^{-14}$, where $A_{\rm gw}$ denotes the standard normalization factor for a power spectral density with power-law exponent $-13/3$. Additional, systematic validation tests are also performed with synthetic data generated independently by adjusting the injected parameters to cover astrophysically plausible ranges. Full posterior distributions are recovered and tested for accuracy. The state-space procedure is memory-light and evaluates the likelihood for a standard-sized PTA dataset in $\lesssim 10^{-1}$ s without optimization on a standard central processing unit., Comment: 10 pages, 4 figures + appendices. Accepted for publication in MNRAS; v2 corrected references

Pulsar timing arrays can detect continuous nanohertz gravitational waves emitted by individual supermassive black hole binaries. The data analysis procedure can be formulated within a time-domain, state-space framework, in which the radio timing observations are related to a temporal sequence of latent states, namely the intrinsic pulsar spin frequency. The achromatic wandering of the pulsar spin frequency is tracked using a Kalman filter concurrently with the pulse frequency modulation induced by a gravitational wave from a single source. The modulation is the sum of terms proportional to the gravitational wave strain at the Earth and at every pulsar in the array. Here we generalize previous state-space formulations of the pulsar timing array problem to include the pulsar terms; that is, we copy the pulsar terms from traditional, non-state-space analyses over to the state-space framework. The performance of the generalized Kalman filter is tested using astrophysically representative software injections in Gaussian measurement noise. It is shown that including the pulsar terms corrects for previously identified biases in the parameter estimates (especially the sky position of the source) which also arise in traditional matched-filter analyses that exclude the pulsar terms. Additionally, including the pulsar terms decreases the minimum detectable strain by $14\%$. Overall, the study verifies that the pulsar terms do not raise any special extra impediments for the state-space framework, beyond those studied in traditional analyses. The inspiral-driven evolution of the wave frequency at the Earth and at the retarded time at every pulsar in the array is also investigated., Comment: 24 pages, 13 figures. Accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2409.14613

Continuous nanohertz gravitational waves from individual supermassive black hole binaries may be detectable with pulsar timing arrays. A novel search strategy is developed, wherein intrinsic achromatic spin wandering is tracked simultaneously with the modulation induced by a single gravitational wave source in the pulse times of arrival. A two-step inference procedure is applied within a state-space framework, such that the modulation is tracked with a Kalman filter, which then provides a likelihood for nested sampling. The procedure estimates the static parameters in the problem, such as the sky position of the source, without fitting for ensemble-averaged statistics such as the power spectral density of the timing noise, and therefore complements traditional parameter estimation methods. It also returns the Bayes factor relating a model with a single gravitational wave source to one without, complementing traditional detection methods. It is shown via astrophysically representative software injections in Gaussian measurement noise that the procedure distinguishes a gravitational wave from pure noise down to a characteristic wave strain of $h_0 \approx 2 \times 10^{-15}$. Full posterior distributions of model parameters are recovered and tested for accuracy. There is a bias of $\approx 0.3$ rad in the marginalised one-dimensional posterior for the orbital inclination $\iota$, introduced by dropping the so-called `pulsar terms'. Smaller biases $\lesssim 10 \%$ are also observed in other static parameters., Comment: 26 pages, 11 figures. Accepted for publication in MNRAS

We present results of a follow-up search for continuous gravitational waves (CWs) associated with sub-threshold candidates from the LIGO-Virgo-KAGRA (LVK) All-Sky All-Frequency (ASAF) directed radiometer analysis, using Advanced LIGO data from the third observing run (O3). Each ASAF candidate corresponds to a $1/32\,\,\mathrm{Hz}$ frequency band and ${\sim}13\,\,\mathrm{deg}^2$ sky pixel. Assuming they represent possible CW sources, we analyze all $515$ ASAF candidates using a semi-coherent, $\mathcal{F}$-statistic-based matched filter search. The search algorithm incorporates a hidden Markov model (HMM), expanding the signal model to allow frequency spin-wandering, as well as unmodeled frequency evolution of less than $10^{-5}$ Hz per day that is not captured by the searched range of $\pm10^{-9}\,\,\mathrm{Hz/s}$ in frequency derivative. Significance thresholds with a $5\%$ probability of false alarm per ASAF candidate are determined empirically by searching detector noise at various off-target sky positions. We obtain $14$ outliers surviving a set of vetoes designed to eliminate instrumental artifacts. Upon further investigation, these outliers are deemed unlikely to represent astrophysical signals. We estimate the sensitivity of our search to both isolated and binary sources with orbital period greater than one year by recovering simulated signals added to detector data. The minimum detectable strain amplitude at $95\%$ confidence for isolated (long-period binary) sources is $h_0^{95\%} = 8.8\times 10^{-26}$ ($9.4\times 10^{-26}$) at a frequency of $222.6\,\,\mathrm{Hz}$. While this study focuses on ASAF sub-threshold candidates, the method presented could be applied to follow up candidates from future all-sky CW searches, complementing currently existing methods., Comment: 20 pages, 12 figures. Published in PRD

We report the radio and high-energy properties of a new outburst from the radio-loud magnetar 1E 1547.0$-$5408. Following the detection of a short burst from the source with Swift-BAT on 2022 April 7, observations by NICER detected an increased flux peaking at $(6.0 \pm 0.4) \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ in the soft X-ray band, falling to the baseline level of $1.7\times10^{-11}$ erg s$^{-1}$ cm$^{-2}$ over a 17-day period. Joint spectroscopic measurements by NICER and NuSTAR indicated no change in the hard non-thermal tail despite the prominent increase in soft X-rays. Observations at radio wavelengths with Murriyang, the 64-m Parkes radio telescope, revealed that the persistent radio emission from the magnetar disappeared at least 22 days prior to the initial Swift-BAT detection and was re-detected two weeks later. Such behavior is unprecedented in a radio-loud magnetar, and may point to an unnoticed slow rise in the high-energy activity prior to the detected short-bursts. Finally, our combined radio and X-ray timing revealed the outburst coincided with a spin-up glitch, where the spin-frequency and spin-down rate increased by $0.2 \pm 0.1$ $\mu$Hz and $(-2.4 \pm 0.1) \times 10^{-12}$ s$^{-2}$ respectively. A linear increase in spin-down rate of $(-2.0 \pm 0.1) \times 10^{-19}$ s$^{-3}$ was also observed over 147 d of post-outburst timing. Our results suggest that the outburst may have been associated with a reconfiguration of the quasi-polar field lines, likely signalling a changing twist, accompanied by spatially broader heating of the surface and a brief quenching of the radio signal, yet without any measurable impact on the hard X-ray properties., Comment: Accepted to ApJ, 18 pages, 7 figures. Fixed typo in abstract and discussion

We present a search for continuous gravitational wave signals from an unidentified pulsar potentially powering HESS~J1427-608, a spatially unresolved TeV point source detected by the High Energy Stereoscopic System (HESS). The search uses a semi-coherent algorithm, which combines the maximum likelihood $\mathcal{F}$~statistic with a hidden Markov model to efficiently detect and track quasi-monochromatic signals that wander randomly in frequency. It uses data from the second observing run of the Advanced Laser Interferometer Gravitational-Wave Observatory. Multi-wavelength observations of the HESS source are combined with the proprieties of the population of TeV-bright pulsar wind nebulae to constrain the search parameters. We find no evidence of gravitational wave emission from this target. We set upper limits on the characteristic wave strain $h_0^{95\%}$ (for circularly polarized signals) at $95\%$ confidence level in sample sub-bands and interpolate it to estimate the sensitivity in the full band. We find $h_0^{95\%} = 1.3\times 10^{-25}$ near 185~Hz. The implied constraints on the ellipticity and \textit{r}-mode amplitude reach $\epsilon\leq 10^{-5}$ and $\alpha \leq 10^{-3}$ at 200~Hz, respectively., Comment: 16 pages, 7 figures

Following up large numbers of candidates in continuous gravitational wave searches presents a challenge, particularly in regard to computational power and the time required to manually scrutinize each of the candidates. It is important to design and test good follow-up procedures that are safe (i.e., minimize false dismissals) and computationally efficient across many search configurations. We investigate two follow-up procedures, or "vetoes," both of which exploit the Doppler modulation predicted in astrophysical signals. In particular, we introduce the concept of using an effective point spread function as part of our veto criteria. We take advantage of a well-established semicoherent search algorithm based on a hidden Markov model to study various search configurations and to generalize the veto criteria by considering the overall veto performance in terms of efficiency and safety. The results can serve as a guideline for follow-up studies in future continuous gravitational wave searches using a hidden Markov model algorithm. The results also apply qualitatively to other semicoherent search algorithms., Comment: 24 pages, 15 figures

Radio pulsar glitches probe far-from-equilibrium processes involving stress accumulation and relaxation in neutron star interiors. Previous studies of glitch rates have focused on individual pulsars with as many recorded glitches as possible. In this work we analyze glitch rates using all available data including objects that have glitched never or once. We assume the glitch rate follows a homogeneous Poisson process, and therefore exclude pulsars which exhibit quasiperiodic glitching behavior. Calculating relevant Bayes factors shows that a model in which the glitch rate $\lambda$ scales as a power of the characteristic age $\tau$ is preferred over models which depend arbitrarily on powers of the spin frequency $\nu$ and/or its time derivative $\dot{\nu}$. For $\lambda = A (\tau/\tau\vref)^{-\gamma}$, where $\tau_{\rm ref}=1\ {\rm yr}$ is a reference time, the posterior distributions are unimodal with $A=\ModelAAglitch\ \rm{yr}^{-1}$, and $\gamma=\ModelAgammaglitch$. Importantly, the data exclude with 99\% confidence the possibility $\gamma=1$ canvassed in the literature. When objects with zero recorded glitches are included, the age-based rate law is still preferred and the posteriors change to give $A=\ModelAAall\ \rm{yr}^{-1}$, and $\gamma=\ModelAgammaall$. The updated estimates still support increased glitch activity for younger pulsars, while demonstrating that the large number of objects with zero glitches contain important statistical information about the rate, provided that they are part of the same population as opposed to a disjoint population which never glitches for some unknown physical reason., Comment: 18 pages, 11 figures. Accepted for publication in MNRAS

The $\mathcal{F}$-statistic is a detection statistic used widely in searches for continuous gravitational waves with terrestrial, long-baseline interferometers. A new implementation of the $\mathcal{F}$-statistic is presented which accelerates the existing "resampling" algorithm using graphics processing units (GPUs). The new implementation runs between 10 and 100 times faster than the existing implementation on central processing units without sacrificing numerical accuracy. The utility of the GPU implementation is demonstrated on a pilot narrowband search for four newly discovered millisecond pulsars in the globular cluster Omega Centauri using data from the second Laser Interferometer Gravitational-Wave Observatory observing run. The computational cost is $17.2$ GPU-hours using the new implementation, compared to 1092 core-hours with the existing implementation., Comment: Accepted for publication in Classical and Quantum Gravity. 19 pages, 4 figures

We report on a timing programme of 74 young pulsars that have been observed by the Parkes 64-m radio telescope over the past decade. Using modern Bayesian timing techniques, we have measured the properties of 124 glitches in 52 of these pulsars, of which 74 are new. We demonstrate that the glitch sample is complete to fractional increases in spin-frequency greater than $\Delta\nu^{90\%}_{g}/\nu \approx 8.1 \times 10^{-9}$. We measure values of the braking index, $n$, in 33 pulsars. In most of these pulsars, their rotational evolution is dominated by episodes of spin-down with $n > 10$, punctuated by step changes in the spin-down rate at the time of a large glitch. The step changes are such that, averaged over the glitches, the long-term $n$ is small. We find a near one-to-one relationship between the inter-glitch value of $n$ and the change in spin-down of the previous glitch divided by the inter-glitch time interval. We discuss the results in the context of a range of physical models., Comment: 24 pages, 9 figures and 9 tables. Published in MNRAS. (Fixed typo in the abstract and conclusion)

In certain pulsar timing experiments, where observations are scheduled approximately periodically (e.g. daily), timing models with significantly different frequencies (including but not limited to glitch models with different frequency increments) return near-equivalent timing residuals. The average scheduling aperiodicity divided by the phase error due to time-of-arrival uncertainties is a useful indicator of when the degeneracy is important. Synthetic data are used to explore the effect of this degeneracy systematically. It is found that phase-coherent tempo2 or temponest-based approaches are biased sometimes toward reporting small glitch sizes regardless of the true glitch size. Local estimates of the spin frequency alleviate this bias. A hidden Markov model is free from bias towards small glitches and announces explicitly the existence of multiple glitch solutions but sometimes fails to recover the correct glitch size. Two glitches in the UTMOST public data release are re-assessed, one in PSR J1709$-$4429 at MJD 58178 and the other in PSR J1452$-$6036 at MJD 58600. The estimated fractional frequency jump in PSR J1709$-$4429 is revised upward from $\Delta f/f = (54.6\pm 1.0) \times 10^{-9}$ to $\Delta f/f = (2432.2 \pm 0.1) \times 10^{-9}$ with the aid of additional data from the Parkes radio telescope. We find that the available UTMOST data for PSR J1452$-$6036 are consistent with $\Delta f/f = 270 \times 10^{-9} + N/(fT)$ with $N = 0,1,2$, where $T \approx 1\,\text{sidereal day}$ is the observation scheduling period. Data from the Parkes radio telescope can be included, and the $N = 0$ case is selected unambiguously with a combined dataset., Comment: Accepted for publication in MNRAS. 16 pages, 12 figures

We pursue a novel strategy towards a first detection of continuous gravitational waves from rapidly-rotating deformed neutron stars. Computational power is focused on a narrow region of signal parameter space selected by a strategically-chosen benchmark. We search data from the 2nd observing run of the LIGO Observatory with an optimised analysis run on graphics processing units. While no continuous waves are detected, the search achieves a sensitivity to gravitational wave strain of $h_0 = 1.01{\times}10^{-25}$ at 90% confidence, 24% to 69% better than past searches of the same parameter space. Constraints on neutron star deformity are within theoretical maxima, thus a detection by this search was not inconceivable., Comment: 8 pages, 4 figures; accepted for publication in Physical Review D; v2 corrects number of segments stated in Section IV

Isolated neutron stars are prime targets for continuous-wave (CW) searches by ground-based gravitational$-$wave interferometers. Results are presented from a CW search targeting ten pulsars. The search uses a semicoherent algorithm, which combines the maximum-likelihood $\mathcal{F}$-statistic with a hidden Markov model (HMM) to efficiently detect and track quasi$-$monochromatic signals which wander randomly in frequency. The targets, which are associated with TeV sources detected by the High Energy Stereoscopic System (H.E.S.S.), are chosen to test for gravitational radiation from young, energetic pulsars with strong $\mathrm{\gamma}$-ray emission, and take maximum advantage of the frequency tracking capabilities of HMM compared to other CW search algorithms. The search uses data from the second observing run of the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). It scans 1$-$Hz sub-bands around $f_*$, 4$f_*$/3, and 2$f_*$, where $f_*$ denotes the star's rotation frequency, in order to accommodate a physically plausible frequency mismatch between the electromagnetic and gravitational-wave emission. The 24 sub-bands searched in this study return 5,256 candidates above the Gaussian threshold with a false alarm probability of 1$\%$ per sub-band per target. Only 12 candidates survive the three data quality vetoes which are applied to separate non$-$Gaussian artifacts from true astrophysical signals. CW searches using the data from subsequent observing runs will clarify the status of the remaining candidates., Comment: 28 pages, 22 figures; accepted for publication in Physical Review D

Low mass X-ray binaries are prime targets for continuous gravitational wave searches by ground-based interferometers. Results are presented from a search for five low-mass X-ray binaries whose spin frequencies and orbital elements are measured accurately from X-ray pulsations: HETE J1900.1-2455, IGR J00291+5934, SAX J1808.4-3658, XTE J0929-314, and XTE J1814-338. Data are analysed from Observing Run 2 of the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). The search algorithm uses a hidden Markov model to track spin wandering, the $\mathcal{J}$-statistic maximum likelihood matched filter to track orbital phase, and a suite of five vetoes to reject artefacts from non-Gaussian noise. The search yields a number of low-significance, above threshold candidates consistent with the selected false-alarm probability. The candidates will be followed up in subsequent observing runs., Comment: 24 pages, 11 figures; accepted for publication in Physical Review D

We calculate the baryon asymmetry of the Universe in the Z3-invariant Next-to-Minimal Supersymmetric Standard Model where the interactions of the singlino provide the necessary source of charge and parity violation. Using the closed time path formalism, we derive and solve transport equations for the cases where the singlet acquires a vacuum expectation value (VEV) before and during the electroweak phase transition. We perform a detailed scan to show how the baryon asymmetry varies throughout the relevant parameter space. Our results show that the case where the singlet acquires a VEV during the electroweak phase transition typically generates a larger baryon asymmetry, although we expect that the case where the singlet acquires a VEV first is far more common for any model in which parameters unify at a high scale. Finally, we examine the dependence of the baryon asymmetry on the three-body interactions involving gauge singlets., Comment: 24 pages, version submitted to the journal