2026-07-03 | | Total: 30
We construct a perturbative scalar-tensor solution describing a central inhomogeneity embedded in an evolving cosmological background, with the aim of studying black hole persistence through a nonsingular bounce. Scalar-tensor gravity provides a natural framework for realizing bouncing cosmologies, while the inclusion of a localized inhomogeneity makes the field equations substantially more difficult to solve. We therefore adopt a perturbative scheme, with perturbative parameter $ε$, in which the leading-order equations are solved by a spatially flat bouncing FLRW spacetime sourced by a radiation perfect fluid. At next order, a central inhomogeneity is introduced through a generalized McVittie geometry, with the perturbations encoded in the corresponding first-order metric and scalar-field functions. We first allow an anisotropic fluid with radial and tangential pressures, whose diagonal components solve the diagonal field equations. The field equations are solved as a series expansion up to $\mathcal{O}(η^4)$ near the bounce at $η=0$. The resulting perfect fluid solution contains three arbitrary functions which are constrained by requiring the spacetime to asymptote to FLRW as $r\to\infty$. With suitable initial conditions preserving the parabolic structure of the bounce, the integration constant $d_0$ emerges as the true perturbative parameter: all perturbations vanish as $d_0\to0$. Finally, we find a small evolving horizon, $r_h\sim d_0$, which we interpret as the horizon of the central inhomogeneity. Its persistence through the bounce supports the interpretation of a black hole surviving the cosmological transition, and its evolution is not symmetric about $η=0$.
Recently, a distinct non--quasinormal mode component of black hole binary merger radiation, named the direct wave, has been identified. The frequency and damping time of the direct wave have been associated with properties of the remnant horizon. This has led to direct-wave based analysis of GW250114, including a test of Hawking's area law. However, as we demonstrate here using numerical relativity strain data, the direct wave frequency is not correlated with the horizon frequency or surface gravity, other than an incidental crossing around $χ_f \approx 0.7$ corresponding closely to the remnant spin of GW250114. We show that while the instantaneous frequency of the direct wave is quasi-stable, the damping time shows significant evolution and therefore a single damped sinusoid model, containing a fixed damping time, is not appropriate. We further show that an evolving frequency model based on horizon properties also does not model the direct wave component for systems with large remnant spins. We demonstrate that testing Hawking's area law with a horizon frequency based on the direct wave interpretation will lead to apparent violations of Hawking's area law when no violation actually occurs. Our results therefore indicate that the direct wave is not a reliable probe of the remnant horizon's properties.
The investigation of nonclassical correlations in curved spacetimes offers key insights into the intersection of quantum information theory and gravitational physics. This paper studies two nonlocal correlation measures, non local advantage of quantum coherence (NAQC) and Bell nonlocality (BN) in a $d$-dimensional spherically symmetric Einstein-Gauss-Bonnet (EGB) black hole spacetime. We consider two observers (Alice and Rob) initially sharing a maximally entangled Bell state: Alice freely falls into the black hole (inertial Kruskal frame), while Rob accelerates outside the horizon (non-inertial Schwarzschild-like frame). The Unruh-Hawking effect modifies Rob's field modes, requiring Bogoliubov transformations to relate the two frames. We derive the mixed bipartite density matrix for fermionic fields and analytical expressions for NAQC and BN, which depend on Hawking temperature (itself governed by $α$, $d$, and $r_h$). Our results show both correlations degrade monotonically with increasing Hawking temperature, confirm the NAQC-BN hierarchical relationship persists in EGB spacetime, and highlighting the impact of high curvature corrections on quantum resources.
In this work, we investigate boson star models within the framework of teleparallel gravity with non-minimal coupling, and obtain static, spherically symmetric solutions for both the ground state and excited states. The results indicate that the energy density of the excited-state solutions can become negative. For these solutions, the four commonly used energy conditions are no longer satisfied. In contrast, for all the ground-state solutions we have studied, the energy density remains positive and all four energy conditions are consistently satisfied. Moreover, considering the importance of astrophysical observations, the gravitational-wave signals from Extreme-Mass-Ratio Inspirals (EMRIs) composed of these boson stars are investigated. Our results reveal that the frequency-domain characteristic strain of these waveforms falls within the detectability range of LISA, which can provide potential evidence for distinguishing compact astrophysical objects.
We investigate static and spherically symmetric neutron star solutions endowed with primary scalar hair in a subfamily of Degenerate-Higher-Order-Scalar-Tensor (DHOST) theories of gravity. By solving the modified Tolman-Oppenheimer-Volkoff (TOV) equations, we construct equilibrium configurations for polytropic and realistic equations of state and analyse the impact of the scalar hair on the stellar structure. We examine the resulting metric and scalar field profiles as well as the mass-radius relation, showing deviations from the predictions of General Relativity (GR). Positive scalar charges lead to more compact stars than in GR and, above a critical threshold, to singularities. Observations could therefore put stringent constraints on the parameters characterising the beyond-GR effects in these theories and their potential scalar hair.
We investigate Lyapunov exponents as dynamical probes of black hole phase transitions in dilatonic Reissner--Nordström--AdS black holes within Einstein--Maxwell--dilaton theory. The thermodynamic quantities and the Lyapunov exponent of charged probe particles were analyzed in both the Einstein and string frames, thus providing a direct comparison between the thermodynamic phase structure of the black hole and that captured by the Lyapunov exponent. Thermodynamic quantities, including the Hawking temperature and Wald entropy, remained constant under conformal frame transformations, yielding identical phase structures in the two frames. In contrast, the Lyapunov exponent exhibited non-trivial frame dependence for massive probe particles due to dilaton coupling, while no frame dependence was found in the massless limit. Numerical analysis revealed that the phase structure features captured by the Lyapunov exponent, including characteristic cusp behavior and transition points, were independent of the choice of frame, despite the Lyapunov exponent itself being frame-dependent. Therefore, the Lyapunov exponent exhibited frame-dependent values, while the critical structure it captures remained constant across conformal frames.
We present an exact solution to the Einstein-Maxwell equations that describes a static black hole coexisting with anisotropic matter immersed in an external magnetic field, obtained via the Harrison transformation. Our findings reveal that an increase in the anisotropic matter parameter systematically suppresses the local chaotic behavior, as indicated by a reduction in the Lyapunov exponent. Conversely, variations in the external magnetic field lead to qualitative changes in global chaotic behavior. This is analyzed through Poincaré sections, which demonstrate transitions between regular and chaotic trajectories resulting from the nonlinear gravitational-magnetic interactions. These factors play distinct yet complementary roles in shaping chaotic particle dynamics around black holes. This study would offer a new theoretical framework for exploring non-integrable particle motion within magnetized black hole spacetimes and for probing a black hole at the galactic center, where magnetic fields may arise from plasma effects surrounding astrophysical black holes.
We constructed rotating black holes with equal angular momentum in five dimensional Einstein-Bumblebee gravity with and without cosmological constant. Their thermodynamic properties are examined via two distinct methods: the Wald formalism and the Komar integral. Notably, the conserved charges, including mass, angular momentum, and entropy, computed from these two approaches differ by a constant prefactor that is solely determined by the Bumblebee coupling. Subsequently, we apply the Kerr/CFT correspondence to derive the microscopic entropy of these black holes and find that it precisely reproduces the entropy in Komar-integral version, rather than the Wald entropy.
We study the cosmological dynamics of interacting dark energy and dark matter in Classical Einstein Gravity and Loop Quantum Cosmology. Two dark matter scenarios are considered: superfluid dark matter described by a generalized cubic equation of state and the standard pressureless fluid. The dark energy component is modeled using both a generalized nonlinear equation of state and a constant equation of state. We examine two phenomenological interaction terms, $Q=α\dotρ_m$ and $Q=β\dotρ_d$, which govern the energy transfer between the dark sectors. In classical gravity, the pressureless matter model exhibits stable late-time attractors, whereas the superfluid dark matter model admits only saddle and non-hyperbolic critical points. Extending the analysis to Loop Quantum Cosmology, quantum geometric corrections replace the Big Bang singularity with a nonsingular quantum bounce, and significantly modify the phase-space dynamics. As a result, the stable attractors of the classical pressureless matter model disappear, and all interacting models possess only saddle and non-hyperbolic critical points. These findings highlight the significant influence of both dark matter properties and quantum gravitational effects on the asymptotic evolution of interacting dark-sectors.
We investigate the dynamics of charged test particles and quasi-periodic oscillations around a Reissner--Nordström-like black hole in Kalb--Ramond (KR) gravity in the presence of an external magnetic test field. The KR background introduces a Lorentz-violating parameter $\ell$, which modifies the spacetime geometry, horizon structure, circular orbits, and characteristic frequencies of particle motion. In contrast to the standard Wald-type prescription, the magnetic-field configuration is constructed from the source-free Maxwell equation on the charged KR background, allowing the magnetic profile to be consistently adapted to the modified geometry. We derive the equations of motion, the effective potential, the conditions for circular orbits, and the orbital and radial epicyclic frequencies of charged particles. The results show that the black-hole charge $Q/M$, the KR parameter $\ell$, the specific particle charge $ε$, and the magnetic coupling $β=bM$ jointly affect the innermost stable circular orbit (ISCO) and the quasi-periodic oscillation (QPO) frequencies. We then apply the obtained frequencies to the relativistic precession model, where the upper QPO frequency is identified with the orbital frequency and the lower one with the periastron-precession frequency. Using the observed twin-peak QPO data of GRO J1655--40, XTE J1550--564, and M82 X-1, we perform a Markov chain Monte Carlo analysis to constrain the model parameters. The obtained posterior constraints indicate that the charged KR black-hole model with an external magnetic field can consistently reproduce the observed QPO pairs within the adopted parameter ranges. These findings suggest that QPO observations may serve as a useful phenomenological tool for probing Lorentz-violating black-hole geometries and electromagnetic effects in strong-gravity environments.
Asymptotic zero Arnowitt-Deser-Misner (ADM) mass wormholes, such as the zero-mass traversable Ellis-Bronnikov wormhole, are of great interest for astrophysical applications such as in the galactic microlensing. However, when considered individually, they are unstable to small perturbations. On the other hand, there is a possibility that they can be stable as an interior partner of a traversable glued wormhole obtained by suitably gluing the interior to the observationally relevant massive exterior spacetimes across spherically symmetric thin shells. Although the exterior spacetime has non-zero ADM mass, massless interior partner remains massless sharing the stability of the glued wormhole. The dynamics of the thin-shell then demarcates the stability regions of the glued wormhole that we wish to study here by employing the novel concepts of thin-shell "mass" and of "external force" constraints discovered by Garcia, Lobo and Visser. We shall consider two classes, where the zero ADM mass interior are Kim-Lee wormholes glued to the exterior Schwarzschild vacuum and Reissner-Nordström spacetime respectively. It turns out that the stability regions in both cases are almost similar although the two interior Kim-Lee spacetimes are physically very different, one scalar charged and the other electrically charged. The conditions under which the stability of glued wormholes could be achieved are analyzed in detail.
In the presence of very intense electromagnetic fields, the refractive index of vacuum is modified such that light velocity is less than $c$ and ultrarelativistic charged particles can be faster than light and can induce Cherenkov radiation. We present the comparison of the Cherenkov radiation produced by the Euler-Heisenberg theory for critical and supercritical magnetic fields. We also make the comparison between the Cherenkov and the synchrotron radiation produced by the charged particles.
The Laser Interferometer Space Antenna (LISA) is a space-based gravitational wave observatory that consists of three spacecraft in a near-equilateral triangular formation. The spacecraft orbits are typically assumed to be perfectly known in LISA data analysis studies, but in reality, the orbit determination process introduces uncertainties in the spacecraft positions and velocities. In this work, we investigate how these uncertainties propagate into the LISA detector output and the impact of neglecting the spacecraft velocities. We quantify these errors in the knowledge of the LISA response using mismatches and discuss the implications for gravitational wave data analysis. We find that spacecraft orbit uncertainties impact the LISA response knowledge at high frequencies with worst mismatch below $10^{-7}$. The effect of neglecting the spacecraft velocities is largest at frequencies around $10^{-4}$ Hz with mismatches of order $10^{-4}$. For a galactic binary with frequency $10^{-4}$ Hz and SNR=200 observed for one year, we find that neglecting the spacecraft velocities in the response leads to less than 1-$σ$ biases in the parameter estimates. This work provides the first characterization of how errors in the LISA gravitational wave response propagate from gravitational wave strain through detector output to estimated parameters.
A hypersurface anchored variational extension of General Relativity is formulated in which the Einstein-Hilbert action is supplemented by a diffeomorphism invariant functional supported on an embedded spacelike hypersurface whose embedding is varied independently of the spacetime metric. The resulting Euler-Lagrange system consists of the Einstein equations with a localized distributional contribution together with an anchoring equation determining admissible embeddings. For the admissible class of hypersurface functionals considered here, the bulk field equations retain the standard second order principal structure away from the hypersurface and no additional propagating bulk gravitational degrees of freedom are introduced. Under ellipticity and invertibility assumptions, local persistence and linear stability of anchoring hypersurfaces follow from standard implicit function and elliptic estimates. The anchoring condition is generically inequivalent to local slicing gauge conditions and therefore defines a genuine variational restriction rather than a coordinate choice. In the canonical formulation, the momentum constraints retain their standard form, whereas the Hamiltonian constraint acquires a hypersurface supported term. The corresponding smeared generators close in the weak sense: the Dirac algebra is recovered in the bulk, with deviations confined to localized distributional surface contributions. The construction therefore defines a constrained sector of the classical solution space of General Relativity consisting of spacetimes that admit at least one embedded hypersurface satisfying the anchoring equation. In homogeneous cosmology, the localized term induces matching conditions across the anchoring surface, allowing finite transitions in the expansion rate while preserving standard evolution away from the transition hypersurface.
We present a covariant mechanism in which a smooth change of metric signature, from a Euclidean to a Lorentzian regime, drives a finite interval of accelerated expansion. The transition, encoded by a scalar interpolator along a timelike congruence, occurs on a codimension-one hypersurface where the continued metric is degenerate but curvature invariants remain finite, so the surface is curvature-regular. Using this covariant continuation, we rewrite the Einstein tensor of the continued metric as a localized, interpolator-dependent effective source for the post-transition Lorentzian branch, yielding a purely geometric stress tensor supported near the crossing. In the Lorentzian regime, we derive a model-independent, local criterion for acceleration: inflation persists while the interpolator's slope exceeds a critical value fixed by the extrinsic curvature and the spatial Ricci curvature on the initial hypersurface, and ends when this inequality is first saturated. Standard smooth profiles (tanh, generalized logistic, and power-law/arctan) admit closed-form expressions for the proper-time duration of the accelerated epoch, showing that, for fixed geometric data, the profile shape controls this duration. The construction provides a non-singular, inflaton-free route from a regular Euclidean origin to an early Lorentzian phase of accelerated expansion, in a manner compatible with no--boundary--type boundary conditions.
We investigate the epistemic opacity of computer simulations and machine learning methods in the context of black hole imaging. We argue that there are forms of opacity-including opacity resulting from the use of machine learning-which do not need to affect the reliability of an inference when it is seen as a part of a broader inferential framework. We propose conditions under which that can plausibly be the case, and discuss how opaque methods can be useful in the context of the (next generation) Event Horizon Telescope. However, we also argue that at least one problematic form of opacity is currently present in black hole imaging: GRMHD models of Sagittarius A* are opaque. This form of opacity signals the limitations of current understanding of the models of this source, and constrains the potential uses of ML models in future observations.
In our recent work, arXiv:2604.11474 [hep-th], we have shown that effective detection of gravitons is possible using an array of charged harmonic oscillators in a dynamical electromagnetic field. Using the interaction Hamiltonian of the identical model, we find out that a systematic way of population inversion of the gravitons is possible in ultra-cold atomic systems. We find out that the exponential growth depends strictly on the number of bosons in the system as well as their inherent squeezing of the matter wave packets. A coherent source of gravitons may lead directly to an unavoidable evidence on the existence of gravitons and based on this analysis we propose an experimental proposal for generating true graviton laser.
Binary (and trinary) radio pulsars are natural laboratories in space for understanding gravity in the strong field regime, with many unique and precise tests carried out so far, including the most precise tests of the strong equivalence principle and of the radiative properties of gravity. The Square Kilometre Array (SKA) telescope, with its high sensitivity in the Southern Hemisphere, will vastly improve the timing precision of recycled pulsars, allowing for a deeper search of potential deviations from general relativity (GR) in currently known systems. A Galactic census of pulsars will, in addition, will yield the discovery of dozens of relativistic pulsar systems, including potentially pulsar -- black hole binaries, which can be used to test the cosmic censorship hypothesis and the ``no-hair'' theorem. Aspects of gravitation to be explored include tests of strong equivalence principles, gravitational dipole radiation, extra field components of gravitation, gravitomagnetism, and spacetime symmetries. In this chapter, we describe the kinds of gravity tests possible with binary pulsar and outline the features and abilities that SKA must possess to best contribute to this science.
Typicality arguments for emergent spacetime rely on the Lubkin-Page bounds, which show that generic quantum states have vanishing correlations between subsystems. These bounds assume a tensor-product Hilbert space (a Type~I von~Neumann algebra), but the observable algebras in quantum field theory and quantum gravity are generically Type~II or Type~III, raising the question of whether the bounds survive. We prove that they do for all Type~II von~Neumann factors. For the hyperfinite Type~II$_1$ factor with a tripartite decomposition $R \cong A \otimes B \otimes E$, the mutual information between subsystems $A$ and $B$ vanishes as $O((d_A d_B / d_E)^2)$ in finite-dimensional approximations, provided $d_A d_B \leq d_E$ (Theorem~1). For Type~II$_\infty$ factors, including the gravitational algebras constructed via the crossed-product method by Witten and by Chandrasekaran, Longo, Penington, and Witten, the bound acquires an additional exponential suppression controlled by the Bekenstein-Hawking entropy (Theorem~2). We identify the obstructions to extending the result to Type~III factors and discuss the open question of whether the commutant of the observable algebra can serve as a natural thermal bath that tightens the bound further.
We develop a semi-analytic framework for stationary, axisymmetric GRMHD jets that efficiently generates resolved polarized images from the near-horizon region out to $\sim 10^5\,r_g$ across a broad parameter space, enabling rapid exploration of how gravity and magnetohydrodynamic flows imprint scale-dependent signatures on jet morphology and polarization. We identify a new scale-dependent separation in polarimetric diagnostics. Outside the photon ring, plasma loading strongly modifies the polarization-angle profile of the integrated jet-layer emission through inertia-driven winding of the magnetic field. At large image-plane radii, the polarization angle follows a power-law in radius, with an index determined by the jet collimation profile. Near the horizon, in contrast, jets converge to a universal polarization pattern controlled solely by black hole spin. This convergence is hierarchical: differences in velocity and magnetic-field structure are erased first, whereas collimation-dependent differences persist to smaller radii, thereby allowing these effects to be disentangled. These results establish a largely achromatic polarimetric diagnostic that connects GRMHD jet dynamics to resolved image structure, with direct implications for high-resolution polarimetry and for constraining black hole spin and jet formation.
Formation of primordial black holes (PBHs) requires a significant enhancement of curvature perturbations. This mechanism leaves a twofold gravitational-wave (GW) signature: a \emph{low-frequency} stochastic background of scalar-induced GWs (SIGWs) and a distinct \emph{high-frequency} signal from subsequent PBH binary mergers. We leverage this shared origin to establish a consistent, \emph{model-independent} connection between these two observables for a monochromatic PBH mass function. Using PBH abundance constraints on the primordial curvature power spectrum, we evaluate the stochastic SIGW background for spherical and ellipsoidal collapse models, demonstrating that the ellipsoidal scenario yields a significantly stronger signal. Furthermore, we analyze the GW signal from PBH binary mergers and find a direct correspondence between the SIGW frequency and the innermost stable circular orbit (ISCO) frequency of the binaries. Because GW emission is nearly maximal near the ISCO, we additionally show that the peak of the full merger GW spectrum relates to the ISCO frequency via $f_{\text{peak}} = 1.79 \, f_{\text{ISCO}}$, a relation that is independent of the binary masses. Remarkably, this unified framework connects these distinct GW channels, enabling the same primordial fluctuations to be probed across widely separated frequency bands.
With the advancements in clock timing technology, increasingly smaller time differences can be distinguished. Therefore, it is critical to investigate the fractional frequency shift of clocks at different locations on Earth. In this paper, we study it systematically under the influence of a subtle lunar tidal potential based on a new method. Our calculations in the geocentric Fermi frame show that when two clocks are located at the same latitude, the longitude difference changes the fractional frequency shift between them. A similar phenomenon occurs when there is a difference in latitude between two clocks on the ground at the same longitude. Interestingly, when the Moon's longitude changes, the phase and amplitude of the lunar tidal fractional frequency shift between two clocks with the same longitude difference will change, while the change in the Moon's latitude only affects the amplitude of the fractional frequency shift of these two clocks. Our results provide useful information for the calibration and synchronization of clocks on Earth.
The nonseparability of post-Newtonian (PN) Hamiltonian systems typically necessitates the use of computationally expensive implicit integrators. Recent research overcomes this limitation by embedding the dynamics into a doubled phase space, which enables the development of explicit symplectic methods. However, existing specially designed explicit integrators suffer from order reduction for high-order methods when the time stepsize is small, i.e., $h <\varepsilon^3$. In this paper, we propose a novel extension and splitting approach for the doubled Hamiltonian, under which specially designed explicit symplectic integrators can be constructed. It is shown that the proposed integrators achieve genuine high-order convergence without order reduction and take advantage of the small PN parameter $\varepsilon$. Numerical results from simulations with 2PN spinning binaries demonstrate superior long-term conservation of invariants and significantly higher computational efficiency compared to both implicit methods and existing explicit splitting techniques.
The magneto-rotational instability can generate strong, turbulent substructure within magnetised shear flows. The efficacy of the mechanism as a function of microphysical aspects of the fluid, such as stratification and diffusivity, has been explored extensively. One aspect that has not been studied thus far, however, is whether the instability can also operate in solids. Motivated by the possibility that solid regions within planets or degenerate stars may rotate differentially with respect to liquid or gaseous layers during some phase of their life, we examine the extent to which elasticity suppresses the instability. A simplified, plane-parallel analysis reveals that only in cases where the flow is strongly sheared, such that the magnetic tension that would result from the instability in a liquid exceeds the shear modulus of the elastic cavity, can magnetic growth occur. In the context of dynamical tides in binary neutron-star mergers, this implies that the magnetic field can be amplified in the crust prior to coalescence only if the star boasts a spin frequency of $\gtrsim 300$Hz. If viscous heating weakens the crystalline structure prior to resonance, the required spin frequency is reduced.
In this work, we focus our analysis of Fractional-Dimension Gravity (FDG) on our home galaxy, the Milky Way (MW), by using the latest Gaia DR3 data as well as previous rotation curve (RC) data for this galaxy. FDG is an alternative gravitational model (previously known as Newtonian Fractional-Dimension Gravity - NFDG) which does not require the dark matter (DM) paradigm. The MW is studied here with the methods of FDG and its observed rotation curves are successfully reproduced by using a variable fractional dimension $D\left (R\right)$, following previous studies of several other galaxies which were analyzed with the same methodology. An alternative dimension function $D_{m}\left(R \right)$, based on the mass-dimension field equation, was also used and yielded less accurate fits to the experimental data. In addition, we also considered possible implications of the FDG metric, based on the presence of additional weights, on the structure of Special Relativity (SR) for spacetimes with fractional dimension. One notable outcome of this analysis is the possibility of an effective superluminal motion in galactic regions where the space dimension is $D<3$. Although this result is very speculative, it opens interesting new perspectives for possible interstellar travel in our galaxy.