General Relativity and Quantum Cosmology (gr-qc)

Abstract: In this work, we present a non-GR full waveform for general parametrization of axisymmetric black holes by extending our previous PSI model. Our model comprises two main components: an inspiral part obtained by using phenomenological method in frequency-domain and a ringdown part derived from quasinormal modes associated with photon motion. For quantitatively revealing the influence of the deviation from Kerr black holes on the waveforms, we specify our model to the bumpy black holes, which are typical examples of non-GR black holes. The results show that the deviation from the Kerr quadrupole moment could be measured in a high accuracy. The new waveform model can be directly used to test black holes for the LIGO-Virgo-KAGRA observations, the third generation detectors and space-borne interferometers.

Abstract: The activation temperature $T$ in the de Sitter environment is twice larger than the Gibbons-Hawking temperature, related to the cosmological horizon. We consider the activation temperature as the local temperature of the de Sitter vacuum, and construct the local thermodynamics of the de Sitter state. This thermodynamics includes also the gravitational coupling $K$ and the scalar Riemann curvature ${\cal R}$ as the thermodynamically conjugate variables. These variables modify the thermodynamics of the Gibbs-Duhem relation in the de Sitter state. The free energy density is proportional to $-T^2$, which is similar to that in the non-relativistic Fermi liquids and in relativistic matter with equation of state $w=1$. The local entropy is proportional to the local temperature, while the total entropy inside the cosmological horizon is $A/4G$, where $A$ is the area of the horizon. This entropy is usually interpreted as the entropy of the cosmological horizon. We also consider the possible application of the de Sitter thermodynamics to the Schwarzschild-de Sitter black hole and to black and white holes with the de Sitter cores.

Abstract: The hyperbolic encounters of two massive objects are characterized by the emission of a gravitational wave burst, with most of the energy released during the closest approach (near the periastron). The detection of such events, different from the well-known inspiral emission, would be an interesting discovery and provide complementary information to observations of binary mergers of black holes and neutron stars in the observable Universe, shedding light, for instance, on the clustering properties of black holes and providing valuable hints on their formation scenario. Here, we analyse the dynamics of such phenomena in the simplest case where two compact objects follow unbound/hyperbolic orbits. Moreover, we explore the effects of orbital precession on the gravitational wave emission, since the precession encodes certain general relativistic effects between two bodies. We also provide templates for the strain of gravitational waves and the power spectrum for the emission, and analytical expressions for the memory effect associated with such signals.

Abstract: We generelize the Hamiltonian picture of General Relativity coupled to classical matter, known as geometrodynamics, to the case where such matter is described by a Quantum Field Theory in Curved Spacetime, but gravity is still described by a classical metric tensor field over a spatial hypersurface and its associated momentum. Thus, in our approach there is no non-dynamic background structure, apart from the manifold of events, and the gravitational and quantum degrees of freedom have their dynamics inextricably coupled. Given the Hamiltonian nature of the framework, there is no need to search for a consistent quantum stress-energy tensor, but instead we work with the generators of hypersurface deformations over the manifold of quantum states. The construction relies heavily on the differential geometry of a fibration of the set of quantum states over the set of gravitational variables, and the introduction of a notion of quantum connection. The most remarkable physical implications of the construction are norm conservation of the quantum state (even if the total dynamics are non-unitary), the clear identification of the hybrid conserved quantities and the description of a dynamical backreaction of quantum matter on geometry and vice versa, which shall modify the physical properties the gravitational field would have in the absence of backreaction.

Abstract: We investigate the process of the Schwinger effect by axial coupling in the natural single-field inflation model in two parts. First we consider the Schwinger effect when the conformal invariance of Maxwell action should be broken by axial coupling $ I(\phi)F_{\mu\nu}\tilde{F}^{\mu\nu} $ with the inflaton field by identifying the standard horizon scale $ k=aH $ at the very beginning of inflation for additional boundary term and use several values of coupling constant $ \chi_{1} $ and estimate electric and magnetic energy densities and energy density of produced charged particles due to the Schwinger effect.We find that for both coupling functions the energy density of the produced charged particles due to the Schwinger effect is so high and spoils inflaton field.In fact the strong coupling or back-reaction occurs because the energy density of produced charged particles is exceeding of inflaton field.We use two coupling functions to break conformal invariance of maxwell action.The simplest coupling function $ I\left(\phi\right)=\chi_{1}\frac{\phi}{M_{p}} $ and a curvature based coupling function $ I\left(\phi\right)= 12\chi_{1}e^{\left(\sqrt{\frac{2}{3}}\frac{\phi}{M_{p}}\right)}\left[\frac{1}{3M_{p}^{2}}\left(4V\left(\phi\right)\right)+\frac{\sqrt{2}}{\sqrt{3}M_{p}}\left(\frac{dV}{d\phi}\right)\right] $ where $V\left(\phi\right) $ is the potential of natural inflation. In second part , in oder to avoid strong back-reaction problem we identify the horizon scale $ k_{H}=aH|\zeta| , \zeta=\frac{{I}^{\prime}\left(\phi\right)\dot{\phi}}{H} $ in which a given Fourier begins to become tachyonically unstable.The effect of this scale is reducing the value of coupling constant $ \chi_{1} $ and weakening the back-reaction problem but in both cases strong coupling or strong back-reaction exists and the Schwinger effect is impossible.

Abstract: We study the vacuum zero point energy associated to a scalar field with an arbitrary mass and conformal coupling in a dS background. Employing dimensional regularization scheme, we calculate the regularized zero point energy density, pressure and the trace of the energy momentum tensor. It is shown that the classical relation $\langle T \rangle =-4 \langle \rho \rangle$ for the vacuum stress energy tensor receives anomalous quantum correction which depends on the mass and the conformal coupling while the relation $\langle \rho \rangle = - \langle P \rangle$ does hold. We calculate the density contrast associated to the vacuum zero point energy and show that $\delta \rho \sim \langle \rho \rangle$ indicating an inhomogeneous and non-perturbative distribution of the zero point energy. Finally, we calculate the skewness associated to the distribution of the zero point energy and pressure and show that they are highly non-Gaussian.

Abstract: Spinning neutron stars (NSs) will emit continuous gravitational waves (GWs) that carry a wealth of information about the compact object. If such a signal is detected, it will provide us with new insight into the physical properties of the matter under extreme conditions. According to binary population synthesis simulations, future space-based GW detectors, such as LISA and TianQin, can potentially detect some double NSs in tight binaries with orbital periods shorter than 10 minutes. Targeted searches for continuous GWs from the spinning NS in such a binary system identified by LISA/TianQin will be possible with the proposed next-generation ground-based GW observatories, such as Cosmic Explorer and Einstein Telescope. Searching for continuous GWs from such a tight binary system requires highly accurate waveform templates that account for the interaction of the NS with its companion. In this spirit, we derive analytic approximate GWs emitted by a triaxial non-aligned NS in a binary system in which the effects of spin-orbit coupling have been incorporated. The difference with the widely used waveform for the isolated NS is estimated and the parameter estimation accuracy of the signals using Cosmic Explorer is calculated. For a typical tight double NS system with a 6~min orbital period, the angular frequency correction of the spinning NS in this binary due to spin precession is $\sim 10^{-6}~{\rm Hz}$, which is in the same order of magnitude as the angular frequency of orbital precession. The fitting factor between the waveforms with and without spin precession will drop to less than 0.97 after a few days ($\sim 10^5~{\rm s}$). We find that spin-orbit coupling has the potential to improve the accuracy of parameter estimation, especially for the binary inclination angle and spin precession cone opening angle, by up to 3 orders of magnitude.

Abstract: Strong evidence of the existence of the Stochastic Gravitational-Wave Background (SGWB) has been reported by the NANOGrav, PPTA, EPTA and CPTA collaborations. The Bayesian posteriors of the Gravitational-Wave Background (GWB) amplitude and spectrum are compatible with current astrophysical predictions for the GWB from the population of supermassive black hole binaries (SMBHBs). In this paper, we discuss the corrections arising from the extra scalar or vector radiation to the characteristic dimensionless strain in PTA experiments and explore the possibility to detect charges surrounding massive black holes, which could give rise to SGWB with vector or scalar polarizations. The parametrized frequency-dependent characteristic dimensionless strain is used to take a Bayesian analysis and the Bayes factor is also computed for charged and neutral SMBHBs. The Bayesian posterior of GWB tensor amplitude is $\log_{10} A_T=-14.85^{+0.26}_{-0.38}$ and spectral exponent $\alpha=-0.60^{+0.32}_{-0.36}$. The Bayesian posterior for vector or scalar amplitude $A_{V, S}$ is nearly flat and there is nearly no constraint from the current observation data. The Bayesian factor is $0.71$ far less than 100, so the current observation can not support the existence of the charged SMBHB.

Abstract: We demonstrate that the asymptotic oscillatory tails of massive gravitons, present in both massive theories of gravity and effectively in extra-dimensional scenarios, could potentially contribute to gravitational waves with very long wavelengths. However, their impact on recent pulsar timing array observations is expected to be relatively small, predominantly consisting of radiation emitted by black holes in our region of the Milky Way.

Abstract: Next-generation ground-based gravitational-wave observatories such as the Einstein Telescope and Cosmic Explorer will detect $O(10^{5}-10^{6})$ signals from compact binary coalescences every year, the exact number depending on uncertainties in the binary merger rate. Several overlapping signals will be present in band at any given time, generating a confusion noise background. We study how this confusion noise affects constraints on possible deviations from general relativity induced by modified gravity and environmental effects. Confusion noise impacts only the signals that last longer in band. Even for a "golden" GW170817-like signal, the constraints broaden by a factor in the range $[10\%,40\%]$ $([70\%,110\%])$ for the fiducial (highest) value of the local binary neutron star merger rate. Our ability to test general relativity or constrain environmental effects will be limited by systematic errors, and not by confusion noise.

Abstract: The entropy produced by the Unruh radiation is estimated and compared to the entropy of a Schwarzschild black hole. We simulate a spherical system of mass $ M $ by set of Unruh horizons and estimate the total entropy of the outgoing radiation. Dependence on mass and spin of the emitted particles is taken into account. The obtained results can be easily extended to any other intrinsic degrees of freedom of outgoing particles. The ratio of Unruh entropy to the Schwarzschild black hole entropy is derived in exact analytical form. For large black holes this ratio is highly sensitive to quantum numbers of emitted quanta, e.g., spin $ s $, for which it varies from $ 0\% $ for $ s = 0 $ to $ 19\% $ for $ s = 5/2 $.

Abstract: The description of gravity in the form of an embedding theory is based on the hypothesis that our space-time is a four-dimensional surface in a flat ten-dimensional space. The choice of standard Einstein-Hilbert action leads in this case to more general field equations than Einstein's equations. By writing them in the form of Einstein's equations with the contribution of additional fictitious matter, one can try to interpret this matter as dark matter. In order to study the behavior of this fictitious matter near the centers of real galaxies, we discuss an analytical method of obtaining corresponding matter density profiles. This method is based on the consideration of the distribution function of particles over all possible trajectories and allows us to estimate the type (core or cusp) of the emerging density profile.

Abstract: We compute analytically differential invariants for accelerating, rotating and charged black holes with a cosmological constant $\Lambda$. In particular, we compute in closed form novel explicit algebraic expressions for curvature invariants constructed from covariant derivatives of the Riemann and Weyl tensors, such as the Karlhede and the Lake-Abdelqader invariants, for the Kerr-Newman-(anti-)de Sitter and accelerating Kerr-Newman-(anti-)de Sitter black hole spacetimes. We explicitly show that some of the computed curvature invariants are vanishing on the event and Cauchy horizons and/or the ergosurface of the accelerating, charged and rotating black holes with a non-zero cosmological constant. Therefore they can serve as possible detectors of the event horizon and ergosurface for such black hole metrics which belong to the most general type D solution of the Einstein-Maxwell equations with a cosmological constant.

Abstract: Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produce consistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using the IllinoisGRMHD code that is part of an open-source community software for simulation in relativistic astrophysics, the Einstein Toolkit. We verify that numerical results of simulating a single isolated neutron star with IllinoisGRMHD can be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find that IllinoisGRMHD reproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version of IlliinoisGRMHD reproduces these results once we account for bug fixes which has occurred since the original publication

Abstract: Gravitational wave signals from extreme mass ratio inspirals are a key target for space-based gravitational wave detectors. These systems are typically modeled as a distributionally-forced Teukolsky equation, where the smaller black hole is treated as a Dirac delta distribution. Time-domain solvers often use regularization approaches that approximate the Dirac distribution that often introduce small length scales and are a source of systematic error, especially near the smaller black hole. We describe a multi-domain discontinuous Galerkin method for solving the distributionally-forced Teukolsky equation that describes scalar fields evolving on a Kerr spacetime. To handle the Dirac delta, we expand the solution in spherical harmonics and recast the sourced Teukolsky equation as a first-order, one-dimensional symmetric hyperbolic system. This allows us to derive the method's numerical flux to correctly account for the Dirac delta. As a result, our method achieves global spectral accuracy even at the source's location. To connect the near field to future null infinity, we use the hyperboloidal layer method, allowing us to supply outer boundary conditions and providing direct access to the far-field waveform. We document several numerical experiments where we test our method, including convergence tests against exact solutions, energy luminosities for circular orbits, the scheme's superconvergence properties at future null infinity, and the late-time tail behavior of the scalar field. We also compare two systems that arise from different choices of the first-order reduction variables, finding that certain choices are numerically problematic in practice. The methods developed here may be beneficial when computing gravitational self-force effects, where the regularization procedure has been developed for the spherical harmonic modes and high accuracy is needed at the Dirac delta's location.