BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew Melatos

Andrew Melatos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.SRphysics.soc-phastro-ph.IMmath-phmath.MPastro-phastro-ph.COastro-ph.GAq-fin.GNcond-mat.quant-gascond-mat.stat-mechhep-phphysics.ed-phphysics.flu-dynq-fin.TR

Catalog footprint

What is connected

42works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

42 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
preprint2016arXiv

Ernst formulation of axisymmetric fields in $f(R)$ gravity: applications to neutron stars and gravitational waves

The Ernst formulation of the Einstein equations is generalised to accommodate $f(R)$ theories of gravity. It is shown that, as in general relativity, the axisymmetric $f(R)$ field equations for a vacuum spacetime that is either stationary or cylindrically symmetric reduce to a single, non-linear differential equation for a complex-valued scalar function. As a worked example, we apply the generalised Ernst equations to derive a $f(R)$ generalisation of the Zipoy-Voorhees metric, which may be used to describe the gravitational field outside of an ellipsoidal neutron star. We also apply the theory to investigate the phase speed of large-amplitude gravitational waves in $f(R)$ gravity in the context of soliton-like solutions that display shock-wave behaviour across the causal boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detecting gravitational waves from mountains on neutron stars in the Advanced Detector Era

Rapidly rotating Neutron Stars (NSs) in Low Mass X-ray Binaries (LMXBs) are thought to be interesting sources of Gravitational Waves (GWs) for current and next generation ground based detectors, such as Advanced LIGO and the Einstein Telescope. The main reason is that many of the NS in these systems appear to be spinning well below their Keplerian breakup frequency, and it has been suggested that torques associated with GW emission may be setting the observed spin period. This assumption has been used extensively in the literature to assess the strength of the likely gravitational wave signal. There is now, however, a significant amount of theoretical and observation work that suggests that this may not be the case, and that GW emission is unlikely to be setting the spin equilibrium period in many systems. In this paper we take a different starting point and predict the GW signal strength for two physical mechanisms that are likely to be at work in LMXBs: crustal mountains due to thermal asymmetries and magnetically confined mountains. We find that thermal crustal mountains in transient LMXBs are unlikely to lead to detectable GW emission, while persistent systems are good candidates for detection by Advanced LIGO and by the Einstein Telescope. Detection prospects are pessimistic for the magnetic mountain case, unless the NS has a buried magnetic field of $B\approx 10^{12}$ G, well above the typically inferred exterior dipole fields of these objects. Nevertheless, if a system were to be detected by a GW observatory, cyclotron resonant scattering features in the X-ray emission could be used to distinguish between the two different scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electromagnetic spin down of a transient accreting millisecond pulsar during quiescence

The measured spin-down rates in quiescence of the transient accreting millisecond pulsars IGR J00291+5934, XTE J1751-305, SAX J1808.4-3658, and Swift J1756.9-2508 have been used to estimate the magnetic moments of these objects assuming standard magnetic dipole braking. It is shown that this approach leads to an overestimate, if the amount of residual accretion is enough to distort the magnetosphere away from a force-free configuration, through magnetospheric mass loading or crushing, so that the lever arm of the braking torque migrates inside the light cylinder. We derive an alternative spin-down formula and calculate the residual accretion rates where the formula is applicable. As a demonstration, we apply the alternative spin-down formula to produce updated magnetic moment estimates for the four objects above. We note that, based on current uncertain observations of quiescent accretion rates, magnetospheric mass loading and crushing are neither firmly indicated nor ruled out in these four objects. Because quiescent accretion rates are not measured directly (only upper limits are placed), it is impossible to be confident without more data whether the thresholds for magnetospheric mass loading or crushing are reached or not.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interpreting the AXP 1E 2259+586 antiglitch as a change in internal magnetization

The sudden spin-down event ('anti-glitch') observed in AXP 1E 2259+586 on 2012 April 21 was arguably caused by a decay of its internal toroidal magnetic field component, which turns a stable prolate configuration into an unstable one. We refine previous models of this process by modelling the star's magnetic field self-consistently as a 'twisted torus' configuration in non-barotropic equilibrium (which allows us to explore a greater range of equilibrium configurations). It is shown that, if the star's magnetic field is purely dipolar, the change in the toroidal field strength required to produce an anti-glitch of the observed size can be ~ 10 times larger than previously calculated. If the star has a quadrupolar magnetic field component, then an anti-glitch of similar magnitude can be produced via a decay of the quadrupole component, in addition to a decay of the toroidal component. We show that, if the quadrupole component decays, the minimum initial toroidal field strength and the change in toroidal field strength needed to produce the observed anti-glitch are lower than in the pure dipole twisted torus. In addition, we predict the maximum anti-glitch sizes, assuming that they are caused by a change in ellipticity, in four glitching magnetars and discuss the implications for energetics of accompanying X-ray bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetars: Properties, Origin and Evolution

Magnetars are neutron stars in which a strong magnetic field is the main energy source. About two dozens of magnetars, plus several candidates, are currently known in our Galaxy and in the Magellanic Clouds. They appear as highly variable X-ray sources and, in some cases, also as radio and/or optical pulsars. Their spin periods (2-12 s) and spin-down rates (~10^{-13}-10^{-10} s/s) indicate external dipole fields of ~10^{13-15} G, and there is evidence that even stronger magnetic fields are present inside the star and in non-dipolar magnetospheric components. Here we review the observed properties of the persistent emission from magnetars, discuss the main models proposed to explain the origin of their magnetic field and present recent developments in the study of their evolution and connection with other classes of neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetic Field Generation in Stars

Enormous progress has been made on observing stellar magnetism in stars from the main sequence through to compact objects. Recent data have thrown into sharper relief the vexed question of the origin of stellar magnetic fields, which remains one of the main unanswered questions in astrophysics. In this chapter we review recent work in this area of research. In particular, we look at the fossil field hypothesis which links magnetism in compact stars to magnetism in main sequence and pre-main sequence stars and we consider why its feasibility has now been questioned particularly in the context of highly magnetic white dwarfs. We also review the fossil versus dynamo debate in the context of neutron stars and the roles played by key physical processes such as buoyancy, helicity, and superfluid turbulence,in the generation and stability of neutron star fields. Independent information on the internal magnetic field of neutron stars will come from future gravitational wave detections. Thus we maybe at the dawn of a new era of exciting discoveries in compact star magnetism driven by the opening of a new, non-electromagnetic observational window. We also review recent advances in the theory and computation of magnetohydrodynamic turbulence as it applies to stellar magnetism and dynamo theory. These advances offer insight into the action of stellar dynamos as well as processes whichcontrol the diffusive magnetic flux transport in stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Models of Pulsar Glitches

Radio pulsars provide us with some of the most stable clocks in the universe. Nevertheless several pulsars exhibit sudden spin-up events, known as glitches. More than forty years after their first discovery, the exact origin of these phenomena is still open to debate. It is generally thought that they an observational manifestation of a superfluid component in the stellar interior and provide an insight into the dynamics of matter at extreme densities. In recent years there have been several advances on both the theoretical and observational side, that have provided significant steps forward in our understanding of neutron star interior dynamics and possible glitch mechanisms. In this article we review the main glitch models that have been proposed and discuss our understanding, in the light of current observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutron star deformation due to poloidal-toroidal magnetic fields of arbitrary multipole order: a new analytic approach

A recipe is presented to construct an analytic, self-consistent model of a non-barotropic neutron star with a poloidal-toroidal field of arbitrary multipole order, whose toroidal component is confined in a torus around the neutral curve inside the star, as in numerical simulations of twisted tori. The recipe takes advantage of magnetic-field-aligned coordinates to ensure continuity of the mass density at the surface of the torus. The density perturbation and ellipticity of such a star are calculated in general and for the special case of a mixed dipole-quadrupole field as a worked example. The calculation generalises previous work restricted to dipolar, poloidal-toroidal and multipolar, poloidal-only configurations. The results are applied, as an example, to magnetars whose observations (e.g., spectral features and pulse modulation) indicate that the internal magnetic fields may be at least one order of magnitude stronger than the external fields, as inferred from their spin downs, and are not purely dipolar.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pinned vortex hopping in a neutron star crust

The motion of superfluid vortices in a neutron star crust is at the heart of most theories of pulsar glitches. Pinning of vortices to ions can decouple the superfluid from the crust and create a reservoir of angular momentum. Sudden large scale unpinning can lead to an observable glitch. In this paper we investigate the scattering of a free vortex off a pinning potential and calculate its mean free path, in order to assess whether unpinned vortices can skip multiple pinning sites and come close enough to their neighbours to trigger avalanches, or whether they simply hop from one pinning site to another giving rise to a more gradual creep. We find that there is a significant range of parameter space in which avalanches can be triggered, thus supporting the hypothesis that they may lie at the origin of pulsar glitches. For realistic values of the pinning force and superfluid drag parameters we find that avalanches are more likely in the higher density regions of the crust where pinning is stronger. Physical differences in stellar parameters, such as mass and temperature, may lead to a switch between creep-like motion and avalanches, explaining the different characteristics of glitching pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pulsar timing noise and the minimum observation time to detect gravitational waves with pulsar timing arrays

The sensitivity of pulsar timing arrays to gravitational waves is, at some level, limited by timing noise. Red timing noise - the stochastic wandering of pulse arrival times with a red spectrum - is prevalent in slow-spinning pulsars and has been identified in many millisecond pulsars. Phenomenological models of timing noise, such as from superfluid turbulence, suggest that the timing noise spectrum plateaus below some critical frequency, $f_c$, potentially aiding the hunt for gravitational waves. We examine this effect for individual pulsars by calculating minimum observation times, $T_{\rm min}(f_c)$, over which the gravitational wave signal becomes larger than the timing noise plateau. We do this in two ways: 1) in a model-independent manner, and 2) by using the superfluid turbulence model for timing noise as an example to illustrate how neutron star parameters can be constrained. We show that the superfluid turbulence model can reproduce the data qualitatively from a number of pulsars observed as part of the Parkes Pulsar Timing Array. We further show how a value of $f_c$, derived either through observations or theory, can be related to $T_{\rm min}$. This provides a diagnostic whereby the usefulness of timing array pulsars for gravitational-wave detection can be quantified.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing modified gravity and no-hair relations for the Kerr-Newman metric through quasi-periodic oscillations of galactic microquasars

We construct multipole moments for stationary, asymptotically flat, spacetime solutions to higher-order curvature theories of gravity. The moments are defined using $3+1$ techniques involving timelike Killing vector constructions as in the classic papers by Geroch and Hansen. Using the fact that the Kerr-Newman metric is a vacuum solution to a particular class of $f(R)$ theories of gravity, we compute all its moments, and find that they admit recurrence relations similar to those for the Kerr solution in general relativity. It has been proposed previously that modelling the measured frequencies of quasi-periodic oscillations from galactic microquasars enables experimental tests of the no-hair theorem. We explore the possibility that, even if the no-hair relation is found to break down in the context of general relativity, there may be an $f(R)$ counterpart that is preserved. We apply the results to the microquasars GRS $1915$+$105$ and GRO J$1655$-$40$ using the diskoseismology and kinematic resonance models, and constrain the spins and `charges' [which are not really electric charges in the $f(R)$ context] of their black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cyclotron line signatures of thermal and magnetic mountains from accreting neutron stars

Cyclotron resonance scattering features (CRSFs) in the X-ray spectrum of an accreting neutron star are modified differently by accretion mounds sustained by magnetic and thermocompositional gradients. It is shown that one can discriminate, in principle, between mounds of different physical origins by studying how the line energy, width, and depth of a CRSF depend on the orientation of the neutron star, accreted mass, surface temperature distribution, and equation of state. CRSF signatures including gravitational light bending are computed for both phase-resolved and phase-averaged spectra on the basis of self-consistent Grad-Shafranov mound equilibria satisfying a global flux-freezing constraint. The prospects of multimessenger X-ray and gravitational-wave observations with future instruments are canvassed briefly.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational waves from rapidly rotating neutron stars

Rapidly rotating neutron stars in Low Mass X-ray Binaries have been proposed as an interesting source of gravitational waves. In this chapter we present estimates of the gravitational wave emission for various scenarios, given the (electromagnetically) observed characteristics of these systems. First of all we focus on the r-mode instability and show that a 'minimal' neutron star model (which does not incorporate exotica in the core, dynamically important magnetic fields or superfluid degrees of freedom), is not consistent with observations. We then present estimates of both thermally induced and magnetically sustained mountains in the crust. In general magnetic mountains are likely to be detectable only if the buried magnetic field of the star is of the order of $B\approx 10^{12}$ G. In the thermal mountain case we find that gravitational wave emission from persistent systems may be detected by ground based interferometers. Finally we re-asses the idea that gravitational wave emission may be balancing the accretion torque in these systems, and show that in most cases the disc/magnetosphere interaction can account for the observed spin periods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Implementation of the frequency-modulated sideband search method for gravitational waves from low mass X-ray binaries

We describe the practical implementation of the sideband search, a search for periodic gravitational waves from neutron stars in binary systems. The orbital motion of the source in its binary system causes frequency-modulation in the combination of matched filters known as the $\mathcal{F}$-statistic. The sideband search is based on the incoherent summation of these frequency-modulated $\mathcal{F}$-statistic sidebands. It provides a new detection statistic for sources in binary systems, called the $\mathcal{C}$-statistic. The search is well suited to low-mass X-ray binaries, the brightest of which, called Sco X-1, is an ideal target candidate. For sources like Sco X-1, with well constrained orbital parameters, a slight variation on the search is possible. The extra orbital information can be used to approximately demodulate the data from the binary orbital motion in the coherent stage, before incoherently summing the now reduced number of sidebands. We investigate this approach and show that it improves the sensitivity of the standard Sco X-1 directed sideband search. Prior information on the neutron star inclination and gravitational wave polarization can also be used to improve upper limit sensitivity. We estimate the sensitivity of a Sco X-1 directed sideband search on 10 days of LIGO data and show that it can beat previous upper limits in current LIGO data, with a possibility of constraining theoretical upper limits using future advanced instruments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutron star deformation due to multipolar magnetic fields

Certain multi-wavelength observations of neutron stars, such as intermittent radio emissions from rotation-powered pulsars beyond the pair-cascade death line, the pulse profile of the magnetar SGR 1900+14 after its 1998 August 27 giant flare, and X-ray spectral features of PSR J0821-4300 and SGR 0418+5729, suggest that the magnetic fields of non-accreting neutron stars are not purely dipolar and may contain higher-order multipoles. Here, we calculate the ellipticity of a non-barotropic neutron star with (i) a quadrupole poloidal-toroidal field, and (ii) a purely poloidal field containing arbitrary multipoles, deriving the relation between the ellipticity and the multipole amplitudes. We present, as a worked example, a purely poloidal field comprising dipole, quadrupole, and octupole components. We show the correlation between field energy and ellipticity for each multipole, that the l=4 multipole has the lowest energy, and that l=5 has the lowest ellipticity. We show how a mixed multipolar field creates an observationally testable mismatch between the principal axes of inertia (to be inferred from gravitational wave data) and the magnetic inclination angle. Strong quadrupole and octupole components (with amplitudes ~ 10^2 times higher than the dipole) in SGR 0418+5729 still yield ellipticity ~ 10^{-8}, consistent with current gravitational wave upper limits. The existence of higher multipoles in fast-rotating objects (e.g., newborn magnetars) has interesting implications for the braking law and hence phase tracking during coherent gravitational wave searches.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pulsar timing noise from superfluid turbulence

Shear-driven turbulence in the superfluid interior of a neutron star exerts a fluctuating torque on the rigid crust, causing the rotational phase to walk randomly. The phase fluctuation spectrum is calculated analytically for incompressible Kolmogorov turbulence and is found to be red; the half-power point is set by the observed spin-down rate, the crust-superfluid lag, and the dynamical response time of the superfluid. Preliminary limits are placed on the latter quantities using selected time- and frequency-domain data. It is found that measurements of the normalization and slope of the power spectrum are reproduced for reasonable choices of the turbulence parameters. The results point preferentially to the neutron star interior containing a turbulent superfluid rather than a turbulent Navier-Stokes fluid. The implications for gravitational wave detection by pulsar timing arrays are discussed briefly.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Realistic neutron star constraints on bosonic asymmetric dark matter

It has been argued that the existence of old neutron stars excludes the possibility of non-annihilating light bosonic dark matter, such as that arising in asymmetric dark matter scenarios. If non-annihilating dark matter is captured by neutron stars, the density will eventually become sufficient for black hole formation. However, the dynamics of collapse is highly sensitive to dark-matter self-interactions. Repulsive self-interactions, even if extremely weak, can prevent black hole formation. We argue that self-interactions will necessarily be present, and estimate their strength in representative models. We also consider co-annihilation of dark matter with nucleons, which arises naturally in many asymmetric dark matter models, and which again acts to prevent black hole formation. We demonstrate how the excluded region of the dark-matter parameter space shrinks as the strength of such interactions is increased, and conclude that neutron star observations do not exclude most realistic bosonic asymmetric dark matter models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stochastic gravitational wave background from hydrodynamic turbulence in differentially rotating neutron stars

Hydrodynamic turbulence driven by crust-core differential rotation imposes a fundamental noise floor on gravitational wave observations of neutron stars. The gravitational wave emission peaks at the Kolmogorov decoherence frequency which, for reasonable values of the crust-core shear, ΔΩ, occurs near the most sensitive part of the frequency band for ground-based, long-baseline interferometers. We calculate the energy density spectrum of the stochastic gravitational wave background from a cosmological population of turbulent neutron stars generalising previous calculations for individual sources. The spectrum resembles a piecewise power law, Ω_{gw}(ν)=Ω_αν^α, with α=-1 and 7 above and below the decoherence frequency respectively, and its normalisation scales as Ω_α\propto(ΔΩ)^{7}. Non-detection of a stochastic signal by Initial LIGO implies an upper limit on ΔΩand hence by implication on the internal relaxation time-scale for the crust and core to come into co-rotation, τ_{d}=ΔΩ/\dotΩ, where \dotΩ is the observed electromagnetic spin-down rate, with τ_{d}\lesssim 10^{7} yr for accreting millisecond pulsars and τ_{d}\lesssim 10^{5} yr for radio-loud pulsars. Target limits on τ_{d} are also estimated for future detectors, namely Advanced LIGO and the Einstein Telescope, and are found to be astrophysically interesting.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Tilted torus magnetic fields in neutron stars and their gravitational wave signatures

Gravitational-wave diagnostics are developed for discriminating between varieties of mixed poloidal-toroidal magnetic fields in neutron stars, with particular emphasis on differentially rotating protoneutron stars. It is shown that tilted torus magnetic fields, defined as the sum of an internal/external poloidal component, whose axis of symmetry is tilted with respect to the rotation axis, and an internal toroidal component, whose axis of symmetry is aligned with the rotation axis, deform the star triaxially, unlike twisted torus fields, which deform the star biaxially. Utilizing an analytic tilted torus example, we show that these two topologies can be distinguished by their gravitational wave spectrum and polarization phase portraits. For example, the relative amplitudes and frequencies of the spectral peaks allows one to infer the relative strengths of the toroidal and poloidal components of the field, and the magnetic inclination angle. Finally, we show how a tilted torus field arises naturally from magnetohydrodynamic simulations of differentially rotating neutron stars, and how the gravitational wave spectrum evolves as the internal toroidal field winds up. These results point to the sorts of experiments that may become possible once gravitational wave interferometers detect core-collapse supernovae routinely.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Knock-on processes in superfluid vortex avalanches and pulsar glitch statistics

A framework is presented for a statistical theory of neutron star glitches, motivated by the results emerging from recent Gross-Pitaevskii simulations of pinned, decelerating quantum condensates. It is shown that the observed glitch size distributions cannot be reproduced if superfluid vortices unpin independently via a Poisson process; the central limit theorem yields a narrow Gaussian for the size distribution, instead of the broad, power-law tail observed. This conclusion is not altered fundamentally when a range of pinning potentials is included, which leads to excavation of the potential distribution of occupied sites, vortex accumulation at strong pinning sites, and hence the occasional, abnormally large glitch. Knock-on processes are therefore needed to make the unpinning rate of a vortex conditional on the pinning state of its near and/or remote neighbours, so that the Gaussian size distributions resulting generically from the central limit theorem are avoided. At least two knock-on processes, nearest- neighbour proximity knock-on and remote acoustic knock-on, are clearly evident in the Gross-Pitaevskii simulation output. It is shown that scale-invariant (i.e. power-law) vortex avalanches occur when knock-on is included, provided that two specific relations hold between the temperature and spin-down torque. This fine tuning is unlikely in an astronomical setting, leaving the overall problem partly unsolved. A state-dependent Poisson formalism is presented which will form the basis of future studies in this area.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Loan and nonloan flows in the Australian interbank network

High-value transactions between Australian banks are settled in the Reserve Bank Information and Transfer System (RITS) administered by the Reserve Bank of Australia. RITS operates on a real-time gross settlement (RTGS) basis and settles payments sourced from the SWIFT, the Austraclear, and the interbank transactions entered directly into RITS. In this paper, we analyse a dataset received from the Reserve Bank of Australia that includes all interbank transactions settled in RITS on an RTGS basis during five consecutive weekdays from 19 February 2007 inclusive, a week of relatively quiescent market conditions. The source, destination, and value of each transaction are known, which allows us to separate overnight loans from other transactions (nonloans) and reconstruct monetary flows between banks for every day in our sample. We conduct a novel analysis of the flow stability and examine the connection between loan and nonloan flows. Our aim is to understand the underlying causal mechanism connecting loan and nonloan flows. We find that the imbalances in the banks' exchange settlement funds resulting from the daily flows of nonloan transactions are almost exactly counterbalanced by the flows of overnight loans. The correlation coefficient between loan and nonloan imbalances is about -0.9 on most days. Some flows that persist over two consecutive days can be highly variable, but overall the flows are moderately stable in value. The nonloan network is characterised by a large fraction of persistent flows, whereas only half of the flows persist over any two consecutive days in the loan network. Moreover, we observe an unusual degree of coherence between persistent loan flow values on Tuesday and Wednesday. We probe static topological properties of the Australian interbank network and find them consistent with those observed in other countries.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Designing a cross-correlation search for continuous-wave gravitational radiation from a neutron star in the supernova remnant SNR 1987A

A strategy is devised for a semi-coherent cross-correlation search for a young neutron star in the supernova remnant SNR 1987A, using science data from the Initial LIGO and/or Virgo detectors. An astrophysical model for the gravitational wave phase is introduced which describes the star's spin down in terms of its magnetic field strength $B$ and ellipticity $ε$, instead of its frequency derivatives. The model accurately tracks the gravitational wave phase from a rapidly decelerating neutron star under the restrictive but computationally unavoidable assumption of constant braking index, an issue which has hindered previous searches for such young objects. The theoretical sensitivity is calculated and compared to the indirect, age-based wave strain upper limit. The age-based limit lies above the detection threshold in the frequency band 75\,Hz $\lesssim ν\lesssim 450$\,Hz. The semi-coherent phase metric is also calculated and used to estimate the optimal search template spacing for the search. The range of search parameters that can be covered given our computational resources ($\sim 10^9$ templates) is also estimated. For Initial LIGO sensitivity, in the frequency band between 50\,Hz and 500\,Hz, in the absence of a detected signal, we should be able to set limits of $B \gtrsim 10^{11}$\,G and $ε\lesssim 10^{-4}$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational wave emission from a magnetically deformed non-barotropic neutron star

A strong candidate for a source of gravitational waves is a highly magnetised, rapidly rotating neutron star (magnetar) deformed by internal magnetic stresses. We calculate the mass quadrupole moment by perturbing a zeroth-order hydrostatic equilibrium by an axisymmetric magnetic field with a \emph{linked poloidal-toroidal structure}. In this work, we do \emph{not} require the model star to obey a barotropic equation of state (as a realistic neutron star is not barotropic), allowing us to explore the hydromagnetic equilibria with fewer constraints. We derive the relation between the ratio of poloidal-to-total field energy $Λ$ and ellipticity $ε$ and briefly compare our results to those obtained using the barotropic assumption. Then, we present some examples of how our results can be applied to astrophysical contexts. First, we show how our formulae, in conjunction with current gravitational wave (non-)detections of the Crab pulsar and the Cassiopeia A central compact object (Cas A CCO), can be used to constrain the strength of the internal toroidal fields of those objects. We find that, for the Crab pulsar (whose canonical equatorial dipole field strength, inferred from spin down, is $4\times 10^8$ T) to emit detectable gravitational radiation, the neutron star must have a strong toroidal field component, with maximum internal toroidal field strength $B_{\mathrm{tm}}=7\times 10^{12}$ T; for gravitational waves to be detected from the Cas A CCO at 300 Hz, $B_{\mathrm{tm}}\sim 10^{13}$ T, whereas detection at 100 Hz would require $B_{\mathrm{tm}}\sim 10^{14}$ T. Using our results, we also show how the gravitational wave signal emitted by a magnetar immediately after its birth (assuming it is born rapidly rotating, with $Λ\lesssim 0.2$) makes such a newborn magnetar a stronger candidate for gravitational wave detection than, for example, an SGR giant flare.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Instability-driven interfacial dynamo in protoneutron stars

The existence of a tachocline in the Sun has been proven by helioseismology. It is unknown whether a similar shear layer, widely regarded as the seat of magnetic dynamo action, also exists in a protoneutron star. Sudden jumps in magnetic diffusivity $η$ and turbulent vorticity $α$, for example at the interface between the neutron-finger and convective zones, are known to be capable of enhancing mean-field dynamo effects in a protoneutron star. Here we apply the well-known, plane-parallel, MacGregor-Charbonneau analysis of the Solar interfacial dynamo to the protoneutron star problem and calculate the growth rate analytically under a range of conditions. It is shown that, like the Solar dynamo, it is impossible to achieve self-sustained growth if the discontinuities in $α$, $η$, and shear are coincident and the magnetic diffusivity is isotropic. In contrast, when the jumps in $η$ and $α$ are situated away from the shear layer, self-sustained growth is possible for $P\lesssim 49.8$ ms (if the velocity shear is located at $0.3R$) or $P\lesssim 83.6$ ms (if the velocity shear is located at $0.6R$). This translates into stronger shear and/or $α$-effect than in the Sun. Self-sustained growth is also possible if the magnetic diffusivity if anisotropic, through the ${\bfΩ}\times{\bf{J}}$ effect, even when the $α$, $η$, and shear discontinuities are coincident.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Interpreting superfluid spin up through the response of the container

A recipe is presented for interpreting non-invasively the transport processes at work during relaxation of a cylindrical, superfluid-filled vessel, after it is accelerated impulsively and then allowed to respond to the viscous torque exerted by the contained fluid. The recipe exploits a recently published analytic solution for Ekman pumping in a two-component superfluid, which treats the back-reaction self-consistently in arbitrary geometry for the first time. The applicability of the recipe to He II, 3He, 3He-4He mixtures and Bose-Einstein condensates is assessed, and the effects of turbulence discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Updated gravitational-wave upper limits on the internal magnetic field strength of recycled pulsars

Recent calculations of the hydromagnetic deformation of a stratified, non-barotropic neutron star are generalized to describe objects with superconducting interiors, whose magnetic permeability μis much smaller than the vacuum value μ_0. It is found that the star remains oblate if the poloidal magnetic field energy is \gtrsim 40% of total magnetic field energy, that the toroidal field is confined to a torus which shrinks as μdecreases, and that the deformation is much larger (by a factor \sim μ_0/μ) than in a non-superconducting object. The results are applied to the latest direct and indirect upper limits on gravitational-wave emission from Laser Interferometer Gravitational Wave Observatory (LIGO) and radio pulse timing (spin-down) observations of 81 millisecond pulsars, to show how one can use these observations to infer the internal field strength. It is found that the indirect spin-down limits already imply astrophysically interesting constraints on the poloidal-toroidal field ratio and diamagnetic shielding factor (by which accretion reduces the observable external magnetic field, e.g. by burial). These constraints will improve following gravitational-wave detections, with implications for accretion-driven magnetic field evolution in recycled pulsars and the hydromagnetic stability of these objects' interiors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A note on the theory of fast money flow dynamics

The gauge theory of arbitrage was introduced by Ilinski in [arXiv:hep-th/9710148] and applied to fast money flows in [arXiv:cond-mat/9902044]. The theory of fast money flow dynamics attempts to model the evolution of currency exchange rates and stock prices on short, e.g.\ intra-day, time scales. It has been used to explain some of the heuristic trading rules, known as technical analysis, that are used by professional traders in the equity and foreign exchange markets. A critique of some of the underlying assumptions of the gauge theory of arbitrage was presented by Sornette in [arXiv:cond-mat/9804045]. In this paper, we present a critique of the theory of fast money flow dynamics, which was not examined by Sornette. We demonstrate that the choice of the input parameters used in [arXiv:cond-mat/9902044] results in sinusoidal oscillations of the exchange rate, in conflict with the results presented in [arXiv:cond-mat/9902044]. We also find that the dynamics predicted by the theory are generally unstable in most realistic situations, with the exchange rate tending to zero or infinity exponentially.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational-wave spin-down and stalling lower limits on the electrical resistivity of the accreted mountain in a millisecond pulsar

The electrical resistivity of the accreted mountain in a millisecond pulsar is limited by the observed spin-down rate of binary radio millisecond pulsars (BRMSPs) and the spins and X-ray fluxes of accreting millisecond pulsars (AMSPs). We find $η\ge 10^{-28}\,\mathrm{s}\, (τ_\mathrm{SD}/1\,\mathrm{Gyr})^{-0.8}$ (where $τ_\mathrm{SD}$ is the spin-down age) for BRMSPs and $η\ge 10^{-25}\,\mathrm{s}\,(\dot{M}_\mathrm{a}/\dot{M}_\mathrm{E})^{0.6}$ (where $\dot{M}_\mathrm{a}$ and $\dot{M}_\mathrm{E}$ are the actual and Eddington accretion rates) for AMSPs. These limits are inferred assuming that the mountain attains a steady state, where matter diffuses resistively across magnetic flux surfaces but is replenished at an equal rate by infalling material. The mountain then relaxes further resistively after accretion ceases. The BRMSP spin-down limit approaches the theoretical electron-impurity resistivity at temperatures $\ga 10^5$ K for an impurity concentration of $\sim 0.1$, while the AMSP stalling limit falls two orders of magnitude below the theoretical electron-phonon resistivity for temperatures above $10^8$ K. Hence BRMSP observations are already challenging theoretical resistivity calculations in a useful way. Next-generation gravitational-wave interferometers will constrain $η$ at a level that will be competitive with electromagnetic observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Laplace transform analysis of a multiplicative asset transfer model

We analyze a simple asset transfer model in which the transfer amount is a fixed fraction $f$ of the giver's wealth. The model is analyzed in a new way by Laplace transforming the master equation, solving it analytically and numerically for the steady-state distribution, and exploring the solutions for various values of $f\in(0,1)$. The Laplace transform analysis is superior to agent-based simulations as it does not depend on the number of agents, enabling us to study entropy and inequality in regimes that are costly to address with simulations. We demonstrate that Boltzmann entropy is not a suitable (e.g. non-monotonic) measure of disorder in a multiplicative asset transfer system and suggest an asymmetric stochastic process that is equivalent to the asset transfer model.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Gravitational Radiation from Accreting Millisecond Pulsars

It is widely assumed that the observed reduction of the magnetic field of millisecond pulsars can be connected to the accretion phase during which the pulsar is spun up by mass accretion from a companion. A wide variety of reduction mechanisms have been proposed, including the burial of the field by a magnetic mountain, formed when the accreted matter is confined to the poles by the tension of the stellar magnetic field. A magnetic mountain effectively screens the magnetic dipole moment. On the other hand, observational data suggests that accreting neutron stars are sources of gravitational waves, and magnetic mountains are a natural source of a time-dependent quadrupole moment. We show that the emission is sufficiently strong to be detectable by current and next generation long-baseline interferometers. Preliminary results from fully three-dimensional magnetohydrodynamic (MHD) simulations are presented. We find that the initial axisymmetric state relaxes into a nearly axisymmetric configuration via toroidal magnetic modes. A substantial quadrupole moment is still present in the final state, which is stable (in ideal MHD) yet highly distorted.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Collapsed and Extended Cold Dark Matter Haloes in Softened N-Body Gravity

The statistical mechanics of N cold dark matter (CDM) particles interacting via a softened gravitational potential is reviewed in the microcanonical ensemble and mean-field limit. A phase diagram for the system is computed as a function of the total energy E and gravitational softening length epsilon. For softened systems, two stable phases exist: a collapsed phase, whose radial density profile rho(r) is a central Dirac cusp, and an extended phase, for which rho(r) has a central core and rho(r) ~ r^{-2.2} at large r. It is shown that many N-body simulations of CDM haloes in the literature inadvertently sample the collapsed phase only, even though this phase is unstable when there is zero softening. Consequently, there is no immediate reason to expect agreement between simulated and observed profiles unless the gravitational potential is appreciably softened in nature.

0 citations0 reviewsNo code linkedOpen access