Source author record

Ik Siong Heng

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

17 published item(s)

preprint2022arXiv

Utilizing Gaussian mixture models in all-sky searches for short-duration gravitational wave bursts

Coherent WaveBurst is a generic, multidetector gravitational wave burst search based on the excess power approach. The coherent WaveBurst algorithm currently employed in the all-sky short-duration gravitational wave burst search uses a conditional approach on selected attributes in the multidimensional event attribute space to distinguish between noisy events from that of astrophysical origin. We have been developing a supervised machine learning approach based on the Gaussian mixture modeling to model the attribute space for signals as well as noise events to enhance the probability of burst detection [Gayathri et al.Phys. Rev. D 102, 104023 (2020)]. We further extend the GMM approach to the all-sky short-duration coherent WaveBurst search as a postprocessing step on events from the first half of the third observing run (O3a). We show an improvement in sensitivity to generic gravitational wave burst signal morphologies as well as the astrophysical source such as core-collapse supernova models due to the application of our Gaussian mixture model approach to coherent WaveBurst triggers. The Gaussian mixture model method recovers the gravitational wave signals from massive compact binary coalescences identified by coherent WaveBurst targeted for binary black holes in GWTC-2, with better significance than the all-sky coherent WaveBurst search. No additional significant gravitational wave bursts are observed.

preprint2021arXiv

Bayesian parameter estimation using conditional variational autoencoders for gravitational-wave astronomy

Gravitational wave (GW) detection is now commonplace and as the sensitivity of the global network of GW detectors improves, we will observe $\mathcal{O}(100)$s of transient GW events per year. The current methods used to estimate their source parameters employ optimally sensitive but computationally costly Bayesian inference approaches where typical analyses have taken between 6 hours and 5 days. For binary neutron star and neutron star black hole systems prompt counterpart electromagnetic (EM) signatures are expected on timescales of 1 second -- 1 minute and the current fastest method for alerting EM follow-up observers, can provide estimates in $\mathcal{O}(1)$ minute, on a limited range of key source parameters. Here we show that a conditional variational autoencoder pre-trained on binary black hole signals can return Bayesian posterior probability estimates. The training procedure need only be performed once for a given prior parameter space and the resulting trained machine can then generate samples describing the posterior distribution $\sim 6$ orders of magnitude faster than existing techniques.

preprint2021arXiv

Exploring the sky localization and early warning capabilities of third generation gravitational wave detectors in three-detector network configurations

This work characterises the sky localization and early warning performance of networks of third generation gravitational wave detectors, consisting of different combinations of detectors with either the Einstein Telescope or Cosmic Explorer configuration in sites in North America, Europe and Australia. Using a Fisher matrix method which includes the effect of earth rotation, we estimate the sky localization uncertainty for $1.4\text{M}\odot$-$1.4\text{M}\odot$ binary neutron star mergers at distances $40\text{Mpc}$, $200\text{Mpc}$, $400\text{Mpc}$, $800\text{Mpc}$, $1600\text{Mpc}$, and an assumed astrophysical population up to redshift of 2 to characterize its performance for binary neutron star observations. We find that, for binary neutron star mergers at $200\text{Mpc}$ and a network consisting of the Einstein Telescope, Cosmic Explorer and an extra Einstein Telescope-like detector in Australia(2ET1CE), the upper limit of the size of the 90% credible region for the best localized 90% signals is $0.25\text{deg}^2$. For the simulated astrophysical distribution, this upper limit is $91.79\text{deg}^2$. If the Einstein Telescope-like detector in Australia is replaced with a Cosmic Explorer-like detector(1ET2CE), for $200\text{Mpc}$ case, the upper limit is $0.18\text{deg}^2$, while for astrophysical distribution, it is $56.77\text{deg}^2$. We note that the 1ET2CE network can detect 7.2% more of the simulated astrophysical population than the 2ET1CE network. In terms of early warning performance, we find that a network of 2ET1CE and 1ET2CE networks can both provide early warnings of the order of 1 hour prior to merger with sky localization uncertainties of 30 square degrees or less. Our study concludes that the 1ET2CE network is a good compromise between binary neutron stars detection rate, sky localization and early warning capabilities.

preprint2020arXiv

A Precessing Numerical Relativity Waveform Surrogate Model for Binary Black Holes: A Gaussian Process Regression Approach

Gravitational wave astrophysics relies heavily on the use of matched filtering both to detect signals in noisy data from detectors, and to perform parameter estimation on those signals. Matched filtering relies upon prior knowledge of the signals expected to be produced by a range of astrophysical systems, such as binary black holes. These waveform signals can be computed using numerical relativity techniques, where the Einstein field equations are solved numerically, and the signal is extracted from the simulation. Numerical relativity simulations are, however, computationally expensive, leading to the need for a surrogate model which can predict waveform signals in regions of the physical parameter space which have not been probed directly by simulation. We present a method for producing such a surrogate using Gaussian process regression which is trained directly on waveforms generated by numerical relativity. This model returns not just a single interpolated value for the waveform at a new point, but a full posterior probability distribution on the predicted value. This model is therefore an ideal component in a Bayesian analysis framework, through which the uncertainty in the interpolation can be taken into account when performing parameter estimation of signals.

preprint2020arXiv

Comparing Short Gamma-Ray Burst Jet Structure Models

A structured gamma-ray burst jet could explain the dimness of the prompt emission observed from GRB$\,170817$A but the exact form of this structure is still ambiguous. However, with the promise of future joint gravitational wave and gamma-ray burst observations, we shall be able to examine populations of binary neutron star mergers rather than a case-by-case basis. We present an analysis that considers gravitational wave triggered binary neutron star events both with and without short gamma-ray burst counterparts assuming that events without a counterpart were observed off-axis. This allows for Bayes factors to be calculated to compare different jet structure models. We perform model comparison between a Gaussian and power-law apparent jet structure on simulated data to demonstrate that the correct model can be distinguished with a log Bayes factor of $>5$ after less than 100 events. Constraints on the apparent structure jet model parameters are also made. After 25(100) events the angular width of the core of a power-law jet structure can be constrained within a $90\%$ credible interval of width $ \sim9.1(4.4)^{\circ} $, and the outer beaming angle to be within $\sim19.9(8.5)^{\circ}$. Similarly we show the width of a Gaussian jet structure to be constrained to $\sim2.8(1.6)^{\circ}$.

preprint2020arXiv

Optic disc and fovea localisation in ultra-widefield scanning laser ophthalmoscope images captured in multiple modalities

We propose a convolutional neural network for localising the centres of the optic disc (OD) and fovea in ultra-wide field of view scanning laser ophthalmoscope (UWFoV-SLO) images of the retina. Images captured in both reflectance and autofluorescence (AF) modes, and central pole and eyesteered gazes, were used. The method achieved an OD localisation accuracy of 99.4% within one OD radius, and fovea localisation accuracy of 99.1% within one OD radius on a test set comprising of 1790 images. The performance of fovea localisation in AF images was comparable to the variation between human annotators at this task. The laterality of the image (whether the image is of the left or right eye) was inferred from the OD and fovea coordinates with an accuracy of 99.9%

preprint2019arXiv

Detection and Classification of Supernova Gravitational Waves Signals: A Deep Learning Approach

We demonstrate the application of a convolutional neural network to the gravitational wave signals from core collapse supernovae. Using simulated time series of gravitational wave detectors, we show that based on the explosion mechanisms, a convolutional neural network can be used to detect and classify the gravitational wave signals buried in noise. For the waveforms used in the training of the convolutional neural network, our results suggest that a network of advanced LIGO, advanced VIRGO and KAGRA, or a network of LIGO A+, advanced VIRGO and KAGRA is likely to detect a magnetorotational core collapse supernovae within the Large and Small Magellanic Clouds, or a Galactic event if the explosion mechanism is the neutrino-driven mechanism. By testing the convolutional neural network with waveforms not used for training, we show that the true alarm probabilities are 52% and 83% at 60 kpc for waveforms R3E1AC and R4E1FC L. For waveforms s20 and SFHx at 10 kpc, the true alarm probabilities are 70% and 93% respectively. All at false alarm probability equal to 10%.

preprint2016arXiv

Inferring the core-collapse supernova explosion mechanism with gravitational waves

A detection of a core-collapse supernova (CCSN) gravitational-wave (GW) signal with an Advanced LIGO and Virgo detector network may allow us to measure astrophysical parameters of the dying massive star. GWs are emitted from deep inside the core and, as such, they are direct probes of the CCSN explosion mechanism. In this study we show how we can determine the CCSN explosion mechanism from a GW supernova detection using a combination of principal component analysis and Bayesian model selection. We use simulations of GW signals from CCSN exploding via neutrino-driven convection and rapidly-rotating core collapse. Previous studies have shown that the explosion mechanism can be determined using one LIGO detector and simulated Gaussian noise. As real GW detector noise is both non-stationary and non-Gaussian we use real detector noise from a network of detectors with a sensitivity altered to match the advanced detectors design sensitivity. For the first time we carry out a careful selection of the number of principal components to enhance our model selection capabilities. We show that with an advanced detector network we can determine if the CCSN explosion mechanism is neutrino-driven convection for sources in our Galaxy and rapidly-rotating core collapse for sources out to the Large Magellanic Cloud.

preprint2015arXiv

Classification methods for noise transients in advanced gravitational-wave detectors

Noise of non-astrophysical origin will contaminate science data taken by the Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) and Advanced Virgo gravitational-wave detectors. Prompt characterization of instrumental and environmental noise transients will be critical for improving the sensitivity of the advanced detectors in the upcoming science runs. During the science runs of the initial gravitational-wave detectors, noise transients were manually classified by visually examining the time-frequency scan of each event. Here, we present three new algorithms designed for the automatic classification of noise transients in advanced detectors. Two of these algorithms are based on Principal Component Analysis. They are Principal Component Analysis for Transients (PCAT), and an adaptation of LALInference Burst (LIB). The third algorithm is a combination of an event generator called Wavelet Detection Filter (WDF) and machine learning techniques for classification. We test these algorithms on simulated data sets, and we show their ability to automatically classify transients by frequency, SNR and waveform morphology.

preprint2015arXiv

Global Optimization for Future Gravitational Wave Detectors' Sites

We consider the optimal site selection of future generations of gravitational wave detectors. Previously, Raffai et al. optimized a 2-detector network with a combined figure of merit. This optimization was extended to networks with more than two detectors in a limited way by first fixing the parameters of all other component detectors. In this work we now present a more general optimization that allows the locations of all detectors to be simultaneously chosen. We follow the definition of Raffai et al. on the metric that defines the suitability of a certain detector network. Given the locations of the component detectors in the network, we compute a measure of the network's ability to distinguish the polarization, constrain the sky localization and reconstruct the parameters of a gravitational wave source. We further define the `flexibility index' for a possible site location, by counting the number of multi-detector networks with a sufficiently high Figure of Merit that include that site location. We confirm the conclusion of Raffai et al., that in terms of flexibility index as defined in this work, Australia hosts the best candidate site to build a future generation gravitational wave detector. This conclusion is valid for either a 3-detector network or a 5-detector network. For a 3-detector network site locations in Northern Europe display a comparable flexibility index to sites in Australia. However for a 5-detector network, Australia is found to be a clearly better candidate than any other location.

preprint2014arXiv

A Bayesian approach to multi-messenger astronomy: Identification of gravitational-wave host galaxies

We present a general framework for incorporating astrophysical information into Bayesian parameter estimation techniques used by gravitational wave data analysis to facilitate multi-messenger astronomy. Since the progenitors of transient gravitational wave events, such as compact binary coalescences, are likely to be associated with a host galaxy, improvements to the source sky location estimates through the use of host galaxy information are explored. To demonstrate how host galaxy properties can be included, we simulate a population of compact binary coalescences and show that for ~8.5% of simulations with in 200Mpc, the top ten most likely galaxies account for a ~50% of the total probability of hosting a gravitational wave source. The true gravitational wave source host galaxy is in the top ten galaxy candidates ~10% of the time. Furthermore, we show that by including host galaxy information, a better estimate of the inclination angle of a compact binary gravitational wave source can be obtained. We also demonstrate the flexibility of our method by incorporating the use of either B or K band into our analysis.

preprint2014arXiv

Efficient Exploration of Multi-Modal Posterior Distributions

The Markov Chain Monte Carlo (MCMC) algorithm is a widely recognised as an efficient method for sampling a specified posterior distribution. However, when the posterior is multi-modal, conventional MCMC algorithms either tend to become stuck in one local mode, become non-Markovian or require an excessively long time to explore the global properties of the distribution. We propose a novel variant of MCMC, mixed MCMC, which exploits a specially designed proposal density to allow the generation candidate points from any of a number of different modes. This new method is efficient by design, and is strictly Markovian. We present our method and apply it to a toy model inference problem to demonstrate its validity.

preprint2014arXiv

Enhancing gravitational wave astronomy with galaxy catalogues

Joint gravitational wave (GW) and electromagnetic (EM) observations, as a key research direction in multi-messenger astronomy, will provide deep insight into the astrophysics of a vast range of astronomical phenomena. Uncertainties in the source sky location estimate from gravitational wave observations mean follow-up observatories must scan large portions of the sky for a potential companion signal. A general frame of joint GW-EM observations is presented by a multi-messenger observational triangle. Using a Bayesian approach to multi-messenger astronomy, we investigate the use of galaxy catalogue and host galaxy information to reduce the sky region over which follow-up observatories must scan, as well as study its use for improving the inclination angle estimates for coalescing binary compact objects. We demonstrate our method using a simulated neutron stars inspiral signal injected into simulated Advanced detectors noise and estimate the injected signal sky location and inclination angle using the Gravitational Wave Galaxy Catalogue. In this case study, the top three candidates in rank have $72\%$, $15\%$ and $8\%$ posterior probability of being the host galaxy, receptively. The standard deviation of cosine inclination angle (0.001) of the neutron stars binary using gravitational wave-galaxy information is much smaller than that (0.02) using only gravitational wave posterior samples.

preprint2014arXiv

Glitch or Anti-Glitch: a Bayesian View

The sudden spin-down in the rotation of magnetar 1E 2259+586 observed by Archibald et al. (2013) was a rare event. However this particular event, referred to as an anti-glitch, was followed by another event which Archibald et al. (2013) suggested could either be a conventional glitch or another anti-glitch. Although there is no accompanied radiation activity or pulse profile change, there is decisive evidence for the existence of the second timing event, judging from the timing data. We apply Bayesian Model Selection to quantitatively determine which of these possibilities better explains the observed data. We show that the observed data strongly supports the presence of two successive anti-glitches with a Bayes Factor, often called the odds ratio, greater than 40. Furthermore, we show that the second anti-gtlich has an associated frequency change $Δν$ of $-8.2 \times 10^{-8}$ Hz. We discuss the implications of these results for possible physical mechanisms behind this anti-glitch.

preprint2014arXiv

Investigating Binary Black Hole Mergers with Principal Component Analysis

Despite recent progress in numerical simulations of the coalescence of binary black hole systems, highly asymmetric spinning systems and the construction of accurate physical templates remain challenging and computationally expensive. We explore the feasibility of a prompt and robust test of whether the signals exhibit evidence for generic features that can educate new simulations. We form catalogs of numerical relativity waveforms with distinct physical effects and compute the relative probability that a gravitational wave signal belongs to each catalog. We introduce an algorithm designed to perform this task for coalescence signals using principal component analysis of waveform catalogs and Bayesian model selection and demonstrate its effectiveness.

preprint2013arXiv

Optimal Networks of Future Gravitational-Wave Telescopes

We aim to find the optimal site locations for a hypothetical network of 1-3 triangular gravitational-wave telescopes. We define the following N-telescope figures of merit (FoMs) and construct three corresponding metrics: (a) capability of reconstructing the signal polarization; (b) accuracy in source localization; and (c) accuracy in reconstructing the parameters of a standard binary source. We also define a combined metric that takes into account the three FoMs with practically equal weight. After constructing a geomap of possible telescope sites, we give the optimal 2-telescope networks for the four FoMs separately in example cases where the location of the first telescope has been predetermined. We found that based on the combined metric, placing the first telescope to Australia provides the most options for optimal site selection when extending the network with a second instrument. We suggest geographical regions where a potential second and third telescope could be placed to get optimal network performance in terms of our FoMs. Additionally, we use a similar approach to find the optimal location and orientation for the proposed LIGO-India detector within a five-detector network with Advanced LIGO (Hanford), Advanced LIGO (Livingston), Advanced Virgo, and KAGRA. We found that the FoMs do not change greatly in sites within India, though the network can suffer a significant loss in reconstructing signal polarizations if the orientation angle of an L-shaped LIGO-India is not set to the optimal value of ~58.2deg(+k*90deg) (measured counterclockwise from East to the bisector of the arms).

preprint2013arXiv

The Transient Gravitational-Wave Sky

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.