Gravitational waves from merging compact binaries are sensitive to line-of-sight acceleration (LOSA) induced by a massive companion in their vicinity. Interestingly, the leading-order contributions of LOSA and residual orbital eccentricity to the Fourier phase of the inspiral waveform exhibit similar frequency dependence, raising the possibility that a small eccentricity could mimic LOSA effects in transient GW events such as GW190814. We perform Bayesian inference using the IMRPhenomXPHM waveform family as the baseline LIGO-Virgo-KAGRA waveform model, augmented with leading-order LOSA and residual eccentricity corrections while using 32 seconds of data associated with GW190814. For a LOSA-only analysis, we find no evidence for a non-zero LOSA effect in GW190814, with a Bayes factor relative to the baseline model of approximately 0.22, consistent with the findings of Hendriks et al. and in tension with the claim by Yang et al., who employed only 4 seconds of GW190814 data. In a joint analysis that includes both leading-order LOSA and eccentricity effects, we obtain informative posteriors for both parameters, with representative values a/c approximately -2.8 x 10^{-3} s^{-1} and e_0 approximately 0.11. However, the corresponding Bayes factor relative to the baseline model is approximately 0.64, suggesting that the 32-second data do not provide significant evidence for either LOSA or residual eccentricity in GW190814. Further, our Bayesian runs reveal a strong correlation between the LOSA and eccentricity parameters, indicating a significant degeneracy in their imprint on the inspiral phase. This finding is consistent with our theoretical arguments and most likely explains the non-zero parameter estimates obtained in the joint analysis. less

We construct charged wormhole solutions with an even number of asymptotically flat regions in the four-dimensional Einstein-Maxwell-massless phantom scalar system via the Harrison transformation. The solutions are characterized by five parameters: the mass $M$, the electric charge $Q_\mathrm{e}$, the magnetic charge $Q_\mathrm{m}$, the scalar charge $P$ and the number of sheets $2n$. The regularity condition then determines the throat radius. Although the Harrison transformation directly generates the solutions only in the parameter region $Q_{\mathrm{e}}^2 + Q_{\mathrm{m}}^2 < M^2$, we show that regular solutions exist in a wider parameter region beyond this bound. In addition, we introduce a spheroidal coordinate system that covers one complete asymptotically flat region and its adjacent ones, and allows the solution to be expressed in a simple form. less

Ghost-free non-local gravity is investigated with regards to its late-time dynamics. Viable solutions in this model are confronted with the observational data including the Pantheon+ catalogue of Type Ia supernovae, the Dark Energy Spectroscopic Instrument, the measurements of baryon acoustic oscillations and the Hubble parameter estimations $H(z)$. The ghost-free non-local gravity is found to be successful in these tests in comparison to the $Λ$CDM model and can be also comparable with the generalized exponential $F(R)$ gravity scenario. However the model encounters difficulties when the data from the above observations and the cosmic microwave background radiation data are combined. In tests with the whole set of Pantheon+, DESI, $H(z)$ and CMB data, the generalized exponential $F(R)$ model is essentially more successful. This success is related with the dynamical behavior of its effective dark energy equation of state evolving from a phantom to a quintessence phase during the late-time epoch, whereas the ghost-free non-local model demonstrates only a quintessence behavior. Hence the ghost-free non-local gravity scenario is successful only when the Pantheon+, DESI and $H(z)$ data are considered. The generalized exponential $F(R)$ model satisfies the viability conditions and in tests with all observational data including CMB surpasses the $Λ$CDM model in $χ^2$ statistics and also with information criteria. less

We investigate photon-axion conversion in the vicinity of rotating Kerr black holes where strong gravity traps photons on near-circular trajectories, effectively enhancing the path length. We explore the observable signatures of such a conversion near the photon region. The process, driven by ambient magnetic fields, is significantly more efficient around supermassive black holes such as M87*, since the luminosity of photons increases with the mass of the BH. By numerically evaluating photon path lengths (on which the conversion depends), we analyze how key parameters-photon frequency, axion mass, photon-axion coupling, magnetic field strength, plasma density, and black hole spin-affect the conversion probability and the resultant dimming of photon spectral luminosity. We find that the conversion is most efficient at high frequencies (X-rays and gamma rays), while the frequency window associated with efficient conversion widens with an increase in the photon-axion coupling and a decrease in the electron density and the axion mass. The magnitude of dimming of the photon spectral luminosity depends primarily on the magnetic field, the photon-axion coupling and the BH spin. Our study reveals that rotating black holes generally exhibit enhanced dimming compared to static ones. Thus, if future telescopes achieving a resolution $\sim 10^{-5}$ arcsec in the X-ray/gamma-ray band detect a dimming of the photon spectral luminosity, then they can provide interesting constraints on the axion mass and its coupling with photons. less

We study the dynamics of a binary orbiting a supermassive black hole (SMBH), focusing on both binary scattering in unbound orbits and eccentric von Zeipel-Lidov-Kozai (vZLK) oscillations in bound orbits. The motion is described in a local inertial frame in Kerr spacetime, where tidal effects are encoded in the Riemann curvature. For unbound (parabolic and hyperbolic) orbits, we identify four scattering regimes-adiabatic, tidally affected, chaotic, and disruptive-depending on the binary semi-major axis. As the binary becomes softer, tidal interactions near periapsis lead to strong eccentricity excitation, large changes in the orbital parameters, and eventually chaotic behavior or tidal disruption, with a sensitive dependence on the argument of periapsis. For eccentric bound (elliptic) orbits, the vZLK mechanism differs qualitatively from the standard one, although the $z$-component of the angular momentum in the local inertial frame remains approximately conserved. The evolution proceeds on a dynamical timescale and exhibits step-like changes driven by repeated periapsis passages, which can be interpreted as a sequence of scattering events. We refer to this behavior as scattering-type vZLK oscillations. The rotation of the SMBH also modifies the oscillation profiles, although its effect is less significant than the dependence on the initial orbital parameters. These results suggest a unified picture of periapsis-driven tidal dynamics in galactic nuclei. less

Building on Padmanabhan's entropy functional, originally introduced to derive Einstein's equations and highlight the emergent nature of gravity, we demonstrate its robustness in a broader context. Using the same entropy density, we show that the Brown-York (BY) tensor arises naturally as the projection of the canonical momentum conjugate to the normal vectors on the relevant hypersurface, thereby providing a common construction applicable to both timelike and null hypersurfaces. This perspective also offers insight into the structural differences of the null BY tensor, including its non-symmetric character. We further extend the analysis to scalar-tensor theories, showing that the entropy-based formulation reproduces the expected equations of motion along with the corresponding BY tensor, and, clarifies its non-conservation in the presence of additional scalar field which is non-minimally coupled. Our results provide a coherent variational interpretation of quasi-local gravitational quantities and reveal a common underlying structure linking bulk dynamics and boundary momentum. less

We initiate the study of gravitational-wave lensing in the wave-optics regime within modified gravity. We consider a phenomenological setup in which the gravitational-wave amplitude obeys a curvature-coupled propagation equation. This framework reproduces the standard GR behaviour in the geometric-optics regime, while leading to qualitatively different infrared dynamics. In particular, the usual argument implying that the amplification factor approaches unity in the zero-frequency limit no longer applies. This is due to the persistence of curvature-induced interactions in the infrared, which modify the natural propagation basis itself. As a result, the standard Fresnel treatment ceases to be valid at sufficiently low frequency. The correct infrared regime is instead controlled by an interacting static Green function, with a finite-frequency completion provided by a partial-wave formulation. We show that this structure admits an equivalent distorted-wave interpretation, in which the curvature interaction is absorbed into a dressed reference propagation basis, while the residual lensing effect is encoded in finite-frequency phase shifts. We further demonstrate that these phenomena admit a natural interpretation in the language of scattering amplitudes. Wave-optics lensing can therefore probe propagation-level departures from GR that remain entirely invisible in geometric optics. less

Close hyperbolic encounters between black holes produce distinctive bursts of gravitational radiation with a time-frequency morphology that is qualitatively different from that of quasi-circular inspirals. Expected to arise in dense stellar environments through dynamical interactions, these encounters probe formation channels and mass ranges inaccessible to isolated binary evolution, making them a compelling target for current and next-generation detectors. In this work, we reanalyze \totalevents high-mass events from the LIGO-Virgo-KAGRA catalogs using the hyperbolic configuration of the~\dali~waveform model. We compare these with analyses using the quasi-circular, precessing configuration of the same model, computing Bayes factors to evaluate which description is favored by the data. We find that most events strongly to mildly favor the quasi-circular, precessing scenario, except for GW190521. For this event, we find that the signal is best fit by a dynamical capture waveform, with Bayes factor $\ln \mathcal{B}^{\rm hyp}_{\rm prec}=3.71^{+0.11}_{-0.11}$. We confirm this preference via further analyses with~\dali~in different configurations (quasi-circular, non-precessing; eccentric, non-precessing; and eccentric, precessing), as well as one using the quasi-circular, precessing numerical relativity surrogate model \nrsur. We also highlight the results we obtain for GW231123, another high-mass signal linked to evidence of strong precession, for which we find strong preference for the quasi-circular, precessing scenario, with $\ln \mathcal{B}^{\rm hyp}_{\rm prec}=-15.80^{+0.24}_{-0.24}$. The analysis of mock signals generated with the best fitting waveforms for GW190521 and GW231123 suggest that the former might belong to a region of parameter space where high-mass, bound, precessing signals can be hard to distinguish from dynamical captures in parameter estimation. less

In a companion paper arXiv:2510.27676, we introduced a non-perturbative classical renormalisation group (RG) flow equation as a novel method for treating strongly interacting problems in general relativity, with a prominent application to the two-body problem. While we demonstrated that it reproduces perturbation theory, via the Post-Minkowskian (PM) expansion, and its computational efficiency in reproducing the 1PN Post-Newtonian action, its derivation was heuristic. In this work, we place this flow equation on a firm formal foundation. In particular, we demonstrate that a Legendre transform maps the classical analogue of the Polchinski equation precisely to our classical RG equation. This establishes a duality between equivalent, exact RG equations for the gravitational effective action. The result, combined with the successful applications in arXiv:2510.27676, solidifies the classical RG framework as a powerful and rigorous new approach to the general relativistic two-body problem and gravitational wave physics. less

We study the entropy of static, spherically symmetric black holes in diffeomorphism-invariant theories with nonminimal matter--curvature couplings, using the covariant phase space formalism. For regular bifurcate Killing horizons, the Iyer--Wald construction gives the standard Wald entropy. If a matter field cannot be smoothly extended to the regular bifurcation surface, however, the horizon surface charge variation can contain finite contributions that are not included in the Wald entropy density. In the representative obtained by directly varying the action, and after ordinary work terms are subtracted, we decompose the entropy entering the first law of black hole thermodynamics as \(S_{\mathrm H}=S_{\mathrm W}+S_1+ΔS\). Here \(S_{\mathrm W}\) is the Wald entropy, \(S_1\) is the non-Wald part of the Noether charge, and \(ΔS\) is the remaining integrable part of the horizon surface charge variation. Applying this criterion to Kalb--Ramond, bumblebee, and extended Gauss--Bonnet black holes, we find that the regular Kalb--Ramond branch has \(S_{\mathrm H}=S_{\mathrm W}\), the bumblebee branches yield either \(S_1=0\) with \(ΔS\neq0\) or a cancellation between \(S_1\) and \(ΔS\), and the Weyl-vector extended Gauss--Bonnet examples require both corrections. This gives a direct test of whether the Wald entropy density is sufficient, or whether the full horizon surface charge variation has to be used. less