Researcher profile

Lijing Shao

Lijing Shao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
28works
0followers
11topics
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

28 published item(s)

preprint2026arXiv

An agentic framework for gravitational-wave counterpart association in the multi-messenger era

With the detection of gravitational waves (GWs), multi-messenger astronomy has opened a new window for advancing our understanding of astrophysics, dense matter, gravitation, and cosmology. The GW sources detected to date are from mergers of compact object binaries, which possess the potential to generate detectable electromagnetic (EM) counterparts. Searching for associations between GW signals and their EM counterparts is an essential step toward enabling subsequent multi-messenger studies. In the era of next-generation GW and EM detectors, the rapid increase in the number of events brings not only unprecedented scientific opportunities, but also substantial challenges to the existing data analysis paradigm. To help address these challenges, we develop GW-Eyes, an agentic framework powered by large language models (LLMs). For the first time, GW-Eyes integrates domain-specific tools and autonomously performs counterpart association tasks between GW and candidate EM events. It supports natural language interaction to assist human experts with auxiliary tasks such as catalog management, skymap visualization, and rapid verification. Our framework leverages the complex decision-making capabilities of LLMs and their traceable reasoning processes, offering a new perspective to the multi-messenger astronomy.

preprint2025arXiv

Lightcurve Features of Magnetar-Powered Superluminous Supernovae with Gravitational-Wave Emission and High-Energy Leakage

Superluminous supernovae (SLSNe) are a distinct class of stellar explosions, exhibiting peak luminosities 10-100 times brighter than those of normal SNe. Their extreme luminosities cannot be explained by the radioactive decay of $^{56}\mathrm{Ni}$ and its daughter $^{56}\mathrm{Co}$ alone. Consequently, models invoking newly formed millisecond magnetars have been widely proposed, capable of supplying additional energy through magnetic dipole radiation. For these rapidly rotating magnetars, however, gravitational-wave (GW) emission may also contribute significantly to the spin-down, particularly during their early evolutionary stages. While high-energy photons initially remain trapped within the optically thick ejecta, they will eventually escape as the ejecta becomes transparent during the expansion, thereby influencing the late-time lightcurve. In this work, we adopt an analytical framework to systematically explore the combined effects of GW emission and high-energy leakage on the lightcurve of SLSNe. Compared to scenarios that neglect these processes, we find that for magnetars with initial spin periods of millisecond, the combined influence suppresses early-time luminosities but enhances late-time emission. We further investigate the effects of the neutron-star equation of state to the lightcurve, GW emission efficiency, ejecta mass, and other relevant quantities. Our results highlight the complex interplay between GW-driven spin-down and radiative transport in shaping the observable features of SLSNe, offering new insights into diagnosing the nature of their central engines.

preprint2025arXiv

The double neutron star PSR J1946+2052 I. Masses and tests of general relativity

We conducted high-precision timing of PSR J1946+2052 to determine the masses of the two neutron stars in the system, test general relativity (GR) and assessed the system's potential for future measurement of the moment of inertia of the pulsar. We analysed seven years of timing data from the Arecibo 305-m radio telescope, the Green Bank Telescope (GBT), and the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The data processing accounted for dispersion measure variations and relativistic spin precession-induced profile evolution. We employed both DDFWHE and DDGR binary models to measure the spin parameters, kinematic parameters and orbital parameters. The timing campaign has resulted in the precise measurement of five post-Keplerian parameters, which yield very precise masses for the system and three tests of general relativity. One of these is the second most precise test of the radiative properties of gravity to date: the intrinsic orbital decay, $\dot{P}_{\rm b,int}=-1.8288(16)\times10^{-12}\rm\,s\,s^{-1}$, represents $1.00005(91)$ of the GR prediction, indicating that the theory has passed this stringent test. The other two tests, of the Shapiro delay parameters, have precisions of 6\% and 5\% respectively; this is caused by the moderate orbital inclination of the system, $\sim 74^{\circ}$; the measurements of the Shapiro delay parameters also agree with the GR predictions. Additionally, we analysed the higher-order contributions of $\dotω$, including the Lense-Thirring contribution. Both the second post-Newtonian and the Lense-Thirring contributions are larger than the current uncertainty of $\dotω$ ($δ\dotω=4\times10^{-4}\,\rm deg\,yr^{-1}$), leading to the higher-order correction for the total mass.

preprint2023arXiv

Ultra low-mass and small-radius white dwarfs made of heavy elements

Seven ultra low-mass and small-radius white dwarfs (LSPM J0815+1633, LP 240-30, BD+20 5125B, LP 462-12, WD J1257+5428, 2MASS J13453297+4200437, and SDSS J085557.46+053524.5) have been recently identified with masses ranging from $\sim$0.02 $M_\odot$ to $\sim$0.08 $M_\odot$ and radii from $\sim$ 4270 km to 10670 km. The mass-radius measurements of these white dwarfs pose challenges to traditional white dwarf models assuming they are mostly made of nuclei lighter than $^{56}$Fe. In this work we consider the possibility that those white dwarfs are made of heavier elements. Due to the small charge-to-mass ratios in heavy elements, the electron number density in white dwarf matter is effectively reduced, which reduces the pressure with additional contributions of lattice energy and electron polarization corrections. This consequently leads to white dwarfs with much smaller masses and radii, which coincide with the seven ultra low-mass and small-radius white dwarfs. The corresponding equation of state and matter contents of dense stellar matter with and without reaching the cold-catalyzed ground state are presented, which are obtained using the latest Atomic Mass Evaluation (AME 2020). Further observations are necessary to unveil the actual matter contents in those white dwarfs via, e.g., spectroscopy, asteroseismology, and discoveries of other ultra low-mass and small-radius white dwarfs.

preprint2022arXiv

Closing a spontaneous-scalarization window with binary pulsars

Benefitting from the unequaled precision of the pulsar timing technique, binary pulsars are important testbeds of gravity theories, providing some of the tightest bounds on alternative theories of gravity. One class of well-motivated alternative gravity theories, the scalar-tensor gravity, predict large deviations from general relativity for neutron stars through a nonperturbative phenomenon known as spontaneous scalarization. This effect, which cannot be tested in the Solar System, can now be tightly constrained using the latest results from the timing of a set of 7 binary pulsars (PSRs J0348+0432, J0737$-$3039A, J1012+5307, J1738+0333, J1909$-$3744, J1913+1102, and J2222$-$0137), especially with the updated parameters of PSRs J0737$-$3039A, J1913+1102, and J2222$-$0137. Using new timing results, we constrain the neutron star's effective scalar coupling, which describes how strongly neutron stars couple to the scalar field, to a level of $\left|α_{\rm A}\right| \lesssim 6 \times 10^{-3}$ in a Bayesian analysis. Our analysis is thorough, in the sense that our results apply to all neutron star masses and all reasonable equations of state of dense matters, in the full relevant parameter space. It excludes the possibility of spontaneous scalarization of neutron stars, at least within a class of scalar-tensor gravity theories.

preprint2022arXiv

Cosmology and perturbations in tachyonic massive gravity

As massive gravity and its extensions offer physically well-defined gravitational theories with a nonzero graviton mass, we present a new extension of the de Rham-Gabadadze-Tolley (dRGT)massive gravity, which is tachyonic massive gravity theory. We firstly introduce the new extension of the dRGT massive gravity, constructed by adding a tachyonic term. We then find the cosmological background equations, and present the analysis of self-accelerating solutions. We examine the tensor perturbations to calculate the dispersion relation of gravitational waves (GWs). In a special case, we consider a constant tachyon potential for the tachyon field, and calculate the equations of motion and self-accelerating solutions. Finally, we investigate the background perturbations, which include tensor, vector, and scalar perturbations in this case. We calculate the dispersion relation of GWs in the FLRW cosmology in a tachyonic massive gravity theory. These analyses provide potential inputs to future applications in cosmology and GWs propagation.

preprint2022arXiv

Cosmology of Dirac-Born-Infeld dRGT massive gravity

We introduce the cosmological analysis of the Dirac-Born-Infeld dRGT massive gravity theory which is a new extension of de Rham-Gabadadze-Tolley (dRGT) massive gravity. In this theory, we consider the Dirac-Born-Infeld (DBI) scalar field which is coupled to the graviton field. Moreover, we perform the cosmological background equations, and we demonstrate the self-accelerating background solutions. We show that the theory consists of self-accelerating solutions with an effective cosmological constant. In the following, we exhibit tensor perturbations analyses and achieve the dispersion relation of gravitational waves. We analyze the propagation of gravitational perturbation in the Friedmann-Lemaître-Robertson-Walker cosmology in the DBI dRGT massive gravity. Finally, we present the vector and scalar perturbations to show the stability conditions of the theory.

preprint2022arXiv

Electromagnetic Follow-up Observations of Binary Neutron Star Mergers with Early Warnings from Decihertz Gravitational-wave Observatories

We investigate the prospects of electromagnetic follow-up observations for binary neutron star (BNS) mergers, with the help of early warnings from decihertz gravitational-wave (GW) observatories, B-DECIGO and DO-Optimal. Extending the previous work, we not only give quick assessments of joint short $γ$-ray burst (sGRB) detection rates for different $γ$-ray satellites and BNS population models, but also elaborate on the analyses and results on multi-band kilonova detections for survey telescopes with different limiting magnitudes. During an assumed 4-year mission time for decihertz GW observatories, we find that for the goals of electromagnetic follow-ups, DO-Optimal performs better than B-DECIGO as a whole on the detection rate, and has a larger detectable distance for joint sGRB/kilonova searches. Taking the log-normal population model for BNS mergers and a one-day early-warning time as an example, we discuss the accuracy in localization and timing, as well as the redshift distributions for various synergy observations with electromagnetic facilities and decihertz GW detectors. Based on our analyses, we propose a feasible "wait-for" pattern as a novel detecting mode for future multi-messenger astrophysics.

preprint2022arXiv

Extending the Fisher Information Matrix in Gravitational-wave Data Analysis

The Fisher information matrix (FM) plays an important role in forecasts and inferences in many areas of physics. While giving fast parameter estimation with the Gaussian likelihood approximation in the parameter space, the FM can only give the ellipsoidal posterior contours of parameters and lose the higher-order information beyond Gaussianity. We extend the FM in gravitational-wave (GW) data analysis using the Derivative Approximation for LIkelihoods (DALI), a method to expand the likelihood while keeping it positive definite and normalizable at every order, for more accurate forecasts and inferences. When applied to the two real GW events, GW150914 and GW170817, DALI can reduce the difference between FM approximation and the real posterior by 5 times in the best case. The calculation time of DALI and FM is at the same order of magnitude, while obtaining the real full posterior will take several orders of magnitude longer. Besides more accurate approximations, higher-order correction from DALI provides a fast assessment on the FM analysis and gives suggestions for complex sampling techniques which are computationally intensive. We recommend using the DALI method as an extension to the FM method in GW data analysis to pursue better accuracy while still keeping the speed.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2022arXiv

Realistic Detection and Early Warning of Binary Neutron Stars with Decihertz Gravitational-wave Observatories

We investigated the detection rates and early warning parameters of binary neutron star (BNS) populations with decihertz gravitational-wave observatories in a realistic detecting strategy. Assuming 4 years' operation of B-DECIGO, we based on parameter precision to classify the detectable BNSs into three categories: (a) sources that merge within 1 year, which could be localized with an uncertainty of $ΔΩ\sim 10^{0}$ deg$^2$; (b) sources that merge in 1-4 years, which take up three quarters of the total events and yield the most precise angular resolution with $ΔΩ\sim 10^{-2}$ deg$^2$ and time-of-merger accuracy with $Δt_c\sim 10^{-1}$ s; and (c) sources that do not merge during the 4-yr mission window, which enable possible early warnings, with $ΔΩ\sim 10^{-1}$ deg$^2$ and $Δt_c\sim 10^{0}$ s. Furthermore, we compared the pros and cons of B-DECIGO with the third-generation ground-based detectors, and explored the prospects of detections using 3 other decihertz observatories and 4 BNS population models. In realistic observing scenarios, we found that decihertz detectors could even provide early-warning alerts to a source decades before its merger while their localizations are still as accurate as ground-based facilities. Finally we found a decrease of events when considering the confusion noise, but this could be partially solved by a proper noise subtraction.

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.

preprint2022arXiv

Static spherical vacuum solutions in the bumblebee gravity model

The bumblebee gravity model is a vector-tensor theory of gravitation where the vector field nonminimally couples to the Ricci tensor. By investigating the vacuum field equations with spherical symmetry, we find two families of black-hole (BH) solutions in this model: one has a vanishing radial component of the vector field and the other has a vanishing radial component of the Ricci tensor. When the coupling between the vector field and the Ricci tensor is set to zero, the first family becomes the Reissner-Nordström solution while the second family degenerates to the Schwarzschild solution with the vector field being zero. General numerical solutions in both families are obtained for nonzero coupling between the vector field and the Ricci tensor. Besides BH solutions, we also reveal the existence of solutions that have a nonvanishing $tt$-component of the metric on the supposed event horizon where the $rr$-component of the metric diverges while the curvature scalars are finite. These solutions are not supported by existing observations but present certain properties that are of academic interests. We conclude the study by putting the BH solutions into tests against the Solar-system observations and the images of supermassive BHs.

preprint2022arXiv

The Eccentric and Accelerating Stellar Binary Black Hole Mergers in Galactic Nuclei: Observing in Ground and Space Gravitational Wave Observatories

We study the stellar binary black holes (BBHs) inspiralling/merging in galactic nuclei based on our numerical method GNC. We find that $3-40\%$ of all new born BBHs will finally merge due to various dynamical effects. In a five year&#39;s mission, up to $10^4$, $10^5$, $\sim100$ of BBHs inspiralling/merging in galactic nuclei can be detected with SNR$>8$ in aLIGO, Einstein/DECIGO, TianQin/LISA/TaiJi, respectively. About tens are detectable in both LISA/TaiJi/TianQin and aLIGO. These BBHs have two unique characteristics: (1) Significant eccentricities. $1-3\%$, $2-7\%$, or $30-90\%$ of them is with $e_i>0.1$ when they enter into aLIGO, Einstein, or space observatories, respectively. Such high eccentricities provide a possible explanation for that of GW 190521. Most highly-eccentric BBHs are not detectable in LISA/Tianqin/TaiJi before entering into aLIGO/Einstein as their strain become significant only at $f_{\rm GW}\gtrsim0.1$ Hz. DECIGO become an ideal observatory to detect those events as it can fully cover the rising phase. (2) Up to $2\%$ of BBHs can inspiral/merge at distances $\lesssim10^3 r_{\rm SW}$ from the massive black hole (MBH), with significant accelerations, such that the Doppler phase drift of $\sim10-10^5$ of them can be detectable with SNR$>8$ in space observatories. The energy density of the gravitational wave backgrounds (GWB) contributed by these BBHs deviate from the powerlaw slope of $2/3$ at $f_{\rm GW}\lesssim 1$mHz. The high eccentricity, significant accelerations and different profile of GWB of these sources make them distinguishable, thus interesting for future GW detections and tests of relativities.

preprint2022arXiv

The Photon Ring in M87*

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called &#34;photon rings&#34; -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

preprint2021arXiv

Neutron stars in massive scalar-Gauss-Bonnet gravity: Spherical structure and time-independent perturbations

The class of scalar-tensor theories with the scalar field coupling to the Gauss-Bonnet invariant has drawn great interest since solutions of spontaneous scalarization were found for black holes in these theories. We contribute to the existing literature a detailed study of the spontaneously scalarized neutron stars (NSs) in a typical theory where the coupling function of the scalar field takes the quadratic form and the scalar field is massive. The investigation here includes the spherical solutions of the NSs as well as their perturbative properties, namely the tidal deformability and the moment of inertia, treated in a unified and extendable way under the framework of spherical decomposition. We find that while the mass, the radius, and the moment of inertia of the spontaneously scalarized NSs show very moderate deviations from those of the NSs in general relativity (GR), the tidal deformability exhibits significant differences between the solutions in GR and the solutions of spontaneous scalarization for certain values of the parameters in the scalar-Gauss-Bonnet theory. As a result, the celebrated universal relation between the moment of inertia and the tidal deformability of neutron stars breaks down. With the mass and the tidal deformability of NSs attainable in the gravitational waves from binary NS mergers, the radius measurable using the X-ray satellites, and the moment of inertia accessible via the high-precision pulsar timing techniques, future multi-messenger observations can be contrasted with the theoretical results and provide us necessary information for building up theories beyond GR.

preprint2021arXiv

Signature of Lorentz Violation in Continuous Gravitational-Wave Spectra of Ellipsoidal Neutron Stars

We study effects of Lorentz-invariance violation on the rotation of neutron stars (NSs) in the minimal gravitational Standard-Model Extension framework, and calculate the quadrupole radiation generated by them. Aiming at testing Lorentz invariance with observations of continuous gravitational waves (GWs) from rotating NSs in the future, we compare the GW spectra of a rotating ellipsoidal NS under Lorentz-violating gravity with those of a Lorentz-invariant one. The former are found to possess frequency components higher than the second harmonic, which does not happen for the latter, indicating those higher frequency components to be potential signatures of Lorentz violation in continuous GW spectra of rotating NSs.

preprint2020arXiv

Combined search for anisotropic birefringence in the gravitational-wave transient catalog GWTC-1

The discovery of gravitational waves (GWs) provides an unprecedented arena to test general relativity, including the gravitational Lorentz invariance violation (gLIV). In the propagation of GWs, a generic gLIV leads to anisotropy, dispersion, and birefringence. GW events constrain the anisotropic birefringence particularly well. Kostelecký and Mewes (2016) performed a preliminary analysis for GW150914. We improve their method and extend the analysis systematically to the whole GW transient catalog, GWTC-1. This is the first global analysis of the spacetime anisotropic Lorentzian structure with a catalog of GWs, where multiple events are crucial in breaking the degeneracy among gLIV parameters. With the absence of abnormal propagation, we obtain new limits on 34 coefficients for gLIV in the nonminimal gravity that surpass previous limits by $\sim 10^2$-$10^5$.

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

Multiband Observation of LIGO/Virgo Binary Black Hole Mergers in the Gravitational-wave Transient Catalog GWTC-1

The Advanced LIGO and Virgo detectors opened a new era to study black holes (BHs) in our Universe. A population of stellar-mass binary BHs (BBHs) are discovered to be heavier than previously expected. These heavy BBHs provide us an opportunity to achieve multiband observation with ground-based and space-based gravitational-wave (GW) detectors. In this work, we use BBHs discovered by the LIGO/Virgo Collaboration as indubitable examples, and study in great detail the prospects for multiband observation with GW detectors in the near future. We apply the Fisher matrix to spinning, non-precessing inspiral-merger-ringdown waveforms, while taking the motion of space-based GW detectors fully into account. Our analysis shows that, detectors with decihertz sensitivity are expected to log stellar-mass BBH signals with very large signal-to-noise ratio, and provide accurate parameter estimation, including the sky location and time to coalescence. Furthermore, the combination of multiple detectors will achieve unprecedented measurement of BBH properties. As an explicit example, we present the multiband sensitivity to the generic dipole radiation for BHs, which is vastly important for the equivalence principle in the foundation of gravitation, in particular for those theories that predict curvature-induced scalarization of BHs.

preprint2020arXiv

Neutron Star Structure in the Minimal Gravitational Standard-Model Extension and the Implication to Continuous Gravitational Waves

Tiny violation of Lorentz invariance has been the subject of theoretic study and experimental test for a long time. We use the Standard-Model Extension (SME) framework to investigate the effect of the minimal Lorentz violation on the structure of a neutron star. A set of hydrostatic equations with modifications from Lorentz violation are derived, and then the modifications are isolated and added to the Tolman-Oppenheimer-Volkoff (TOV) equation as the leading-order Lorentz-violation corrections in relativistic systems. A perturbation solution to the leading-order modified TOV equations is found. The quadrupole moments due to the anisotropy in the structure of neutron stars are calculated and used to estimate the quadrupole radiation of a spinning neutron star with the same deformation. The calculation puts forward a new test for Lorentz invariance in the strong-field regime when continuous gravitational waves are observed in the future.

preprint2020arXiv

New Graviton Mass Bound from Binary Pulsars

In Einstein&#39;s general relativity, gravity is mediated by a massless metric field. The extension of general relativity to consistently include a mass for the graviton has profound implications for gravitation and cosmology. Salient features of various massive gravity theories can be captured by Galileon models, the simplest of which is the cubic Galileon. The presence of the Galileon field leads to additional gravitational radiation in binary pulsars where the Vainshtein mechanism is less suppressed than its fifth-force counterpart, which deserves a detailed confrontation with observations. We prudently choose fourteen well-timed binary pulsars, and from their intrinsic orbital decay rates we put a new bound on the graviton mass, $m_g \lesssim 2 \times 10^{-28}\,{\rm eV}/c^2$ at the 95% confidence level, assuming a flat prior on $\ln m_g$. It is equivalent to a bound on the graviton Compton wavelength $λ_g \gtrsim 7 \times 10^{21}\,{\rm m}$. Furthermore, we extensively simulate times of arrival for pulsars in orbit around stellar-mass black holes and the supermassive black hole at the Galactic center, and investigate their prospects in probing the cubic Galileon theory in the near future.

preprint2020arXiv

Tests of conservation laws in post-Newtonian gravity with binary pulsars

General relativity is a fully conservative theory, but there exist other possible metric theories of gravity. We consider non-conservative ones with a parameterized post-Newtonian (PPN) parameter, $ζ_2$. A non-zero $ζ_2$ induces a self-acceleration for the center of mass of an eccentric binary pulsar system, which contributes to the second time derivative of the pulsar spin frequency, $\ddotν$. In our work, using the method in Will (1992), we provide an improved analysis with four well-timed, carefully-chosen binary pulsars. In addition, we extend Will&#39;s method and derive $ζ_2$&#39;s effect on the third time derivative of the spin frequency, $\dddotν$. For PSR B1913+16, the constraint from $\dddotν$ is even tighter than that from $\ddotν$. We combine multiple pulsars with Bayesian inference, and obtain an upper limit, $\left|ζ_{2}\right|<1.3\times10^{-5}$ at 95% confidence level, assuming a flat prior in $\log_{10} \left| ζ_{2}\right|$. It improves the existing bound by a factor of three. Moreover, we propose an analytical timing formalism for $ζ_2$. Our simulated times of arrival with simplified assumptions show binary pulsars&#39; capability in limiting $ζ_{2}$, and useful clues are extracted for real data analysis in future. In particular, we discover that for PSRs B1913+16 and J0737$-$3039A, $\dddotν$ can yield more constraining limits than $\ddotν$.

preprint2020arXiv

The CPT-violating effects on neutrons&#39; gravitational bound state

In this work, the CPT-violating (CPTV) interactions on neutrons&#39; gravitational bound state are studied. With simple analytical solutions, we provide a preliminary investigation on the Lorentz-violation (LV) induced spin precession due to the $\vecσ\cdot\vec{\tilde{b}}(1+gz)$ and $\bar{b}/m_{_I}\vecσ\cdot\hat{\vec{p}}$ couplings, where $\vec{\tilde{b}}$ and $\bar{b}$ represent LV coefficients. The helicity-dependent couplings can induce unusual phase evolutions with position and momentum dependence. As $\vec{\tilde{b}}$ varies with time due to the Earth&#39;s motion, the spin polarization also shows a sidereal time dependence, and it may be enhanced with time for ultra-stable polarized state of neutrons. The inseparability of the spin-momentum coupling of the $\bar{b}$-term can also lead to motional dependent polarization state. With the precisely measured transition frequency between different gravitational bound states, we get a rough bound $|\vec{\tilde{b}}|<3.9\times10^{-3}$GeV for unpolarized neutrons. If the spin-flip transition frequency can reach comparable precision in the future, the bound can be improved to the level of $10^{-24}$GeV. The test of weak equivalence principle with polarized atom may also improve it significantly.

preprint2020arXiv

Triaxially-deformed Freely-precessing Neutron Stars: Continuous electromagnetic and gravitational radiation

The shape of a neutron star (NS) is closely linked to its internal structure and the equation of state of supranuclear matters. A rapidly rotating, asymmetric NS in the Milky Way undergoes free precession, making it a potential source for multimessenger observation. The free precession could manifest in (i) the spectra of continuous gravitational waves (GWs) in the kilohertz band for ground-based GW detectors, and (ii) the timing behavior and pulse-profile characteristics if the NS is monitored as a pulsar with radio and/or X-ray telescopes. We extend previous work and investigate in great detail the free precession of a triaxially deformed NS with analytical and numerical approaches. In particular, its associated continuous GWs and pulse signals are derived. Explicit examples are illustrated for the continuous GWs, as well as timing residuals in both time and frequency domains. These results are ready to be used for future multimessenger observation of triaxially-deformed freely-precessing NSs, in order to extract scientific implication as much as possible.

preprint2019arXiv

Gravitational-wave merging events from the dynamics of stellar mass binary black holes around the massive black hole in a galactic nucleus

We study the dynamical evolution of the stellar mass binary black holes (BBHs) in a galactic nucleus that contains a massive black hole (MBH). For a comprehensive study of their merging events, we consider simultaneously the non-resonant and resonant relaxations of the BBHs, the binary-single encounters of the BBHs with the field stars, the Kozai-Lidov (KL) oscillation and the close encounters between the BBHs and the central MBH, which usually lead to binaries&#39; tidal disruptions. As the BBHs are usually heavier than the background stars, they sink to the center by mass segregation, making the KL oscillation an important effect in merging BBHs. The binary-single encounters can not only lead to softening and ionization of the BBHs, it can also make them hardening, that increases the merging rates significantly. The mergers of BBHs are mainly contributed by galaxies containing MBHs less massive than $10^8 M_\odot$ and the total event rates are likely in orders of $1$--$10$ Gpc$^{-3}$ yr$^{-1}$, depending on the detailed assumptions of the nucleus clusters. About $3-10\%$ of these BBH mergers are with eccentricity $\ge 0.01$ when their gravitational wave oscillating frequencies enter the LIGO band ($10$\,Hz). Our results show that merging the BBHs within galactic nuclei can be an important source of the merging events detected by the Advanced LIGO/Virgo detectors, and they can be distinguished from BBH mergers from the galactic fields and globular clusters when enough events are accumulated.

preprint2019arXiv

Pulsar tests of the gravitational Lorentz violation

Pulsars are precision celestial clocks. When being put in a binary, the ticking conveys the secret of underlying spacetime geometrodynamics. We use pulsars to test if the gravitational interaction possesses a tiny deviation from Einstein&#39;s General Relativity (GR). In the framework of Standard-Model Extension (SME), we systematically search for Lorentz-violating operators cataloged by (a) the minimal couplings of mass dimension 4, (b) the CPT symmetry of mass dimension 5, and (c) the gravitational weak equivalence principle (GWEP) of mass dimension 8. No deviation from GR was found yet.

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.