Next Article in Journal
Scalar Induced Gravitational Waves Review
Next Article in Special Issue
A Comprehensive Study of Bright Fermi-GBM Short Gamma-ray Bursts: I. Multi-Pulse Lightcurves and Multi-Component Spectra
Previous Article in Journal
Multimessenger Astronomy with Neutrinos
Previous Article in Special Issue
A Bayesian Inference Framework for Gamma-ray Burst Afterglow Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multimodal Analysis of Gravitational Wave Signals and Gamma-Ray Bursts from Binary Neutron Star Mergers

1
European Gravitational Observatory (EGO), Via E. Amaldi, 5, 56021 Cascina, Italy
2
Scuola Normale Superiore (SNS), Piazza dei Cavalieri, 7, 56126 Pisa, Italy
3
INFN, Sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
4
INFN, Sezione di Roma Tor Vergata, Via della Ricerca Scientifica, 1, 00133 Rome, Italy
5
Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, 00-716 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Universe 2021, 7(11), 394; https://doi.org/10.3390/universe7110394
Submission received: 27 September 2021 / Revised: 18 October 2021 / Accepted: 19 October 2021 / Published: 21 October 2021
(This article belongs to the Special Issue Waiting for GODOT—Present and Future of Multi-Messenger Astronomy)

Abstract

:
A major boost in the understanding of the universe was given by the revelation of the first coalescence event of two neutron stars (GW170817) and the observation of the same event across the entire electromagnetic spectrum. With third-generation gravitational wave detectors and the new astronomical facilities, we expect many multi-messenger events of the same type. We anticipate the need to analyse the data provided to us by such events not only to fulfil the requirements of real-time analysis, but also in order to decipher the event in its entirety through the information emitted in the different messengers using machine learning. We propose a change in the paradigm in the way we do multi-messenger astronomy, simultaneously using the complete information generated by violent phenomena in the Universe. What we propose is the application of a multimodal machine learning approach to characterize these events.

1. Introduction

The detection of gravitational waves (GWs) from the inspiral phase and coalescence of a pair of neutron stars (NS) on 17 August 2017 [1] and the following observations of the event in its electromagnetic (EM) counterparts (see [2] and references therein) marked the beginning of multi-messenger astronomy with GWs.
For the first time, we observed the coalescence of two NSs through GWs and EM radiation across the entire electromagnetic spectrum, thanks to the participation of more than 70 astronomical observatories to the EM follow-up campaign. Multi-messenger astronomy opens up not only new scenarios for the observation of the Universe and new perspectives for the investigations of astronomical objects, but also new challenges for the way of extracting all information that these astrophysical events bring with them. The synergy between the information that only GWs can provide and the concomitant observations through other detectors of the EM and neutrino counterparts can strongly accelerate our knowledge of the Universe. It is clear that multi-messenger astronomy discloses the need for new paradigms for data analysis and introduces new challenges for real-time analysis, and there are many efforts ongoing to face them that involve the use of machine learning techniques (see, e.g., [3,4,5,6,7,8,9,10,11,12]). Multimodal machine learning MMML analysis is efficiently applied in many fields of data analysis for the more inclusive interpretation of events where several modalities are concurrent, such as in a video with audio; images with captions; or images, text, and sound [13]. To our knowledge, these techniques have never been applied in the interpretation of astrophysical data, where signals of different nature can be almost simultaneous. In this paper we introduce, for the first time, a multimodal machine learning analysis applied to astrophysical transient signals such as the case of GWs and gamma-ray bursts (GRBs).
In Section 2, we will describe the importance of multi-messenger astronomy; in Section 3, we will introduce the multimodal analysis and some examples of application; in Section 4, we will describe how we can implement the multimodal analysis for multi-messenger events; and in Section 5, we will report a proof of concept for the the application of MMML to GW and GRB data. In the conclusion, we will discuss other possible applications.

2. Multi-Messenger Observations as a Powerful Tool to Investigate the Extreme Universe

The joint observation of GW170817 and its EM counterparts clearly demonstrated the enormous informative power of multi-messenger astronomy with GWs. For instance, the joint detection of GW170817 and GRB 170817A represents the first direct proof that NS-NS mergers are progenitors of short GRBs [14]. In addition to this, the observed time delay between the GW and the gamma-ray signal (∼1.7 s) allowed us to put constraints on the difference between the speed of gravity and the speed of light, which has been estimated to be between −3 × 10 15 and + 7 × 10 16 times the speed of light [14]. After the joint detection of GW170817 and GRB 170817A and the release of a three-detector, a well-constrained GW skymap has been key for the identification of other EM counterparts [2]: this allowed us to obtain more insights into the physics of the source. For instance, the multi-wavelength EM observations associated with GW170817 allowed us to infer some basic properties of short GRB jets. The temporal evolution of the X-ray [15] and radio [16] light curves, together with the very low gamma-ray luminosity of GRB 170817A, suggested two possible scenarios: an off-axis GRB with a relativistic, structured jet or a “cocoon” emission from the relativistic jet shocking its surrounding non-relativistic material [14]. Subsequent very long baseline interferometry observations have been crucial to discriminate between these two scenarios: they allowed astrophysicists to put constraints on the size of the source and on its displacement, which were found to be consistent with a structured, relativistic jet [17,18].
The detection of the optical/NIR counterpart to GW170817 (AT2017gf), first reported by [19] and later by other teams (see [2] and references therein), allowed for the first time to identify the host galaxy of a GW event and led to the first spectroscopic identification of a kilonova [20,21], thus expanding our knowledge of heavy element nucleosynthesis in the Universe. The joint GW and kilonova observation also allowed us to investigate in more detail the neutron star equation of state (EOS). For instance, [22] found the lower bound on the tidal deformability parameter through the interpretation of the UV/optical/IR counterpart of GW170817 with kilonova models, combined with new numerical relativity results; by combining this result with the constraints obtained with GW data alone, they have shown that both extremely stiff and soft EOSs are tentatively ruled out. In addition, Ref. [23] presented a Bayesian parameter estimation combining information from GW170817, GRB170817A, and AT2017gfo, and with this analysis, they were able to obtain multi-messenger constraints on the EOS and on the binary properties. More recently, [24] performed a joint analysis of GW170817, GRB170817A, and AT2017gfo, together with GW190425, and they combined these with previous measurements of pulsars using X-ray and radio observations, as well as and nuclear-theory computations, to put constraints on the EOS. A deeper knowledge of the EOS is also fundamental to understand which is the outcome of coalescing binary systems and therefore to constrain the short GRB central engine (see, e.g., [25]).
Finally, the estimate of the luminosity distance with GWs, together with the estimate of the redshift obtained from the host galaxy, allowed us to estimate the Hubble constant with a totally new approach, independent from previous measurements [26]. Joint observations of GW and EM are thus a remarkable instrument for unveiling the physics of some of the most extreme phenomena in the Universe.

3. Artificial Intelligence via Multimodal Inputs

Multimodal machine learning (MML) is a multidisciplinary research area that addresses some of the main objectives of artificial intelligence (AI) by incorporating and creating models that can process and link information from multiple modal inputs, differing in the representation (text, images, etc.), dimensions (1D, 2D, etc.), and input data sources. By considering data from multiple modalities, it is possible to take into account the complementary information among them, which in turn leads to more robust predictions that reflect patterns not available when working with individual modalities.
The multimodal approach is already used in a wide variety of artificial intelligence problems, such as use for virtual assistance, image captioning, question answering, etc. One must consider the advantages in obtaining multiple input data characterizing the same event, as well as features extracted in different domain decompositions. In this way, we have the ability to capture otherwise unidentifiable details [13]. In Figure 1, we reported a schematic view of the general idea underlying the multimodal analysis. The input samples contain different kinds of signals and representations which can be encoded in a deep network analysis and concatenated, and the later stage can be used for classification/regression analysis.
Machine learning (ML) and deep learning (DL) techniques, which have changed the way data are processed in recent years, have already been implemented in the gravitational wave astronomy community [27], thanks to the computational resources available to us in recent years and the implementation of algorithms that allow the effective use of graphical processing units (GPUs). Recently, DL has even been successfully applied to multimodal machine learning problems, with the aim of learning useful joint representations in data fusion applications [28]. In [29], we introduced a first approach to MML, based on the merging of deep learning pipeline outputs on two different kinds of inputs: time series and images in the scheme, which is called late fusion multimodal [30].
The challenge will also be to apply these techniques to data from different instruments with different outputs that characterize a multi-messenger event. We want to show how MML could be used to process a 1D strain as well as 2D spectrograms from GW detectors with sparse light-curves collected by astronomical telescopes in order to infer astrophysical information from the common sources. At the same time, even the ability to caption GW data with an associated GRB event could help in quickly identifying source parameters.

4. From Multi-Messenger Observations to Multimodal Analysis

In the next years, second-generation GW interferometers (Advanced LIGO [31], Advanced Virgo [32], and KAGRA [33,34]) will take data with increased sensitivity, and third-generation GW detectors (such as the Einstein Telescope [35]) will become operative; furthermore, many new telescopes will start taking data (e.g., CTA [36] and LSST [37]). We therefore expect an increase in the data rates and in the data complexity. In order to maximize the scientific return of future multi-messenger observations, there is a need to develop new approaches to analyse large streams of EM, GW, and neutrino data, taking into account the differences in instrument sensitivities, spatial and temporal coverage, data formats, etc.; furthermore, novel tools are needed to combine information in an efficient way from multi-messenger observations, allowing us to infer the properties of the astrophysical sources and their environment. For instance, joint GW and EM observations can be used to put more stringent constraints on the tidal deformability parameter (Λ) and therefore on the EOS of NS (see, e.g., [22,38]; see also Section 2). To do this, we need to perform a GW parameter estimation with an accurate family of template waveforms, as well as a detailed comparison between the EM observations and the existing kilonova theoretical models; information obtained from the two messengers should then be combined in a consistent way to put constraints on the EOS. Such constraints, together with the estimates of the masses of the two NSs as obtained with GWs, could also help us to obtain insight into the outcome of the coalescence (that can be an NS or a BH) [24]. In the case of a coincident short GRB observation, the additional characterization of the X-ray afterglow emission will allow us to eventually probe the GRB-magnetar model, based on which magnetars are the GRB central engine; such a model has revealed itself very successful in reproducing the observed properties for the sub-class of short GRBs showing an X-ray plateau [39] and/or an extended emission [40,41], but a direct proof of this connection is still missing.
In this work, we propose a new paradigm for analysing the multi-messenger data we will collect with next generation instruments.
In Figure 2, we sketched an example of multimodal analysis for multi-messenger events. We can approach the information of the diverse messengers through a dedicated pipeline, representing them in the best format for feature extraction. We can analyse these representations using the best suited machine learning workflow to maximize the capability of prediction and, at a final stage, combine the output of the extracted features in a multimodal model. We are depicting here a future vision, where, in an open access environment, we can analyse shared data using shared software on cloud systems [42].

5. Application to Astrophysical Sources: The Case of Binaries of Compact Objects

We decided to test our idea on a set of simulated short GRB light curves and associated GW events, with a focus on NS-NS mergers. Specifically, we divided this task into three steps. (1) We simulated a sample of NS-NS merging systems populating the universe volume that can be explored with next-generation GW interferometers. Specifically, we assigned to each NS-NS system a luminosity distance that is randomly extracted from a uniform distribution in the range between 1 and 500 Mpc to cover a realistic range of the matched filter signal-to-noise ratio (SNR), varying from 4 to 20. We then assumed that both components of the binary systems have a mass distribution equal to a uniform distribution between 1 and 2.5 M; the two distributions were assumed to be uncorrelated. The inclination angle ( θ i ) of the systems was chosen taking into account that the GRB emission is collimated, and these sources are typically detected when they are on-axis, i.e., with the jet pointing towards the observer. According to the few estimates currently available, the jet opening angle θ j ranges from 3 ° to 8 ° (see, e.g., [43] and references therein); we therefore restricted the range of possible values of θ i to take this observation into account. (2) Following the approach presented in [44], we assumed that all the NS-NS mergers are associated with a short GRB, and we simulated the associated high-energy afterglow light curves that could be observed by the LAT instrument onboard the Fermi satellite [45], using GRB 090510 [46] as a template. (3) We simulated the GW signal associated with the NS-NS mergers using the TaylorF2 waveform model [47]. We simulated the noise data for a GW detector such as the Einstein Telescope [48], where we injected the NS-NS merger GW signals. The simulations were performed using the pyCBC library [49].
We propose an MML pipeline consisting of two convolutional neural networks (CNN) concatenated at the output in order to estimate the redshift of the GRB and GW sources. To convert luminosity distances in redshifts, we use the cosmological parameters in [50]. A 2D CNN takes the time-frequency image representation of the GW signal as input. The images are built from the detector strain time series containing the injected NS-NS GW signals. As a first step, the strain is whitened in the time domain by means of an auto-regressive (AR) model [51] to remove the stationary noise component. The 60 s long segments containing the chirp signals are then converted into a time-frequency representation based on the continuous wavelet transform, using Morlet wavelets (using ssqueezepy library (July, 2021)1) (see [52,53,54]). Finally, the images are compressed to dimensions in pixels of 128 × 256 . For GRB lightcurves, we used the data in the time domain. The length of the GRB simulated data was kept up to 1000 points, where most of the information was encoded2. Examples of the chosen representations are shown in Figure 3: on the left side the image for the GW signal, and on the right side, the time domain light curve data.
The CNN processing GW data consists of 5 layers with the following number of filters: 64, 32, 16, 16, and 32, and kernels of dimensions ( 3 , 3 ). After every convolutional layer, we applied maxpooling with kernels ( 2 , 2 ). The 1D CNN processing the light curves consists of 3 layers with the following number of filters: 80, 40, and 40, and kernels 5, 3, and 3. Additionally, in this case, after every convolutional layer, we applied maxpooling with kernel 2. After the last layers, both CNNs were flattened and concatenated as one. The concatenated output is fed to a fully connected layer which outputs the prediction for the source redshift. All layers have a ReLU activation function, with the exception of the final layer with a linear activation. The MML model was trained on the training data set using the mean squared error (MSE) loss function and Adam optimizer with a learning rate of α = 0.001 and a learning rate decay of 0.066667. The model is summarised in Figure 4. For training and evaluation, the dataset was divided with the following scheme: 70 % training set, 20 % test set, and 10 % validation set. Training was carried out with a batch size of 16 samples for 100 epochs. The total number of trainable parameters is 61,881. As a pre-processing step, we applied minmax scaling to the inputs. Shuffling is also applied.
In Figure 5, we report the mean squared error loss evolution over the training epochs for both the train and validation sets. It can be observed that the MSE converged well over the 100 epochs algorithm and that there is no overfitting.
In Figure 6, we show the comparison between the predicted redshifts and the true redshifts of the simulated sources on the test data set. It can be seen that for the lowest values of redshifts, the predictions are in good agreement with the real values, while for the highest redshifts, the scatter around the line y = x increases, and the predicted values are typically lower than the true ones. This is partially due to the fact that for the most distant sources, the GW signal has a low SNR value. Additionally, we used a limited dataset, and this could have affected our results. However, with this work, we only wanted to investigate the feasibility of multi-modal ML analysis; a more detailed study, including a larger simulated dataset and possibly real EM and GW data, will be presented elsewhere.
In Figure 7, we show the histogram of the the relative difference between the predicted redshifts and the true redshifts. It can be seen that the histogram has a peak at zero, meaning that for the majority of the simulated sources, the algorithm allowed for correctly estimating the redshift.

6. Outlook and Perspective

Our team’s innovative multimodal machine learning-based approach consists of analysing each multi-messenger event simultaneously in its entirety through the information we can gather from GWs or electromagnetic or neutrino emissions, which will be the inputs for a single MMML pipeline. This allows us to increase the information we can extract, thanks to the concomitant analysis of every piece of information. Just as it is easier for us to understand language when we associate sound to video, or to understand people’s moods by analysing text and images, we will apply multimodal analysis techniques to astrophysical transient events.
This is a new paradigm for analysing the data we will collect in future astro and particle experiments. We tested the idea on a data set we created simulating both GW and GRB emissions, just to set up a proof of concept. We built an MML pipeline taking as input the images obtained through a time frequency representation of whitened GW data for the Einstein Telescope detector and the simulated GRB light curves for a Fermi-like detector. We want to underline that the work here presented is just an example of a basic application of multimodal analysis for multi-messenger events, based only on GW strain data and gamma-ray GRB light curves, with a simple neural network architecture. The results, even if preliminary, are very encouraging, and we plan to continue with our approach on real data or including other input messengers. We will implement a more sophisticated and comprehensive analysis by adding data in different formats, at different levels (from raw data to high-level data), and also related to other messengers (e.g., neutrinos and photons at other wavelengths). For instance, we want to include in our framework raw data from Imaging Atmospheric Cherenkov Telescopes (IACTs) such as CTA: 2D images that represent the tracks left in the telescope’s camera by showers of particles that can be induced, for example, by very high-energy (VHE, E > 100 GeV) gamma rays emitted by GRBs. At the same time, we plan to optimize the analysis workflows by using feature extraction and engineering, data balance, and augmentation for less represented classes of events and by working with deeper networks. It is worth emphasizing that this approach we are proposing represents a new challenge for data scientists and astrophysicists, even in anticipation of the larger number of events we will be able to detect with third-generation GW detectors.

Author Contributions

Conceptualization, E.C.; Data curation, E.C., B.P., A.I. and F.M.; Formal analysis, E.C., B.P., A.I. and F.M.; Methodology, E.C., B.P., A.I. and F.M.; Writing—original draft, E.C., B.P., A.I. and F.M. All authors have read and agreed to the published version of the manuscript.

Funding

FM is supported by the Polish National Science Centre grants no. 2016/22/E/ST9/00037, 2017/26/M/ST9/00978 and 2020/37/N/ST9/02151, Polish National Agency for Academic Exchange grant PPN/IWA/2019/1/00157.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

This article/publication is based upon work from COST Action CA17137, supported by COST (European Cooperation in Science and Technology). We acknowledge support from the ESCAPE project with grant no. GA:824064.

Conflicts of Interest

The authors declare no conflict of interest.

Notes

1
2
We used a time resolution of 0.01 s, so 1000 points correspond to a time interval of 10 s.

References

  1. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [Green Version]
  2. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; Affeldt, C.; et al. Multi-messenger Observations of a Binary Neutron Star Merger. Astrophys. J. 2017, 848, L12. [Google Scholar] [CrossRef]
  3. Nordin, J.; Brinnel, V.; van Santen, J.; Bulla, M.; Feindt, U.; Franckowiak, A.; Fremling, C.; Gal-Yam, A.; Giomi, M.; Kowalski, M.; et al. Transient processing and analysis using AMPEL: Alert management, photometry, and evaluation of light curves. Astron. Astrophys. 2019, 631, A147. [Google Scholar] [CrossRef]
  4. Chatterjee, D.; Narayan, G.; Aleo, P.D.; Malanchev, K.; Muthukrishna, D. El-CID: A filter for Gravitational-wave Electromagnetic Counterpart Identification. arXiv 2021, arXiv:2108.04166. [Google Scholar]
  5. Stachie, C.; Coughlin, M.W.; Christensen, N.; Muthukrishna, D. Using machine learning for transient classification in searches for gravitational-wave counterparts. Mon. Not. R. Astron. Soc. 2020, 497, 1320–1331. [Google Scholar] [CrossRef]
  6. Dax, M.; Green, S.R.; Gair, J.; Macke, J.H.; Buonanno, A.; Schölkopf, B. Real-time gravitational-wave science with neural posterior estimation. arXiv 2021, arXiv:2106.12594. [Google Scholar]
  7. Krastev, P.G.; Gill, K.; Villar, V.A.; Berger, E. Detection and parameter estimation of gravitational waves from binary neutron-star mergers in real LIGO data using deep learning. Phys. Lett. B 2021, 815, 136161. [Google Scholar] [CrossRef]
  8. Singh, S.; Singh, A.; Prajapati, A.; Pathak, K.N. Deep learning for estimating parameters of gravitational waves. Mon. Not. R. Astron. Soc. 2021, 508, 1358–1370. [Google Scholar] [CrossRef]
  9. Breschi, M.; Gamba, R.; Bernuzzi, S. Bayesian inference of multimessenger astrophysical data: Methods and applications to gravitational waves. Phys. Rev. D 2021, 104, 042001. [Google Scholar] [CrossRef]
  10. Krastev, P.G. Real-time detection of gravitational waves from binary neutron stars using artificial neural networks. Phys. Lett. B 2020, 803, 135330. [Google Scholar] [CrossRef]
  11. Sadeh, I. Data-driven Detection of Multimessenger Transients. Astrophys. J. Lett. 2020, 894, L25. [Google Scholar] [CrossRef]
  12. Green, S.R.; Gair, J. Complete parameter inference for GW150914 using deep learning. arXiv 2020, arXiv:2008.03312. [Google Scholar]
  13. Baltrusaitis, T.; Ahuja, C.; Morency, L.P. Multimodal Machine Learning: A Survey and Taxonomy. IEEE Trans. Pattern Anal. Mach. Intell. 2019, 41, 423–443. [Google Scholar] [CrossRef] [Green Version]
  14. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. Gravitational Waves and Gamma-Rays from a Binary Neutron Star Merger: GW170817 and GRB 170817A. Astrophys. J. Lett. 2017, 848, L13. [Google Scholar] [CrossRef]
  15. Troja, E.; Piro, L.; van Eerten, H.; Wollaeger, R.T.; Im, M.; Fox, O.D.; Butler, N.R.; Cenko, S.B.; Sakamoto, T.; Fryer, C.L.; et al. The X-ray counterpart to the gravitational-wave event GW170817. Nature 2017, 551, 71–74. [Google Scholar] [CrossRef] [Green Version]
  16. Hallinan, G.; Corsi, A.; Mooley, K.P.; Hotokezaka, K.; Nakar, E.; Kasliwal, M.M.; Kaplan, D.L.; Frail, D.A.; Myers, S.T.; Murphy, T.; et al. A radio counterpart to a neutron star merger. Science 2017, 358, 1579–1583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Mooley, K.P.; Deller, A.T.; Gottlieb, O.; Nakar, E.; Hallinan, G.; Bourke, S.; Frail, D.A.; Horesh, A.; Corsi, A.; Hotokezaka, K. Superluminal motion of a relativistic jet in the neutron-star merger GW170817. Nature 2018, 561, 355–359. [Google Scholar] [CrossRef]
  18. Ghirlanda, G.; Salafia, O.S.; Paragi, Z.; Giroletti, M.; Yang, J.; Marcote, B.; Blanchard, J.; Agudo, I.; An, T.; Bernardini, M.G.; et al. Compact radio emission indicates a structured jet was produced by a binary neutron star merger. Science 2019, 363, 968–971. [Google Scholar] [CrossRef] [Green Version]
  19. Coulter, D.A.; Kilpatrick, C.D.; Siebert, M.R.; Foley, R.J.; Shappee, B.J.; Drout, M.R.; Simon, J.S.; Piro, A.L.; Rest, A.; One-Meter Two-Hemisphere (1M2H) Collaboration. LIGO/Virgo G298048: Potential optical counterpart discovered by Swope telescope. GRB Coord. Netw. 2017, 21529, 1. [Google Scholar]
  20. Pian, E.; D’Avanzo, P.; Benetti, S.; Branchesi, M.; Brocato, E.; Campana, S.; Cappellaro, E.; Covino, S.; D’Elia, V.; Fynbo, J.P.U.; et al. Spectroscopic identification of r-process nucleosynthesis in a double neutron-star merger. Nature 2017, 551, 67–70. [Google Scholar] [CrossRef]
  21. Smartt, S.J.; Chen, T.-W.; Jerkstrand, A.; Coughlin, M.; Kankare, E.; Sim, S.A.; Fraser, M.; Inserra, C.; Maguire, K.; Chambers, K.C.; et al. A kilonova as the electromagnetic counterpart to a gravitational-wave source. Nature 2017, 551, 75–79. [Google Scholar] [CrossRef] [PubMed]
  22. Radice, D.; Perego, A.; Zappa, F.; Bernuzzi, S. GW170817: Joint Constraint on the Neutron Star Equation of State from Multimessenger Observations. Astrophys. J. Lett. 2018, 852, L29. [Google Scholar] [CrossRef] [Green Version]
  23. Coughlin, M.W.; Dietrich, T.; Margalit, B.; Metzger, B.D. Multimessenger Bayesian parameter inference of a binary neutron star merger. Mon. Not. R. Astron. Soc. 2019, 489, L91–L96. [Google Scholar] [CrossRef]
  24. Dietrich, T.; Coughlin, M.W.; Pang, P.T.H.; Bulla, M.; Heinzel, J.; Issa, L.; Tews, I.; Antier, S. Multimessenger constraints on the neutron-star equation of state and the Hubble constant. Science 2020, 370, 1450–1453. [Google Scholar] [CrossRef]
  25. Patricelli, B.; Bernardini, M.G. Can we constrain the aftermath of binary neutron star mergers with short gamma-ray bursts? Mon. Not. R. Astron. Soc. 2020, 499, L96–L100. [Google Scholar] [CrossRef]
  26. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. A gravitational-wave standard siren measurement of the Hubble constant. Nature 2017, 551, 85–88. [Google Scholar] [CrossRef] [Green Version]
  27. Cuoco, E.; Powell, J.; Cavaglià, M.; Ackley, K.; Bejger, M.; Chatterjee, C.; Coughlin, M.; Coughlin, S.; Easter, P.; Essick, R. Enhancing Gravitational-Wave Science with Machine Learning. Mach. Learn. Sci. Technol. 2020, 2, 011002. [Google Scholar] [CrossRef]
  28. Ngiam, J.; Khosla, A.; Kim, M.; Nam, J.; Lee, H.; Ng, A.Y. Multimodal Deep Learning. In Proceedings of the 28th International Conference on International Conference on Machine Learning, Bellevue, WA, USA, 28 June–2 July 2011; Omnipress: Madison, WI, USA, 2011; pp. 689–696. [Google Scholar]
  29. Iess, A.; Cuoco, E.; Morawski, F.; Powell, J. Core-Collapse supernova gravitational-wave search and deep learning classification. Mach. Learn. Sci. Technol. 2020, 1, 025014. [Google Scholar] [CrossRef]
  30. Pandeya, Y.; Lee, J. Deep learning-based late fusion of multimodal information for emotion classification of music video. Multimed. Tools Appl. 2021, 80, 2887–2905. [Google Scholar] [CrossRef]
  31. LIGO Scientific Collaboration; Aasi, J.; Abbott, B.P.; Abbott, R.; Abbott, T.; Abernathy, M.R.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; et al. Advanced LIGO. Class. Quantum Gravity 2015, 32, 074001. [Google Scholar] [CrossRef]
  32. Acernese1, F.; Agathos, M.; Agatsuma, K.; Aisa, D.; Allemandou, N.; Allocca, A.; Amarni, J.; Astone, P.; Balestri, G.; Ballardin, G.; et al. Advanced Virgo: A second-generation interferometric gravitational wave detector. Class. Quantum Gravity 2015, 32, 024001. [Google Scholar] [CrossRef] [Green Version]
  33. Somiya, K. Detector configuration of KAGRA-the Japanese cryogenic gravitational-wave detector. Class. Quantum Gravity 2012, 29, 124007. [Google Scholar] [CrossRef] [Green Version]
  34. Aso, Y.; Michimura, Y.; Somiya, K.; Ando, M.; Miyakawa, O.; Sekiguchi, T.; Tatsumi, D.; Yamamoto, H. Interferometer design of the KAGRA gravitational wave detector. Phys. Rev. D 2013, 88, 043007. [Google Scholar] [CrossRef] [Green Version]
  35. Andersson, N. Whispers from the Edge of Physics. J. Astrophys. Astron. 2017, 38, 58. [Google Scholar] [CrossRef] [Green Version]
  36. Acharya, B.S.; Actis, M.; Aghajani, T.; Agnetta, G.; Aguilar, J.; Aharonian, F.; Ajello, M.; Akhperjanian, A.; Alcubierre, M.; Aleksić, J.; et al. Introducing the CTA concept. Astropart. Phys. 2013, 43, 3–18. [Google Scholar] [CrossRef] [Green Version]
  37. Ivezić, Ž.; Kahn, S.M.; Tyson, J.A.; Abel, B.; Acosta, E.; Allsman, R.; Alonso, D.; AlSayyad, Y.; Anderson, S.F.; Andrew, J.; et al. LSST: From Science Drivers to Reference Design and Anticipated Data Products. Astrophys. J. 2019, 873, 111. [Google Scholar] [CrossRef]
  38. Coughlin, M.W.; Dietrich, T.; Doctor, Z.; Kasen, D.; Coughlin, S.; Jerkstrand, A.; Leloudas, G.; McBrien, O.; Metzger, B.D.; O’Shaughnessy, R.; et al. Constraints on the neutron star equation of state from AT2017gfo using radiative transfer simulations. Mon. Not. R. Astron. Soc. 2018, 480, 3871–3878. [Google Scholar] [CrossRef] [Green Version]
  39. Rowlinson, A.; O’Brien, P.T.; Metzger, B.D.; Tanvir, N.R.; Levan, A.J. Signatures of magnetar central engines in short GRB light curves. Mon. Not. R. Astron. Soc. 2013, 430, 1061–1087. [Google Scholar] [CrossRef] [Green Version]
  40. Metzger, B.D.; Quataert, E.; Thompson, T.A. Short-duration gamma-ray bursts with extended emission from protomagnetar spin-down. Mon. Not. R. Astron. Soc. 2008, 385, 1455–1460. [Google Scholar] [CrossRef] [Green Version]
  41. Bucciantini, N.; Metzger, B.D.; Thompson, T.A.; Quataert, E. Short gamma-ray bursts with extended emission from magnetar birth: Jet formation and collimation. Mon. Not. R. Astron. Soc. 2012, 419, 1537–1545. [Google Scholar] [CrossRef] [Green Version]
  42. Allen, M.G.; Lamanna, G.; Espinal, X.; Graf, K.; van Haarlem, M.; Serjeant, S.; Bird, I.; Cuoco, E.; Wagh, J. ESCAPE–addressing Open Science challenges. arXiv 2020, arXiv:astro-ph.IM/2012.11534. [Google Scholar]
  43. Fong, W.; Berger, E.; Metzger, B.D.; Margutti, R.; Chornock, R.; Migliori, G.; Foley, R.J.; Zauderer, B.A.; Lunnan, R.; Laskar, T.; et al. Short GRB 130603B: Discovery of a Jet Break in the Optical and Radio Afterglows, and a Mysterious Late-time X-Ray Excess. Astrophys. J. 2014, 780, 118. [Google Scholar] [CrossRef] [Green Version]
  44. Patricelli, B.; Razzano, M.; Cella, G.; Fidecaro, F.; Pian, E.; Branchesi, M.; Stamerra, A. Prospects for joint observations of gravitational waves and gamma rays from merging neutron star binaries. J. Cosmol. Astropart. Phys. 2016, 2016, 056. [Google Scholar] [CrossRef] [Green Version]
  45. Atwood, W.B.; Abdo, A.A.; Ackermann, M.; Althouse, W.; Anderson, B.; Axelsson, M.; Baldini, L.; Ballet, J.; Band, D.L.; Barbiellini, G.; et al. The Large Area Telescope on the Fermi Gamma-Ray Space Telescope Mission. Astrophys. J. 2009, 697, 1071–1102. [Google Scholar] [CrossRef] [Green Version]
  46. Ackermann, M.; Asano, K.; Atwood, W.B.; Axelsson, M.; Baldini, L.; Ballet, J.; Barbiellini, G.; Baring, M.G.; Bastieri, D.; Bechtol, K.; et al. Fermi Observations of GRB 090510: A Short-Hard Gamma-ray Burst with an Additional, Hard Power-law Component from 10 keV TO GeV Energies. Astrophys. J. 2010, 716, 1178–1190. [Google Scholar] [CrossRef] [Green Version]
  47. Buonanno, A.; Iyer, B.R.; Ochsner, E.; Pan, Y.; Sathyaprakash, B.S. Comparison of post-Newtonian templates for compact binary inspiral signals in gravitational-wave detectors. Phys. Rev. D 2009, 80, 084043. [Google Scholar] [CrossRef] [Green Version]
  48. Hild, S.; Abernathy, M.; Acernese, F.; Amaro-Seoane, P.; Andersson, N.; Arun, K.; Barone, F.; Barr, B.; Barsuglia, M.; Beker, M.; et al. Sensitivity studies for third-generation gravitational wave observatories. Class. Quantum Gravity 2011, 28, 094013. [Google Scholar] [CrossRef]
  49. Nitz, A.; Harry, I.; Brown, D.; Biwer, C.M.; Willis, J.; Canton, T.D.; Capano, C.; Pekowsky, L.; Dent, T.; Williamson, A.R.; et al. gwastro/pycbc: PyCBC Release v1.16.10. 2020. Available online: https://zenodo.org/record/4063644#.YXDgmHkRXIU (accessed on 18 October 2021). [CrossRef]
  50. Ade, P.A.; Aghanim, N.; Arnaud, M.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Banday, A.J.; Barreiro, R.B.; Bartlett, J.G.; Bartolo, N.; et al. Planck 2015 results. XIII. Cosmological parameters. Astron. Astrophys. 2016, 594, A13. [Google Scholar] [CrossRef] [Green Version]
  51. Cuoco, E.; Calamai, G.; Fabbroni, L.; Losurdo, G.; Mazzoni, M.; Stanga, R.; Vetrano, F. On-line power spectra identification and whitening for the noise in interferometric gravitational wave detectors. Class. Quantum Gravity 2001, 18, 1727. [Google Scholar] [CrossRef] [Green Version]
  52. Thakur, G.; Wu, H.T. Synchrosqueezing-based Recovery of Instantaneous Frequency from Nonuniform Samples. SIAM J. Math. Anal. 2011, 43, 2078. [Google Scholar] [CrossRef] [Green Version]
  53. Thakur, G.; Brevdo, E.; Fučkar, N.S.; Wu, H.T. The Synchrosqueezing algorithm for time-varying spectral analysis: Robustness properties and new paleoclimate applications. Signal Process. 2013, 93, 1079. [Google Scholar] [CrossRef] [Green Version]
  54. Daubechies, I.; Wu, H.T. Synchrosqueezed Wavelet Transforms: A Tool for Empirical ModeDecomposition. Appl. Comput. Harmon. Anal. 2011, 30, 243. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Simple multimodal machine learning analysis workflow example.
Figure 1. Simple multimodal machine learning analysis workflow example.
Universe 07 00394 g001
Figure 2. An astrophysical phenomenon such as core-collapse supernova (CCSN) or NS-NS or BH-NS coalescence (multi-messenger events) can manifest itself through different signal types, such as gravitational waves, gamma-rays, X-rays, optical and radio emission, and neutrinos. The different modalities have their own representations in different domains. By using DL and ML models, we can use the extracted features to perform model prediction at an early stage. At a later stage, we can use the predicted features by combining them in the global MML model.
Figure 2. An astrophysical phenomenon such as core-collapse supernova (CCSN) or NS-NS or BH-NS coalescence (multi-messenger events) can manifest itself through different signal types, such as gravitational waves, gamma-rays, X-rays, optical and radio emission, and neutrinos. The different modalities have their own representations in different domains. By using DL and ML models, we can use the extracted features to perform model prediction at an early stage. At a later stage, we can use the predicted features by combining them in the global MML model.
Universe 07 00394 g002
Figure 3. Sample binary neutron star with component masses 2.4 M and 1.8 M used as input to the ML model, at a distance of 283 Mpc. The gravitational wave signal (left) and the GRB light curve (right) both contribute to the prediction. We turned off the x/y label (time-frequency information) to feed the CCN-2D network with images.
Figure 3. Sample binary neutron star with component masses 2.4 M and 1.8 M used as input to the ML model, at a distance of 283 Mpc. The gravitational wave signal (left) and the GRB light curve (right) both contribute to the prediction. We turned off the x/y label (time-frequency information) to feed the CCN-2D network with images.
Universe 07 00394 g003
Figure 4. Illustration of the multimodal machine learning model adopted in the analysis to compute redshifts of joint GW and GRB sources. The network reads two different types of data, images, and time series to tackle the regression problem.
Figure 4. Illustration of the multimodal machine learning model adopted in the analysis to compute redshifts of joint GW and GRB sources. The network reads two different types of data, images, and time series to tackle the regression problem.
Universe 07 00394 g004
Figure 5. Training (in blue) and validation (in orange) MSE loss over epochs.
Figure 5. Training (in blue) and validation (in orange) MSE loss over epochs.
Universe 07 00394 g005
Figure 6. Predicted redshifts versus the true redshifts of the simulated sources (blue points); the red solid line represents the expectations in the case that the predicted redshifts exactly match the true redshifts.
Figure 6. Predicted redshifts versus the true redshifts of the simulated sources (blue points); the red solid line represents the expectations in the case that the predicted redshifts exactly match the true redshifts.
Universe 07 00394 g006
Figure 7. Histogram of the relative error in the estimated distance with respect to true values.
Figure 7. Histogram of the relative error in the estimated distance with respect to true values.
Universe 07 00394 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cuoco, E.; Patricelli, B.; Iess, A.; Morawski, F. Multimodal Analysis of Gravitational Wave Signals and Gamma-Ray Bursts from Binary Neutron Star Mergers. Universe 2021, 7, 394. https://doi.org/10.3390/universe7110394

AMA Style

Cuoco E, Patricelli B, Iess A, Morawski F. Multimodal Analysis of Gravitational Wave Signals and Gamma-Ray Bursts from Binary Neutron Star Mergers. Universe. 2021; 7(11):394. https://doi.org/10.3390/universe7110394

Chicago/Turabian Style

Cuoco, Elena, Barbara Patricelli, Alberto Iess, and Filip Morawski. 2021. "Multimodal Analysis of Gravitational Wave Signals and Gamma-Ray Bursts from Binary Neutron Star Mergers" Universe 7, no. 11: 394. https://doi.org/10.3390/universe7110394

APA Style

Cuoco, E., Patricelli, B., Iess, A., & Morawski, F. (2021). Multimodal Analysis of Gravitational Wave Signals and Gamma-Ray Bursts from Binary Neutron Star Mergers. Universe, 7(11), 394. https://doi.org/10.3390/universe7110394

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop