BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Andrew Melatos

Andrew Melatos contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcastro-ph.IMastro-ph.GAphysics.ed-phphysics.soc-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2026arXiv

Discovering pulsars in compact binaries with a hidden Markov model

Discovering radio pulsars in compact binaries, whose orbital periods $P_{\rm b}$ satisfy $P_{\rm b} \lesssim 1 \, \rm{day}$, is computationally challenging, because the time-dependent pulse frequency $f_{\rm p}(t)$ is strongly Doppler modulated by the binary motion. Here we present a new, fast, semi-coherent detection scheme based on a hidden Markov model (HMM) combined with a maximum likelihood matched filter, the Schuster periodogram. The HMM scheme complements traditional acceleration searches by dividing $f_{\rm p}(t)$ into piecewise-constant blocks and tracking the block-to-block evolution efficiently using dynamic programming. Monte Carlo simulations show that the new method can detect compact binaries with flux densities $S \geq 0.50 \, \rm{mJy}$ and orbital periods $P_{\rm b} \geq 0.012 \, \rm{day}$ under observing conditions (e.g.\ cadence) typical of radio pulsar surveys, with and without impulsive, narrowband radio frequency interference. The new method is fast; it employs the classic Viterbi algorithm to solve the HMM recursively. The central processing unit run time scales nominally as $T_{\rm run} \approx 2.8 \, N_B (N_T/10^2) (N_Q \ln N_Q/10^4 \ln 10^4) \, {\rm s}$ for $N_B$ subbands, $N_T$ coherent segments, and $N_Q$ frequency bins.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An estimate of the stochastic gravitational wave background from the MassiveBlackII simulation

A population of super-massive black hole binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of $10^{-9}$--$10^{-7}$ Hz. Detection of this signal is a current observational goal and so predictions of its characteristics are of significant interest. In this work we use super-massive black hole binary mergers from the MassiveBlackII simulation to estimate the characteristic strain of the stochastic background. We examine both a gravitational wave driven model of binary evolution and a model which also includes the effects of stellar scattering and a circumbinary gas disk. Results are consistent with PTA upper limits and similar to estimates in the literature. The characteristic strain at a reference frequency of $1 yr^{-1}$ is found to be $A_{yr^{-1}} = 6.9 \times 10^{-16}$ and $A_{yr^{-1}} = 6.4 \times 10^{-16}$ in the gravitational-wave driven and stellar scattering/gas disk cases, respectively. Using the latter approach, our models show that the SGWB is mildly suppressed compared to the purely gravitational wave driven model as frequency decreases inside the PTA frequency band.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An updated glitch rate law inferred from radio pulsars

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 $λ$ scales as a power of the characteristic age $τ$ is preferred over models which depend arbitrarily on powers of the spin frequency $ν$ and/or its time derivative $\dotν$. For $λ= A (τ/τ\vref)^{-γ}$, where $τ_{\rm ref}=1\ {\rm yr}$ is a reference time, the posterior distributions are unimodal with $A=\ModelAAglitch\ \rm{yr}^{-1}$, and $γ=\ModelAgammaglitch$. Importantly, the data exclude with 99\% confidence the possibility $γ=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 $γ=\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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Continuous gravitational waves in the lab: recovering audio signals with a table-top optical microphone

Gravitational-wave observatories around the world are searching for continuous waves: persistent signals from sources such as spinning neutron stars. These searches use sophisticated statistical techniques to look for weak signals in noisy data. In this paper, we demonstrate these techniques using a table-top model gravitational-wave detector: a Michelson interferometer where sound is used as an analog for gravitational waves. Using signal processing techniques from continuous-wave searches, we demonstrate the recovery of tones with constant and wandering frequencies. We also explore the use of the interferometer as a teaching tool for educators in physics and electrical engineering by using it as an "optical microphone" to capture music and speech. A range of filtering techniques used to recover signals from noisy data are detailed in the Supplementary Material. Here, we present highlights of our results using a combined notch plus Wiener filter and the statistical log minimum mean-square error (logMMSE) estimator. Using these techniques, we easily recover recordings of simple chords and drums, but complex music and speech are more challenging. This demonstration can be used by educators in undergraduate laboratories and can be adapted for communicating gravitational-wave and signal-processing topics to non-specialist audiences.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deep exploration for continuous gravitational waves at 171--172 Hz in LIGO second observing run data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Discerning media bias within a network of political allies and opponents: the idealized example of a biased coin

Perceptions of political bias in the media are formed directly, through the independent consumption of the published outputs of a media organization, and indirectly, through observing the collective responses of political allies and opponents to the same published outputs. A network of Bayesian learners is constructed to model this system, in which the bias perceived by each agent obeys a probability density function, which is updated according to Bayes's theorem given data about the published outputs and the beliefs of the agent's political allies and opponents. The Bayesian framework allows for uncertain beliefs, multimodal probability distribution functions, and antagonistic interactions with opponents, not just cooperation with allies. Numerical simulations are performed to test the idealized example of inferring the bias of a coin. It is found that some agents converge on the wrong conclusion faster than others converge on the right conclusion under a surprisingly broad range of conditions, when antagonistic interactions are present which "lock out" some agents from the truth, e.g. in Barabási-Albert networks. It is also found that structurally unbalanced networks routinely experience turbulent nonconvergence, where some agents fail to achieve a steady-state belief, e.g. when they are allies of two agents who are opponents themselves. The subtle phenomenon of long-term intermittency is also explored.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Graphics processing unit implementation of the F-statistic for continuous gravitational wave searches

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Parameter estimation of a two-component neutron star model with spin wandering

It is an open challenge to estimate systematically the physical parameters of neutron star interiors from pulsar timing data while separating spin wandering intrinsic to the pulsar (achromatic timing noise) from measurement noise and chromatic timing noise (due to propagation effects). In this paper we formulate the classic two-component, crust-superfluid model of neutron star interiors as a noise-driven, linear dynamical system and use a state-space-based expectation-maximization method to estimate the system parameters using gravitational-wave and electromagnetic timing data. Monte Carlo simulations show that we can accurately estimate all six parameters of the two-component model provided that electromagnetic measurements of the crust angular velocity, and gravitational-wave measurements of the core angular velocity, are both available. When only electromagnetic data are available we can recover the overall relaxation time-scale, the ensemble-averaged spin-down rate, and the strength of the white-noise torque on the crust. However, the estimates of the secular torques on the two components and white noise torque on the superfluid are biased significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ground Based Gravitational Wave Astronomy in the Asian Region

The current gravitational wave detectors have identified a surprising population of heavy stellar mass black holes, and an even larger population of coalescing neutron stars. The first observations have led to many dramatic discoveries and the confirmation of general relativity in very strong gravitational fields. The future of gravitational wave astronomy looks bright, especially if additional detectors with greater sensitivity, broader bandwidth, and better global coverage can be implemented. The first discoveries add impetus to gravitational wave detectors designed to detect in the nHz, mHz and kHz frequency bands. This paper reviews the century-long struggle that led to the recent discoveries, and reports on designs and possibilities for future detectors. The benefits of future detectors in the Asian region are discussed, including analysis of the benefits of a detector located in Australia.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Long-term statistics of pulsar glitches triggered by a Brownian stress accumulation process

A microphysics-agnostic meta-model of rotational glitches in rotation-powered pulsars is developed, wherein the globally averaged internal stress accumulates as a Brownian process between glitches, and a glitch is triggered once a critical threshold is surmounted. Precise, falsifiable predictions are made regarding long-term event statistics in individual pulsars. For example, the Spearman cross-correlation coefficient between the size of a glitch and the waiting time until the next glitch should exceed 0.25 in all pulsars. Among the six pulsars with the most recorded glitches, PSR J0537$-$6910 and PSR J0835$-$4510 are consistent with the predictions of the meta-model, while PSR J1740$-$3015 and PSR J0631$+$1036 are not. PSR J0534$+$2200 and PSR J1341$-$6220 are only consistent with the meta-model, if there exists an undetected population of small glitches with small waiting times, which we do not resolve. The results are compared with a state-dependent Poisson process, another microphysics-agnostic meta-model in the literature. The results are also applied briefly to recent pulse-to-pulse observations of PSRJ0835$-$4510, which appear to reveal evidence for a negative fluctuation in rotation frequency just prior to the 2016 glitch.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Search for gravitational waves from five low mass X-ray binaries in the second Advanced LIGO observing run with an improved hidden Markov model

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.

0 citations0 reviewsNo code linkedOpen access