Showing 1–50 of 91 results for author: Chua, A

SPIRIT: Low Power Seizure Prediction using Unsupervised Online-Learning and Zoom Analog Frontends

Authors: Aviral Pandey, Adelson Chua, Ryan Kaveh, Justin Doong, Rikky Muller

Abstract: Early prediction of seizures and timely interventions are vital for improving patients' quality of life. While seizure prediction has been shown in software-based implementations, to enable timely warnings of upcoming seizures, prediction must be done on an edge device to reduce latency. Ideally, such devices must also be low-power and track long-term drifts to minimize maintenance from the user.… ▽ More Early prediction of seizures and timely interventions are vital for improving patients' quality of life. While seizure prediction has been shown in software-based implementations, to enable timely warnings of upcoming seizures, prediction must be done on an edge device to reduce latency. Ideally, such devices must also be low-power and track long-term drifts to minimize maintenance from the user. This work presents SPIRIT: Stochastic-gradient-descent-based Predictor with Integrated Retraining and In situ accuracy Tuning. SPIRIT is a complete system-on-a-chip (SoC) integrating an unsupervised online-learning seizure prediction classifier with eight 14.4 uW, 0.057 mm2, 90.5 dB dynamic range, Zoom Analog Frontends. SPIRIT achieves, on average, 97.5%/96.2% sensitivity/specificity respectively, predicting seizures an average of 8.4 minutes before they occur. Through its online learning algorithm, prediction accuracy improves by up to 15%, and prediction times extend by up to 7x, without any external intervention. Its classifier consumes 17.2 uW and occupies 0.14 mm2, the lowest reported for a prediction classifier by >134x in power and >5x in area. SPIRIT is also at least 5.6x more energy efficient than the state-of-the-art. △ Less

Submitted 7 September, 2024; originally announced September 2024.

Assessment of scattered and leakage radiation from ultra-portable digital chest X-ray systems: An independent study

Authors: Leonie E Paulis, Roald S Schnerr, Jarred Halton, Zhi Zhen Qin, Arlene Chua

Abstract: Ultraportable X-ray devices are ideal for TB screening in resource-limited settings. Unfortunately, guidelines on the radiation safety of these devices are lacking. The aim of the study was to determine the radiation dose by scattered and leakage radiation of ultraportable X-ray devices to provide a basis for these guidelines. Radiation dose measurements were performed with four ultraportable X-ra… ▽ More Ultraportable X-ray devices are ideal for TB screening in resource-limited settings. Unfortunately, guidelines on the radiation safety of these devices are lacking. The aim of the study was to determine the radiation dose by scattered and leakage radiation of ultraportable X-ray devices to provide a basis for these guidelines. Radiation dose measurements were performed with four ultraportable X-ray devices that meet the WHO-IAEA criteria. An anthropomorphic thorax phantom was positioned in posterior-anterior orientation in a clinically representative X-ray setup. X-ray exams were acquired with the following scan parameters: 90kV, 2.5mAs and maximum mAs, 1m and 1.8m source-to-skin distance. The entrance skin dose was measured at the center of the phantom. The scattered radiation dose was measured at 1m from the phantom as a function of scatter angle. Leakage radiation was measured at 0.5m from the X-ray tube with collimators closed and covered with additional lead. From the scatter measurements, the 'safe distance bunny' was determined, which was the minimum distance to the phantom to stay below the international dose limit to the public (1mSv per year) at a given workload. The longest distance (related to highest scatter dose) was observed behind the edges of the detector and back towards the X-ray tube, whereas the shortest distance (related to lowest scatter dose) was observed to the sides of the phantom. For the radiographer position, the total radiation dose by scattered and leakage radiation was determined in various scenarios. In most cases, the total radiation dose of ultraportable X-ray devices can be kept below 1mSv per year by employing basic radiation safety rules:1. Reduce time in the field, 2. Increase distance to the X-ray source, and 3. Use shielding measures (lead apron). Ultraportable X-ray devices can be safely used for TB screening when using adequate precautions. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Probing fundamental physics with Extreme Mass Ratio Inspirals: a full Bayesian inference for scalar charge

Authors: Lorenzo Speri, Susanna Barsanti, Andrea Maselli, Thomas P. Sotiriou, Niels Warburton, Maarten van de Meent, Alvin J. K. Chua, Ollie Burke, Jonathan Gair

Abstract: Extreme Mass Ratio Inspirals (EMRIs) are key sources for the future space-based gravitational wave detector LISA, and are considered promising probes of fundamental physics. Here, we present the first complete Bayesian analysis of EMRI signals in theories with an additional massless scalar, which could arise in an extension of General Relativity or of the Standard Model of Particle Physics. We dev… ▽ More Extreme Mass Ratio Inspirals (EMRIs) are key sources for the future space-based gravitational wave detector LISA, and are considered promising probes of fundamental physics. Here, we present the first complete Bayesian analysis of EMRI signals in theories with an additional massless scalar, which could arise in an extension of General Relativity or of the Standard Model of Particle Physics. We develop a waveform model accurate at adiabatic order for equatorial eccentric orbits around spinning black holes. Using full Bayesian inference, we forecast LISA's ability to probe the presence of new fundamental fields with EMRI observations. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Repeating Nuclear Transients as Candidate Electromagnetic Counterparts of LISA Extreme Mass Ratio Inspirals

Authors: Shubham Kejriwal, Vojtech Witzany, Michal Zajacek, Dheeraj R. Pasham, Alvin J. K. Chua

Abstract: Extreme-mass-ratio inspirals (EMRIs) are one of the primary targets for the recently adopted millihertz gravitational-wave (GW) observatory LISA. Some previous studies have argued that a fraction of all EMRIs form in matter-rich environments, and can potentially explain the dozens of soft X-ray band ($\sim 10^{-1} \rm keV$), low-frequency ($\sim 0.1$ mHz) periodic phenomena known as quasi-periodic… ▽ More Extreme-mass-ratio inspirals (EMRIs) are one of the primary targets for the recently adopted millihertz gravitational-wave (GW) observatory LISA. Some previous studies have argued that a fraction of all EMRIs form in matter-rich environments, and can potentially explain the dozens of soft X-ray band ($\sim 10^{-1} \rm keV$), low-frequency ($\sim 0.1$ mHz) periodic phenomena known as quasi-periodic eruptions (QPEs) and quasi-periodic oscillations (QPOs). Here, using a representative EMRI population retrofitted with cutoffs on LISA-band SNRs and luminosity distances to account for the sensitivity of current instruments, we estimate the mean frequency band in which QPEs and QPOs originating from detectable LISA EMRIs may be emitting an X-ray signal ``today'' (i.e., in 2024) to be $0.46 \pm 0.22$ mHz. We also model the well-known QPO source, RE J1034+396, which falls in this frequency band, as an EMRI assuming its primary black hole mass to be $10^6-10^7 M_\odot$. Through a prior-predictive analysis, we estimate the orbiting compact object's mass to be $46^{+ 10}_{-40} M_\odot$ and the source's LISA-band SNR as $\approx 14$, highlighting it as a candidate multi-messenger EMRI target. We also highlight the role of current and near-future X-ray and UV observatories in enabling multi-messenger observations of EMRIs in conjunction with LISA, and conclude with a discussion of caveats of the current analysis, such as the exclusion of eccentricity and inclination from the model, and the measurability of sub-solar mass compact object EMRIs. △ Less

Submitted 1 July, 2024; v1 submitted 1 April, 2024; originally announced April 2024.

Comments: (Before-proofs-accepted) 15 + 1 pages, 10 + 1 figures

Kazhdan-Lusztig map and Langlands duality

Authors: Anlong Chua

Abstract: Let $G$ be a connected reductive group over $\mathbb{C}$ with Weyl group $W$. Following a suggestion of Bezrukavnikov, we define a map from two-sided cells to conjugacy classes in $W$ using the geometry of the affine flag variety. This is an affine analog of the classical story of two-sided cells of $W$, special nilpotent orbits and special representations of $W$. The proof involves studying Sprin… ▽ More Let $G$ be a connected reductive group over $\mathbb{C}$ with Weyl group $W$. Following a suggestion of Bezrukavnikov, we define a map from two-sided cells to conjugacy classes in $W$ using the geometry of the affine flag variety. This is an affine analog of the classical story of two-sided cells of $W$, special nilpotent orbits and special representations of $W$. The proof involves studying Springer representations appearing in the cohomology of affine Springer fibers. This relies on tools developed by Yun in his papers on Global Springer theory. Our result provides evidence for a conjecture of Lusztig on strata in a connected reductive group. △ Less

Submitted 13 March, 2024; v1 submitted 11 March, 2024; originally announced March 2024.

Comments: Corrected typos

Generalizing Geometric Nonwindowed Scattering Transforms on Compact Riemannian Manifolds

Authors: Albert Chua, Yang Yang

Abstract: Let $\mathcal{M}$ be a compact, smooth, $n$-dimensional Riemannian manifold without boundary. In this paper, we generalize nonwindowed geometric scattering transforms, which we formulate as $\mathbf{L}^q(\mathcal{M})$ norms of a cascade of geometric wavelet transforms and modulus operators. We then provide weighted measures for these operators, prove that these operators are well-defined under spe… ▽ More Let $\mathcal{M}$ be a compact, smooth, $n$-dimensional Riemannian manifold without boundary. In this paper, we generalize nonwindowed geometric scattering transforms, which we formulate as $\mathbf{L}^q(\mathcal{M})$ norms of a cascade of geometric wavelet transforms and modulus operators. We then provide weighted measures for these operators, prove that these operators are well-defined under specific conditions on the manifold, invariant to the action of isometries, and stable to diffeomorphisms for $λ$-bandlimited functions. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: Initial draft. Comments welcome!

LISA Definition Study Report

Authors: Monica Colpi, Karsten Danzmann, Martin Hewitson, Kelly Holley-Bockelmann, Philippe Jetzer, Gijs Nelemans, Antoine Petiteau, David Shoemaker, Carlos Sopuerta, Robin Stebbins, Nial Tanvir, Henry Ward, William Joseph Weber, Ira Thorpe, Anna Daurskikh, Atul Deep, Ignacio Fernández Núñez, César García Marirrodriga, Martin Gehler, Jean-Philippe Halain, Oliver Jennrich, Uwe Lammers, Jonan Larrañaga, Maike Lieser, Nora Lützgendorf , et al. (86 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the e… ▽ More The Laser Interferometer Space Antenna (LISA) is the first scientific endeavour to detect and study gravitational waves from space. LISA will survey the sky for Gravitational Waves in the 0.1 mHz to 1 Hz frequency band which will enable the study of a vast number of objects ranging from Galactic binaries and stellar mass black holes in the Milky Way, to distant massive black-hole mergers and the expansion of the Universe. This definition study report, or Red Book, presents a summary of the very large body of work that has been undertaken on the LISA mission over the LISA definition phase. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 155 pages, with executive summary and table of contents

Impact of Correlations on the Modeling and Inference of Beyond Vacuum-GR Effects in Extreme-Mass-Ratio Inspirals

Authors: Shubham Kejriwal, Lorenzo Speri, Alvin J. K. Chua

Abstract: In gravitational-wave astronomy, extreme-mass-ratio-inspiral (EMRI) sources for the upcoming LISA observatory have the potential to serve as high-precision probes of astrophysical environments in galactic nuclei, and of potential deviations from general relativity (GR). Such ``beyond vacuum-GR'' effects are often modeled as perturbations to the evolution of vacuum EMRIs under GR. Previous studies… ▽ More In gravitational-wave astronomy, extreme-mass-ratio-inspiral (EMRI) sources for the upcoming LISA observatory have the potential to serve as high-precision probes of astrophysical environments in galactic nuclei, and of potential deviations from general relativity (GR). Such ``beyond vacuum-GR'' effects are often modeled as perturbations to the evolution of vacuum EMRIs under GR. Previous studies have reported unprecedented constraints on these effects by examining the inference of one effect at a time. However, a more realistic analysis would require the simultaneous inference of multiple beyond vacuum-GR effects. The parameters describing such effects are generally significantly correlated with each other and the vacuum EMRI parameters. We explicitly show how these correlations remain even if any modeled effect is absent in the actual signal, and how they cause inference bias when any effect in the signal is absent in the analysis model. This worsens the overall measurability of the whole parameter set, challenging the constraints found by previous studies, and posing a general problem for the modeling and inference of beyond vacuum-GR effects in EMRIs. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 6 + 3 pages, 3 figures

Waveform Modelling for the Laser Interferometer Space Antenna

Authors: LISA Consortium Waveform Working Group, Niayesh Afshordi, Sarp Akçay, Pau Amaro Seoane, Andrea Antonelli, Josu C. Aurrekoetxea, Leor Barack, Enrico Barausse, Robert Benkel, Laura Bernard, Sebastiano Bernuzzi, Emanuele Berti, Matteo Bonetti, Béatrice Bonga, Gabriele Bozzola, Richard Brito, Alessandra Buonanno, Alejandro Cárdenas-Avendaño, Marc Casals, David F. Chernoff, Alvin J. K. Chua, Katy Clough, Marta Colleoni, Mekhi Dhesi, Adrien Druart , et al. (121 additional authors not shown)

Abstract: LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic… ▽ More LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome. △ Less

Submitted 20 December, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

Calibrating approximate Bayesian credible intervals of gravitational-wave parameters

Authors: Ruiting Mao, Jeong Eun Lee, Ollie Burke, Alvin J. K. Chua, Matthew C. Edwards, Renate Meyer

Abstract: Approximations are commonly employed in realistic applications of scientific Bayesian inference, often due to convenience if not necessity. In the field of gravitational-wave (GW) data analysis, fast-to-evaluate but approximate waveform models of astrophysical GW signals are sometimes used in lieu of more accurate models to infer properties of a true GW signal buried within detector noise. In addi… ▽ More Approximations are commonly employed in realistic applications of scientific Bayesian inference, often due to convenience if not necessity. In the field of gravitational-wave (GW) data analysis, fast-to-evaluate but approximate waveform models of astrophysical GW signals are sometimes used in lieu of more accurate models to infer properties of a true GW signal buried within detector noise. In addition, a Fisher-information-based normal approximation to the posterior distribution can also be used to conduct inference in bulk, without the need for extensive numerical calculations such as Markov chain Monte Carlo (MCMC) simulations. Such approximations can generally lead to an inaccurate posterior distribution with poor statistical coverage of the true posterior. In this article, we present a novel calibration procedure that calibrates the credible sets for a family of approximate posterior distributions, to ensure coverage of the true posterior at a level specified by the analyst. Tools such as autoencoders and artificial neural networks are used within our calibration model to compress the data (for efficiency) and to perform tasks such as logistic regression. As a proof of principle, we demonstrate our formalism on the GW signal from a high-mass binary black hole merger, a promising source for the near-future space-based GW observatory LISA. △ Less

Submitted 30 January, 2024; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 24 pages, 12 figures

Journal ref: Phys. Rev. D 109, 083002 (2024)

Fast and Fourier: Extreme Mass Ratio Inspiral Waveforms in the Frequency Domain

Authors: Lorenzo Speri, Michael L. Katz, Alvin J. K. Chua, Scott A. Hughes, Niels Warburton, Jonathan E. Thompson, Christian E. A. Chapman-Bird, Jonathan R. Gair

Abstract: Extreme Mass Ratio Inspirals (EMRIs) are one of the key sources for future space-based gravitational wave interferometers. Measurements of EMRI gravitational waves are expected to determine the characteristics of their sources with sub-percent precision. However, their waveform generation is challenging due to the long duration of the signal and the high harmonic content. Here, we present the firs… ▽ More Extreme Mass Ratio Inspirals (EMRIs) are one of the key sources for future space-based gravitational wave interferometers. Measurements of EMRI gravitational waves are expected to determine the characteristics of their sources with sub-percent precision. However, their waveform generation is challenging due to the long duration of the signal and the high harmonic content. Here, we present the first ready-to-use Schwarzschild eccentric EMRI waveform implementation in the frequency domain for use with either graphics processing units (GPUs) or central processing units (CPUs). We present the overall waveform implementation and test the accuracy and performance of the frequency domain waveforms against the time domain implementation. On GPUs, the frequency domain waveform takes in median $0.044$ seconds to generate and is twice as fast to compute as its time domain counterpart when considering massive black hole masses $\geq 2 \times 10^6 \,{\rm M_\odot}$ and initial eccentricities $e_0 > 0.2$. On CPUs, the median waveform evaluation time is $5$ seconds, and it is five times faster in the frequency domain than in the time domain. Using a sparser frequency array can further speed up the waveform generation, reaching up to $ 0.3$ seconds. This enables us to perform, for the first time, EMRI parameter inference with fully relativistic waveforms on CPUs. Future EMRI models which encompass wider source characteristics (particularly black hole spin and generic orbit geometries) will require significantly more harmonics. Frequency-domain models will be essential analysis tools for these astrophysically realistic and important signals. △ Less

Submitted 15 January, 2024; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: 23 pages, 6 figures

Journal ref: Front. Appl. Math. Stat. 9 (2023)

Improving the scalability of Gaussian-process error marginalization in gravitational-wave inference

Authors: Miaoxin Liu, Xiao-Dong Li, Alvin J. K. Chua

Abstract: The accuracy of Bayesian inference can be negatively affected by the use of inaccurate forward models. In the case of gravitational-wave inference, accurate but computationally expensive waveform models are sometimes substituted with faster but approximate ones. The model error introduced by this substitution can be mitigated in various ways, one of which is by interpolating and marginalizing over… ▽ More The accuracy of Bayesian inference can be negatively affected by the use of inaccurate forward models. In the case of gravitational-wave inference, accurate but computationally expensive waveform models are sometimes substituted with faster but approximate ones. The model error introduced by this substitution can be mitigated in various ways, one of which is by interpolating and marginalizing over the error using Gaussian process regression. However, the use of Gaussian process regression is limited by the curse of dimensionality, which makes it less effective for analyzing higher-dimensional parameter spaces and longer signal durations. In this work, to address this limitation, we focus on gravitational-wave signals from extreme-mass-ratio inspirals as an example, and propose several significant improvements to the base method: an improved prescription for constructing the training set, GPU-accelerated training algorithms, and a new likelihood that better adapts the base method to the presence of detector noise. Our results suggest that the new method is more viable for the analysis of realistic gravitational-wave data. △ Less

Submitted 28 July, 2023; v1 submitted 14 July, 2023; originally announced July 2023.

Journal ref: Phys. Rev. D 108, 103027 (2023)

Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: II. Posterior predictive distributions and pseudo Bayes factors

Authors: Patrick M. Meyers, Katerina Chatziioannou, Michele Vallisneri, Alvin J. K. Chua

Abstract: The detection of nanoHertz gravitational waves through pulsar timing arrays hinges on identifying a common stochastic process affecting all pulsars in a correlated way across the sky. In the presence of other deterministic and stochastic processes affecting the time-of-arrival of pulses, a detection claim must be accompanied by a detailed assessment of the various physical or phenomenological mode… ▽ More The detection of nanoHertz gravitational waves through pulsar timing arrays hinges on identifying a common stochastic process affecting all pulsars in a correlated way across the sky. In the presence of other deterministic and stochastic processes affecting the time-of-arrival of pulses, a detection claim must be accompanied by a detailed assessment of the various physical or phenomenological models used to describe the data. In this study, we propose posterior predictive checks as a model-checking tool that relies on the predictive performance of the models with regards to new data. We derive and study predictive checks based on different components of the models, namely the Fourier coefficients of the stochastic process, the correlation pattern, and the timing residuals. We assess the ability of our checks to identify model misspecification in simulated datasets. We find that they can accurately flag a stochastic process spectral shape that deviates from the common power-law model as well as a stochastic process that does not display the expected angular correlation pattern. Posterior predictive likelihoods derived under different assumptions about the correlation pattern can further be used to establish detection significance. In the era of nanoHertz gravitational wave detection from different pulsar-timing datasets, such tests represent an essential tool in assessing data consistency and supporting astrophysical inference. △ Less

Submitted 12 February, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 18 pages, 9 figures

Posterior predictive checking for gravitational-wave detection with pulsar timing arrays: I. The optimal statistic

Authors: Michele Vallisneri, Patrick M. Meyers, Katerina Chatziioannou, Alvin J. K. Chua

Abstract: A gravitational-wave background can be detected in pulsar-timing-array data as Hellings--Downs correlations among the timing residuals measured for different pulsars. The optimal statistic implements this concept as a classical null-hypothesis statistical test: a null model with no correlations can be rejected if the observed value of the statistic is very unlikely under that model. To address the… ▽ More A gravitational-wave background can be detected in pulsar-timing-array data as Hellings--Downs correlations among the timing residuals measured for different pulsars. The optimal statistic implements this concept as a classical null-hypothesis statistical test: a null model with no correlations can be rejected if the observed value of the statistic is very unlikely under that model. To address the dependence of the statistic on the uncertain pulsar noise parameters, the pulsar-timing-array community has adopted a hybrid classical--Bayesian scheme (Vigeland et al. 2018) in which the posterior distribution of the noise parameters induces a posterior distribution for the statistic. In this article we propose a rigorous interpretation of the hybrid scheme as an instance of posterior predictive checking, and we introduce a new summary statistic (the Bayesian signal-to-noise ratio) that should be used to accurately quantify the statistical significance of an observation instead of the mean posterior signal-to-noise ratio, which does not support such a direct interpretation. In addition to falsifying the no-correlation hypothesis, the Bayesian signal-to-noise ratio can also provide evidence supporting the presence of Hellings--Downs correlations. We demonstrate our proposal with simulated datasets based on NANOGrav's 12.5-yr data release. We also establish a relation between the posterior distribution of the statistic and the Bayes factor in favor of correlations, thus calibrating the Bayes factor in terms of hypothesis-testing significance. △ Less

Submitted 12 February, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 12 pages, 8 figures

ChatGPT or academic scientist? Distinguishing authorship with over 99% accuracy using off-the-shelf machine learning tools

Authors: Heather Desaire, Aleesa E. Chua, Madeline Isom, Romana Jarosova, David Hua

Abstract: ChatGPT has enabled access to AI-generated writing for the masses, and within just a few months, this product has disrupted the knowledge economy, initiating a culture shift in the way people work, learn, and write. The need to discriminate human writing from AI is now both critical and urgent, particularly in domains like higher education and academic writing, where AI had not been a significant… ▽ More ChatGPT has enabled access to AI-generated writing for the masses, and within just a few months, this product has disrupted the knowledge economy, initiating a culture shift in the way people work, learn, and write. The need to discriminate human writing from AI is now both critical and urgent, particularly in domains like higher education and academic writing, where AI had not been a significant threat or contributor to authorship. Addressing this need, we developed a method for discriminating text generated by ChatGPT from (human) academic scientists, relying on prevalent and accessible supervised classification methods. We focused on how a particular group of humans, academic scientists, write differently than ChatGPT, and this targeted approach led to the discovery of new features for discriminating (these) humans from AI; as examples, scientists write long paragraphs and have a penchant for equivocal language, frequently using words like but, however, and although. With a set of 20 features, including the aforementioned ones and others, we built a model that assigned the author, as human or AI, at well over 99% accuracy, resulting in 20 times fewer misclassified documents compared to the field-leading approach. This strategy for discriminating a particular set of humans writing from AI could be further adapted and developed by others with basic skills in supervised classification, enabling access to many highly accurate and targeted models for detecting AI usage in academic writing and beyond. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Long Range Constraints for Neural Texture Synthesis Using Sliced Wasserstein Loss

Authors: Liping Yin, Albert Chua

Abstract: In the past decade, exemplar-based texture synthesis algorithms have seen strong gains in performance by matching statistics of deep convolutional neural networks. However, these algorithms require regularization terms or user-added spatial tags to capture long range constraints in images. Having access to a user-added spatial tag for all situations is not always feasible, and regularization terms… ▽ More In the past decade, exemplar-based texture synthesis algorithms have seen strong gains in performance by matching statistics of deep convolutional neural networks. However, these algorithms require regularization terms or user-added spatial tags to capture long range constraints in images. Having access to a user-added spatial tag for all situations is not always feasible, and regularization terms can be difficult to tune. Thus, we propose a new set of statistics for texture synthesis based on Sliced Wasserstein Loss, create a multi-scale method to synthesize textures without a user-added spatial tag, study the ability of our proposed method to capture long range constraints, and compare our results to other optimization-based, single texture synthesis algorithms. △ Less

Submitted 9 February, 2024; v1 submitted 20 November, 2022; originally announced November 2022.

Comments: Added extra ablation studies

On Generalizations of the Nonwindowed Scattering Transform

Authors: Albert Chua, Matthew Hirn, Anna Little

Abstract: In this paper, we generalize finite depth wavelet scattering transforms, which we formulate as $\Lb^q(\mathbb{R}^n)$ norms of a cascade of continuous wavelet transforms (or dyadic wavelet transforms) and contractive nonlinearities. We then provide norms for these operators, prove that these operators are well-defined, and are Lipschitz continuous to the action of $C^2$ diffeomorphisms in specific… ▽ More In this paper, we generalize finite depth wavelet scattering transforms, which we formulate as $\Lb^q(\mathbb{R}^n)$ norms of a cascade of continuous wavelet transforms (or dyadic wavelet transforms) and contractive nonlinearities. We then provide norms for these operators, prove that these operators are well-defined, and are Lipschitz continuous to the action of $C^2$ diffeomorphisms in specific cases. Lastly, we extend our results to formulate an operator invariant to the action of rotations $R \in \text{SO}(n)$ and an operator that is equivariant to the action of rotations of $R \in \text{SO}(n)$. △ Less

Submitted 13 September, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: Corrected small typos throughout. Diffeomorphism stability definition in introduction was changed. The mapping Phi isn't necessarily from V to V, which is reflected in the results we found

A one-stop function for gravitational-wave detection, identification and inference

Authors: Alvin J. K. Chua

Abstract: I define here a novel function on a modeled space of gravitational-wave signals, before studying its properties as a statistic for detection, as an objective function for identification, and as an effective likelihood function for inference. The main motivation behind this work is the open data-analysis problem for signals from extreme-mass-ratio inspirals, which is severely hindered by the presen… ▽ More I define here a novel function on a modeled space of gravitational-wave signals, before studying its properties as a statistic for detection, as an objective function for identification, and as an effective likelihood function for inference. The main motivation behind this work is the open data-analysis problem for signals from extreme-mass-ratio inspirals, which is severely hindered by the presence of strong non-local parameter degeneracy in the signal space. I demonstrate the utility of the proposed function for the analysis of such signals, and suggest various possible directions for future research. △ Less

Submitted 1 November, 2022; v1 submitted 17 May, 2022; originally announced May 2022.

Comments: Published version

Journal ref: Phys. Rev. D 106, 104051 (2022)

Modeling transient resonances in extreme-mass-ratio inspirals

Authors: Priti Gupta, Lorenzo Speri, Béatrice Bonga, Alvin J. K. Chua, Takahiro Tanaka

Abstract: Extreme-mass-ratio inspirals are one of the most exciting and promising target sources for space-based interferometers (such as LISA, TianQin). The observation of their emitted gravitational waves will offer stringent tests on general theory of relativity, and provide a wealth of information about the dense environment in galactic centers. To unlock such potential, it is necessary to correctly cha… ▽ More Extreme-mass-ratio inspirals are one of the most exciting and promising target sources for space-based interferometers (such as LISA, TianQin). The observation of their emitted gravitational waves will offer stringent tests on general theory of relativity, and provide a wealth of information about the dense environment in galactic centers. To unlock such potential, it is necessary to correctly characterize EMRI signals. However, resonances are a phenomena that occurs in EMRI systems and can impact parameter inference, and therefore the science outcome, if not properly modeled. Here, we explore how to model resonances and develop an efficient implementation. Our previous work has demonstrated that tidal resonances induced by the tidal field of a nearby astrophysical object alters the orbital evolution, leading to a significant dephasing across observable parameter space. Here, we extensively explore a more generic model for the tidal perturber with additional resonance combinations, to study the dependence of resonance strength on the intrinsic orbital and tidal parameters. To analyze the resonant signals, accurate templates that correctly incorporate the effects of the tidal field are required. The evolution through resonances is obtained using a step function, whose amplitude is calculated using an analytic interpolation of the resonance jumps. We benchmark this procedure by comparing our approximate method to a numerical evolution. We find that there is no significant error caused by this simplified prescription, as far as the astronomically reasonable range in the parameter space is concerned. Further, we use Fisher matrices to study both the measurement precision of parameters and the systematic bias due to inaccurate modeling. Modeling of self-force resonances can also be carried out using the implementation presented in this study, which will be crucial for EMRI waveform modeling. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: 18 pages, 14 figures

Journal ref: Phys. Rev. D 106, 104001 (2022)

Assessing the data-analysis impact of LISA orbit approximations using a GPU-accelerated response model

Authors: Michael L. Katz, Jean-Baptiste Bayle, Alvin J. K. Chua, Michele Vallisneri

Abstract: The analysis of gravitational wave (GW) datasets is based on the comparison of measured time series with theoretical templates of the detector's response to a variety of source parameters. For LISA, the main scientific observables will be the so-called time-delay interferometry (TDI) combinations, which suppress the otherwise overwhelming laser noise. Computing the TDI response to GW involves proj… ▽ More The analysis of gravitational wave (GW) datasets is based on the comparison of measured time series with theoretical templates of the detector's response to a variety of source parameters. For LISA, the main scientific observables will be the so-called time-delay interferometry (TDI) combinations, which suppress the otherwise overwhelming laser noise. Computing the TDI response to GW involves projecting the GW polarizations onto the LISA constellation arms, and then combining projections delayed by a multiple of the light propagation time along the arms. Both computations are difficult to perform efficiently for generic LISA orbits and GW signals. Various approximations are currently used in practice, e.g., assuming constant and equal armlengths, which yields analytical TDI expressions. In this article, we present 'fastlisaresponse', a new efficient GPU-accelerated code that implements the generic TDI response to GWs in the time domain. We use it to characterize the parameter-estimation bias incurred by analyzing loud Galactic-binary signals using the equal-armlength approximation. We conclude that equal-armlength parameter-estimation codes should be upgraded to the generic response if they are to achieve optimal accuracy for high (but reasonable) SNR sources within the actual LISA data. △ Less

Submitted 1 September, 2022; v1 submitted 13 April, 2022; originally announced April 2022.

Comments: 15 pages, 7 figures, 2 tables

Journal ref: Phys. Rev. D 106, 103001 (2022)

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

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

Authors: Soichiro Isoyama, Ryuichi Fujita, Alvin J. K. Chua, Hiroyuki Nakano, Adam Pound, Norichika Sago

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 12 June, 2022; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: 9+5 pages,4+1 figures,5PN-e10 Teukolsky data are all available at the Black Hole Perturbation Club: https://sites.google.com/view/bhpc1996/data v2:updated references,expanded supplemental discussion of initial phases.v3:contains minor corrections and more references.Matches version published in PRL

Report number: YITP-21-114,KUNS-2901,OCU-PHYS-551,AP-GR-175

Journal ref: Phys. Rev. Lett. 128, 231101 (2022)

SOUL: An Energy-Efficient Unsupervised Online Learning Seizure Detection Classifier

Authors: Adelson Chua, Michael I. Jordan, Rikky Muller

Abstract: Implantable devices that record neural activity and detect seizures have been adopted to issue warnings or trigger neurostimulation to suppress epileptic seizures. Typical seizure detection systems rely on high-accuracy offline-trained machine learning classifiers that require manual retraining when seizure patterns change over long periods of time. For an implantable seizure detection system, a l… ▽ More Implantable devices that record neural activity and detect seizures have been adopted to issue warnings or trigger neurostimulation to suppress epileptic seizures. Typical seizure detection systems rely on high-accuracy offline-trained machine learning classifiers that require manual retraining when seizure patterns change over long periods of time. For an implantable seizure detection system, a low power, at-the-edge, online learning algorithm can be employed to dynamically adapt to the neural signal drifts, thereby maintaining high accuracy without external intervention. This work proposes SOUL: Stochastic-gradient-descent-based Online Unsupervised Logistic regression classifier. After an initial offline training phase, continuous online unsupervised classifier updates are applied in situ, which improves sensitivity in patients with drifting seizure features. SOUL was tested on two human electroencephalography (EEG) datasets: the CHB-MIT scalp EEG dataset, and a long (>100 hours) NeuroVista intracranial EEG dataset. It was able to achieve an average sensitivity of 97.5% and 97.9% for the two datasets respectively, at >95% specificity. Sensitivity improved by at most 8.2% on long-term data when compared to a typical seizure detection classifier. SOUL was fabricated in TSMC's 28 nm process occupying 0.1 mm2 and achieves 1.5 nJ/classification energy efficiency, which is at least 24x more efficient than state-of-the-art. △ Less

Submitted 17 May, 2022; v1 submitted 1 October, 2021; originally announced October 2021.

Comments: Copyright 2022 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works

Journal ref: IEEE Journal of Solid-State Circuits, 2022

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

Authors: Alvin J. K. Chua, Curt J. Cutler

Abstract: 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… ▽ More 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. △ Less

Submitted 31 December, 2022; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 39 pages, 22 figures. This paper is long - but mainly because it is written in an expository/pedagogical manner, is thoroughly signposted, and contains multiple large figures

Journal ref: Phys. Rev. D 106, 124046 (2022)

The Effect of Mission Duration on LISA Science Objectives

Authors: Pau Amaro Seoane, Manuel Arca Sedda, Stanislav Babak, Christopher P. L. Berry, Emanuele Berti, Gianfranco Bertone, Diego Blas, Tamara Bogdanović, Matteo Bonetti, Katelyn Breivik, Richard Brito, Robert Caldwell, Pedro R. Capelo, Chiara Caprini, Vitor Cardoso, Zack Carson, Hsin-Yu Chen, Alvin J. K. Chua, Irina Dvorkin, Zoltan Haiman, Lavinia Heisenberg, Maximiliano Isi, Nikolaos Karnesis, Bradley J. Kavanagh, Tyson B. Littenberg , et al. (16 additional authors not shown)

Abstract: 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 asses… ▽ More 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. △ Less

Submitted 12 January, 2022; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 50 pages, 19 figures, 5 tables. Matches version published in GERG

FastEMRIWaveforms: New tools for millihertz gravitational-wave data analysis

Authors: Michael L. Katz, Alvin J. K. Chua, Lorenzo Speri, Niels Warburton, Scott A. Hughes

Abstract: 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 w… ▽ More 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. △ Less

Submitted 27 September, 2021; v1 submitted 9 April, 2021; originally announced April 2021.

Comments: 26 pages, 12 Figures, FastEMRIWaveforms Package: bhptoolkit.org/FastEMRIWaveforms/

Journal ref: Phys. Rev. D 104, 064047 (2021)

Importance of tidal resonances in extreme-mass-ratio inspirals

Authors: Priti Gupta, Béatrice Bonga, Alvin J. K. Chua, Takahiro Tanaka

Abstract: Extreme mass ratio inspirals (EMRIs) will be important sources for future space-based gravitational-wave detectors. In recent work, tidal resonances in binary orbital evolution induced by the tidal field of nearby stars or black holes have been identified as being potentially significant in the context of extreme mass-ratio inspirals. These resonances occur when the three orbital frequencies descr… ▽ More Extreme mass ratio inspirals (EMRIs) will be important sources for future space-based gravitational-wave detectors. In recent work, tidal resonances in binary orbital evolution induced by the tidal field of nearby stars or black holes have been identified as being potentially significant in the context of extreme mass-ratio inspirals. These resonances occur when the three orbital frequencies describing the orbit are commensurate. During the resonance, the orbital parameters of the small body experience a jump leading to a shift in the phase of the gravitational waveform. In this paper, we treat the tidal perturber as stationary and restricted to the equatorial plane, and present a first study of how common and important such resonances are over the entire orbital parameter space. We find that a large proportion of inspirals encounter a low-order resonance in the observationally important regime. While the instantaneous effect of a tidal resonance is small, its effect on the accumulated phase of the gravitational waveform of an EMRI system can be significant due to its many cycles in band; we estimate that the effect is detectable for a significant fraction of sources. We also provide fitting formulae for the induced change in the constants of motion of the orbit due to the tidal resonance for several low-order resonances. △ Less

Submitted 27 July, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 19 pages, 14 figures

Journal ref: Phys. Rev. D 104, 044056 (2021)

Adiabatic waveforms for extreme mass-ratio inspirals via multivoice decomposition in time and frequency

Authors: Scott A. Hughes, Niels Warburton, Gaurav Khanna, Alvin J. K. Chua, Michael L. Katz

Abstract: We compute adiabatic waveforms for extreme mass-ratio inspirals (EMRIs) by "stitching" together a long inspiral waveform from a sequence of waveform snapshots, each of which corresponds to a particular geodesic orbit. We show that the complicated total waveform can be regarded as a sum of "voices." Each voice evolves in a simple way on long timescales, a property which can be exploited to efficien… ▽ More We compute adiabatic waveforms for extreme mass-ratio inspirals (EMRIs) by "stitching" together a long inspiral waveform from a sequence of waveform snapshots, each of which corresponds to a particular geodesic orbit. We show that the complicated total waveform can be regarded as a sum of "voices." Each voice evolves in a simple way on long timescales, a property which can be exploited to efficiently produce waveform models that faithfully encode the properties of EMRI systems. We look at examples for a range of different orbital geometries: spherical orbits, equatorial eccentric orbits, and one example of generic (inclined and eccentric) orbits. To our knowledge, this is the first calculation of a generic EMRI waveform that uses strong-field radiation reaction. We examine waveforms in both the time and frequency domains. Although EMRIs evolve slowly enough that the stationary phase approximation (SPA) to the Fourier transform is valid, the SPA calculation must be done to higher order for some voices, since their instantaneous frequency can change from chirping forward ($\dot f > 0$) to chirping backward ($\dot f < 0$). The approach we develop can eventually be extended to more complete EMRI waveform models, for example to include effects neglected by the adiabatic approximation such as the conservative self force and spin-curvature coupling. △ Less

Submitted 28 March, 2024; v1 submitted 4 February, 2021; originally announced February 2021.

Comments: 31 pages, 17 figures. Update: corrected a sign error in a term in Eq. (B5). This factor was correct in the related code, it was purely a typo in the manuscript. An erratum to PRD has been submitted. Many thanks to Angelica Albertini for bringing this error to our attention. Update 2: corrected a missing factor of $m$ in Eq. (3.30)

Journal ref: Phys. Rev. D 103, 104014 (2021)

Cyclic orthogonal convolutions for long-range integration of features

Authors: Federica Freddi, Jezabel R Garcia, Michael Bromberg, Sepehr Jalali, Da-Shan Shiu, Alvin Chua, Alberto Bernacchia

Abstract: In Convolutional Neural Networks (CNNs) information flows across a small neighbourhood of each pixel of an image, preventing long-range integration of features before reaching deep layers in the network. We propose a novel architecture that allows flexible information flow between features $z$ and locations $(x,y)$ across the entire image with a small number of layers. This architecture uses a cyc… ▽ More In Convolutional Neural Networks (CNNs) information flows across a small neighbourhood of each pixel of an image, preventing long-range integration of features before reaching deep layers in the network. We propose a novel architecture that allows flexible information flow between features $z$ and locations $(x,y)$ across the entire image with a small number of layers. This architecture uses a cycle of three orthogonal convolutions, not only in $(x,y)$ coordinates, but also in $(x,z)$ and $(y,z)$ coordinates. We stack a sequence of such cycles to obtain our deep network, named CycleNet. As this only requires a permutation of the axes of a standard convolution, its performance can be directly compared to a CNN. Our model obtains competitive results at image classification on CIFAR-10 and ImageNet datasets, when compared to CNNs of similar size. We hypothesise that long-range integration favours recognition of objects by shape rather than texture, and we show that CycleNet transfers better than CNNs to stylised images. On the Pathfinder challenge, where integration of distant features is crucial, CycleNet outperforms CNNs by a large margin. We also show that even when employing a small convolutional kernel, the size of receptive fields of CycleNet reaches its maximum after one cycle, while conventional CNNs require a large number of layers. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 11 pages, 5 figures

Exciting modes due to the aberration of gravitational waves: Measurability for extreme-mass-ratio inspirals

Authors: Alejandro Torres-Orjuela, Pau Amaro Seoane, Zeyuan Xuan, Alvin J. K. Chua, María J. B. Rosell, Xian Chen

Abstract: Gravitational waves from a source moving relative to us can suffer from special-relativistic effects such as aberration. The required velocities for these to be significant are on the order of $1000\,\textrm{km s}^{-1}$. This value corresponds to the velocity dispersion that one finds in clusters of galaxies. Hence, we expect a large number of gravitational-wave sources to have such effects imprin… ▽ More Gravitational waves from a source moving relative to us can suffer from special-relativistic effects such as aberration. The required velocities for these to be significant are on the order of $1000\,\textrm{km s}^{-1}$. This value corresponds to the velocity dispersion that one finds in clusters of galaxies. Hence, we expect a large number of gravitational-wave sources to have such effects imprinted in their signals. In particular, the signal from a moving source will have its higher modes excited, i.e., $(3,3)$ and beyond. We derive expressions describing this effect, and study its measurability for the specific case of a circular, non-spinning extreme-mass-ratio inspiral. We find that the excitation of higher modes by a peculiar velocity of $1000\,\textrm{km\,s}^{-1}$ is detectable for such inspirals with signal-to-noise ratios of $\gtrsim20$. Using a Fisher matrix analysis, we show that the velocity of the source can be measured to a precision of just a few percent for a signal-to-noise ratio of 100. If the motion of the source is ignored parameter estimates could be biased, e.g., the estimated masses of the components through a Doppler shift. Conversely, by including this effect in waveform models, we could measure the velocity dispersion of clusters of galaxies at distances inaccessible to light. △ Less

Submitted 11 July, 2021; v1 submitted 29 October, 2020; originally announced October 2020.

Comments: Accepted PRL

Journal ref: Phys. Rev. Lett. 127, 041102 (2021)

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

Authors: Alvin J. K. Chua, Michael L. Katz, Niels Warburton, Scott A. Hughes

Abstract: 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 w… ▽ More 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. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Lett. 126, 051102 (2021)

On parametric tests of relativity with false degrees of freedom

Authors: Alvin J. K. Chua, Michele Vallisneri

Abstract: 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 relativi… ▽ More 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. △ Less

Submitted 16 June, 2020; originally announced June 2020.

Comments: 4 pages, 2 figures

GPU-accelerated massive black hole binary parameter estimation with LISA

Authors: Michael L. Katz, Sylvain Marsat, Alvin J. K. Chua, Stanislav Babak, Shane L. Larson

Abstract: 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… ▽ More 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. △ Less

Submitted 13 July, 2020; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: 27 pages, 13 figures

Journal ref: Phys. Rev. D 102, 023033 (2020)

Efficient attention guided 5G power amplifier digital predistortion

Authors: Alexandru Cioba, Alvin Chua, Da-shan Shiu, Ting-Hsun Kuo, Chia-Sheng Peng

Abstract: We investigate neural network (NN) assisted techniques for compensating the non-linear behaviour and the memory effect of a 5G PA through digital predistortion (DPD). Traditionally, the most prevalent compensation technique computes the compensation element using a Memory Polynomial Model (MPM). Various neural network proposals have been shown to improve on this performance. However, thus far they… ▽ More We investigate neural network (NN) assisted techniques for compensating the non-linear behaviour and the memory effect of a 5G PA through digital predistortion (DPD). Traditionally, the most prevalent compensation technique computes the compensation element using a Memory Polynomial Model (MPM). Various neural network proposals have been shown to improve on this performance. However, thus far they mostly come with prohibitive training or inference costs for real world implementations. In this paper, we propose a DPD architecture that builds upon the practical MPM formulation governed by neural attention. Our approach enables a set of MPM DPD components to individually learn to target different regions of the data space, combining their outputs for a superior overall compensation. Our method produces similar performance to that of higher capacity NN models with minimal complexity. Finally, we view our approach as a framework that can be extended to a wide variety of local compensator types. △ Less

Submitted 30 March, 2020; originally announced March 2020.

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

Authors: Alvin J. K. Chua, Natalia Korsakova, Christopher J. Moore, Jonathan R. Gair, Stanislav Babak

Abstract: 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… ▽ More 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. △ Less

Submitted 13 February, 2020; v1 submitted 24 December, 2019; originally announced December 2019.

Comments: Published version

Journal ref: Phys. Rev. D 101, 044027 (2020)

Learning Bayesian posteriors with neural networks for gravitational-wave inference

Authors: Alvin J. K. Chua, Michele Vallisneri

Abstract: 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 detec… ▽ More 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. △ Less

Submitted 29 January, 2020; v1 submitted 12 September, 2019; originally announced September 2019.

Comments: (Superior-to-)published version; source code and trained networks available at https://github.com/vallis/truebayes

Journal ref: Phys. Rev. Lett. 124, 041102 (2020)

The unique potential of extreme mass-ratio inspirals for gravitational-wave astronomy

Authors: Christopher P. L. Berry, Scott A. Hughes, Carlos F. Sopuerta, Alvin J. K. Chua, Anna Heffernan, Kelly Holley-Bockelmann, Deyan P. Mihaylov, M. Coleman Miller, Alberto Sesana

Abstract: The inspiral of a stellar-mass compact object into a massive ($\sim 10^{4}$-$10^{7} M_{\odot}$) black hole produces an intricate gravitational-wave signal. Due to the extreme-mass ratios involved, these systems complete $\sim 10^{4}$-$10^{5}$ orbits, most of them in the strong-field region of the massive black hole, emitting in the frequency range $\sim10^{-4}-1~$Hz. This makes them prime sources… ▽ More The inspiral of a stellar-mass compact object into a massive ($\sim 10^{4}$-$10^{7} M_{\odot}$) black hole produces an intricate gravitational-wave signal. Due to the extreme-mass ratios involved, these systems complete $\sim 10^{4}$-$10^{5}$ orbits, most of them in the strong-field region of the massive black hole, emitting in the frequency range $\sim10^{-4}-1~$Hz. This makes them prime sources for the space-based observatory LISA (Laser Interferometer Space Antenna). LISA observations will enable high-precision measurements of the physical characteristics of these extreme-mass-ratio inspirals (EMRIs): redshifted masses, massive black hole spin and orbital eccentricity can be determined with fractional errors $\sim 10^{-4}$-$10^{-6}$, the luminosity distance with better than $\sim 10\%$ precision, and the sky localization to within a few square degrees. EMRIs will provide valuable information about stellar dynamics in galactic nuclei, as well as precise data about massive black hole populations, including the distribution of masses and spins. They will enable percent-level measurements of the multipolar structure of massive black holes, precisely testing the strong-gravity properties of their spacetimes. EMRIs may also provide cosmographical data regarding the expansion of the Universe if inferred source locations can be correlated with galaxy catalogs. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 13 pages, 1 figure; white paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics); Section III is arXiv-only

Sampling from manifold-restricted distributions using tangent bundle projections

Authors: Alvin J. K. Chua

Abstract: 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… ▽ More 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. △ Less

Submitted 16 October, 2019; v1 submitted 13 November, 2018; originally announced November 2018.

Comments: Published version; Python implementation available at https://github.com/alvincjk/sampling-manifold-restricted-gaussians

Journal ref: Stat. Comput. 30, 587 (2020)

Reduced-order modeling with artificial neurons for gravitational-wave inference

Authors: Alvin J. K. Chua, Chad R. Galley, Michele Vallisneri

Abstract: Gravitational-wave data analysis is rapidly absorbing techniques from deep learning, with a focus on convolutional networks and related methods that treat noisy time series as images. We pursue an alternative approach, in which waveforms are first represented as weighted sums over reduced bases (reduced-order modeling); we then train artificial neural networks to map gravitational-wave source para… ▽ More Gravitational-wave data analysis is rapidly absorbing techniques from deep learning, with a focus on convolutional networks and related methods that treat noisy time series as images. We pursue an alternative approach, in which waveforms are first represented as weighted sums over reduced bases (reduced-order modeling); we then train artificial neural networks to map gravitational-wave source parameters into basis coefficients. Statistical inference proceeds directly in coefficient space, where it is theoretically straightforward and computationally efficient. The neural networks also provide analytic waveform derivatives, which are useful for gradient-based sampling schemes. We demonstrate fast and accurate coefficient interpolation for the case of a four-dimensional binary-inspiral waveform family, and discuss promising applications of our framework in parameter estimation. △ Less

Submitted 30 May, 2019; v1 submitted 13 November, 2018; originally announced November 2018.

Comments: Published version

Journal ref: Phys. Rev. Lett. 122, 211101 (2019)

Tests of General Relativity with GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy , et al. (1123 additional authors not shown)

Abstract: The recent discovery by Advanced LIGO and Advanced Virgo of a gravitational wave signal from a binary neutron star inspiral has enabled tests of general relativity (GR) with this new type of source. This source, for the first time, permits tests of strong-field dynamics of compact binaries in presence of matter. In this paper, we place constraints on the dipole radiation and possible deviations fr… ▽ More The recent discovery by Advanced LIGO and Advanced Virgo of a gravitational wave signal from a binary neutron star inspiral has enabled tests of general relativity (GR) with this new type of source. This source, for the first time, permits tests of strong-field dynamics of compact binaries in presence of matter. In this paper, we place constraints on the dipole radiation and possible deviations from GR in the post-Newtonian coefficients that govern the inspiral regime. Bounds on modified dispersion of gravitational waves are obtained; in combination with information from the observed electromagnetic counterpart we can also constrain effects due to large extra dimensions. Finally, the polarization content of the gravitational wave signal is studied. The results of all tests performed here show good agreement with GR. △ Less

Submitted 29 July, 2019; v1 submitted 1 November, 2018; originally announced November 2018.

Comments: 15 pages, 4 figures. Matches journal submission

Report number: LIGO-P1800059

Journal ref: Phys. Rev. Lett. 123, 011102 (2019)

A Fermi Gamma-ray Burst Monitor Search for Electromagnetic Signals Coincident with Gravitational-Wave Candidates in Advanced LIGO's First Observing Run

Authors: The Fermi Gamma-ray Burst Monitor Team, The LIGO Scientific Collaboration, the Virgo Collaboration, :, E. Burns, A. Goldstein, C. M. Hui, L. Blackburn, M. S. Briggs, V. Connaughton, R. Hamburg, D. Kocevski, P. Veres, C. A. Wilson-Hodge, E. Bissaldi, W. H. Cleveland, M. M. Giles, B. Mailyan, C. A. Meegan, W. A. Paciesas, S. Poolakkil, R. D. Preece, J. L. Racusin, O. J. Roberts, A. von Kienlin , et al. (1139 additional authors not shown)

Abstract: We present a search for prompt gamma-ray counterparts to compact binary coalescence gravitational wave (GW) candidates from Advanced LIGO's first observing run (O1). As demonstrated by the multimessenger observations of GW170817/GRB 170817A, electromagnetic and GW observations provide complementary information about the astrophysical source and, in the case of weaker candidates, may strengthen the… ▽ More We present a search for prompt gamma-ray counterparts to compact binary coalescence gravitational wave (GW) candidates from Advanced LIGO's first observing run (O1). As demonstrated by the multimessenger observations of GW170817/GRB 170817A, electromagnetic and GW observations provide complementary information about the astrophysical source and, in the case of weaker candidates, may strengthen the case for an astrophysical origin. Here we investigate low-significance GW candidates from the O1 compact-binary coalescence searches using the Fermi Gamma-ray Burst Monitor (GBM), leveraging its all-sky and broad energy coverage. Candidates are ranked and compared to background to measure significance. Those with false alarm rates of less than 10^-5 Hz (about one per day) are used as the search sample for gamma-ray follow-up. No GW candidates were found to be coincident with gamma-ray transients independently identified by blind searches of the GBM data. In addition, GW candidate event times were followed up by a separate targeted search of GBM data. Among the resulting GBM events, the two with lowest false alarm rates were the gamma-ray transient GW150914-GBM presented in Connaughton et al. (2016) and a solar flare in chance coincidence with a GW candidate. △ Less

Submitted 18 November, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: Published in ApJ

Search for gravitational waves from a long-lived remnant of the binary neutron star merger GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy , et al. (1116 additional authors not shown)

Abstract: One unanswered question about the binary neutron star coalescence GW170817 is the nature of its post-merger remnant. A previous search for post-merger gravitational waves targeted high-frequency signals from a possible neutron star remnant with a maximum signal duration of 500 s. Here we revisit the neutron star remnant scenario with a focus on longer signal durations up until the end of the Secon… ▽ More One unanswered question about the binary neutron star coalescence GW170817 is the nature of its post-merger remnant. A previous search for post-merger gravitational waves targeted high-frequency signals from a possible neutron star remnant with a maximum signal duration of 500 s. Here we revisit the neutron star remnant scenario with a focus on longer signal durations up until the end of the Second Advanced LIGO-Virgo Observing run, 8.5 days after the coalescence of GW170817. The main physical scenario for such emission is the power-law spindown of a massive magnetar-like remnant. We use four independent search algorithms with varying degrees of restrictiveness on the signal waveformand different ways of dealing with noise artefacts. In agreement with theoretical estimates, we find no significant signal candidates. Through simulated signals, we quantify that with the current detector sensitivity, nowhere in the studied parameter space are we sensitive to a signal from more than 1 Mpc away, compared to the actual distance of 40 Mpc. This study however serves as a prototype for post-merger analyses in future observing runs with expected higher sensitivity. △ Less

Submitted 4 October, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: main paper: 9 pages and 3 figures; total with appendices: 24 pages and 9 figures. Full UL tables available as MRT ancillary files

Report number: LIGO-P1800195

Journal ref: The Astrophysical Journal 875:160 (2019)

Constraining the p-mode--g-mode tidal instability with GW170817

Authors: The LIGO Scientific Collaboration, The Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy , et al. (1114 additional authors not shown)

Abstract: We analyze the impact of a proposed tidal instability coupling $p$-modes and $g$-modes within neutron stars on GW170817. This non-resonant instability transfers energy from the orbit of the binary to internal modes of the stars, accelerating the gravitational-wave driven inspiral. We model the impact of this instability on the phasing of the gravitational wave signal using three parameters per sta… ▽ More We analyze the impact of a proposed tidal instability coupling $p$-modes and $g$-modes within neutron stars on GW170817. This non-resonant instability transfers energy from the orbit of the binary to internal modes of the stars, accelerating the gravitational-wave driven inspiral. We model the impact of this instability on the phasing of the gravitational wave signal using three parameters per star: an overall amplitude, a saturation frequency, and a spectral index. Incorporating these additional parameters, we compute the Bayes Factor ($\ln B^{pg}_{!pg}$) comparing our $p$-$g$ model to a standard one. We find that the observed signal is consistent with waveform models that neglect $p$-$g$ effects, with $\ln B^{pg}_{!pg} = 0.03^{+0.70}_{-0.58}$ (maximum a posteriori and 90% credible region). By injecting simulated signals that do not include $p$-$g$ effects and recovering them with the $p$-$g$ model, we show that there is a $\simeq 50\%$ probability of obtaining similar $\ln B^{pg}_{!pg}$ even when $p$-$g$ effects are absent. We find that the $p$-$g$ amplitude for 1.4 $M_\odot$ neutron stars is constrained to $\lesssim \text{few}\times10^{-7}$, with maxima a posteriori near $\sim 10^{-7}$ and $p$-$g$ saturation frequency $\sim 70\, \mathrm{Hz}$. This suggests that there are less than a few hundred excited modes, assuming they all saturate by wave breaking. For comparison, theoretical upper bounds suggest a $p$-$g$ amplitude $\lesssim 10^{-6}$ and $\lesssim 10^{3}$ modes saturating by wave breaking. Thus, the measured constraints only rule out extreme values of the $p$-$g$ parameters. They also imply that the instability dissipates $\lesssim 10^{51}\, \mathrm{ergs}$ over the entire inspiral, i.e., less than a few percent of the energy radiated as gravitational waves. △ Less

Submitted 19 September, 2018; v1 submitted 26 August, 2018; originally announced August 2018.

Comments: 7 pages, 2 figures

Journal ref: Phys. Rev. Lett. 122, 061104 (2019)

Search for sub-solar mass ultracompact binaries in Advanced LIGO's first observing run

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato , et al. (1113 additional authors not shown)

Abstract: We present the first Advanced LIGO and Advanced Virgo search for ultracompact binary systems with component masses between 0.2 $M_\odot$ - 1.0 $M_\odot$ using data taken between September 12, 2015 and January 19, 2016. We find no viable gravitational wave candidates. Our null result constrains the coalescence rate of monochromatic (delta function) distributions of non-spinning (0.2 $M_\odot$, 0.2… ▽ More We present the first Advanced LIGO and Advanced Virgo search for ultracompact binary systems with component masses between 0.2 $M_\odot$ - 1.0 $M_\odot$ using data taken between September 12, 2015 and January 19, 2016. We find no viable gravitational wave candidates. Our null result constrains the coalescence rate of monochromatic (delta function) distributions of non-spinning (0.2 $M_\odot$, 0.2 $M_\odot$) ultracompact binaries to be less than $1.0 \times 10^6 \text{Gpc}^{-3} \text{yr}^{-1}$ and the coalescence rate of a similar distribution of (1.0 $M_\odot$, 1.0 $M_\odot$) ultracompact binaries to be less than $1.9 \times 10^4 \text{Gpc}^{-3} \text{yr}^{-1}$ (at 90 percent confidence). Neither black holes nor neutron stars are expected to form below ~ 1 solar mass through conventional stellar evolution, though it has been proposed that similarly low mass black holes could be formed primordially through density fluctuations in the early universe. Under a particular primordial black hole binary formation scenario, we constrain monochromatic primordial black hole populations of 0.2 $M_\odot$ to be less than $33\%$ of the total dark matter density and monochromatic populations of 1.0 $M_\odot$ to be less than $5\%$ of the dark matter density. The latter strengthens the presently placed bounds from micro-lensing surveys of MAssive Compact Halo Objects (MACHOs) provided by the MACHO and EROS collaborations. △ Less

Submitted 15 August, 2018; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: 12 pages, 3 figures

Report number: LIGO-DCC-P1800158-v12

Journal ref: Phys. Rev. Lett. 121, 231103 (2018)

GW170817: Measurements of Neutron Star Radii and Equation of State

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy , et al. (1127 additional authors not shown)

Abstract: On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of th… ▽ More On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function $p(ρ)$ of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as $R_1=10.8^{+2.0}_{-1.7}$ km for the heavier star and $R_2= 10.7^{+2.1}_{-1.5}$ km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than $1.97 \,M_\odot$ as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain $R_1= 11.9^{+1.4}_{-1.4}$ km and $R_2= 11.9^{+1.4}_{-1.4}$ km at the 90% credible level. Finally, we obtain constraints on $p(ρ)$ at supranuclear densities, with pressure at twice nuclear saturation density measured at $3.5^{+2.7}_{-1.7}\times 10^{34} \,\mathrm{dyn}/\mathrm{cm}^{2}$ at the 90% level. △ Less

Submitted 15 October, 2018; v1 submitted 29 May, 2018; originally announced May 2018.

Comments: 10 pages, 3 figures; v2 matches published version; data associated with the figures can be found at https://dcc.ligo.org/LIGO-P1800115/public

Report number: LIGO-P1800115

Journal ref: Phys. Rev. Lett. 121, 161101 (2018)

Properties of the binary neutron star merger GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy , et al. (1126 additional authors not shown)

Abstract: On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary's properties, including component masses, spin… ▽ More On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary's properties, including component masses, spins, and tidal parameters, using the known source location, improved modeling, and recalibrated Virgo data. We extend the range of gravitational-wave frequencies considered down to 23 Hz, compared to 30 Hz in the initial analysis. We also compare results inferred using several signal models, which are more accurate and incorporate additional physical effects as compared to the initial analysis. We improve the localization of the gravitational-wave source to a 90% credible region of $16~\mathrm{deg}^2$. We find tighter constraints on the masses, spins, and tidal parameters, and continue to find no evidence for nonzero component spins. The component masses are inferred to lie between 1.00 and 1.89 $M_\odot$ when allowing for large component spins, and to lie between 1.16 and 1.60 $M_\odot$ (with a total mass $2.73^{+0.04}_{-0.01} \, M_\odot$) when the spins are restricted to be within the range observed in Galactic binary neutron stars. Under minimal assumptions about the nature of the compact objects, our constraints for the tidal deformability parameter $\tilde Λ$ are $(0,630)$ when we allow for large component spins, and $300^{+420}_{-230}$ (using a 90% highest posterior density interval) when restricting the magnitude of the component spins, ruling out several equation-of-state models at the 90% credible level. Finally, with LIGO and GEO600 data, we use a Bayesian analysis to place upper limits on the amplitude and spectral energy density of a possible post-merger signal. (Abridged) △ Less

Submitted 22 January, 2019; v1 submitted 29 May, 2018; originally announced May 2018.

Comments: 32 pages, 15 figures. Version 3 is the final published version; results are unchanged compared to previous versions. Data behind the figures, including posterior samples, are available at dcc.ligo.org/LIGO-P1800061/public. The gravitational wave strain data for this event are available at the LIGO Open Science Center, gw-openscience.org/events/GW170817

Journal ref: Phys. Rev. X 9, 011001 (2019)

Towards a framework for testing general relativity with extreme-mass-ratio-inspiral observations

Authors: Alvin J. K. Chua, Sonke Hee, Will J. Handley, Edward Higson, Christopher J. Moore, Jonathan R. Gair, Michael P. Hobson, Anthony N. Lasenby

Abstract: Extreme-mass-ratio-inspiral observations from future space-based gravitational-wave detectors such as LISA will enable strong-field tests of general relativity with unprecedented precision, but at prohibitive computational cost if existing statistical techniques are used. In one such test that is currently employed for LIGO black-hole binary mergers, generic deviations from relativity are represen… ▽ More Extreme-mass-ratio-inspiral observations from future space-based gravitational-wave detectors such as LISA will enable strong-field tests of general relativity with unprecedented precision, but at prohibitive computational cost if existing statistical techniques are used. In one such test that is currently employed for LIGO black-hole binary mergers, generic deviations from relativity are represented by $N$ deformation parameters in a generalised waveform model; the Bayesian evidence for each of its $2^N$ combinatorial submodels is then combined into a posterior odds ratio for modified gravity over relativity in a null-hypothesis test. We adapt and apply this test to a generalised model for extreme-mass-ratio inspirals constructed on deformed black-hole spacetimes, and focus our investigation on how computational efficiency can be increased through an evidence-free method of model selection. This method is akin to the algorithm known as product-space Markov chain Monte Carlo, but uses nested sampling and improved error estimates from a rethreading technique. We perform benchmarking and robustness checks for the method, and find order-of-magnitude computational gains over regular nested sampling in the case of synthetic data generated from the null model. △ Less

Submitted 25 April, 2018; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: Published version

Journal ref: Mon. Notices Royal Astron. Soc. 478, 28 (2018)

A Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca , et al. (1075 additional authors not shown)

Abstract: The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational w… ▽ More The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually-unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically-polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy-densities of tensor, vector, and scalar modes at 95% credibility to $Ω^T_0 < 5.6 \times 10^{-8}$, $Ω^V_0 < 6.4\times 10^{-8}$, and $Ω^S_0 < 1.1\times 10^{-7}$ at a reference frequency $f_0 = 25$ Hz. △ Less

Submitted 2 October, 2019; v1 submitted 27 February, 2018; originally announced February 2018.

Comments: Minor updates to match published version

Journal ref: Phys. Rev. Lett. 120, 201102 (2018)

Full Band All-sky Search for Periodic Gravitational Waves in the O1 LIGO Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Afrough, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, B. Allen, G. Allen, A. Allocca, P. A. Altin , et al. (1077 additional authors not shown)

Abstract: We report on a new all-sky search for periodic gravitational waves in the frequency band 475-2000 Hz and with a frequency time derivative in the range of [-1.0e-8, +1e-9] Hz/s. Potential signals could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run O1. No gravitational wave… ▽ More We report on a new all-sky search for periodic gravitational waves in the frequency band 475-2000 Hz and with a frequency time derivative in the range of [-1.0e-8, +1e-9] Hz/s. Potential signals could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGO's first observational run O1. No gravitational wave signals were observed, and upper limits were placed on their strengths. For completeness, results from the separately published low frequency search 20-475 Hz are included as well. Our lowest upper limit on worst-case (linearly polarized) strain amplitude h_0 is 4e-25 near 170 Hz, while at the high end of our frequency range we achieve a worst-case upper limit of 1.3e-24. For a circularly polarized source (most favorable orientation), the smallest upper limit obtained is ~1.5e-25. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Journal ref: Phys. Rev. D 97, 102003 (2018)

Coprime Subdegrees of Twisted Wreath Permutation Groups

Authors: Alexander Y. Chua, Michael Giudici, Luke Morgan

Abstract: Dolfi, Guralnick, Praeger and Spiga asked if there exist infinitely many primitive groups of twisted wreath type with nontrivial coprime subdegrees. Here we settle this question in the affirmative. We construct infinite families of primitive twisted wreath permutation groups with nontrivial coprime subdegrees. In particular, we define a primitive twisted wreath group $G(m,q)$ constructed from the… ▽ More Dolfi, Guralnick, Praeger and Spiga asked if there exist infinitely many primitive groups of twisted wreath type with nontrivial coprime subdegrees. Here we settle this question in the affirmative. We construct infinite families of primitive twisted wreath permutation groups with nontrivial coprime subdegrees. In particular, we define a primitive twisted wreath group $G(m,q)$ constructed from the nonabelian simple group $\text{PSL}(2,q)$ and a primitive permutation group of diagonal type with socle $\text{PSL}(2,q)^m$, and determine all values of $m$ and $q$ for which $G(m,q)$ has nontrivial coprime subdegrees. In the case where $m=2$ and $q\notin\{7,11,29\}$ we obtain a full classification of all pairs of nontrivial coprime subdegrees. △ Less

Submitted 8 January, 2018; originally announced January 2018.

MSC Class: 20B15

Journal ref: Proceedings of the Edinburgh Mathematical Society 62 (2019) 1137-1162