Researcher profile

Zhoujian Cao

Zhoujian Cao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

18 published item(s)

preprint2026arXiv

Constraining Lorentz and parity violations in gravity with multiband gravitational wave observations

This study evaluates the capability of future multi-band observations of gravitational waves emitted from binary black hole coalescences, utilizing joint third-generation ground-based (CE, ET) and space-based (LISA, Taiji, TianQin) detector networks, to constrain parity and Lorentz symmetry violations in the gravitational sector. We model these effects through a parameterized waveform framework that incorporates a set of parameters that quantify potential deviations from general relativity. The frequency-dependence of their effects is described by power-law indices $β$ (i.e., $β_{\bar ν}$, $β_{\bar μ}$, $β_ν$, and $β_μ$). By analyzing events such as a high-signal noise ratio (SNR) &#34;golden event&#34; like GW250114 and a massive binary system like GW231123 (total mass $190-265 M_\odot$) using two networks of ground- and space-based detectors, we demonstrate that multi-band observations can significantly improve the current constraints on Lorentz and parity violations by several order of magnitude, for both high-frequency ($β> 0$) and low-frequency ($β< 0$) modifications. Our Bayesian analysis reveals that while the exceptional SNR of the GW250114-like event yields superior constraints for high-frequency modifications ($β> 0$), the massive nature of GW231123 provides more stringent limits for low-frequency effects ($β< 0$). This work highlights the critical value of future multi-band gravitational wave astronomy for conducting precision tests of general relativity across diverse binary populations.

preprint2026arXiv

Impact of particle production during inflation on the CMB detection

This work focuses on particle production described by a nonminimally coupled model during inflation. In this model, three parameters determine the characteristic frequency and strength of the induced gravitational waves (GWs). Considering the impact of particle production on inflation, we identify the parameter values that generate the strongest GWs without violating the slow-roll mechanism at the CMB scale. However, even with such extreme parameters, the power spectrum of induced GWs is only about $0.3\%$ of that of vacuum GWs. This contribution remains insignificant when identifying the primary source of the detected CMB B-mode polarization. Furthermore, when our results are integrated with the constraints driven by P+ACT+LB+BK18, the contribution of induced GWs at CMB scales becomes negligible. In contrast, their impact on the scalar spectral index $n_s$ proves significant. For a range of parameter values, the Starobinsky inflation model yields predictions for $n_s$ that are consistent with the measurements obtained from P+ACT+LB+BK18.

preprint2026arXiv

Multiband Gravitational Wave Detection Prospects for M31 UCXB-1 System in Low and Middle Frequency Band

The recent discovery of M31 UCXB-1, the first extragalactic ultracompact X-ray binary (UCXB) with an orbital period of $T_{\rm orb} \sim 465$ s, presents a unique laboratory for studying close binary evolution and an unprecedented target for continuous gravitational wave (GW) searches. Its identification as a strong candidate black hole-white dwarf (BH-WD) system, combined with its exceptionally short period and high X-ray luminosity, suggests it may be one of the most vital low-frequency GW sources in M31. In this \textit{Letter}, we investigate the detectability of its GW signal for future space-borne detectors in multiband GW detection. We find that while its signal-to-noise ratio (S/N) for low-frequency detectors remains marginal for high-confidence detection, middle-frequency detectors such as DECIGO and BBO are far more promising, potentially achieving S/N $\varrho>8$ within reasonable observational duration. With a primary mass of only $m_1 > 5.4M_\odot$ (or $6.6M_\odot$), the network of all low and middle frequency detector (or BBO alone) is capable of detecting GW from this system with a $\varrho > 8$, during 10-year observation. Furthermore, orbital eccentricity can enhance the GW strain at higher harmonics, further improving its detectability, especially for middle-frequency detectors. This study establishes M31 UCXB-1 as a key prototype of short-period UCXBs, cementing its role as a cornerstone for multiband, multi-messenger astrophysics and a vital bridge between X-ray astronomy and the future GW era.

preprint2025arXiv

The transverse-traceless gauge and the gauge problem of second order gravitational waves

The gauge problem arises in the second order gravitational waves due to the mode mixing. Here, we introduce the transverse-traceless (TT) gauge to cosmological backgrounds, and find that if we choose the TT gauge at first order, the second order tensor mode would be gauge invariant. Analogous to the Ricci flat spacetime, the vacuum condition is the key to guarantee the existence of the TT gauge on cosmological backgrounds. When we have the vacuum condition, the Poisson gauge, the uniform curvature gauge, the synchronous gauge and the total matter gauge are all equivalent to the TT gauge. Once the vacuum condition is approximately satisfied, the Poisson gauge would reduce to the TT gauge at the same order of approximation. With the sub-horizon limit, the vacuum condition could be obtained approximately, and the Poisson gauge, the uniform curvature gauge and the synchronous gauge are all approximated TT gauge. Our findings explain several existing results in the literature and indicate that the proposed TT gauge is useful to discuss higher order gravitational waves.

preprint2024arXiv

Detection of the Permanent Strain Offset Component of Gravitational-Wave Memory in Black Hole Mergers

We propose a novel approach to detecting the elusive gravitational-wave memory predicted by general relativity to accompany black hole mergers: direct measurement of the permanent space-time strain offset. Compared to previous techniques modeling and disentangling both the &#34;chirp&#34; and memory signals, this approach has several advantages: it targets the feature of the signal carrying nearly all its Shannon information, has great simplicity, circumvents the need for precise modeling of the time evolution of all components of the gravitational wave signal, and uses only data largely free of the more complicated chirp signal. The frequency spectrum of the predicted memory signal is roughly similar to that of the chirp signal. However its inclusion of lower frequencies, where noise and data calibration are problematic, makes detection difficult but not impossible. We applied this novel analysis, implemented with a template-like algorithm, to a selection of 67 observations of 41 black hole mergers in the LIGO/Virgo Gravitational Wave Transient Catalog. Statistical significance was assessed by analyzing many time-shifted intervals. The result: a few possible detections ($2σ-4σ$) and many upper limits. The probability that a random ensemble of 67 strain time series, with the same noise but no memory signals, will yield a particular figure-of-merit computed for the actual data is approximately 0.1. Several validation checks proved useless, partly due to large measurement and theoretical uncertainties, so these results should be viewed with reservation. Appendices contain MatLab code for various operations, including an algorithm for the complex Fourier transform of arbitrarily spaced data.

preprint2022arXiv

Direct measurement of the distribution of dark matter with strongly lensed gravitational waves

In this Letter, we present a new idea of probing the distribution of dark matter exhibiting elastic and velocity-independent self-interactions. These interactions might be revealed in multiple measurements of strongly lensed gravitational waves, which can be observationally explored to determine the strength of self-scatterings. Specifically, each individual galactic-scale strong-lensing system whose source is a coalescing compact binary emitting gravitational waves will provide a model-independent measurement of the shear viscosity of dark matter along the line of sight. These individual measurements could be a probe of large-scale distribution of dark matter and its properties. Our results indicate that with 10-1000 strongly lensed gravitational waves from ET and DECIGO, robust constraints on the large-scale distribution of self-interacting dark matter might be produced. More stringent limits on the dark matter scattering cross-section per unit mass ($σ_χ/m_χ$) relevant to galaxy and cluster scales are also expected, compared with the conservative estimates obtained in the electromagnetic domain. Finally, we discuss the effectiveness of our method in the context of self-interacting dark matter particle physics.

preprint2022arXiv

Ensemble of Deep Convolutional Neural Networks for real-time gravitational wave signal recognition

With the rapid development of deep learning technology, more and more researchers apply it to gravitational wave (GW) data analysis. Previous studies focused on a single deep learning model. In this paper we design an ensemble algorithm combining a set of convolutional neural networks (CNN) for GW signal recognition. The whole ensemble model consists of two sub-ensemble models. Each sub-ensemble model is also an ensemble model of deep learning. The two sub-ensemble models treat data of Hanford and Livinston detectors respectively. Proper voting scheme is adopted to combine the two sub-ensemble models to form the whole ensemble model. We apply this ensemble model to all reported GW events in the first observation and second observation runs (O1/O2) by LIGO-VIRGO Scientific Collaboration. We find that the ensemble algorithm can clearly identify all binary black hole merger events except GW170818. We also apply the ensemble model to one month (August 2017) data of O2. There is no false trigger happens although only O1 data are used for training. Our test results indicate that the ensemble learning algorithms can be used in real-time GW data analysis.

preprint2022arXiv

Population Properties of Gravitational-Wave Neutron Star--Black Hole Mergers

Over the course of the third observing run of LIGO-Virgo-KAGRA Collaboration, several gravitational-wave (GW) neutron star--black hole (NSBH) candidates have been announced. By assuming these candidates are real signals and of astrophysical origins, we analyze the population properties of the mass and spin distributions for GW NSBH mergers. We find that the primary BH mass distribution of NSBH systems, whose shape is consistent with that inferred from the GW binary BH (BBH) primaries, can be well described as a power-law with an index of $α= 4.8^{+4.5}_{-2.8}$ plus a high-mass Gaussian component peaking at $\sim33^{+14}_{-9}\,M_\odot$. The NS mass spectrum could be shaped as a near flat distribution between $\sim1.0-2.1\,M_\odot$. The constrained NS maximum mass agrees with that inferred from NSs in our Galaxy. If GW190814 and GW200210 are NSBH mergers, the posterior results of the NS maximum mass would be always larger than $\sim2.5\,M_\odot$ and significantly deviate from that inferred in the Galactic NSs. The effective inspiral spin and effective precession spin of GW NSBH mergers are measured to potentially have near-zero distributions. The negligible spins for GW NSBH mergers imply that most events in the universe should be plunging events, which supports the standard isolated formation channel of NSBH binaries. More NSBH mergers to be discovered in the fourth observing run would help to more precisely model the population properties of cosmological NSBH mergers.

preprint2022arXiv

Search for the birefringence of gravitational waves with the third observing run of Advanced LIGO-Virgo

Gravitational waves would attain birefringence during their propagation from distant sources to the Earth, when the CPT symmetry is broken. If it was sizeable enough, such birefringence could be measured by the Advanced LIGO, Virgo and KAGRA detector network. In this work, we place constraints on the birefringence of gravitational waves with the third observing run of this network, i.e. two catalogues GWTC-2 and GWTC-3. For the dispersion relation $ω^{2}=k^{2}\pm2ζk^{3}$, our analysis shows the up-to-date strictest limit on the CPT-violating parameter, i.e. \textcolor{black}{{$ζ=4.07^{+5.91}_{-5.79}\times10^{-17}\mathrm{m}$}}, at $68\%$ confidence level. This limit is stricter by $\sim$5 times than the existing one {($\sim 2\times10^{-16}$m)} and stands for the first $\sim$10GeV-scale test of the CPT symmetry in gravitational waves. The results of Bayes factor strongly disfavor the birefringence scenario of gravitational waves.

preprint2022arXiv

Simultaneous bounds on the gravitational dipole radiation and varying gravitational constant from compact binary inspirals

Compact binaries are an important class of gravitational-wave (GW) sources that can be detected by current and future GW observatories. They provide a testbed for general relativity (GR) in the highly dynamical strong-field regime. Here, we use GWs from inspiraling binary neutron stars and binary black holes to investigate dipolar gravitational radiation (DGR) and varying gravitational constant predicted by some alternative theories to GR, such as the scalar-tensor gravity. Within the parametrized post-Einsteinian framework, we introduce the parametrization of these two effects simultaneously into compact binaries&#39; inspiral waveform and perform the Fisher-information-matrix analysis to estimate their simultaneous bounds. In general, the space-based GW detectors can give a tighter limit than ground-based ones. The tightest constraints can reach $σ_B<3\times10^{-11}$ for the DGR parameter $B$ and $σ_{\dot{G}}/G < 7\times10^{-9} \, {\rm yr}^{-1} $ for the varying $G$, when the time to coalescence of the GW event is close to the lifetime of space-based detectors. In addition, we analyze the correlation between these two effects and highlight the importance of considering both effects in order to arrive at more realistic results.

preprint2021arXiv

A higher-multipole gravitational waveform model for an eccentric binary black holes based on the effective-one-body-numerical-relativity formalism

We construct a new factorized waveform including $(l,|m|)=(2,2),(2,1),(3,3),(4,4)$ modes based on effective-one-body (EOB) formalism, which is valid for spinning binary black holes (BBH) in general equatorial orbit. When combined with the dynamics of $\texttt{SEOBNRv4}$, the $(l,|m|)=(2,2)$ mode waveform generated by this new waveform can fit the original $\texttt{SEOBNRv4}$ waveform very well in the case of a quasi-circular orbit. We have calibrated our new waveform model to the Simulating eXtreme Spacetimes (SXS) catalog. The comparison is done for BBH with total mass in $(20,200)M_\odot$ using Advanced LIGO designed sensitivity. For the quasi-circular cases we have compared our $(2,2)$ mode waveforms to the 281 numerical relativity (NR) simulations of BBH along quasi-circular orbits. All of the matching factors are bigger than 98\%. For the elliptical cases, 24 numerical relativity simulations of BBH along an elliptic orbit are used. For each elliptical BBH system, we compare our modeled gravitational polarizations against the NR results for different combinations of the inclination angle, the initial orbit phase and the source localization in the sky. We use the the minimal matching factor respect to the inclination angle, the initial orbit phase and the source localization to quantify the performance of the higher modes waveform. We found that after introducing the high modes, the minimum of the minimal matching factor among the 24 tested elliptical BBHs increases from 90\% to 98\%. Following our previous $\texttt{SEOBNRE}$ waveform model, we call our new waveform model $\texttt{SEOBNREHM}$. Our $\texttt{SEOBNREHM}$ waveform model can match all tested 305 SXS waveforms better than 98\% including highly spinning ($χ=0.99$) BBH, highly eccentric ($e\approx0.15$) BBH and large mass ratio ($q=10$) BBH.

preprint2021arXiv

Extending the frequency reach of pulsar timing array based gravitational wave search without high cadence observations

Gravitational wave (GW) searches using pulsar timing arrays (PTAs) are assumed to be limited by the typical average observational cadence of $1/(2~{\rm weeks})$ for a single pulsar to GW frequencies $\lesssim 4\times 10^{-7}$ Hz. We show that this assumption is incorrect and that a PTA can detect signals with much higher frequencies, which are preserved in the data due to aliasing, by exploiting asynchronous observations from multiple pulsars. This allows an observation strategy that is scalable to future large-scale PTAs containing $O(10^3)$ pulsars, enabled by the Five-hundred meter Aperture Spherical Telescope and the Square Kilometer Array, without requiring a higher per-pulsar observation cadence. We show that higher frequency GW observations, reaching up to $4\times 10^{-4}$ Hz with an SKA-era PTA, have significant astrophysical implications, such as (i) a three orders of magnitude better constraint than current high-cadence observations on GW strain in the $[10,400]$ $μ{\rm Hz}$ band, and (ii) sensitive tests of the no-hair theorem in the mass range of supermassive black hole binaries using their inspiral, merger, and ringdown signals.

preprint2021arXiv

Probing the large-scale structure of the universe through gravitational-wave observations

The improvements in the sensitivity of the gravitational wave (GW) network enable the detection of several large redshift GW sources by third-generation GW detectors. These advancements provide an independent method to probe the large-scale structure of the universe by using the clustering of the binary black holes. The black hole catalogs are complementary to the galaxy catalogs because of large redshifts of GW events, which may imply that binary black holes (BBHs) are a better choice than galaxies to probe the large-scale structure of the universe and cosmic evolution over a large redshift range. To probe the large-scale structure, we used the sky position of the binary black holes observed by third-generation GW detectors to calculate the angular correlation function (ACF) and the bias factor of the population of binary black holes. This method is also statistically significant as 5000 BBHs are simulated. Moreover, for the third-generation GW detectors, we found that the bias factor can be recovered to within 33$\%$ with an observational time of ten years. This method only depends on the GW source-location posteriors; hence, it can be an independent method to reveal the formation mechanisms and origin of the BBH mergers compared to the electromagnetic method.

preprint2020arXiv

Improved deep learning techniques in gravitational-wave data analysis

In recent years, convolutional neural network (CNN) and other deep learning models have been gradually introduced into the area of gravitational-wave (GW) data processing. Compared with the traditional matched-filtering techniques, CNN has significant advantages in efficiency in GW signal detection tasks. In addition, matched-filtering techniques are based on the template bank of the existing theoretical waveform, which makes it difficult to find GW signals beyond theoretical expectation. In this paper, based on the task of GW detection of binary black holes, we introduce the optimization techniques of deep learning, such as batch normalization and dropout, to CNN models. Detailed studies of model performance are carried out. Through this study, we recommend to use batch normalization and dropout techniques in CNN models in GW signal detection tasks. Furthermore, we investigate the generalization ability of CNN models on different parameter ranges of GW signals. We point out that CNN models are robust to the variation of the parameter range of the GW waveform. This is a major advantage of deep learning models over matched-filtering techniques.

preprint2020arXiv

Physical Implications of the Sub-threshold GRB GBM-190816 and its Associated Sub-threshold Gravitational Wave Event

The LIGO-Virgo and Fermi collaborations recently reported a possible joint detection of a sub-threshold gravitational wave (GW) event and a sub-threshold gamma-ray burst (GRB), GBM-190816, that occurred 1.57 s after the merger. We perform an independent analysis of the publicly available data and investigate the physical implications of this potential association. By carefully studying the following properties of GBM-190816 using Fermi/GBM data, including signal-to-noise ratio, duration, f-parameter, spectral properties, energetic properties, and its compliance with some GRB statistical correlations, we confirm that this event is likely a typical short GRB. Assuming its association with the sub-threshold GW event, the inferred luminosity is $1.47_{-1.04}^{+3.40} \times 10^{49}$ erg s$^{-1}$. Based on the available information of the sub-threshold GW event, we infer the mass ratio q of the compact binary as $q=2.26_{-1.43}^{+2.75}$ according to the reported range of luminosity distance. If the heavier compact object has a mass > 3 solar masses, q can be further constrained to $q=2.26_{-0.12}^{+2.75}$. The leading physical scenario invokes an NS-BH merger system with the NS tidally disrupted. Within this scenario, we constrain the physical properties of such a system to produce a GRB. The GW data may also allow an NS-BH system with no tidal disruption of the NS or a BH-BH merger. We apply the charged compact binary coalescence (cCBC) theory (for both a constant charge and an increasing charge for the merging members) to derive the model parameters to account for GBM-190816 and found that the required parameters are extreme. Finally, we argue that the fact that the observed GW-GRB delay time scale is comparable to that of GW170817/GRB 170817A suggests that the GW-GRB time delay of these two cases is mainly defined by the time scale for the jet to propagate to the GRB emission site.

preprint2019arXiv

Detecting the Beaming Effect of Gravitational Waves

The models currently used in the detection of gravitational waves (GWs) either do not consider a relative motion between the center-of-mass of the source and the observer, or usually only consider its effect on the frequencies of GWs. However, it is known for light waves that a relative motion not only changes the frequencies but also the brightness of the source, the latter of which is called the `beaming effect&#39;. Here we investigate such an effect for GWs and find that the observed amplitude of a GW signal, unlike the behavior of light, is not a monotonic function of the relative velocity and responds differently to the two GW polarizations. We attribute the difference to a rotation of the wave-vector, as well as a reorientation of the GW polarizations. We find that even for velocities as small as $0.25\%$ of the speed of light, ignoring the aforementioned beaming effect could induce a systematic error that is larger than the designated calibration accuracy of LIGO. This error could lead to an incorrect estimation of the distance and orbital inclination of a GW source, or result in a spurious signal that appears to be incompatible with general relativity.

preprint2019arXiv

Gravitational wave signal recognition of O1 data by deep learning

Deep learning method develops very fast as a tool for data analysis these years. Such a technique is quite promising to treat gravitational wave detection data. There are many works already in the literature which used deep learning technique to process simulated gravitational wave data. In this paper we apply deep learning to LIGO O1 data. In order to improve the weak signal recognition we adjust the convolutional neural network (CNN) a little bit. Our adjusted convolutional neural network admits comparable accuracy and efficiency of signal recognition as other deep learning works published in the literature. Based on our adjusted CNN, we can clearly recognize the eleven confirmed gravitational wave events included in O1 and O2. And more we find about 2000 gravitational wave triggers in O1 data.

preprint2019arXiv

Validating the Effective-One-Body Numerical-Relativity Waveform Models for Spin-aligned Binary Black Holes along Eccentric Orbits

Effective-one-body (EOB) numerical-relativity (NR) waveform models for spin-aligned binary black holes (BBHs), known as the SEOBNR waveform models, are based on the EOB theoretical framework and NR simulations. SEOBNR models have played an important role in the LIGO scientific collaboration (LSC) gravitational wave (GW) data analysis for both signal search and parameter estimation. SEOBNR models for quasi-circular orbits have evolved through version 1 to version 4 by extending their validity domain and including more NR results. Along another direction, we recently extended SEOBNRv1 model to SEOBNRE model which is valid for spin-aligned BBH coalescence along eccentric orbits. In this paper we validate this theoretical waveform model by comparing them against the numerical relativity simulation bank, Simulating eXtreme Spacetimes (SXS) catalog. In total, 278 NR waveforms are investigated which include binaries with large eccentricity; large spin and large mass ratio. Our SEOBNRE can model the NR waveforms quite well. The fitting factor for most of the 278 waveforms is larger than 99\%. It indicates that the SEOBNRE model could be used as template waveforms for eccentric spin-aligned BBH coalescence.