Source author record

Alvin J. K. Chua

Alvin J. K. Chua appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

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

19 published item(s)

preprint2022arXiv

Adiabatic waveforms from extreme-mass-ratio inspirals: an analytical approach

Scientific analysis for the gravitational-wave detector LISA will require theoretical waveforms from extreme-mass-ratio inspirals (EMRIs) that extensively cover all possible orbital and spin configurations around astrophysical Kerr black holes. However, on-the-fly calculations of these waveforms have not yet overcome the high dimensionality of the parameter space. To confront this challenge, we present a user-ready EMRI waveform model for generic (eccentric and inclined) orbits in Kerr spacetime, using an analytical self-force approach. Our model accurately covers all EMRIs with arbitrary inclination and black hole spin, up to modest eccentricity ($\lesssim 0.3$) and separation ($\gtrsim2$--$10M$ from the last stable orbit). In that regime, our waveforms are accurate at the leading `adiabatic' order, and they approximately capture transient self-force resonances that significantly impact the gravitational-wave phase. The model fills an urgent need for extensive waveforms in ongoing data-analysis studies, and its individual components will continue to be useful in future science-adequate waveforms.

preprint2022arXiv

Non-local parameter degeneracy in the intrinsic space of gravitational-wave signals from extreme-mass-ratio inspirals

Extreme-mass-ratio inspirals will be prized sources for the upcoming space-based gravitational-wave observatory LISA. The hunt for these is beset by many open theoretical and computational problems in both source modeling and data analysis. We draw attention here to one of the most poorly understood: the phenomenon of non-local correlations in the space of extreme-mass-ratio-inspiral signals. Such correlations are ubiquitous in the continuum of possible signals (degeneracy), and severely hinder the search for actual signals in LISA data. However, they are unlikely to manifest in a realistic set of putative signals (confusion). We develop an inventory of new analysis tools in order to conduct an extensive qualitative study of degeneracy - its nature, causes, and implications. Previously proposed search strategies for extreme-mass-ratio inspirals are reviewed in the light of our results, and additional guidelines are suggested for the scientific analysis of such sources.

preprint2022arXiv

The Effect of Mission Duration on LISA Science Objectives

The science objectives of the LISA mission have been defined under the implicit assumption of a 4 yr continuous data stream. Based on the performance of LISA Pathfinder, it is now expected that LISA will have a duty cycle of $\approx 0.75$, which would reduce the effective span of usable data to 3 yr. This paper reports the results of a study by the LISA Science Group, which was charged with assessing the additional science return of increasing the mission lifetime. We explore various observational scenarios to assess the impact of mission duration on the main science objectives of the mission. We find that the science investigations most affected by mission duration concern the search for seed black holes at cosmic dawn, as well as the study of stellar-origin black holes and of their formation channels via multi-band and multi-messenger observations. We conclude that an extension to 6 yr of mission operations is recommended.

preprint2021arXiv

FastEMRIWaveforms: New tools for millihertz gravitational-wave data analysis

We present the FastEMRIWaveforms (FEW) package, a collection of tools to build and analyze extreme mass ratio inspiral (EMRI) waveforms. Here, we expand on the Physical Review Letter that introduced the first fast and accurate fully-relativistic EMRI waveform template model. We discuss the construction of the overall framework; constituent modules; and the general methods used to accelerate EMRI waveforms. Because the fully relativistic FEW model waveforms are for now limited to eccentric orbits in the Schwarzschild spacetime, we also introduce an improved Augmented Analytic Kludge (AAK) model that describes generic Kerr inspirals. Both waveform models can be accelerated using graphics processing unit (GPU) hardware. With the GPU-accelerated waveforms in hand, a variety of studies are performed including an analysis of EMRI mode content, template mismatch, and fully Bayesian Markov Chain Monte Carlo-based EMRI parameter estimation. We find relativistic EMRI waveform templates can be generated with fewer harmonic modes ($\sim10-100$) without biasing signal extraction. However, we show for the first time that extraction of a relativistic injection with semi-relativistic amplitudes can lead to strong bias and anomalous structure in the posterior distribution for certain regions of parameter space.

preprint2020arXiv

2 Fast 2 Fiducial: Gaussian processes for the interpolation and marginalization of waveform error in extreme-mass-ratio-inspiral parameter estimation

A number of open problems hinder our present ability to extract scientific information from data that will be gathered by the near-future gravitational-wave mission LISA. Many of these relate to the modeling, detection and characterization of signals from binary inspirals with an extreme component-mass ratio of $\lesssim10^{-4}$. In this paper, we draw attention to the issue of systematic error in parameter estimation due to the use of fast but approximate waveform models; this is found to be relevant for extreme-mass-ratio inspirals even in the case of waveforms with $\gtrsim90\%$ overlap accuracy and moderate ($\gtrsim30$) signal-to-noise ratios. A scheme that uses Gaussian processes to interpolate and marginalize over waveform error is adapted and investigated as a possible precursor solution to this problem. Several new methodological results are obtained, and the viability of the technique is successfully demonstrated on a three-parameter example in the setting of the LISA Data Challenge.

preprint2020arXiv

GPU-accelerated massive black hole binary parameter estimation with LISA

The Laser Interferometer Space Antenna (LISA) is slated for launch in the early 2030s. A main target of the mission is massive black hole binaries that have an expected detection rate of $\sim20$ yr$^{-1}$. We present a parameter estimation analysis for a variety of massive black hole binaries. This analysis is performed with a graphics processing unit (GPU) implementation comprising the phenomhm waveform with higher-order harmonic modes and aligned spins; a fast frequency-domain LISA detector response function; and a GPU-native likelihood computation. The computational performance achieved with the GPU is shown to be 500 times greater than with a similar CPU implementation, which allows us to analyze full noise-infused injections at a realistic Fourier bin width for the LISA mission in a tractable and efficient amount of time. With these fast likelihood computations, we study the effect of adding aligned spins to an analysis with higher-order modes by testing different configurations of spins in the injection, as well as the effect of varied and fixed spins during sampling. Within these tests, we examine three different binaries with varying mass ratios, redshifts, sky locations, and detector-frame total masses ranging over three orders of magnitude. We discuss varied correlations between the total masses and mass ratios; unique spin posteriors for the larger mass binaries; and the constraints on parameters when fixing spins during sampling, allowing us to compare to previous analyses that did not include aligned spins.

preprint2020arXiv

Learning Bayesian posteriors with neural networks for gravitational-wave inference

We seek to achieve the Holy Grail of Bayesian inference for gravitational-wave astronomy: using deep-learning techniques to instantly produce the posterior $p(θ|D)$ for the source parameters $θ$, given the detector data $D$. To do so, we train a deep neural network to take as input a signal + noise data set (drawn from the astrophysical source-parameter prior and the sampling distribution of detector noise), and to output a parametrized approximation of the corresponding posterior. We rely on a compact representation of the data based on reduced-order modeling, which we generate efficiently using a separate neural-network waveform interpolant [A. J. K. Chua, C. R. Galley & M. Vallisneri, Phys. Rev. Lett. 122, 211101 (2019)]. Our scheme has broad relevance to gravitational-wave applications such as low-latency parameter estimation and characterizing the science returns of future experiments. Source code and trained networks are available online at https://github.com/vallis/truebayes.

preprint2020arXiv

On parametric tests of relativity with false degrees of freedom

General relativity can be tested by comparing the binary-inspiral signals found in LIGO--Virgo data against waveform models that are augmented with artificial degrees of freedom. This approach suffers from a number of logical and practical pitfalls. 1) It is difficult to ascribe meaning to the stringency of the resultant constraints. 2) It is doubtful that the Bayesian model comparison of relativity against these artificial models can offer actual validation for the former. 3) It is unknown to what extent these tests might detect alternative theories of gravity for which there are no computed waveforms; conversely, when waveforms are available, tests that employ them will be superior.

preprint2020arXiv

Rapid generation of fully relativistic extreme-mass-ratio-inspiral waveform templates for LISA data analysis

The future space mission LISA will observe a wealth of gravitational-wave sources at millihertz frequencies. Of these, the extreme-mass-ratio inspirals of compact objects into massive black holes are the only sources that combine the challenges of strong-field complexity with that of long-lived signals. Such signals are found and characterized by comparing them against a large number of accurate waveform templates during data analysis, but the rapid generation of such templates is hindered by computing the $\sim10^3$-$10^5$ harmonic modes in a fully relativistic waveform. We use order-reduction and deep-learning techniques to derive a global fit for these modes, and implement it in a complete waveform framework with hardware acceleration. Our high-fidelity waveforms can be generated in under $1\,\mathrm{s}$, and achieve a mismatch of $\lesssim 5\times 10^{-4}$ against reference waveforms that take $\gtrsim 10^4$ times longer. This marks the first time that analysis-length waveforms with full harmonic content can be produced on timescales useful for direct implementation in LISA analysis algorithms.

preprint2019arXiv

Sampling from manifold-restricted distributions using tangent bundle projections

A common problem in Bayesian inference is the sampling of target probability distributions at sufficient resolution and accuracy to estimate the probability density, and to compute credible regions. Often by construction, many target distributions can be expressed as some higher-dimensional closed-form distribution with parametrically constrained variables, i.e., one that is restricted to a smooth submanifold of Euclidean space. I propose a derivative-based importance sampling framework for such distributions. A base set of $n$ samples from the target distribution is used to map out the tangent bundle of the manifold, and to seed $nm$ additional points that are projected onto the tangent bundle and weighted appropriately. The method essentially acts as an upsampling complement to any standard algorithm. It is designed for the efficient production of approximate high-resolution histograms from manifold-restricted Gaussian distributions, and can provide large computational savings when sampling directly from the target distribution is expensive.

preprint2016arXiv

Augmented kludge waveforms and Gaussian process regression for EMRI data analysis

Extreme-mass-ratio inspirals (EMRIs) will be an important type of astrophysical source for future space-based gravitational-wave detectors. There is a trade-off between accuracy and computational speed for the EMRI waveform templates required in the analysis of data from these detectors. We discuss how the systematic error incurred by using faster templates may be reduced with improved models such as augmented kludge waveforms, and marginalised over with statistical techniques such as Gaussian process regression.

preprint2016arXiv

Fast methods for training Gaussian processes on large data sets

Gaussian process regression (GPR) is a non-parametric Bayesian technique for interpolating or fitting data. The main barrier to further uptake of this powerful tool rests in the computational costs associated with the matrices which arise when dealing with large data sets. Here, we derive some simple results which we have found useful for speeding up the learning stage in the GPR algorithm, and especially for performing Bayesian model comparison between different covariance functions. We apply our techniques to both synthetic and real data and quantify the speed-up relative to using nested sampling to numerically evaluate model evidences.

preprint2016arXiv

Geometric creation of quantum vorticity

We consider superfluidity and quantum vorticity in rotating spacetimes. The system is described by a complex scalar satisfying a nonlinear Klein-Gordon equation. Rotation terms are identified and found to lead to the transfer of angular momentum of the spacetime to the scalar field. The scalar field responds by rotating, physically behaving as a superfluid, through the creation of quantized vortices. We demonstrate the vortex nucleation through numerical simulation.

preprint2016arXiv

Improving gravitational-wave parameter estimation using Gaussian process regression

Folding uncertainty in theoretical models into Bayesian parameter estimation is necessary in order to make reliable inferences. A general means of achieving this is by marginalizing over model uncertainty using a prior distribution constructed using Gaussian process regression (GPR). As an example, we apply this technique to the measurement of chirp mass using (simulated) gravitational-wave signals from binary black holes that could be observed using advanced-era gravitational-wave detectors. Unless properly accounted for, uncertainty in the gravitational-wave templates could be the dominant source of error in studies of these systems. We explain our approach in detail and provide proofs of various features of the method, including the limiting behavior for high signal-to-noise, where systematic model uncertainties dominate over noise errors. We find that the marginalized likelihood constructed via GPR offers a significant improvement in parameter estimation over the standard, uncorrected likelihood both in our simple one-dimensional study, and theoretically in general. We also examine the dependence of the method on the size of training set used in the GPR; on the form of covariance function adopted for the GPR, and on changes to the detector noise power spectral density.

preprint2016arXiv

Tunable compression of template banks for fast gravitational-wave detection and localisation

One strategy for reducing the online computational cost of matched-filter searches for gravitational waves is to introduce a compressed basis for the waveform template bank in a grid-based search. In this paper, we propose and investigate several tunable compression schemes for a general template bank. Through offline compression, such schemes are shown to yield faster detection and localisation of signals, along with moderately improved sensitivity and accuracy over coarsened banks at the same level of computational cost. This is potentially useful for any search involving template banks, and especially in the analysis of data from future space-based detectors such as eLISA, for which online grid searches are difficult due to the long-duration waveforms and large parameter spaces.

preprint2015arXiv

Improved analytic extreme-mass-ratio inspiral model for scoping out eLISA data analysis

The space-based gravitational-wave detector eLISA has been selected as the ESA L3 mission, and the mission design will be finalised by the end of this decade. To prepare for mission formulation over the next few years, several outstanding and urgent questions in data analysis will be addressed using mock data challenges, informed by instrument measurements from the LISA Pathfinder satellite launching at the end of 2015. These data challenges will require accurate and computationally affordable waveform models for anticipated sources such as the extreme-mass-ratio inspirals (EMRIs) of stellar-mass compact objects into massive black holes. Previous data challenges have made use of the well-known analytic EMRI waveforms of Barack and Cutler, which are extremely quick to generate but dephase relative to more accurate waveforms within hours, due to their mismatched radial, polar and azimuthal frequencies. In this paper, we describe an augmented Barack-Cutler model that uses a frequency map to the correct Kerr frequencies, along with updated evolution equations and a simple fit to a more accurate model. The augmented waveforms stay in phase for months and may be generated with virtually no additional computational cost.

preprint2014arXiv

Electromagnetic signatures of far-field gravitational radiation in the 1+3 approach

Gravitational waves from astrophysical sources can interact with background electromagnetic fields, giving rise to distinctive and potentially detectable electromagnetic signatures. In this paper, we study such interactions for far-field gravitational radiation using the 1+3 approach to relativity. Linearised equations for the electromagnetic field on perturbed Minkowski space are derived and solved analytically. The inverse Gertsenshtein conversion of gravitational waves in a static electromagnetic field is rederived, and the resultant electromagnetic radiation is shown to be significant for highly magnetised pulsars in compact binary systems. We also obtain a variety of nonlinear interference effects for interacting gravitational and electromagnetic waves, although wave-wave resonances previously described in the literature are absent when the electric-magnetic self-interaction is taken into account. The fluctuation and amplification of electromagnetic energy flux as the gravitational wave strength increases towards the gravitational-electromagnetic frequency ratio is a possible signature of gravitational radiation from extended astrophysical sources.

preprint2013arXiv

Galactic Escape Speeds in Mirror and Cold Dark Matter Models

The mirror dark matter (MDM) model of Berezhiani et al. has been shown to reproduce observed galactic rotational curves for a variety of spiral galaxies, and has been presented as an alternative to cold dark matter (CDM) models. We investigate possible additional tests involving the properties of stellar orbits, which may be used to discriminate between the two models. We demonstrate that in MDM and CDM models fitted equally well to a galactic rotational curve, one generally expects predictable differences in escape speeds from the disc. The recent radial velocity (RAVE) survey of the Milky Way has pinned down the escape speed from the solar neighbourhood to $v_{esc}=544^{+64}_{-46}$ km s$^{-1}$, placing an additional constraint on dark matter models. We have constructed an MDM model for the Milky Way based on its rotational curve, and find an escape speed that is just consistent with the observed value given the current errors, which lends credence to the viability of the MDM model. The Gaia-ESO spectroscopic survey is expected to lead to an even more precise estimate of the escape speed that will further constrain dark matter models. However, the largest differences in stellar escape speeds between both models are predicted for dark matter dominated dwarf galaxies such as DDO 154, and kinematical studies of such galaxies could prove key in establishing, or abolishing, the validity of the MDM model.

preprint2011arXiv

Interacting classical and quantum particles

We apply Hall and Reginatto's theory of interacting classical and quantum ensembles to harmonically coupled particles, with a view to understanding its experimental implications. This hybrid theory has no free parameters and makes distinctive predictions that should allow it to be experimentally distinguished from quantum mechanics. It also bears on the questions of quantum measurement and quantum gravity.