Source author record

Zhoujian Cao

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

31 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) "golden event" 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 "chirp" 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' 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'. 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.

preprint2016arXiv

Binary Black Hole merger in f(R) theory

In the near future, gravitational wave detection is set to become an important observational tool for astrophysics. It will provide us with an excellent means to distinguish different gravitational theories. In effective form, many gravitational theories can be cast into an f(R) theory. In this article, we study the dynamics and gravitational waveform of an equal-mass binary black hole system in f(R) theory. We reduce the equations of motion in f(R) theory to the Einstein-Klein-Gordon coupled equations. In this form, it is straightforward to modify our existing numerical relativistic codes to simulate binary black hole mergers in f(R) theory. We considered binary black holes surrounded by a shell of scalar field. We solve the initial data numerically using the Olliptic code. The evolution part is calculated using the extended AMSSNCKU code. Both codes were updated and tested to solve the problem of binary black holes in f(R) theory. Our results show that the binary black hole dynamics in f(R) theory is more complex than in general relativity. In particular, the trajectory and merger time are strongly affected. Via the gravitational wave, it is possible to constrain the quadratic part parameter of f(R) theory in the range |a_2|<10^{11} m^2 . In principle, a gravitational wave detector can distinguish between a merger of binary black hole in f(R) theory and the respective merger in general relativity. Moreover, it is possible to use gravitational wave detection to distinguish between f(R) theory and a non self-interacting scalar field model in general relativity.

preprint2016arXiv

Gravitational wave astronomy: the current status

In the centenary year of Einstein's General Theory of Relativity, this paper reviews the current status of gravitational wave astronomy across a spectrum which stretches from attohertz to kilohertz frequencies. Sect. 1 of this paper reviews the historical development of gravitational wave astronomy from Einstein's first prediction to our current understanding the spectrum. It is shown that detection of signals in the audio frequency spectrum can be expected very soon, and that a north-south pair of next generation detectors would provide large scientific benefits. Sect. 2 reviews the theory of gravitational waves and the principles of detection using laser interferometry. The state of the art Advanced LIGO detectors are then described. These detectors have a high chance of detecting the first events in the near future. Sect. 3 reviews the KAGRA detector currently under development in Japan, which will be the first laser interferometer detector to use cryogenic test masses. Sect. 4 of this paper reviews gravitational wave detection in the nanohertz frequency band using the technique of pulsar timing. Sect. 5 reviews the status of gravitational wave detection in the attohertz frequency band, detectable in the polarisation of the cosmic microwave background, and discusses the prospects for detection of primordial waves from the big bang. The techniques described in sects. 1-5 have already placed significant limits on the strength of gravitational wave sources. Sects. 6 and 7 review ambitious plans for future space based gravitational wave detectors in the millihertz frequency band. Sect. 6 presents a roadmap for development of space based gravitational wave detectors by China while sect. 7 discusses a key enabling technology for space interferometry known as time delay interferometry.

preprint2016arXiv

Multi-horizon and Critical Behavior in Gravitational Collapse of Massless Scalar

This paper studies the whole process of gravitational collapse and accretion of a massless scalar field in asymptotically flat spacetime. Two kinds of initial configurations are considered. One is the initial data without black hole, the other contains a black hole. Under suitable initial conditions, we find that multi-horizon will appear, which means that the initial black hole formed by gravitational collapse or existing at the beginning will instantly expand and suddenly grow, rather than grows gradually in the accretion process. A new type of critical behavior is found around the instant expansion. The numerical computation shows that the critical exponents are universal.

preprint2015arXiv

Improved fast-rotating black hole evolution simulations with modified Baumgarte-Shapiro-Shibata-Nakamura formulation

Different formulations of Einstein's equations used in numerical relativity can affect not only the stability but also the accuracy of numerical simulations. In the original Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation, the loss of the angular momentum, $J$, is non-negligible in highly spinning single black hole evolutions. This loss also appears, usually right after the merger, in highly spinning binary black hole simulations, The loss of $J$ may be attributed to some unclear numerical dissipation. Reducing unphysical dissipation is expected to result in more stable and accurate evolutions. In the previous work \cite{yhlc12} we proposed several modifications which are able to prevent black hole evolutions from the unphysical dissipation, and the resulting simulations are more stable than in the traditional BSSN formulation. Specifically, these three modifications (M1, M2, and M3) enhance the effects of stability, hyperbolicity, and dissipation of the formulation. We experiment further in this work with these modifications, and demonstrate that these modifications improve the accuracy and also effectively suppress the loss of $J$, particularly in the black hole simulations with the initially large ratio of $J$ and the square of the ADM mass.

preprint2015arXiv

Pseudotensor applied to Numerical Relativity in Calculating Global Quantities

In this work we apply the Landau-Lifshitz pseudotensor flux formalism to the calculation of the total mass and the total angular momentum during the evolution of a binary black hole system. We also compare its performance with the traditional integrations for the global quantities. It shows that the advantage of the pseudotensor flux formalism is the smoothness of the numerical value of the global quantities, especially of the total angular momentum. Although the convergence behavior of the global quantities with the pseudotensor flux method is only comparable with the ones with the traditional method, the smoothness of its numerical value allows using a larger radius for surface integration to obtain more accurate result.

preprint2014arXiv

Descope of the ALIA mission

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

preprint2012arXiv

Compact binary evolutions with the Z4c formulation

Numerical relativity simulations of compact binaries with the Z4c and BSSNOK formulations are compared. The Z4c formulation is advantageous in every case considered. In simulations of non-vacuum spacetimes the constraint violations due to truncation errors are between one and three orders of magnitude lower in the Z4c evolutions. Improvements are also found in the accuracy of the computed gravitational radiation. For equal-mass irrotational binary neutron star evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of four smaller. The quality of the Z4c numerical data is also demonstrated by a remarkably accurate computation of the ADM mass from surface integrals. For equal-mass non-spinning binary puncture black hole evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of two smaller. In the same evolutions we find that away from the punctures the Hamiltonian constraint violation is reduced by between one and two orders of magnitude. Furthermore, the utility of gravitational radiation controlling, constraint preserving boundary conditions for the Z4c formulation is demonstrated. The evolution of spacetimes containing a single compact object confirm earlier results in spherical symmetry. The boundary conditions avoid spurious and non-convergent effects present in high resolution runs with either formulation with a more naive boundary treatment. We conclude that Z4c is preferable to BSSNOK for the numerical solution of the 3+1 Einstein equations with the puncture gauge.

preprint2012arXiv

Modifications for numerical stability of black hole evolution

We experiment with several new modifications for the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of the Einstein field equations and demonstrate how these modifications affect the stability of numerical black hole evolution calculations. With these modifications, we evolve both single non-rotating excised Kerr-Schild black holes and punctured binary black holes, and obtain accurate and stable simulations.

preprint2012arXiv

Perturbational Treatment of the Gravitational Potential Effect on Binary Black Hole Evolution

Binary black hole (BBH) systems are usually located in the gravitational potential well formed by a massive black hole (BH), which is mostly located in the center of a galaxy. In most existing studies, the BBH systems are treated as isolated systems, while the effect of the background is ignored. The validity of the approximation is based on the belief that the background gravitational field from other sources is extremely weak compared with the strong gravitational field produced by the BBH itself during the evolution, and can be neglected in gravitational wave detection. However, it is still interesting to check how valid this approximation is. In this work, instead of simulating the three-BH problem with a fully relativistic treatment, we use a perturbational scheme to investigate the effect of the background gravitational potential on the evolution of a BBH, especially on the waveform of its gravitational radiation. Four scenarios are considered including the head-on collision and the inspiral-to-merger process of a BBH which is either freefalling towards or circularly orbiting around a third large BH. The head-on collision and the circular inspiral are two limits of all possible configurations. The existence of the background gravitational potential changes the arrival time of the gravitational wavefront of a BH, prolongs the wavelength, and increases the gravitational radiation energy. And most interestingly, the background gravitational potential induces the higher-order modes of the gravitational wave of a BBH. These interesting phenomena can be explained by the gravitational redshift effect and the change of eccentricity of a BBH's orbit from the background gravitational potential.

preprint2011arXiv

Constructing EOB dynamics with numerical energy flux for intermediate-mass-ratio inspirals

A new scheme for computing dynamical evolutions and gravitational radiations for intermediate-mass-ratio inspirals (IMRIs) based on an effective one-body (EOB) dynamics plus Teukolsky perturbation theory is built in this paper. In the EOB framework, the dynamics essentially affects the resulted gravitational waveform for binary compact star system. This dynamics includes two parts. One is the conservative part which comes from effective one-body reduction. The other part is the gravitational back reaction which contributes to the shrinking process of the inspiral of binary compact star system. Previous works used analytical waveform to construct this back reaction term. Since the analytical form is based on post-Newtonian expansion, the consistency of this term is always checked by numerical energy flux. Here we directly use numerical energy flux by solving the Teukolsky equation via the frequency-domain method to construct this back reaction term. And the conservative correction to the leading order terms in mass-ratio is included in the deformed-Kerr metric and the EOB Hamiltonian. We try to use this method to simulate not only quasi-circular adiabatic inspiral but also the nonadiabatic plunge phase. For several different spinning black holes, we demonstrate and compare the resulted dynamical evolutions and gravitational waveforms.

preprint2011arXiv

Numerical stability of the Z4c formulation of general relativity

We study numerical stability of different approaches to the discretization of a conformal decomposition of the Z4 formulation of general relativity. We demonstrate that in the linear, constant coefficient regime a novel discretization for tensors is formally numerically stable with a method of lines time-integrator. We then perform a full set of apples with apples tests on the non-linear system, and thus present numerical evidence that both the new and standard discretizations are, in some sense, numerically stable in the non-linear regime. The results of the Z4c numerical tests are compared with those of BSSNOK evolutions. We typically do not employ the Z4c constraint damping scheme and find that in the robust stability and gauge wave tests the Z4c evolutions result in lower constraint violation at the same resolution as the BSSNOK evolutions. In the gauge wave tests we find that the Z4c evolutions maintain the desired convergence factor over many more light-crossing times than the BSSNOK tests. The difference in the remaining tests is marginal.

preprint2010arXiv

Numerical evolution of multiple black holes with accurate initial data

We present numerical evolutions of three equal-mass black holes using the moving puncture approach. We calculate puncture initial data for three black holes solving the constraint equations by means of a high-order multigrid elliptic solver. Using these initial data, we show the results for three black hole evolutions with sixth-order waveform convergence. We compare results obtained with the BAM and AMSS-NCKU codes with previous results. The approximate analytic solution to the Hamiltonian constraint used in previous simulations of three black holes leads to different dynamics and waveforms. We present some numerical experiments showing the evolution of four black holes and the resulting gravitational waveform.