Gravitational Wave Signal Research Articles

Gravitational wave signatures of cogenesis from a burdened PBH

We explore the possibility of explaining the observed dark matter (DM) relic abundance, along with matter-antimatter asymmetry, entirely from the evaporation of primordial black holes (PBH) beyond the semi-classical approximation. We find that, depending on the timing of modification to the semi-classical approximation and the efficiency of the backreaction, it is possible to produce the correct DM abundance for PBHs with masses ≳ 𝒪 (103) g, whereas producing the right amount of baryon asymmetry requires light PBHs with masses ≲ 𝒪 (103) g, satisfying bounds on the PBH mass from the Cosmic Microwave Background and Big Bang Nucleosynthesis. However, in a simplistic scenario, achieving both simultaneously is not feasible, typically because of the stringent Lyman-α constraint on warm dark matter mass. In addition to DM and baryon asymmetry, we also investigate the impact of memory burden on dark radiation, evaporated from PBH, constrained by the effective number of relativistic degrees of freedom Δ N eff. Furthermore, we demonstrate how induced gravitational waves from PBH density fluctuations can provide a window to test the memory-burden effects, thereby placing constraints on either the DM mass scale or the scale of leptogenesis.

An effective model for the tidal disruption of satellites undergoing minor mergers with axisymmetric primaries

According to the hierarchical formation paradigm, galaxies form through mergers of smaller entities and massive black holes (MBHs), if present at their centers, migrate to the nucleus of the newly formed galaxy, where they form binary systems. The formation and evolution of MBH binaries, and in particular their coalescence timescale, is highly relevant for current and future facilities aimed at detecting the gravitational wave signal produced by the MBHs close to coalescence. While most of the studies targeting this process are based on hydrodynamic simulations, the high computational cost makes a complete parameter space exploration prohibitive. Semianalytic approaches represent a valid alternative, but they require ad hoc prescriptions for the mass loss of the merging galaxies in minor mergers due to tidal stripping, which is not commonly considered or is at best modelled assuming very idealised geometries. In this work we propose a novel, effective model for the tidal stripping in axisymmetric potentials, to be implemented in semi-analytic models. We validated our semi-analytic approach against N-body simulations considering different galaxy sizes, inclinations, and eccentricities, finding only a moderate dependence on the orbit eccentricity. In particular, we find that, for almost circular orbits, our model mildly overestimates the mass loss, and this is due to the adjustment of the stellar distribution after the mass is removed. Nonetheless, the model exhibits a very good agreement with simulations in all the considered conditions, and thus represents an extremely powerful addition to semi-analytic calculations.

Bayesian inference of multi-messenger astrophysical data: Joint and coherent inference of gravitational waves and kilonovae

Context. Multi-messenger observations of binary neutron star mergers can provide information on the neutron star’s equation of state (EOS) above the nuclear saturation density by directly constraining the mass-radius diagram. Aims. We present a Bayesian framework for joint and coherent analyses of multi-messenger binary neutron star signals. As a first application, we analyze the gravitational-wave GW170817 and the kilonova (kN) AT2017gfo data. These results are then combined with the most recent X-ray pulsar analyses of PSR J0030+0451 and PSR J0740+6620 to obtain new EOS constraints. Methods. We extend the bajes infrastructure with a joint likelihood for multiple datasets, support for various semi-analytical kN models, and numerical-relativity (NR)-informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. The analysis of GW170817 used the TEOBResumS effective-one-body waveform template to model the gravitational-wave signal. The analysis of AT2017gfo used a baseline multicomponent spherically symmetric model for the kN light curves. Various constraints on the mass-radius diagram and neutron star properties were then obtained by resampling over a set of ten million parameterized EOSs, which was built under minimal assumptions (general relativity and causality). Results. We find that a joint and coherent approach improves the inference of the extrinsic parameters (distance) and, among the intrinsic parameters, the mass ratio. The inclusion of NR-informed relations marks a strong improvement over the case in which an agnostic prior is used on the intrinsic parameters. Comparing Bayes factors, we find that the two observations are better explained by the common source hypothesis only by assuming NR-informed relations. These relations break some of the degeneracies in the employed kN models. The EOS inference folding-in PSR J0952-0607 minimum-maximum mass, PSR J0030+0451 and PSR J0740+6620 data constrains, among other quantities, the neutron star radius to R1.4TOV = 12.30− 0.56+ 0.81 km(R1.4TOV = 13.20− 0.90+ 0.91 km) and the maximum mass to MmaxTOV = 2.28− 0.17+ 0.25M⊙(MmaxTOV = 2.32− 0.19+ 0.30M⊙), where the ST+PDT (PDT-U) analysis of Vinciguerra et al. (2024, ApJ, 961, 62) for PSR J0030+0451 was employed. Hence, the systematics on the PSR J0030+0451 data reduction currently dominate the mass-radius diagram constraints. Conclusions. We conclude that bajes delivers robust analyses in line with other state-of-the-art results in the literature. Strong EOS constraints are provided by pulsars observations, albeit with large systematics in some cases. Current gravitational-wave constraints are compatible with pulsar constraints and can further improve the latter.

Type-I two-Higgs-doublet model and gravitational waves from domain walls bounded by strings

The spontaneous breaking of a U(1) symmetry via an intermediate discrete symmetry may yield a hybrid topological defect of domain walls bounded by cosmic strings. The decay of this defect network leads to a unique gravitational wave signal spanning many orders in observable frequencies, that can be distinguished from signals generated by other sources. We investigate the production of gravitational waves from this mechanism in the context of the type-I two-Higgs-doublet model extended by a U(1)R symmetry, that simultaneously accommodates the seesaw mechanism, anomaly cancellation, and eliminates flavour-changing neutral currents. The gravitational wave spectrum produced by the string-bounded-wall network can be detected for U(1)R breaking scale from 1012 to 1015 GeV in forthcoming interferometers including LISA and Einstein Telescope, with a distinctive f3 slope and inflexion in the frequency range between microhertz and hertz.

Rapid identification of time-frequency domain gravitational wave signals from binary black holes using deep learning* *Supported by the National SKA Program of China (2022SKA0110200, 2022SKA0110203), the National Natural Science Foundation of China (12473001, 11975072, 11875102, 11835009), and the National 111 Project (B16009)

Recent developments in deep learning techniques have provided alternative and complementary approaches to the traditional matched-filtering methods for identifying gravitational wave (GW) signals. The rapid and accurate identification of GW signals is crucial to the advancement of GW physics and multi-messenger astronomy, particularly considering the upcoming fourth and fifth observing runs of LIGO-Virgo-KAGRA. In this study, we used the 2D U-Net algorithm to identify time-frequency domain GW signals from stellar-mass binary black hole (BBH) mergers. We simulated BBH mergers with component masses ranging from 7 to 50 and accounted for the LIGO detector noise. We found that the GW events in the first and second observation runs could all be clearly and rapidly identified. For the third observing run, approximately 80% of the GW events could be identified. In contrast to traditional convolutional neural networks, the U-Net algorithm can output time-frequency domain signal images corresponding to probabilities, providing a more intuitive analysis. In conclusion, the U-Net algorithm can rapidly identify the time-frequency domain GW signals from BBH mergers.

Antineutrino Opacity in Neutron Stars in Models Constrained by Recent Terrestrial Experiments and Astrophysical Observations

Inthe present paper, we investigate neutral current (NC) antineutrino scattering with the constituents of neutron star (NS) matter at zero temperature. The modeling of standard matter in NS is constructed within the framework of both extended relativistic mean-field (E-RMF) and nonrelativistic Korea-IBS-Daegu-SKKU energy density functional (KIDS-EDF) models. In the E-RMF model, we use a new parameter, G3(M), which was constrained by the recent PREX II experiment measurement of neutron distribution in 208Pb, while the KIDS-EDF models are constrained by terrestrial experiments, gravitational-wave signals, and astrophysical observations. Using both realistic and well-constrained matter models, we then calculate the antineutrino differential cross-section (ADCS) and antineutrino mean free path (AMFP) for the interaction between antineutrinos and neutron star (NS) matter constituents using linear response theory. It is found that the AMFP for the KIDS0 and KIDSA models are smaller compared to the SLy4 model and the E-RMF model with the G3(M) parameter. The AMFP result of the Skyrme model with the SLy4 parameter set is found to have a prediction almost similar to that of the E-RMF model with the G3(M) parameter. Contributions of each nucleon to the total AMFP are also presented for the G3(M) model.

Optimized search for a binary black hole merger population in LIGO-Virgo O3 data

Maximizing the number of detections in matched filter searches for compact binary coalescence (CBC) gravitational wave signals requires a model of the source population distribution. In previous searches using the framework, sensitivity to the population of binary black hole (BBH) mergers was improved by restricting the range of filter template mass ratios and use of a simple one-dimensional population model. However, this approach does not make use of our full knowledge of the population and cannot be extended to a full parameter space search. Here, we introduce a new ranking method, based on kernel density estimation with adaptive bandwidth, to accurately model the probability distributions of binary source parameters over a template bank, both for signals and for noise events. We demonstrate this ranking method by conducting a search over LIGO-Virgo O3 data for BBHs with unrestricted mass ratio, using a signal model derived from previous significant detected events. We achieve over 10% increase in sensitive volume for a simple power-law simulated signal population, compared to the previous BBH search. Correspondingly, with the new ranking, eight additional candidate events above an inverse false alarm rate threshold 0.5 yr are identified. Published by the American Physical Society 2024

The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars

Abstract Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio frequency. However, pulse arrival times are also affected by radio-frequency-dependent “chromatic” noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect the achromatic noise properties in each pulsar. The models we compare include existing DM models used by the North American Nanohertz Observatory for Gravitational waves (NANOGrav) and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar data set. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from log 10 A RN = − 14.1 − 0.1 + 0.1 to − 14.7 − 0.5 + 0.3 , and the spectral index broadens from γ RN = 2.6 − 0.4 + 0.5 to γ RN = 3.5 − 0.9 + 1.2 . We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs.

A coordinate-independent formalism for detecting high-frequency gravitational waves

In an external electric or magnetic field, a gravitational wave (GW) may be converted into electromagnetic radiation. We present a coordinate-invariant framework to describe the GW signal in a detector that is based on this effect, such as cavities for axion searches. In this framework, we pay special attention to the definition of manifestly coordinate-independent expressions for the electromagnetic fields that an external observer would detect. A careful assessment of the detector’s perceived motion allows us to treat both its mechanical and its electromagnetic response to the GW consistently. We further introduce well-defined approximations for which this motion may be neglected, and hence provide suggestions on which coordinate frame is suitable to characterise the GW signal in practice. We illustrate our findings in two examples, an infinitesimally thin rod and a spherical electromagnetic cavity.

Examining the effects of dark matter spikes on eccentric intermediate-mass ratio inspirals using N-body simulations

ABSTRACT Recent studies suggest that dark matter (DM) spikes around intermediate-mass black holes could cause observable dephasing in gravitational wave (GW) signals from intermediate-mass ratio inspirals (IMRIs). Previous research primarily used non-self-consistent analytic methods to estimate the impact of DM spikes on eccentric IMRIs. Our study provides the first self-consistent treatment of this phenomenon using N-body simulations, incorporating post-Newtonian effects up to the 2.5 order for accurate and robust results. Contrary to prior works, which posited that the cumulative effect of two-body encounters (dynamical friction; DF) is the primary mechanism for energy dissipation, we reveal that a three-body effect (slingshot mechanism) plays a more significant role in driving the binary system’s energy loss and consequent orbital shrinkage. We find that binaries counter-rotating with respect to the DM spike merge faster, while corotating binaries merge slower, contrary to expectations from the DF theory. Using Fokker–Planck methods, we also assess the presence and detectability of spikes in realistic environments. When interacting with surrounding materials, DM spikes can have shallower slopes and lower densities than previously considered, leading to smaller signals and lower detection prospects via dephasing. Our results suggest that ‘deshifting’ rather than dephasing might be a more optimistic signature, as it is more robust even in low-density environments.

Dissipated correction map method with trapezoidal rule for the simulations of gravitational waves from spinning compact binary

ABSTRACT The correction map method means extended phase-space algorithm with correction map. In our research, we have developed a correction map method, specifically the dissipated correction map method with trapezoidal rule, for numerical simulations of gravitational waves from spinning compact binary systems. This new correction map method, denoted as $CM3$, has shown remarkable performance in various simulation results, such as phase-space distance, dissipated energy error, and gravitational waveform, closely resembling the high-order precision implicit Gaussian algorithm. When compared with the previously used mid-point map which is denoted as $C_2$, the $CM3$ consistently exhibits a closer alignment with the highly accurate Gaussian algorithm in waveform evolution and orbital trajectory analysis. Through detailed comparisons and analyses, it is evident that $CM3$ outperforms other algorithms, including $CM2$ and $C_2$ mentioned in this paper, in terms of accuracy and precision in simulating spinning compact binary systems. The incorporation of the trapezoidal rule and the optimization with a scale factor $\gamma$ have significantly enhanced the performance of $CM3$, making it a promising method for future numerical simulations in astrophysics. With the groundbreaking detection of gravitational waves by the LIGO/VIRGO collaboration, interest in this research domain has soared. Our work contributes valuable insights for the application of matched filtering techniques in the analysis of gravitational wave signals, enhancing the precision and reliability of these detection.

Cosmological dynamics of string theory axion strings

The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu-Goldstone boson of global U(1)PQ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we show that string theory axions behave fundamentally differently from field theory axions in the early Universe. Field theory axions may form axion strings if the PQ phase transition takes place after inflation. In contrast, we show that string theory axions do not generically form axion strings. In special inflationary paradigms, such as D-brane inflation, string theory axion strings may form; however, their tension is parametrically larger than that of field theory axion strings. We then show that such QCD axion strings overproduce the DM abundance for all allowed QCD axion masses and are thus ruled out, except in scenarios with large warping. A loop-hole to this conclusion arises in the axiverse, where an axion string could be composed of multiple different axion mass eigenstates; a heavier eigenstate could collapse the network earlier, allowing for the QCD axion to produce the correct DM abundance and also generating observable gravitational wave signals. Published by the American Physical Society 2024

Dark matter from phase transition generated PBH evaporation with gravitational waves signatures

We study the possibility of generating dark matter (DM) purely from ultralight primordial black hole (PBH) evaporation with the latter being produced from a first order phase transition (FOPT) in the early Universe. If such ultralight PBH leads to an early matter domination, it can give rise to a doubly peaked gravitational wave (GW) spectrum in Hz-kHz ballpark with the low frequency peak generated from PBH density fluctuations being within near future experimental sensitivity. In the subdominant PBH regime, the FOPT generated GW spectrum comes within sensitivity due to absence of entropy dilution. In both the regimes, PBH mass from a few kg can be probed by GW experiments like BBO, ET, CE, UDECIGO etc. while DM mass gets restricted to the superheavy ballpark in the PBH dominance case. Apart from distinct DM mass ranges in the two scenarios, GW observations can differentiate by measuring their distinct spectral shapes. Published by the American Physical Society 2024

Search for gravitational waves using a network of RF cavities

The concept of detecting gravitational waves using RF-cavities in strong magnetic fields has gained considerable interest, using setups currently running axion searches. We propose a novel analysis approach for detecting GHz-regime gravitational waves, potentially from primordial black hole mergers, through synchronous measurements from multiple, distant cavities. While individual cavities may detect gravitational wave signatures, isolating them from noise is challenging since the strain (and thereby power deposition in a cavity) of such mergers is expected to be very small and short-lived. By correlating signals from several geographically separated cavities, we can significantly improve the signal-to-noise ratio and potentially investigate the sources of these waves. A demonstration experiment with a superconducting cavity is currently underway, forming the basis for our data analysis methods and outlining the prospects for the GravNet project. This project should be seen as an effort to bring the axion community together and collaborate in the context of gravitational wave physics.

Gravitational axiverse spectroscopy: seeing the forest for the axions

We consider inflationary models with multiple spectator axions coupled to dark gauge sectors via Chern-Simons (CS) terms. The energy injection into Abelian gauge fields from the axions engenders a multi-peak profile for scalar and tensor spectra. We highlight the constraining power of CMB spectral distortions on the scalar signal and discuss the conditions under which spectator sectors can account for the recently observed stochastic gravitational wave (GW) background in the nHz range. Given the tantalizing prospect of a multi-peak “GW forest” spanning several decades in frequency, we elaborate on possible ultraviolet origins of the spectator models from Type IIB orientifolds. String compactifications generically produce a multitude of axions, the “Axiverse”, from dimensional reduction of p-form gauge fields. The CS coupling of such axions to dark gauge fields in the worldvolume theory of D7-branes can be tuned via multiple brane wrappings and/or quantized gauge field strengths. If string axions coupled to Abelian gauge fields undergo slow-roll during inflation, they produce GW signals with peaked frequency distribution whose magnitude depends on the details of the compactification. We discuss the restrictions on spectator models from consistency and control requirements of the string compactification and thereby motivate models that may live in the string landscape as opposed to the swampland.

Spectral Variance in a Stochastic Gravitational-wave Background from a Binary Population

A population of compact object binaries emitting gravitational waves that are not individually resolvable will form a stochastic gravitational-wave signal. While the expected spectrum over population realizations is well known from Phinney, its higher-order moments have not been fully studied before or computed in the case of arbitrary binary evolution. We calculate analytic scaling relationships as a function of gravitational-wave frequency for the statistical variance, skewness, and kurtosis of a stochastic gravitational-wave signal over population realizations due to finite source effects. If the time derivative of the binary orbital frequency can be expressed as a power law in frequency, we find that these moment quantities also take the form of power-law relationships. We also develop a numerical population synthesis framework against which we compare our analytic results, finding excellent agreement. These new scaling relationships provide physical context to understanding spectral fluctuations in a gravitational-wave background signal and may provide additional information that can aid in explaining the origin of the nanohertz-frequency signal observed by pulsar timing array campaigns.

Gravitational waves from primordial black hole evaporation with large extra dimensions

The spectra of gravitational waves from black hole evaporation generically peak at frequencies of order the Hawking temperature, making this signal ultra-high frequency for primordial black holes evaporating in the early universe. This motivates us to consider small black holes in theories with large extra dimensions, for which the peak frequency can be lowered substantially, since the true bulk Planck scale M * can be much smaller than the effective M Pl. We study the emission of brane-localized gravitons during the Hawking evaporation of ultra-light primordial black holes in the context of theories with large extra dimensions, with the ultimate goal of computing the contribution to the stochastic gravitational wave background. To accurately model black hole evolution, we compute greybody factors for particles of spin-0, 1/2, 1, and 2 emitted on the brane and in the bulk, presuming the majority of emission proceeds during the Schwarzschild phase. We then compute the power spectrum and present day spectral density parameter for brane-localized gravitons contributing to a gravitational wave signal. We find that for an optimal choice of parameters, the peak frequency plateaus in the sub-MHz regime, within range of planned high-frequency gravitational wave detectors, making this scenario a target for detection once their sensitivity exceeds ΔN eff bounds.

Enhancing GWOPS Capabilities for Coordinated Multi-Telescope Detection of Gravitational Wave Electromagnetic Counterparts

The groundbreaking detection of gravitational waves (GWs) has ushered in a new era of astronomical observation, granting us access to cosmic phenomena that are imperceptible to electromagnetic waves. The inherently weak GW signals coupled with the substantial uncertainties in source localization pose significant challenges to the field of astronomy. In this paper, we introduce innovative strategies to enhance the efficiency of observing electromagnetic counterparts to GW events, thereby unlocking further secrets of the cosmos. We present a novel technique for designing observation targets and establishing priorities, progressing from the epicenter to the periphery within the boundaries of the GW error sky region. This method has significantly reduced the average slewing distance of telescopes by 41% compared to traditional methods, thus enhancing observational efficiency. Additionally, we have developed a collaborative observation strategy for telescope networks, allocating observation targets based on the field-of-view (FOV) sizes of individual telescopes. This ensures comprehensive coverage without redundancy, allowing a network of four telescopes to cover a sky area and accumulate observation probability more than four times that of a single telescope operating independently over an equivalent period. Building upon these strategies, we have significantly upgraded GWOPS, the GW Follow-up Observation Planning System developed by the China-VO team, to provide precise observational planning for large FOV (greater than 1 square degree) telescope networks. The system also features a web-based user interface that presents the GW error sky area and observation planning results in a graphical format, significantly improving user interaction and experience. The research presented herein equips astronomers with a robust toolkit, advancing the efficiency of searching for and studying electromagnetic counterparts to GW events, and heralding new frontiers in the research of astrophysics and cosmology.

Gravitational waves from more attractive dark binaries

The detection of gravitational waves (GWs) has led to a deeper understanding of binaries of ordinary astrophysical objects, including neutron stars and black holes. In this work, we point out that binary systems may also exist in a dark sector with astrophysical-mass macroscopic dark matter. These “dark binaries”, when coupled to an additional attractive long-range dark force, may generate a stochastic gravitational wave background (SGWB) with a characteristic spectrum different from ordinary binaries. We find that the SGWB from planet-mass dark binaries is detectable by space- and ground-based GW observatories. The contribution to the SGWB today is smaller from binaries that merge before recombination than after, avoiding constraints on extra radiation degrees of freedom while potentially leaving a detectable GW signal at high frequencies up to tens of GHz.

Fragmentation in Gravitationally Unstable Collapsar Disks and Subsolar Neutron Star Mergers

Although stable neutron stars (NSs) can in principle exist down to masses M ns ≈ 0.1 M ⊙, standard models of stellar core-collapse predict a robust lower limit M ns ≳ 1.2 M ⊙, roughly commensurate with the Chandrasekhar mass M Ch of the progenitor’s iron core (electron fraction Y e ≈ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Y e < 0.5) for which MCh∝Ye2 is reduced to ≲1 M ⊙. Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly rotating stars (“collapsars”), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form subsolar-mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNSs, their subsequent coalescence offers a channel for precipitating subsolar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ∼Hz frequency Doppler modulation of the ssNS-merger gravitational-wave signals due to the binary’s orbital motion in the disk; (2) at least one additional gravitational-wave event (coincident within ≲hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched with r-process elements from the NS merger(s) embedded within the exploding stellar envelope (“kilonovae inside a supernova”).