Showing 1–50 of 374 results for author: Font, A

High-speed reconstruction of long-duration gravitational waves from extreme mass ratio inspirals using sparse dictionary learning

Authors: Charles Badger, José A. Font, Mairi Sakellariadou, Alejandro Torres-Forné

Abstract: Measuring accurate long-duration gravitational waves from extreme mass ratio inspirals (EMRIs) could provide scientifically fruitful knowledge of massive black hole populations and robust tests for general relatively during the LISA mission. However, the immense computational requirements surrounding EMRI data processing and analysis makes their detection and analysis challenging. We further devel… ▽ More Measuring accurate long-duration gravitational waves from extreme mass ratio inspirals (EMRIs) could provide scientifically fruitful knowledge of massive black hole populations and robust tests for general relatively during the LISA mission. However, the immense computational requirements surrounding EMRI data processing and analysis makes their detection and analysis challenging. We further develop and explore a sparse dictionary learning (SDL) algorithm to expeditiously reconstruct EMRI gravitational waveforms lasting as long as 1 year. A suite of year-long EMRI systems are studied to understand the detection and accurate waveform retrieval prospects of the method. We show that full-year EMRIs can be reconstructed within 2 minutes, some with a false alarm rate less than 0.001/yr and with 1.16 day time windows with mismatch as low as 0.06. This provides an encouraging prospect to use the SDL method for long-duration GW searches like that for EMRIs in this study. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Rapid detection of gravitational waves from binary black hole mergers using sparse dictionary learning

Authors: Charles Badger, Rahul Srinivasan, Alejandro Torres-Forné, Marie Anne Bizouard, José A. Font, Mairi Sakellariadou, Astrid Lamberts

Abstract: Current gravitational wave (GW) detection pipelines for compact binary coalescence based on matched-filtering have reported over 90 confident detections during the first three observing runs of the LIGO-Virgo-KAGRA (LVK) detector network. Decreasing the latency of detection, in particular for future detectors anticipated to have high detection rates, remains an ongoing effort. In this paper, we de… ▽ More Current gravitational wave (GW) detection pipelines for compact binary coalescence based on matched-filtering have reported over 90 confident detections during the first three observing runs of the LIGO-Virgo-KAGRA (LVK) detector network. Decreasing the latency of detection, in particular for future detectors anticipated to have high detection rates, remains an ongoing effort. In this paper, we develop and test a sparse dictionary learning (SDL) algorithm for the rapid detection of GWs. We evaluate the algorithms biases and estimate its GW detection rate for an astrophysical population of binary black holes. The SDL algorithm is assessed using both, simulated data injected into the proposed A+ detector sensitivity and real data containing confident detections from the third LVK observing run. We find that our SDL algorithm can reconstruct a single binary black hole signal in less than 1 s. This suggests that SDL could be regarded as a promising approach for rapid, efficient GW detection in future observing runs of ground-based detectors. △ Less

Submitted 29 May, 2024; v1 submitted 27 May, 2024; originally announced May 2024.

Applying machine learning to Galactic Archaeology: how well can we recover the origin of stars in Milky Way-like galaxies?

Authors: Andrea Sante, Andreea S. Font, Sandra Ortega-Martorell, Ivan Olier, Ian G. McCarthy

Abstract: We present several machine learning (ML) models developed to efficiently separate stars formed in-situ in Milky Way-type galaxies from those that were formed externally and later accreted. These models, which include examples from artificial neural networks, decision trees and dimensionality reduction techniques, are trained on a sample of disc-like, Milky Way-mass galaxies drawn from the ARTEMIS… ▽ More We present several machine learning (ML) models developed to efficiently separate stars formed in-situ in Milky Way-type galaxies from those that were formed externally and later accreted. These models, which include examples from artificial neural networks, decision trees and dimensionality reduction techniques, are trained on a sample of disc-like, Milky Way-mass galaxies drawn from the ARTEMIS cosmological hydrodynamical zoom-in simulations. We find that the input parameters which provide an optimal performance for these models consist of a combination of stellar positions, kinematics, chemical abundances ([Fe/H] and [$α$/Fe]) and photometric properties. Models from all categories perform similarly well, with area under the precision-recall curve (PR-AUC) scores of $\simeq 0.6$. Beyond a galactocentric radius of $5$~kpc, models retrieve $>90\%$ of accreted stars, with a sample purity close to $60\%$, however the purity can be increased by adjusting the classification threshold. For one model, we also include host galaxy-specific properties in the training, to account for the variability of accretion histories of the hosts, however this does not lead to an improvement in performance. The ML models can identify accreted stars even in regions heavily dominated by the in-situ component (e.g., in the disc), and perform well on an unseen suite of simulations (the Auriga simulations). The general applicability bodes well for application of such methods on observational data to identify accreted substructures in the Milky Way without the need to resort to selection cuts for minimising the contamination from in-situ stars. △ Less

Submitted 18 June, 2024; v1 submitted 30 April, 2024; originally announced May 2024.

Quasinormal Modes in Modified Gravity using Physics-Informed Neural Networks

Authors: Raimon Luna, Daniela D. Doneva, José A. Font, Jr-Hua Lien, Stoytcho S. Yazadjiev

Abstract: In this paper, we apply a novel approach based on physics-informed neural networks to the computation of quasinormal modes of black hole solutions in modified gravity. In particular, we focus on the case of Einstein-scalar-Gauss-Bonnet theory, with several choices of the coupling function between the scalar field and the Gauss-Bonnet invariant. This type of calculation introduces a number of chall… ▽ More In this paper, we apply a novel approach based on physics-informed neural networks to the computation of quasinormal modes of black hole solutions in modified gravity. In particular, we focus on the case of Einstein-scalar-Gauss-Bonnet theory, with several choices of the coupling function between the scalar field and the Gauss-Bonnet invariant. This type of calculation introduces a number of challenges with respect to the case of General Relativity, mainly due to the extra complexity of the perturbation equations and to the fact that the background solution is known only numerically. The solution of these perturbation equations typically requires sophisticated numerical techniques that are not easy to develop in computational codes. We show that physics-informed neural networks have an accuracy which is comparable to traditional numerical methods in the case of numerical backgrounds, while being very simple to implement. Additionally, the use of GPU parallelization is straightforward thanks to the use of standard machine learning environments. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 13 pages, 9 figures

Numerical relativity surrogate models for exotic compact objects: the case of head-on mergers of equal-mass Proca stars

Authors: Raimon Luna, Miquel Llorens-Monteagudo, Ana Lorenzo-Medina, Juan Calderón Bustillo, Nicolas Sanchis-Gual, Alejandro Torres-Forné, José A. Font, Carlos A. R. Herdeiro, Eugen Radu

Abstract: We present several high-accuracy surrogate models for gravitational-wave signals from equal-mass head-on mergers of Proca stars, computed through the Newman-Penrose scalar $ψ_4$. We also discuss the current state of the model extensions to mergers of Proca stars with different masses, and the particular challenges that these present. The models are divided in two main categories: two-stage and mon… ▽ More We present several high-accuracy surrogate models for gravitational-wave signals from equal-mass head-on mergers of Proca stars, computed through the Newman-Penrose scalar $ψ_4$. We also discuss the current state of the model extensions to mergers of Proca stars with different masses, and the particular challenges that these present. The models are divided in two main categories: two-stage and monolithic. In the two-stage models, a dimensional reduction algorithm is applied to embed the data in a reduced feature space, which is then interpolated in terms of the physical parameters. For the monolithic models, a single neural network is trained to predict the waveform from the input physical parameter. Our model displays mismatches below $10^{-3}$ with respect to the original numerical waveforms. Finally, we demonstrate the usage of our model in full Bayesian parameter inference through the accurate recovery of numerical relativity signals injected in zero-noise, together with the analysis of GW190521. For the latter, we observe excellent agreement with existing results that make use of full numerical relativity. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 16 pages, 14 figures

ARTEMIS emulator: exploring the effect of cosmology and galaxy formation physics on Milky Way-mass haloes and their satellites

Authors: Shaun T. Brown, Azadeh Fattahi, Ian G. McCarthy, Andreea S. Font, Kyle A. Oman, Alexander H. Riley

Abstract: We present the new ARTEMIS emulator suite of high resolution (baryon mass of $2.23 \times 10^{4}$ $h^{-1}$M$_{\odot}$) zoom-in simulations of Milky Way mass systems. Here, three haloes from the original ARTEMIS sample have been rerun multiple times, systematically varying parameters for the stellar feedback model, the density threshold for star formation, the reionisation redshift and the assumed… ▽ More We present the new ARTEMIS emulator suite of high resolution (baryon mass of $2.23 \times 10^{4}$ $h^{-1}$M$_{\odot}$) zoom-in simulations of Milky Way mass systems. Here, three haloes from the original ARTEMIS sample have been rerun multiple times, systematically varying parameters for the stellar feedback model, the density threshold for star formation, the reionisation redshift and the assumed warm dark matter (WDM) particle mass (assuming a thermal relic). From these simulations emulators are trained for a wide range of statistics that allow for fast predictions at combinations of parameters not originally sampled, running in $\sim 1$ms (a factor of $\sim 10^{11}$ faster than the simulations). In this paper we explore the dependence of the central haloes' stellar mass on the varied parameters, finding the stellar feedback parameters to be the most important. When constraining the parameters to match the present-day stellar mass halo mass relation inferred from abundance matching we find that there is a strong degeneracy in the stellar feedback parameters, corresponding to a freedom in formation time of the stellar component for a fixed halo assembly history. We additionally explore the dependence of the satellite stellar mass function, where it is found that variations in stellar feedback, the reionisation redshift and the WDM mass all have a significant effect. The presented emulators are a powerful tool which allows for fundamentally new ways of analysing and interpreting cosmological hydrodynamic simulations. Crucially, allowing their free (subgrid) parameters to be varied and marginalised, leading to more robust constraints and predictions. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 18 pages, 9 figures, submitted to MNRAS

Deep-Learning Classification and Parameter Inference of Rotational Core-Collapse Supernovae

Authors: Solange Nunes, Gabriel Escrig, Osvaldo G. Freitas, José A. Font, Tiago Fernandes, Antonio Onofre, Alejandro Torres-Forné

Abstract: We test Deep-Learning (DL) techniques for the analysis of rotational core-collapse supernovae (CCSN) gravitational-wave (GW) signals by performing classification and parameter inference of the GW strain amplitude ($D \cdot Δh $) and the maximum (peak) frequency $(f_\textrm{peak})$, attained at core bounce. Our datasets are built from a catalog of numerically generated CCSN waveforms assembled by R… ▽ More We test Deep-Learning (DL) techniques for the analysis of rotational core-collapse supernovae (CCSN) gravitational-wave (GW) signals by performing classification and parameter inference of the GW strain amplitude ($D \cdot Δh $) and the maximum (peak) frequency $(f_\textrm{peak})$, attained at core bounce. Our datasets are built from a catalog of numerically generated CCSN waveforms assembled by Richers 2017. Those waveforms are injected into noise from the Advanced LIGO and Advanced Virgo detectors corresponding to the O2 and O3a observing runs. For a signal-to-noise ratio (SNR) above 5, our classification network using time series detects Galactic CCSN GW signals buried in detector noise with a false positive rate (FPR) of 0.10% and a 98% accuracy, being able to detect all signals with SNR>10. The inference of $f_\textrm{peak}$ is more accurate than for $D \cdot Δh $, particularly for our datasets with the shortest time window (0.25 s) and for a minimum SNR=15. From the calibration plots of predicted versus true values of the two parameters, the standard deviation ($σ$) and the slope deviation with respect to the ideal value are computed. We find $σ_{D \cdot Δh}$ = 52.6cm and $σ_{f_\textrm{peak}}$ = 18.3Hz, with respective slope deviations of 11.6% and 8.3%. Our best model is also tested on waveforms from a recent CCSN catalog built by Mitra 2023, different from the one used for the training. For these new waveforms the true values of the two parameters are mostly within the 1$σ$ band around the network's predicted values. Our results show that DL techniques hold promise to infer physical parameters of Galactic rotational CCSN even in the presence of real (non-Gaussian) noise conditions from current GW detectors. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Gravitational-wave imprints of non-convex dynamics in binary neutron star mergers

Authors: Giuseppe Rivieccio, Davide Guerra, Milton Ruiz, José A. Font

Abstract: Explaining gravitational-wave (GW) observations of binary neutron star (BNS) mergers requires an understanding of matter beyond nuclear saturation density. Our current knowledge of the properties of high-density matter relies on electromagnetic and GW observations, nuclear physics experiments, and general relativistic numerical simulations. In this paper we perform numerical-relativity simulations… ▽ More Explaining gravitational-wave (GW) observations of binary neutron star (BNS) mergers requires an understanding of matter beyond nuclear saturation density. Our current knowledge of the properties of high-density matter relies on electromagnetic and GW observations, nuclear physics experiments, and general relativistic numerical simulations. In this paper we perform numerical-relativity simulations of BNS mergers subject to non-convex dynamics, allowing for the appearance of expansive shock waves and compressive rarefactions. Using a phenomenological non-convex equation of state we identify observable imprints on the GW spectra of the remnant. In particular, we find that non-convexity induces a significant shift in the quasi-universal relation between the peak frequency of the dominant mode and the tidal deformability (of order $Δf_{\rm peak}\gtrsim 380\,\rm Hz$) with respect to that of binaries with convex (regular) dynamics. Similar shifts have been reported in the literature, attributed however to first-order phase transitions from nuclear/hadronic matter to deconfined quark matter. We argue that the ultimate origin of the frequency shifts is to be found in the presence of anomalous, non-convex dynamics in the binary remnant. △ Less

Submitted 13 March, 2024; v1 submitted 12 January, 2024; originally announced January 2024.

Comments: 13 pages, 8 figures. Minor changes. Journal version

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

Differences in the properties of disrupted and surviving satellites of Milky-Way-mass galaxies in relation to their host accretion histories

Authors: Salvador E. Grimozzi, Andreea S. Font, María Emilia De Rossi

Abstract: From the chemo-dynamical properties of tidal debris in the Milky Way, it has been inferred that the dwarf satellites that have been disrupted had different chemical abundances from their present-day counterparts of similar mass that survive today, specifically, they had lower [Fe/H] and higher [Mg/Fe]. Here we use the ARTEMIS simulations to study the relation between the chemical abundances of dis… ▽ More From the chemo-dynamical properties of tidal debris in the Milky Way, it has been inferred that the dwarf satellites that have been disrupted had different chemical abundances from their present-day counterparts of similar mass that survive today, specifically, they had lower [Fe/H] and higher [Mg/Fe]. Here we use the ARTEMIS simulations to study the relation between the chemical abundances of disrupted progenitors of MW-mass galaxies and their stellar mass, and the evolution of the stellar mass - metallicity relations (MZR) of this population with redshift. We find that these relations have significant scatter, which correlates with the accretion redshifts ($z_{\rm acc}$) of satellites, and with their cold gas fractions. We investigate the MZRs of dwarf populations accreted at different redshifts and find that they have similar slopes, and also similar with the slope of the MZR of the surviving population ($\approx 0.32$). However, the entire population of disrupted dwarfs displays a steeper MZR, with a slope of $\approx 0.48$, which can be explained by the changes in the mass spectrum of accreted dwarf galaxies with redshift. We find strong relations between the (mass-weighted) $\langle z_{\rm acc} \rangle$ of the disrupted populations and their global chemical abundances ($\langle$[Fe/H]$\rangle$ and $\langle$[Mg/Fe]$\rangle$), which suggests that chemical diagnostics of disrupted dwarfs can be used to infer the types of merger histories of their hosts. For the case of the MW, our simulations predict that the bulk of the disrupted population was accreted at $\langle z_{\rm acc} \rangle \approx 2$, in agreement with other findings. We also find that disrupted satellites form and evolve in denser environments, closer to their hosts, than their present-day counterparts. △ Less

Submitted 24 March, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 22 pages, 16 figures. Accepted to MNRAS

Identifying thermal effects in neutron star merger remnants with model-agnostic waveform reconstructions and third-generation detectors

Authors: Miquel Miravet-Tenés, Davide Guerra, Milton Ruiz, Pablo Cerdá-Durán, José A. Font

Abstract: We explore the prospects for identifying differences in simulated gravitational-wave signals of binary neutron star (BNS) mergers associated with the way thermal effects are incorporated in the numerical-relativity modelling. We consider a hybrid approach in which the equation of state (EoS) comprises a cold, zero temperature, piecewise-polytropic part and a thermal part described by an ideal gas,… ▽ More We explore the prospects for identifying differences in simulated gravitational-wave signals of binary neutron star (BNS) mergers associated with the way thermal effects are incorporated in the numerical-relativity modelling. We consider a hybrid approach in which the equation of state (EoS) comprises a cold, zero temperature, piecewise-polytropic part and a thermal part described by an ideal gas, and a tabulated approach based on self-consistent, microphysical, finite-temperature EoS. We use time-domain waveforms corresponding to BNS merger simulations with four different EoS. Those are injected into Gaussian noise given by the sensitivity of the third-generation detector Einstein Telescope and reconstructed using BayesWave, a Bayesian data-analysis algorithm that recovers the signals through a model-agnostic approach. The two representations of thermal effects result in frequency shifts of the dominant peaks in the spectra of the post-merger signals, for both the quadrupole fundamental mode and the late-time inertial modes. For some of the EoS investigated those differences are large enough to be told apart, especially in the early post-merger phase when the signal amplitude is the loudest. △ Less

Submitted 18 January, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: 15 pages, 8 figures, 3 tables

Determining the core-collapse supernova explosion mechanism with current and future gravitational-wave observatories

Authors: Jade Powell, Alberto Iess, Miquel Llorens-Monteagudo, Martin Obergaulinger, Bernhard Müller, Alejandro Torres-Forné, Elena Cuoco, José A. Font

Abstract: Gravitational waves are emitted from deep within a core-collapse supernova, which may enable us to determine the mechanism of the explosion from a gravitational-wave detection. Previous studies suggested that it is possible to determine if the explosion mechanism is neutrino-driven or magneto-rotationally powered from the gravitational-wave signal. However, long duration magneto-rotational wavefor… ▽ More Gravitational waves are emitted from deep within a core-collapse supernova, which may enable us to determine the mechanism of the explosion from a gravitational-wave detection. Previous studies suggested that it is possible to determine if the explosion mechanism is neutrino-driven or magneto-rotationally powered from the gravitational-wave signal. However, long duration magneto-rotational waveforms, that cover the full explosion phase, were not available during the time of previous studies, and explosions were just assumed to be magneto-rotationally driven if the model was rapidly rotating. Therefore, we perform an updated study using new 3D long-duration magneto-rotational core-collapse supernova waveforms that cover the full explosion phase, injected into noise for the Advanced LIGO, Einstein Telescope and NEMO gravitational-wave detectors. We also include a category for failed explosions in our signal classification results. We then determine the explosion mechanism of the signals using three different methods: Bayesian model selection, dictionary learning, and convolutional neural networks. The three different methods are able to distinguish between neutrino-driven explosions and magneto-rotational explosions, even if the neutrino-driven explosion model is rapidly rotating. However they can only distinguish between the non-exploding and neutrino-driven explosions for signals with a high signal to noise ratio. △ Less

Submitted 28 February, 2024; v1 submitted 29 November, 2023; originally announced November 2023.

Comments: Accepted in Physical Review D

Self-consistent treatment of thermal effects in neutron-star post-mergers: observational implications for third-generation gravitational-wave detectors

Authors: Verónica Villa-Ortega, Ana Lorenzo-Medina, Juan Calderón Bustillo, Milton Ruiz, Davide Guerra, Pablo Cerdá-Duran, José A. Font

Abstract: We assess the impact of accurate, self-consistent modelling of thermal effects in neutron-star merger remnants in the context of third-generation gravitational-wave detectors. This is done through the usage, in Bayesian model selection experiments, of numerical-relativity simulations of binary neutron star (BNS) mergers modelled through: a) nuclear, finite-temperature (or ``tabulated'') equations… ▽ More We assess the impact of accurate, self-consistent modelling of thermal effects in neutron-star merger remnants in the context of third-generation gravitational-wave detectors. This is done through the usage, in Bayesian model selection experiments, of numerical-relativity simulations of binary neutron star (BNS) mergers modelled through: a) nuclear, finite-temperature (or ``tabulated'') equations of state (EoSs), and b) their simplifed piecewise (or ``hybrid'') representation. These cover four different EoSs, namely SLy4, DD2, HShen and LS220. Our analyses make direct use of the Newman-Penrose scalar $ψ_4$ outputted by numerical simulations. Considering a detector network formed by three Cosmic Explorers, we show that differences in the gravitational-wave emission predicted by the two models are detectable with a natural logarithmic Bayes Factor $\log{\cal{B}}\geq 5$ at average distances of $d_L \simeq 50$Mpc, reaching $d_L \simeq 100$Mpc for source inclinations $ι\leq 0.8$, regardless of the EoS. This impact is most pronounced for the HShen EoS. For low inclinations, only the DD2 EoS prevents the detectability of such modelling differences at $d_L \simeq 150$Mpc. Our results suggest that the usage a self-consistent treatment of thermal effects is crucial for third-generation gravitational wave detectors. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: 9 pages, 3 Figures

Numerical evolutions of boson stars in Palatini $f(\mathcal{R})$ gravity

Authors: Andreu Masó-Ferrando, Nicolas Sanchis-Gual, José A. Font, Gonzalo J. Olmo

Abstract: We investigate the time evolution of spherically symmetric boson stars in Palatini $f(\mathcal{R})$ gravity through Numerical Relativity computations. Employing a novel approach that establishes a correspondence between modified gravity with scalar matter and General Relativity with modified scalar matter, we are able to use the techniques of Numerical Relativity to simulate these systems. Specifi… ▽ More We investigate the time evolution of spherically symmetric boson stars in Palatini $f(\mathcal{R})$ gravity through Numerical Relativity computations. Employing a novel approach that establishes a correspondence between modified gravity with scalar matter and General Relativity with modified scalar matter, we are able to use the techniques of Numerical Relativity to simulate these systems. Specifically, we focus on the quadratic theory $f(\mathcal{R})=\mathcal{R}+ξ\mathcal{R}^2$ and compare the obtained solutions with those in General Relativity, exploring both positive and negative values of the coupling parameter $ξ$. Our findings reveal that boson stars in Palatini $f(\mathcal{R})$ gravity exhibit both stable and unstable evolutions. The latter give rise to three distinct scenarios: migration towards a stable configuration, complete dispersion, and gravitational collapse leading to the formation of a baby universe structure. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Taking the Milky Way for a spin: disc formation in the ARTEMIS simulations

Authors: Adam M. Dillamore, Vasily Belokurov, Andrey Kravtsov, Andreea S. Font

Abstract: We investigate the formation (spin-up) of galactic discs in the ARTEMIS simulations of Milky Way-mass galaxies. In almost all galaxies discs spin up at higher [Fe/H] than the Milky Way (MW). Those that contain an analogue of the Gaia Sausage-Enceladus (GSE) spin up at a lower average metallicity than those without. We identify six galaxies with spin-up metallicity similar to that of the MW, which… ▽ More We investigate the formation (spin-up) of galactic discs in the ARTEMIS simulations of Milky Way-mass galaxies. In almost all galaxies discs spin up at higher [Fe/H] than the Milky Way (MW). Those that contain an analogue of the Gaia Sausage-Enceladus (GSE) spin up at a lower average metallicity than those without. We identify six galaxies with spin-up metallicity similar to that of the MW, which form their discs $\sim 8-11$ Gyr ago. Five of these experience a merger similar to the GSE. The spin-up times correlate with the halo masses at early times: galaxies with early spin-up have larger virial masses at a lookback time $t_L=12$ Gyr. The fraction of stars accreted from outside the host galaxy is smaller in galaxies with earlier spin-ups. Accreted fractions small enough to be comparable to the MW are only found in galaxies with the earliest disc formation and large initial virial masses ($M_\mathrm{200c} \approx2\times10^{11}M_\odot$ at $t_L=12$ Gyr). We find that discs form when the halo's virial mass reaches a threshold of $M_\mathrm{200c}\approx(6\pm3)\times10^{11}M_\odot$, independent of the spin-up time. However, the failure to form a disc in other galaxies appears to be instead related to mergers at early times. We also find that discs form when the central potential is not particularly steep. Our results indicate that the MW assembled its mass and formed its disc earlier than the average galaxy of a similar mass. △ Less

Submitted 18 December, 2023; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 9 pages, 7 figures, published in MNRAS

The FLAMINGO project: revisiting the $S_8$ tension and the role of baryonic physics

Authors: Ian G. McCarthy, Jaime Salcido, Joop Schaye, Juliana Kwan, Willem Elbers, Roi Kugel, Matthieu Schaller, John C. Helly, Joey Braspenning, Carlos S. Frenk, Marcel P. van Daalen, Bert Vandenbroucke, Jonah T. Conley, Andreea S. Font, Amol Upadhye

Abstract: A number of recent studies have found evidence for a tension between observations of large-scale structure (LSS) and the predictions of the standard model of cosmology with the cosmological parameters fit to the cosmic microwave background (CMB). The origin of this '$S_8$ tension' remains unclear, but possibilities include new physics beyond the standard model, unaccounted for systematic errors in… ▽ More A number of recent studies have found evidence for a tension between observations of large-scale structure (LSS) and the predictions of the standard model of cosmology with the cosmological parameters fit to the cosmic microwave background (CMB). The origin of this '$S_8$ tension' remains unclear, but possibilities include new physics beyond the standard model, unaccounted for systematic errors in the observational measurements and/or uncertainties in the role that baryons play. Here we carefully examine the latter possibility using the new FLAMINGO suite of large-volume cosmological hydrodynamical simulations. We project the simulations onto observable harmonic space and compare with observational measurements of the power and cross-power spectra of cosmic shear, CMB lensing, and the thermal Sunyaev-Zel'dovich (tSZ) effect. We explore the dependence of the predictions on box size and resolution, cosmological parameters including the neutrino mass, and the efficiency and nature of baryonic 'feedback'. Despite the wide range of astrophysical behaviours simulated, we find that baryonic effects are not sufficiently large to remove the $S_8$ tension. Consistent with recent studies, we find the CMB lensing power spectrum is in excellent agreement with the standard model, whilst the cosmic shear power spectrum, tSZ effect power spectrum, and the cross-spectra between shear, CMB lensing, and the tSZ effect are all in varying degrees of tension with the CMB-specified standard model. These results suggest that some mechanism is required to slow the growth of fluctuations at late times and/or on non-linear scales, but that it is unlikely that baryon physics is driving this modification. △ Less

Submitted 9 October, 2023; v1 submitted 14 September, 2023; originally announced September 2023.

Comments: 26 pages, 12 figures, MNRAS, accepted with minor revisions

Revisiting the $I$-Love-$Q$ relations for superfluid neutron stars

Authors: Eneko Aranguren, José A. Font, Nicolas Sanchis-Gual, Raül Vera

Abstract: We study the tidal problem and the resulting $I$-Love-$Q$ approximate universal relations for rotating superfluid neutron stars in the Hartle-Thorne formalism. Superfluid stars are described in this work by means of a two-fluid model consisting of superfluid neutrons and all other charged constituents. We employ a stationary and axisymmetric perturbation scheme to second order around a static and… ▽ More We study the tidal problem and the resulting $I$-Love-$Q$ approximate universal relations for rotating superfluid neutron stars in the Hartle-Thorne formalism. Superfluid stars are described in this work by means of a two-fluid model consisting of superfluid neutrons and all other charged constituents. We employ a stationary and axisymmetric perturbation scheme to second order around a static and spherically symmetric background. Recently, we used this scheme to study isolated rotating superfluid stars. In this paper it is applied to analyze the axially symmetric sector of the tidal problem in a binary system. We show that a consistent use of perturbative matching theory amends the original two-fluid formalism for the tidal problem to account for the possible non-zero value of the energy density at the boundary of the star. This is exemplified by building numerically different stellar models spanning three equations of state. Significant departures from universality are found when the correct matching relations are not taken into account. We also present an augmented set of universal relations for superfluid neutron stars which includes the contribution to the total mass of the star at second order, $δM$. Therefore, our results complete the set of universal relations for rotating superfluid stars, generalizing our previous findings in the perfect fluid case. △ Less

Submitted 22 November, 2023; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: 18 pages, 3 figures. Reference added, typos and minor corrections, and some discussion expanded, to agree with the published version

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Assessment of a new sub-grid model for magnetohydrodynamical turbulence. II. Kelvin-Helmholtz instability

Authors: Miquel Miravet-Tenés, Pablo Cerdá-Durán, Martin Obergaulinger, José A. Font

Abstract: The modelling of astrophysical systems such as binary neutron star mergers or the formation of magnetars from the collapse of massive stars involves the numerical evolution of magnetised fluids at extremely large Reynolds numbers. This is a major challenge for (unresolved) direct numerical simulations which may struggle to resolve highly dynamical features as, e.g. turbulence, magnetic field ampli… ▽ More The modelling of astrophysical systems such as binary neutron star mergers or the formation of magnetars from the collapse of massive stars involves the numerical evolution of magnetised fluids at extremely large Reynolds numbers. This is a major challenge for (unresolved) direct numerical simulations which may struggle to resolve highly dynamical features as, e.g. turbulence, magnetic field amplification, or the transport of angular momentum. Sub-grid models offer a means to overcome those difficulties. In a recent paper we presented MInIT, an MHD-instability-induced-turbulence mean-field, sub-grid model based on the modelling of the turbulent (Maxwell, Reynolds, and Faraday) stress tensors. While in our previous work MInIT was assessed within the framework of the magnetorotational instability, in this paper we further evaluate the model in the context of the Kelvin-Helmholtz instability (KHI). The main difference with other sub-grid models (as e.g. the alpha-viscosity model or the gradient model) is that in MInIT we track independently the turbulent energy density at sub-grid scales, which is used, via a simple closure relation, to compute the different turbulent stresses relevant for the dynamics. The free coefficients of the model are calibrated using well resolved box simulations of magnetic turbulence generated by the KHI. We test the model against these simulations and show that it yields order-of-magnitude accurate predictions for the evolution of the turbulent Reynolds and Maxwell stresses. △ Less

Submitted 19 October, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: 12 pages, 13 figures Accepted for publication in MNRAS

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Bayesian inference from gravitational waves in fast-rotating, core-collapse supernovae

Authors: Carlos Pastor-Marcos, Pablo Cerdá-Durán, Daniel Walker, Alejandro Torres-Forné, Ernazar Abdikamalov, Sherwood Richers, José Antonio Font

Abstract: Core-collapse supernovae (CCSNe) are prime candidates for gravitational-wave detectors. The analysis of their complex waveforms can potentially provide information on the physical processes operating during the collapse of the iron cores of massive stars. In this work we analyze the early-bounce rapidly rotating CCSN signals reported in the waveform catalog of Richers et al 2017, which comprises o… ▽ More Core-collapse supernovae (CCSNe) are prime candidates for gravitational-wave detectors. The analysis of their complex waveforms can potentially provide information on the physical processes operating during the collapse of the iron cores of massive stars. In this work we analyze the early-bounce rapidly rotating CCSN signals reported in the waveform catalog of Richers et al 2017, which comprises over 1800 axisymmetric simulations extending up to about 10~ms of post-bounce evolution. It was previously established that for a large range of progenitors, the amplitude of the bounce signal, $Δh$, is proportional to the ratio of rotational-kinetic energy to potential energy, T/|W|, and the peak frequency, $f_{\rm peak}$, is proportional to the square root of the central rest-mass density. In this work, we exploit these relations to suggest that it could be possible to use such waveforms to infer protoneutron star properties from a future gravitational wave observation, if the distance and inclination are well known. Our approach relies on the ability to describe a subset of the waveforms in the early post-bounce phase in a simple form depending only on two parameters, $Δh$ and $f_{\rm peak}$. We use this template to perform a Bayesian inference analysis of waveform injections in Gaussian colored noise for a network of three gravitational wave detectors formed by Advanced LIGO and Advanced Virgo. We show that, for a galactic event, it is possible to recover the peak frequency and amplitude with an accuracy better than 10% for about 80% and 60% of the signals, respectively, given known distance and inclination angle. However, inference on waveforms from outside the Richers catalog is not reliable, indicating a need for carefully verified waveforms of the first 10 ms after bounce of rapidly rotating supernovae of different progenitors with agreement between different codes. △ Less

Submitted 26 February, 2024; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 23 pages, 17 figures, accepted for publication in PRD

Comparison of neural network architectures for feature extraction from binary black hole merger waveforms

Authors: Osvaldo Gramaxo Freitas, Juan Calderón Bustillo, José A. Font, Solange Nunes, Antonio Onofre, Alejandro Torres-Forné

Abstract: We evaluate several neural-network architectures, both convolutional and recurrent, for gravitational-wave time-series feature extraction by performing point parameter estimation on noisy waveforms from binary-black-hole mergers. We build datasets of 100,000 elements for each of four different waveform models (or approximants) in order to test how approximant choice affects feature extraction. Our… ▽ More We evaluate several neural-network architectures, both convolutional and recurrent, for gravitational-wave time-series feature extraction by performing point parameter estimation on noisy waveforms from binary-black-hole mergers. We build datasets of 100,000 elements for each of four different waveform models (or approximants) in order to test how approximant choice affects feature extraction. Our choices include \texttt{SEOBNRv4P} and \texttt{IMRPhenomPv3}, which contain only the dominant quadrupole emission mode, alongside \texttt{IMRPhenomPv3HM} and \texttt{NRHybSur3dq8}, which also account for high-order modes. Each dataset element is injected into detector noise corresponding to the third observing run of the LIGO-Virgo-KAGRA (LVK) collaboration. We identify the Temporal Convolutional Network (TCN) architecture as the overall best performer in terms of training and validation losses and absence of overfitting to data. Comparison of results between datasets shows that the choice of waveform approximant for the creation of a dataset conditions the feature extraction ability of a trained network. Hence, care should be taken when building a dataset for the training of neural networks, as certain approximants may result in better network convergence of evaluation metrics. However, this performance does not necessarily translate to data which is more faithful to numerical relativity simulations. We also apply this network on actual signals from LVK runs, finding that its feature-extracting performance can be effective on real data. △ Less

Submitted 1 August, 2023; v1 submitted 31 July, 2023; originally announced July 2023.

Comments: Fixed Fig. 7 not being rendered

Towards numerical-relativity informed effective-one-body waveforms for dynamical capture black hole binaries

Authors: Tomas Andrade, Juan Trenado, Simone Albanesi, Rossella Gamba, Sebastiano Bernuzzi, Alessandro Nagar, Juan Calderon-Bustillo, Nicolas Sanchis-Gual, Jose A. Font, William Cook, Boris Daszuta, Francesco Zappa, David Radice

Abstract: Dynamical captures of black holes may take place in dense stellar media due to the emission of gravitational radiation during a close passage. Detection of such events requires detailed modelling, since their phenomenology qualitatively differs from that of quasi-circular binaries. Very few models can deliver such waveforms, and none includes information from Numerical Relativity (NR) simulations… ▽ More Dynamical captures of black holes may take place in dense stellar media due to the emission of gravitational radiation during a close passage. Detection of such events requires detailed modelling, since their phenomenology qualitatively differs from that of quasi-circular binaries. Very few models can deliver such waveforms, and none includes information from Numerical Relativity (NR) simulations of non quasi-circular coalescences. In this study we present a first step towards a fully NR-informed Effective One Body (EOB) model of dynamical captures. We perform 14 new simulations of single and double encounter mergers, and use this data to inform the merger-ringdown model of the TEOBResumS-Dali approximant. We keep the initial energy approximately fixed to the binary mass, and vary the mass-rescaled, dimensionless angular momentum in the range $(0.6, 1.1)$, the mass ratio in $(1, 2.15)$ and aligned dimensionless spins in $(-0.5, 0.5)$. We find that the model is able to match NR to $97%$, improving previous performances, without the need of modifying the base-line template. Upon NR informing the model, this improves to $99%$ with the exception of one outlier corresponding to a direct plunge. The maximum EOBNR phase difference at merger for the uninformed model is of $0.15$ radians, which is reduced to $0.1$ radians after the NR information is introduced. We outline the steps towards a fully informed EOB model of dynamical captures, and discuss future improvements. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 21 pages, 15 figures

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the Third Observing Run of Advanced LIGO and Advanced Virgo

Authors: Marek J. Szczepańczyk, Yanyan Zheng, Javier M. Antelis, Michael Benjamin, Marie-Anne Bizouard, Alejandro Casallas-Lagos, Pablo Cerdá-Durán, Derek Davis, Dorota Gondek-Rosińska, Sergey Klimenko, Claudia Moreno, Martin Obergaulinger, Jade Powell, Dymetris Ramirez, Brad Ratto, Colter Richarson, Abhinav Rijal, Amber L. Stuver, Paweł Szewczyk, Gabriele Vedovato, Michele Zanolin, Imre Bartos, Shubhagata Bhaumik, Tomasz Bulik, Marco Drago , et al. (13 additional authors not shown)

Abstract: We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed optically within 30 Mpc during the third observing run of Advanced LIGO and Advanced Virgo. No gravitational wave associated with a core-collapse supernova has been identified. We then report the detection efficiency for a variety of possible gravitational-wave emissions. For ne… ▽ More We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed optically within 30 Mpc during the third observing run of Advanced LIGO and Advanced Virgo. No gravitational wave associated with a core-collapse supernova has been identified. We then report the detection efficiency for a variety of possible gravitational-wave emissions. For neutrino-driven explosions, the distance at which we reach 50% detection efficiency is up to 8.9 kpc, while more energetic magnetorotationally-driven explosions are detectable at larger distances. The distance reaches for selected models of the black hole formation, and quantum chromodynamics phase transition are also provided. We then constrain the core-collapse supernova engine across a wide frequency range from 50 Hz to 2 kHz. The upper limits on gravitational-wave energy and luminosity emission are at low frequencies down to $10^{-4}\,M_\odot c^2$ and $6 \times 10^{-4}\,M_\odot c^2$/s, respectively. The upper limits on the proto-neutron star ellipticity are down to 3 at high frequencies. Finally, by combining the results obtained with the data from the first and second observing runs of LIGO and Virgo, we improve the constraints of the parameter spaces of the extreme emission models. Specifically, the proto-neutron star ellipticities for the long-lasting bar mode model are down to 1 for long emission (1 s) at high frequency. △ Less

Submitted 3 July, 2024; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 20 pages, 11 figures; https://dcc.ligo.org/LIGO-P2200361

SP(k) -- A hydrodynamical simulation-based model for the impact of baryon physics on the non-linear matter power spectrum

Authors: Jaime Salcido, Ian G. McCarthy, Juliana Kwan, Amol Upadhye, Andreea S. Font

Abstract: Upcoming large-scale structure surveys will measure the matter power spectrum to approximately percent level accuracy with the aim of searching for evidence for new physics beyond the standard model of cosmology. In order to avoid biasing our conclusions, the theoretical predictions need to be at least as accurate as the measurements for a given choice of cosmological parameters. However, recent t… ▽ More Upcoming large-scale structure surveys will measure the matter power spectrum to approximately percent level accuracy with the aim of searching for evidence for new physics beyond the standard model of cosmology. In order to avoid biasing our conclusions, the theoretical predictions need to be at least as accurate as the measurements for a given choice of cosmological parameters. However, recent theoretical work has shown that complex physical processes associated with galaxy formation (particularly energetic feedback processes associated with stars and especially supermassive black holes) can alter the predictions by many times larger than the required accuracy. Here we present $\texttt{SP(k)}$, a model for the effects of baryon physics on the non-linear matter power spectrum based on a new large suite of hydrodynamical simulations. Specifically, the ANTILLES suite consists of 400 simulations spanning a very wide range of the "feedback landscape" and show that the effects of baryons on the matter power spectrum can be understood at approaching the percent level in terms of the mean baryon fraction of haloes, at scales of up to $k \lesssim 10 \, h \, $Mpc$^{-1}$ and redshifts up to $z=3$. For the range of scales and redshifts that will be probed by forthcoming cosmic shear measurements, most of the effects are driven by galaxy group-mass haloes ($M \sim 10^{13-14}$ M$_\odot$). We present a simple Python implementation of our model, available at $\href{https://github.com/jemme07/pyspk}{\mathrm{https{:}//github.com/jemme07/pyspk}}$, which can be used to incorporate baryon effects in standard gravity-only predictions, allowing for marginalisation over baryon physics within cosmological pipelines. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: Refereed version accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS). 16 pages, 11 figures

Gradient-Annihilated PINNs for Solving Riemann Problems: Application to Relativistic Hydrodynamics

Authors: Antonio Ferrer-Sánchez, José D. Martín-Guerrero, Roberto Ruiz de Austri, Alejandro Torres-Forné, José A. Font

Abstract: We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that requires the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting funct… ▽ More We present a novel methodology based on Physics-Informed Neural Networks (PINNs) for solving systems of partial differential equations admitting discontinuous solutions. Our method, called Gradient-Annihilated PINNs (GA-PINNs), introduces a modified loss function that requires the model to partially ignore high-gradients in the physical variables, achieved by introducing a suitable weighting function. The method relies on a set of hyperparameters that control how gradients are treated in the physical loss and how the activation functions of the neural model are dynamically accounted for. The performance of our GA-PINN model is demonstrated by solving Riemann problems in special relativistic hydrodynamics, extending earlier studies with PINNs in the context of the classical Euler equations. The solutions obtained with our GA-PINN model correctly describe the propagation speeds of discontinuities and sharply capture the associated jumps. We use the relative $l^{2}$ error to compare our results with the exact solution of special relativistic Riemann problems, used as the reference ``ground truth'', and with the error obtained with a second-order, central, shock-capturing scheme. In all problems investigated, the accuracy reached by our GA-PINN model is comparable to that obtained with a shock-capturing scheme and significantly higher than that achieved by a baseline PINN algorithm. An additional benefit worth stressing is that our PINN-based approach sidesteps the costly recovery of the primitive variables from the state vector of conserved ones, a well-known drawback of grid-based solutions of the relativistic hydrodynamics equations. Due to its inherent generality and its ability to handle steep gradients, the GA-PINN method discussed could be a valuable tool to model relativistic flows in astrophysics and particle physics, characterized by the prevalence of discontinuous solutions. △ Less

Submitted 19 May, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 25 pages, 16 figures

Birth of baby universes from gravitational collapse in a modified-gravity scenario

Authors: Andreu Masó-Ferrando, Nicolas Sanchis-Gual, José A. Font, Gonzalo J. Olmo

Abstract: We consider equilibrium models of spherical boson stars in Palatini $f(\mathcal{R})=\mathcal{R}+ξ\mathcal{R}^2$ gravity and study their collapse when perturbed. The Einstein-Klein-Gordon system is solved using a recently established correspondence in an Einstein frame representation. We find that, in that frame, the endpoint is a nonrotating black hole surrounded by a quasi-stationary cloud of sca… ▽ More We consider equilibrium models of spherical boson stars in Palatini $f(\mathcal{R})=\mathcal{R}+ξ\mathcal{R}^2$ gravity and study their collapse when perturbed. The Einstein-Klein-Gordon system is solved using a recently established correspondence in an Einstein frame representation. We find that, in that frame, the endpoint is a nonrotating black hole surrounded by a quasi-stationary cloud of scalar field. However, the dynamics in the $f(\mathcal{R})$ frame is dramatically different. The innermost region of the collapsing object exhibits the formation of a finite-size, exponentially-expanding $\textit{ baby universe}$ connected with the outer (parent) universe via a minimal area surface (a throat or umbilical cord). Our simulations indicate that this surface is at all times hidden inside a horizon, causally disconnecting the baby universe from observers above the horizon. The implications of our findings in other areas of gravitational physics are also discussed. △ Less

Submitted 7 July, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

Magnetized tori with magnetic polarization around Kerr black holes: Variable angular momentum discs

Authors: Sergio Gimeno-Soler, Oscar. M. Pimentel, Fabio D. Lora-Clavijo, Alejandro Cruz-Osorio, José A. Font

Abstract: Analytical models of magnetized, geometrically thick disks are relevant to understand the physical conditions of plasma around compact objects and to explore its emitting properties. This has become increasingly important in recent years in the light of the Event Horizon Telescope observations of Sgr A$^*$ and M87. Models of thick disks around black holes usually consider constant angular momentum… ▽ More Analytical models of magnetized, geometrically thick disks are relevant to understand the physical conditions of plasma around compact objects and to explore its emitting properties. This has become increasingly important in recent years in the light of the Event Horizon Telescope observations of Sgr A$^*$ and M87. Models of thick disks around black holes usually consider constant angular momentum distributions and do not take into account the magnetic response of the fluid to applied magnetic fields. We present a generalization of our previous work on stationary models of magnetized accretion disks with magnetic polarization (Pimentel et al. 2018). This extension is achieved by accounting for non-constant specific angular momentum profiles, done through a two-parameter ansatz for those distributions. We build a large number of new equilibrium solutions of thick disks with magnetic polarization around Kerr black holes, selecting suitable parameter values within the intrinsically substantial parameter space of the models. We study the morphology and the physical properties of those solutions, finding qualitative changes with respect to the constant angular momentum tori solutions (Pimentel et al. 2018). However, the dependences found on the angular momentum distribution or on the black hole spin do not seem to be strong. Some of the new solutions, however, exhibit a local maximum of the magnetization function, absent in standard magnetized tori. Due to the enhanced development of the magneto-rotational instability as a result of magnetic susceptibility, those models might be particularly well-suited to investigate jet formation through general-relativistic MHD simulations. The new equilibrium solutions reported here can be used as initial data in numerical codes to assess the impact of magnetic susceptibility in the dynamics and observational properties of black hole-thick disk systems. △ Less

Submitted 8 July, 2024; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 19 pages, 10 figures. Accepted for publication in Physical Review D

MSC Class: 83C57; 76W05; 85A30; 83C55

Convolutional Neural Networks for the classification of glitches in gravitational-wave data streams

Authors: Tiago S. Fernandes, Samuel J. Vieira, Antonio Onofre, Juan Calderón Bustillo, Alejandro Torres-Forné, José A. Font

Abstract: We investigate the use of Convolutional Neural Networks (including the modern ConvNeXt network family) to classify transient noise signals (i.e.~glitches) and gravitational waves in data from the Advanced LIGO detectors. First, we use models with a supervised learning approach, both trained from scratch using the Gravity Spy dataset and employing transfer learning by fine-tuning pre-trained models… ▽ More We investigate the use of Convolutional Neural Networks (including the modern ConvNeXt network family) to classify transient noise signals (i.e.~glitches) and gravitational waves in data from the Advanced LIGO detectors. First, we use models with a supervised learning approach, both trained from scratch using the Gravity Spy dataset and employing transfer learning by fine-tuning pre-trained models in this dataset. Second, we also explore a self-supervised approach, pre-training models with automatically generated pseudo-labels. Our findings are very close to existing results for the same dataset, reaching values for the F1 score of 97.18% (94.15%) for the best supervised (self-supervised) model. We further test the models using actual gravitational-wave signals from LIGO-Virgo's O3 run. Although trained using data from previous runs (O1 and O2), the models show good performance, in particular when using transfer learning. We find that transfer learning improves the scores without the need for any training on real signals apart from the less than 50 chirp examples from hardware injections present in the Gravity Spy dataset. This motivates the use of transfer learning not only for glitch classification but also for signal classification. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: 15 pages, 14 figures

Prospects for the inference of inertial modes from hypermassive neutron stars with future gravitational-wave detectors

Authors: Miquel Miravet-Tenés, Florencia L. Castillo, Roberto De Pietri, Pablo Cerdá-Durán, José A. Font

Abstract: Some recent, long-term numerical simulations of binary neutron star mergers have shown that the long-lived remnants produced in such mergers might be affected by convective instabilities. Those would trigger the excitation of inertial modes, providing a potential method to improve our understanding of the rotational and thermal properties of neutron stars through the analysis of the modes' imprint… ▽ More Some recent, long-term numerical simulations of binary neutron star mergers have shown that the long-lived remnants produced in such mergers might be affected by convective instabilities. Those would trigger the excitation of inertial modes, providing a potential method to improve our understanding of the rotational and thermal properties of neutron stars through the analysis of the modes' imprint in the late post-merger gravitational-wave signal. In this paper we assess the detectability of those modes by injecting numerically generated post-merger waveforms into colored Gaussian noise of second-generation and future detectors. Signals are recovered using BayesWave, a Bayesian data-analysis algorithm that reconstructs them through a morphology-independent approach using series of sine-Gaussian wavelets. Our study reveals that current interferometers (i.e. the Handford-Livingston-Virgo network) recover the peak frequency of inertial modes only if the merger occurs at distances of up to 1 Mpc. For future detectors such as the Einstein Telescope, the range of detection increases by about a factor 10. △ Less

Submitted 10 May, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Open data from the third observing run of LIGO, Virgo, KAGRA and GEO

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Al-Jodah, C. Alléné, A. Allocca , et al. (1719 additional authors not shown)

Abstract: The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasti… ▽ More The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in April of 2019 and lasting six months, O3b starting in November of 2019 and lasting five months, and O3GK starting in April of 2020 and lasting 2 weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org. The main dataset, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 27 pages, 3 figures

Report number: LIGO-P2200316

Inference of proto-neutron star properties in core-collapse supernovae from a gravitational-wave detector network

Authors: T. Bruel, M-A. Bizouard, M. Obergaulinger, P. Maturana-Russel, A. Torres-Forné, P. Cerdá-Durán, N. Christensen, J. A. Font, R. Meyer

Abstract: The next Galactic core-collapse supernova (CCSN) will be a unique opportunity to study within a fully multi-messenger approach the explosion mechanism responsible for the formation of neutron stars and stellar-mass black holes. State-of-the-art numerical simulations of those events reveal the complexity of the gravitational-wave emission which is highly stochastic. This challenges the possibility… ▽ More The next Galactic core-collapse supernova (CCSN) will be a unique opportunity to study within a fully multi-messenger approach the explosion mechanism responsible for the formation of neutron stars and stellar-mass black holes. State-of-the-art numerical simulations of those events reveal the complexity of the gravitational-wave emission which is highly stochastic. This challenges the possibility to infer the properties of the compact remnant and of its progenitor using the information encoded in the waveforms. In this paper we take further steps in a program we recently initiated to overcome those difficulties. In particular we show how oscillation modes of the proto-neutron star, highly visible in the gravitational-wave signal, can be used to reconstruct the time evolution of their physical properties. Extending our previous work where only the information from a single detector was used we here describe a new data-analysis pipeline that coherently combines gravitational-wave detectors' data and infers the time evolution of a combination of the mass and radius of the compact remnant. The performance of the method is estimated employing waveforms from 2D and 3D CCSN simulations covering a progenitor mass range between 11$\mathrm{M_{\odot}}$\, and 40$\mathrm{M_{\odot}}$\, and different equations of state for both a network of up to five second-generation detectors and the proposed third-generation detectors Einstein Telescope and Cosmic Explorer. Our study shows that it will be possible to infer PNS properties for CCSN events occurring in the vicinity of the Milky Way, up to the Large Magellanic Cloud, with the current generation of gravitational-wave detectors. △ Less

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

Journal ref: Phys. Rev. D 107, 083029 (2023)

Bondi-Hoyle-Lyttleton accretion onto a rotating black hole with ultralight scalar hair

Authors: Alejandro Cruz-Osorio, Luciano Rezzolla, Fabio Duvan Lora-Clavijo, José Antonio Font, Carlos Herdeiro, Eugen Radu

Abstract: We present a numerical study of relativistic Bondi-Hoyle-Lyttleton (BHL) accretion onto an asymptotically flat black hole with synchronized hair. The hair is sourced by an ultralight, complex scalar field, minimally coupled to Einstein's gravity. Our simulations consider a supersonic flow parametrized by the asymptotic values of the fluid quantities and a sample of hairy black holes with different… ▽ More We present a numerical study of relativistic Bondi-Hoyle-Lyttleton (BHL) accretion onto an asymptotically flat black hole with synchronized hair. The hair is sourced by an ultralight, complex scalar field, minimally coupled to Einstein's gravity. Our simulations consider a supersonic flow parametrized by the asymptotic values of the fluid quantities and a sample of hairy black holes with different masses, angular momenta, and amount of scalar hair. For all models, steady-state BHL accretion solutions are attained that are characterized by the presence of a shock-cone and a stagnation point downstream. For the models of the sample with the largest component of scalar field, the shock-cone envelops fully the black hole, transitioning into a bow-shock, and the stagnation points move further away downstream. Analytical expressions for the mass accretion rates are obtained after fitting the numerical results, which can be used to analyze black-hole formation scenarios in the presence of ultralight scalar fields. The formation of a shock-cone leads to regions where sound waves can be trapped and resonant oscillations excited. We measure the frequencies of such quasi-periodic oscillations and point out a possible association with quasi-periodic oscillations in the X-ray light curve of Sgr~A* and microquasars. △ Less

Submitted 16 June, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: 24 pages, 3 tables, 10 figures. Minor changes

Journal ref: Journal of Cosmology and Astroparticle Physics 2023 (08), 057

Solving the Teukolsky equation with physics-informed neural networks

Authors: Raimon Luna, Juan Calderón Bustillo, Juan José Seoane Martínez, Alejandro Torres-Forné, José A. Font

Abstract: We use physics-informed neural networks (PINNs) to compute the first quasi-normal modes of the Kerr geometry via the Teukolsky equation. This technique allows us to extract the complex frequencies and separation constants of the equation without the need for sophisticated numerical techniques, and with an almost immediate implementation under the \texttt{PyTorch} framework. We are able to compute… ▽ More We use physics-informed neural networks (PINNs) to compute the first quasi-normal modes of the Kerr geometry via the Teukolsky equation. This technique allows us to extract the complex frequencies and separation constants of the equation without the need for sophisticated numerical techniques, and with an almost immediate implementation under the \texttt{PyTorch} framework. We are able to compute the oscillation frequencies and damping times for arbitrary black hole spins and masses, with accuracy typically below the percentual level as compared to the accepted values in the literature. We find that PINN-computed quasi-normal modes are indistinguishable from those obtained through existing methods at signal-to-noise ratios (SNRs) larger than 100, making the former reliable for gravitational-wave data analysis in the mid term, before the arrival of third-generation detectors like LISA or the Einstein Telescope, where SNRs of ${\cal O}(1000)$ might be achieved. △ Less

Submitted 8 April, 2024; v1 submitted 12 December, 2022; originally announced December 2022.

Comments: 12 pages, 7 figures. v2: Matches published version, minor typos corrected

Journal ref: Phys.Rev.D 107 (2023) 6, 064025

A revised formalism for slowly-rotating superfluid neutron stars in general relativity

Authors: Eneko Aranguren, José A. Font, Nicolas Sanchis-Gual, Raül Vera

Abstract: We discuss slowly-rotating, general relativistic, superfluid neutron stars in the Hartle-Thorne formulation. The composition of the stars is described by a simple two-fluid model which accounts for superfluid neutrons and all other constituents. We apply a perturbed matching framework to derive a new formalism for slowly-rotating superfluid neutron stars, valid up to second-order perturbation theo… ▽ More We discuss slowly-rotating, general relativistic, superfluid neutron stars in the Hartle-Thorne formulation. The composition of the stars is described by a simple two-fluid model which accounts for superfluid neutrons and all other constituents. We apply a perturbed matching framework to derive a new formalism for slowly-rotating superfluid neutron stars, valid up to second-order perturbation theory, building on the original formulation reported by Andersson and Comer in 2001. The present study constitutes an extension of previous work in the single-fluid case where it was shown that the Hartle-Thorne formalism needs to be amended since it does not provide the correct results when the energy density does not vanish at the surface of the star. We discuss in detail the corrections that need to be applied to the original two-fluid formalism in order to account for non vanishing energy densities at the boundary. In the process, we also find a correction needed in the computation of the deformation of the stellar surface in the original two-fluid model in all cases (irrespective of the value of the energy density at the surface). The discrepancies found between the two formalisms are illustrated by building numerical stellar models, focusing on the comparison in the calculation of the stellar mass, the deformation of the star, and in the Kepler limit of rotation. In particular, using a toy-model equation of state for which the energy density does not vanish at the boundary of the star we demonstrate that the corrections to the formalism we find impact the structure of slowly-rotating superfluid neutron stars in a significant way. △ Less

Submitted 20 March, 2023; v1 submitted 5 December, 2022; originally announced December 2022.

Comments: 24 pages, 16 figures: v2. minor corrections and changes to match the published version

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

A correlation between accreted stellar kinematics and dark matter halo spin in the ARTEMIS simulations

Authors: Adam M. Dillamore, Vasily Belokurov, N. Wyn Evans, Andreea S. Font

Abstract: We report a correlation between the presence of a Gaia-Sausage-Enceladus (GSE) analogue and dark matter halo spin in the ARTEMIS simulations of Milky Way-like galaxies. The haloes which contain a large population of accreted stars on highly radial orbits (like the GSE) have lower spin on average than their counterparts with more isotropic stellar velocity distributions. The median modified spin pa… ▽ More We report a correlation between the presence of a Gaia-Sausage-Enceladus (GSE) analogue and dark matter halo spin in the ARTEMIS simulations of Milky Way-like galaxies. The haloes which contain a large population of accreted stars on highly radial orbits (like the GSE) have lower spin on average than their counterparts with more isotropic stellar velocity distributions. The median modified spin parameters $λ^\prime$ differ by a factor of $\sim1.7$ at the present-day, with a similar value when the haloes far from virial equilibrium are removed. We also show that accreted stars make up a smaller proportion of the stellar populations in haloes containing a GSE analogue, and are stripped from satellites with stellar masses typically $\sim4$ times smaller. Our findings suggest that the higher spin of DM haloes without a GSE-like feature is due to mergers with large satellites of stellar mass $\sim10^{10}M_\odot$, which do not result in prominent radially anisotropic features like the GSE. △ Less

Submitted 5 December, 2022; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 6 pages, 5 figures, accepted for publication by MNRAS Letters

Virgo Detector Characterization and Data Quality: tools

Authors: F. Acernese, M. Agathos, A. Ain, S. Albanesi, A. Allocca, A. Amato, T. Andrade, N. Andres, M. Andrés-Carcasona, T. Andrić, S. Ansoldi, S. Antier, T. Apostolatos, E. Z. Appavuravther, M. Arène, N. Arnaud, M. Assiduo, S. Assis de Souza Melo, P. Astone, F. Aubin, S. Babak, F. Badaracco, M. K. M. Bader, S. Bagnasco, J. Baird , et al. (469 additional authors not shown)

Abstract: Detector characterization and data quality studies -- collectively referred to as {\em DetChar} activities in this article -- are paramount to the scientific exploitation of the joint dataset collected by the LIGO-Virgo-KAGRA global network of ground-based gravitational-wave (GW) detectors. They take place during each phase of the operation of the instruments (upgrade, tuning and optimization, dat… ▽ More Detector characterization and data quality studies -- collectively referred to as {\em DetChar} activities in this article -- are paramount to the scientific exploitation of the joint dataset collected by the LIGO-Virgo-KAGRA global network of ground-based gravitational-wave (GW) detectors. They take place during each phase of the operation of the instruments (upgrade, tuning and optimization, data taking), are required at all steps of the dataflow (from data acquisition to the final list of GW events) and operate at various latencies (from near real-time to vet the public alerts to offline analyses). This work requires a wide set of tools which have been developed over the years to fulfill the requirements of the various DetChar studies: data access and bookkeeping; global monitoring of the instruments and of the different steps of the data processing; studies of the global properties of the noise at the detector outputs; identification and follow-up of noise peculiar features (whether they be transient or continuously present in the data); quick processing of the public alerts. The present article reviews all the tools used by the Virgo DetChar group during the third LIGO-Virgo Observation Run (O3, from April 2019 to March 2020), mainly to analyse the Virgo data acquired at EGO. Concurrently, a companion article focuses on the results achieved by the DetChar group during the O3 run using these tools. △ Less

Submitted 25 March, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Comments: 44 pages, 16 figures. New version, resubmitted to Class. and Quantum Grav. This is the "Tools" part of preprint arXiv:2205.01555 [gr-qc] which has been split into two companion articles: one about the tools and methods, the other about the analyses of the O3 Virgo data

Virgo Detector Characterization and Data Quality: results from the O3 run

Authors: F. Acernese, M. Agathos, A. Ain, S. Albanesi, A. Allocca, A. Amato, T. Andrade, N. Andres, M. Andrés-Carcasona, T. Andrić, S. Ansoldi, S. Antier, T. Apostolatos, E. Z. Appavuravther, M. Arène, N. Arnaud, M. Assiduo, S. Assis de Souza Melo, P. Astone, F. Aubin, S. Babak, F. Badaracco, M. K. M. Bader, S. Bagnasco, J. Baird , et al. (469 additional authors not shown)

Abstract: The Advanced Virgo detector has contributed with its data to the rapid growth of the number of detected gravitational-wave (GW) signals in the past few years, alongside the two Advanced LIGO instruments. First during the last month of the Observation Run 2 (O2) in August 2017 (with, most notably, the compact binary mergers GW170814 and GW170817), and then during the full Observation Run 3 (O3): an… ▽ More The Advanced Virgo detector has contributed with its data to the rapid growth of the number of detected gravitational-wave (GW) signals in the past few years, alongside the two Advanced LIGO instruments. First during the last month of the Observation Run 2 (O2) in August 2017 (with, most notably, the compact binary mergers GW170814 and GW170817), and then during the full Observation Run 3 (O3): an 11-months data taking period, between April 2019 and March 2020, that led to the addition of about 80 events to the catalog of transient GW sources maintained by LIGO, Virgo and now KAGRA. These discoveries and the manifold exploitation of the detected waveforms require an accurate characterization of the quality of the data, such as continuous study and monitoring of the detector noise sources. These activities, collectively named {\em detector characterization and data quality} or {\em DetChar}, span the whole workflow of the Virgo data, from the instrument front-end hardware to the final analyses. They are described in details in the following article, with a focus on the results achieved by the Virgo DetChar group during the O3 run. Concurrently, a companion article describes the tools that have been used by the Virgo DetChar group to perform this work. △ Less

Submitted 25 March, 2023; v1 submitted 14 October, 2022; originally announced October 2022.

Comments: 57 pages, 18 figures. New version, resubmitted to Class. and Quantum Grav. This is the "Results" part of preprint arXiv:2205.01555 [gr-qc] which has been split into two companion articles: one about the tools and methods, the other about the analyses of the O3 Virgo data

Search for gravitational-wave transients associated with magnetar bursts in Advanced LIGO and Advanced Virgo data from the third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1645 additional authors not shown)

Abstract: Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bu… ▽ More Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bursts come from two magnetars, SGR 1935$+$2154 and Swift J1818.0$-$1607. We also include three other electromagnetic burst events detected by Fermi GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper bounds on the root-sum-square of the integrated gravitational-wave strain that reach $2.2 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at 100 Hz for the short-duration search and $8.7 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at $450$ Hz for the long-duration search, given a detection efficiency of 50%. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to $1.8 \times 10^{-22}$ $/\sqrt{\text{Hz}}$. Using the estimated distance to each magnetar, we derive upper bounds on the emitted gravitational-wave energy of $3.2 \times 10^{43}$ erg ($7.3 \times 10^{43}$ erg) for SGR 1935$+$2154 and $8.2 \times 10^{42}$ erg ($2.8 \times 10^{43}$ erg) for Swift J1818.0$-$1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935$+$2154 with available fluence information. The lowest of these ratios is $3 \times 10^3$. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 30 pages with appendices, 5 figures, 10 tables

Report number: LIGO-P2100387

GOGREEN: a critical assessment of environmental trends in cosmological hydrodynamical simulations at z ~ 1

Authors: Egidijus Kukstas, Michael L. Balogh, Ian G. McCarthy, Yannick M. Bahe, Gabriella De Lucia, Pascale Jablonka, Benedetta Vulcani, Devontae C. Baxter, Andrea Biviano, Pierluigi Cerulo, Jeffrey C. Chan, M. C. Cooper, Ricardo Demarco, Alexis Finoguenov, Andreea S. Font, Chris Lidman, Justin Marchioni, Sean McGee, Adam Muzzin, Julie Nantais, Lyndsay Old, Irene Pintos-Castro, Bianca Poggianti, Andrew M. M. Reeves, Gregory Rudnick , et al. (6 additional authors not shown)

Abstract: Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been… ▽ More Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been relatively few comparisons at higher redshifts to date. Here we confront three state-of-the-art suites of simulations (BAHAMAS+MACSIS, EAGLE+Hydrangea, IllustrisTNG) with state-of-the-art observations of the field and cluster environments from the COSMOS/UltraVISTA and GOGREEN surveys, respectively, at z ~ 1 to assess the realism of the simulations and gain insight into the evolution of environmental quenching. We show that while the simulations generally reproduce the stellar content and the stellar mass functions of quiescent and star-forming galaxies in the field, all the simulations struggle to capture the observed quenching of satellites in the cluster environment, in that they are overly efficient at quenching low-mass satellites. Furthermore, two of the suites do not sufficiently quench the highest-mass galaxies in clusters, perhaps a result of insufficient feedback from AGN. The origin of the discrepancy at low stellar masses (Mstar <~ 1E10 Msun), which is present in all the simulations in spite of large differences in resolution, feedback implementations, and hydrodynamical solvers, is unclear. The next generation of simulations, which will push to significantly higher resolution and also include explicit modelling of the cold interstellar medium, may help to shed light on the low-mass tension. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 20 pages, 10 figures, resubmitted to MNRAS following minor revisions

Machine-Learning Love: classifying the equation of state of neutron stars with Transformers

Authors: Gonçalo Gonçalves, Márcio Ferreira, João Aveiro, Antonio Onofre, Felipe F. Freitas, Constança Providência, José A. Font

Abstract: The use of the Audio Spectrogram Transformer (AST) model for gravitational-wave data analysis is investigated. The AST machine-learning model is a convolution-free classifier that captures long-range global dependencies through a purely attention-based mechanism. In this paper a model is applied to a simulated dataset of inspiral gravitational wave signals from binary neutron star coalescences, bu… ▽ More The use of the Audio Spectrogram Transformer (AST) model for gravitational-wave data analysis is investigated. The AST machine-learning model is a convolution-free classifier that captures long-range global dependencies through a purely attention-based mechanism. In this paper a model is applied to a simulated dataset of inspiral gravitational wave signals from binary neutron star coalescences, built from five distinct, cold equations of state (EOS) of nuclear matter. From the analysis of the mass dependence of the tidal deformability parameter for each EOS class it is shown that the AST model achieves a promising performance in correctly classifying the EOS purely from the gravitational wave signals, especially when the component masses of the binary system are in the range $[1,1.5]M_{\odot}$. Furthermore, the generalization ability of the model is investigated by using gravitational-wave signals from a new EOS not used during the training of the model, achieving fairly satisfactory results. Overall, the results, obtained using the simplified setup of noise-free waveforms, show that the AST model, once trained, might allow for the instantaneous inference of the cold nuclear matter EOS directly from the inspiral gravitational-wave signals produced in binary neutron star coalescences. △ Less

Submitted 15 October, 2022; originally announced October 2022.

Comments: 11 pages, 11 figures

Dictionary learning: a novel approach to detecting binary black holes in the presence of Galactic noise with LISA

Authors: Charles Badger, Katarina Martinovic, Alejandro Torres-Forné, Mairi Sakellariadou, José A. Font

Abstract: The noise produced by the inspiral of millions of white dwarf binaries in the Milky Way may pose a threat to one of the main goals of the space-based LISA mission: the detection of massive black hole binary mergers. We present a novel study for reconstruction of merger waveforms in the presence of Galactic confusion noise using dictionary learning. We discuss the limitations of untangling signals… ▽ More The noise produced by the inspiral of millions of white dwarf binaries in the Milky Way may pose a threat to one of the main goals of the space-based LISA mission: the detection of massive black hole binary mergers. We present a novel study for reconstruction of merger waveforms in the presence of Galactic confusion noise using dictionary learning. We discuss the limitations of untangling signals from binaries with total mass from $10^2 M_{\odot}$ to $10^4 M_{\odot}$. Our method proves extremely successful for binaries with total mass greater than $\sim 3\times 10^3$ $ M_{\odot}$ up to redshift 3 in conservative scenarios, and up to redshift 7.5 in optimistic scenarios. In addition, consistently good waveform reconstruction of merger events is found if the signal-to-noise ratio is approximately 5 or greater. △ Less

Submitted 14 October, 2022; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: 6 pages, 3 figures

Assessment of a new sub-grid model for magneto-hydrodynamical turbulence. I. Magnetorotational instability

Authors: Miquel Miravet-Tenés, Pablo Cerdá-Durán, Martin Obergaulinger, José A. Font

Abstract: Insufficient numerical resolution of grid-based, direct numerical simulations (DNS) hampers the development of instabilitydriven turbulence at small (unresolved) scales. As an alternative to DNS, sub-grid models can potentially reproduce the effects of turbulence at small scales in terms of the resolved scales, and hence can capture physical effects with less computational resources. We present a… ▽ More Insufficient numerical resolution of grid-based, direct numerical simulations (DNS) hampers the development of instabilitydriven turbulence at small (unresolved) scales. As an alternative to DNS, sub-grid models can potentially reproduce the effects of turbulence at small scales in terms of the resolved scales, and hence can capture physical effects with less computational resources. We present a new sub-grid model, the MHD-instability-induced-turbulence (MInIT) mean-field model. MInIT is a physically motivated model based on the evolution of the turbulent (Maxwell, Reynolds, and Faraday) stress tensors and their relation with the turbulent energy densities of the magneto-rotational (MRI) and parasitic instabilities, modeled with two partial differential evolution equations with stiff source terms. Their solution allows obtaining the turbulent stress tensors through the constant coefficients that link them to the energy densities. The model is assessed using data from MRI in-box DNS and applying a filtering operation to compare the filtered data with that from the model. Using the $L_2$-norm as the metric for the comparison, we find less than one order-of-magnitude difference between the two sets of data. No dependence on filter size or length scale of unresolved scales is found, as opposed to results using the gradient model (which we also use to contrast our model) in which the $L_2$-norm of some of the stresses increases with filter size. We conclude that MInIT can help DNS by properly capturing small-scale turbulent stresses which has potential implications on the dynamics of highly-magnetized rotating compact objects, such as those formed during binary neutron star mergers. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: 20 pages, 15 figures

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

Impact of the wave-like nature of Proca stars on their gravitational-wave emission

Authors: Nicolas Sanchis-Gual, Juan Calderón Bustillo, Carlos Herdeiro, Eugen Radu, José A. Font, Samson H. W. Leong, Alejandro Torres-Forné

Abstract: We present a systematic study of the dynamics and gravitational-wave emission of head-on collisions of spinning vector boson stars, known as Proca stars. To this aim we build a catalogue of about 800 numerical-relativity simulations of such systems. We find that the wave-like nature of bosonic stars has a large impact on the gravitational-wave emission. In particular, we show that the initial rela… ▽ More We present a systematic study of the dynamics and gravitational-wave emission of head-on collisions of spinning vector boson stars, known as Proca stars. To this aim we build a catalogue of about 800 numerical-relativity simulations of such systems. We find that the wave-like nature of bosonic stars has a large impact on the gravitational-wave emission. In particular, we show that the initial relative phase $Δε=ε_1-ε_2$ of the two complex fields forming the stars (or equivalently, the relative phase at merger) strongly impacts both the emitted gravitational-wave energy and the corresponding mode structure. This leads to a non-monotonic dependence of the emission on the frequency of the secondary star $ω_2$, for fixed frequency $ω_1$ of the primary. This phenomenology, which has not been found for the case of black-hole mergers, reflects the distinct ability of the Proca field to interact with itself in both constructive and destructive manners. We postulate this may serve as a smoking gun to shed light on the possible existence of these objects. △ Less

Submitted 22 March, 2023; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: Updated version

Journal ref: Phys. Rev. D 106, 124011 - Published 6 December 2022

Identification of Binary Neutron Star Mergers in Gravitational-Wave Data Using YOLO One-Shot Object Detection

Authors: João Aveiro, Felipe F. Freitas, Márcio Ferreira, Antonio Onofre, Constança Providência, Gonçalo Gonçalves, José A. Font

Abstract: We demonstrate the application of the YOLOv5 model, a general purpose convolution-based single-shot object detection model, in the task of detecting binary neutron star (BNS) coalescence events from gravitational-wave data of current generation interferometer detectors. We also present a thorough explanation of the synthetic data generation and preparation tasks based on approximant waveform model… ▽ More We demonstrate the application of the YOLOv5 model, a general purpose convolution-based single-shot object detection model, in the task of detecting binary neutron star (BNS) coalescence events from gravitational-wave data of current generation interferometer detectors. We also present a thorough explanation of the synthetic data generation and preparation tasks based on approximant waveform models used for the model training, validation and testing steps. Using this approach, we achieve mean average precision ($\text{mAP}_{[0.50]}$) values of 0.945 for a single class validation dataset and as high as 0.978 for test datasets. Moreover, the trained model is successful in identifying the GW170817 event in the LIGO H1 detector data. The identification of this event is also possible for the LIGO L1 detector data with an additional pre-processing step, without the need of removing the large glitch in the final stages of the inspiral. The detection of the GW190425 event is less successful, which attests to performance degradation with the signal-to-noise ratio. Our study indicates that the YOLOv5 model is an interesting approach for first-stage detection alarm pipelines and, when integrated in more complex pipelines, for real-time inference of physical source parameters. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 11 pages, 9 figures

Searching for vector boson-star mergers within LIGO-Virgo intermediate-mass black-hole merger candidates

Authors: Juan Calderon Bustillo, Nicolas Sanchis-Gual, Samson H. W. Leong, Koustav Chandra, Alejandro Torres-Forne, Jose A. Font, Carlos Herdeiro, Eugen Radu, Isaac C. F. Wong, T. G. F. Li

Abstract: We present the first systematic search for exotic compact mergers in Advanced LIGO and Virgo events. We compare the short gravitational-wave signals GW190521, GW190426$\_$190642, GW200220$\_$061928 and the trigger 200114$\_$020818 (or S200114f) to a new catalogue of 759 numerical simulations of head-on mergers of horizonless exotic compact objects known as Proca stars, interpreted as self-gravitat… ▽ More We present the first systematic search for exotic compact mergers in Advanced LIGO and Virgo events. We compare the short gravitational-wave signals GW190521, GW190426$\_$190642, GW200220$\_$061928 and the trigger 200114$\_$020818 (or S200114f) to a new catalogue of 759 numerical simulations of head-on mergers of horizonless exotic compact objects known as Proca stars, interpreted as self-gravitating lumps of (fuzzy) dark matter sourced by an ultralight (vector) bosonic particle. The Proca-star merger hypothesis is strongly rejected with respect to the black hole merger one by GW190426, weakly rejected by GW200220 and weakly favoured by GW190521 and S200114f. GW190521 and GW200220 yield highly consistent boson masses of $μ_{\rm B} = 8.69^{+0.61}_{-0.75}\times10^{-13}$ eV and $μ_{\rm B} = 9.13^{+1.18}_{-1.30}\times10^{-13}$ eV at the $90\%$ credible level. We conduct a preliminary population study of the compact binaries behind these events. Excluding (including) S200114f as a real event, and ignoring boson-mass consistencies across events, we estimate a fraction of Proca-star mergers of $ζ= 0.27^{+0.43}_{-0.25} \ (0.39^{+0.38}_{-0.33})$. We discuss the impact of boson-mass consistency across events in such estimates. Our results maintain GW190521 as a Proca-star merger candidate and pave the way towards population studies considering exotic compact objects. △ Less

Submitted 16 October, 2023; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: Version accepted in Phys.Rev.D. 27 pages, 16 Figures, 4 Appendixes

Impact of ultralight bosonic dark matter on the dynamical bar-mode instability of rotating neutron stars

Authors: Fabrizio Di Giovanni, Nicolás Sanchis-Gual, Davide Guerra, Miquel Miravet-Tenés, Pablo Cerdá-Durán, José Antonio Font

Abstract: We investigate the effects ultralight bosonic field dark matter may have on the dynamics of unstable differentially-rotating neutron stars prone to the bar-mode instability. To this aim we perform numerical simulations in general relativity of rotating neutron stars accreting an initial spherically symmetric bosonic field cloud, solving the Einstein-(complex, massive) Klein-Gordon-Euler and the Ei… ▽ More We investigate the effects ultralight bosonic field dark matter may have on the dynamics of unstable differentially-rotating neutron stars prone to the bar-mode instability. To this aim we perform numerical simulations in general relativity of rotating neutron stars accreting an initial spherically symmetric bosonic field cloud, solving the Einstein-(complex, massive) Klein-Gordon-Euler and the Einstein-(complex) Proca-Euler systems. We find that the presence of the bosonic field can critically modify the development of the bar-mode instability of neutron stars, depending on the total mass of the bosonic field and on the boson particle mass. In some cases, the accreting bosonic field can even quench the dominant $\ell=m=2$ mode of the bar-deformation by dynamically forming a mixed (fermion-boson) star that retains part of the angular momentum of the original neutron star. However, the mixed star undergoes the development of a mixed bar that leads to significant gravitational-wave emission, substantially different to that of the isolated neutron star. Our results indicate that dark-matter accretion in neutron stars could change the frequency of the expected emission of the bar-mode instability, which would have an important impact on ongoing searches for continuous gravitational waves. △ Less

Submitted 18 July, 2022; v1 submitted 2 June, 2022; originally announced June 2022.

Comments: 18 pages, 8 figures