We ask whether two classical time-machine geometries, the Ori (2005) compact-vacuum-core metric and the Ahmed (2018) four-dimensional generalisation of Misner space, remain admissible exact solutions when the gravitational sector is enlarged to the recently proposed $f(R,\mathcal{L}_{m},Φ,X)$ class, an extension of $f(R,\mathcal{L}_{m})$ that couples curvature, the matter Lagrangian density, a scalar field $Φ$, and its kinetic invariant $X = g^{μν}\nabla_μΦ\nabla_νΦ$. Working with the explicit model $f = R + \mathcal{L}_{m} + (λ/2)\,X$ and a vanishing scalar potential, we compute the curvature invariants, the modified field equations, and the effective stress-energy components produced by the harmonic scalar profile $Φ(x,y) = a(x^{2}-y^{2})/2$ in both backgrounds. The Ricci scalar vanishes for the Ori metric and obeys $R = e^{f}(f_{,xx}+f_{,yy})$ for the Ahmed metric; the kinetic invariant takes the explicit forms $X = a^{2}(x^{2}+y^{2})$ and $X = a^{2}e^{f}(x^{2}+y^{2})$, respectively. Both metrics solve the field equations of the modified theory with anisotropic matter sources, and the chronology-violating regions $g_{zz}<0$ (Ori) and $g_{ψψ}<0$ (Ahmed) survive the modification. Energy-density profiles measured by a closed-timelike-curve observer match those measured by a static observer outside the chronology horizon, so the additional scalar degree of freedom in $f(R,\mathcal{L}_{m},Φ,X)$ gravity does not enforce a chronology-protection mechanism in either background. The conclusion mirrors the parallel result for the Li time-machine and supplies a consistency test for scalar-extended modified gravity in non-globally-hyperbolic settings. less

Gravitoelectromagnetism in the Weyl formalism is investigated through an analysis of the consequences of a restricted gauge symmetry acting on the tensor field $A_{μν}$. The propagator associated with the GEM field is explicitly derived and decomposed within the Barnes--Rivers formalism, revealing contributions from the spin-2, spin-1, and scalar spin-0 sectors. By coupling the theory to conserved sources, it is shown that only the spin-2 and scalar sectors contribute to physical processes, while the spin-1 component decouples. The resulting effective propagator can then be written in a compact metric form closely resembling the graviton propagator of linearized General Relativity. The role of gauge fixing is also analyzed by considering both Lorentz-like and de Donder-type gauge conditions. It is shown that the Lorentz-like gauge is consistent with the restricted gauge symmetry of the theory and leads to gauge-independent physical amplitudes, whereas the de Donder gauge introduces residual gauge-dependent scalar contributions, signaling an incompatibility with the underlying symmetry structure. The gauge symmetry is further extended to the fermionic and electromagnetic sectors through diffeomorphism-like transformations. In both cases, conserved currents are derived and shown to coincide with the corresponding energy-momentum tensors, implying that the GEM field couples to matter in the same manner as in linearized gravity. Finally, the associated Ward identities are verified, providing a nontrivial consistency check of the gauge structure and interaction vertices of the theory. less

Gravitational waves (GWs) provide a powerful, theory-independent probe of the dynamical structure of spacetime and the cosmological background. We study linearized GW propagation in k-essence cosmology, where a non-canonical scalar field describes the dark sector. In the high-frequency (short-wavelength) approximation on a Friedmann--Lemaître--Robertson--Walker (FLRW) background, and restricting to the transverse-traceless tensor sector, we derive a modified evolution equation for tensor perturbations. The GW speed remains strictly luminal, consistent with multimessenger bounds such as GW170817, but the interaction with the background field $\barφ$ induces a time-dependent effective mass-like term $m_{\rm eff}$. This background-induced mass modifies the dispersion relation without introducing additional propagating degrees of freedom, leading to a cumulative, frequency-dependent phase shift in the waveform over cosmological distances. We show that $m_{\rm eff}$ is uniquely determined by background cosmological parameters and can be written as a redshift-dependent function, $m_{\rm eff}(z)$, directly linking GW observables to scalar-field dynamics, while the GW luminosity distance remains identical to its electromagnetic counterpart, preserving standard-siren consistency. We test the scenario through a joint Bayesian analysis that combines cosmic chronometers (CC), BAO, Pantheon+SH0ES, and standard-siren data from GWTC-2.1/3/4. The reconstruction is consistent with current constraints and reproduces the late-time expansion history, while the evolution of $m_{\rm eff}(z)$ offers a new mechanism that may help alleviate the $H_0$ tension. less

Gravitational waves admit a Stokes decomposition into intensity ($I$), circular polarization ($V$), and linear polarization ($Q$, $U$), analogous to Cosmic Microwave Background (CMB) polarimetry. We implement a full-Stokes maximum-likelihood SGWB map-making analysis for ground-based detector networks, promoting the standard cross-correlation data products used in existing pipelines to a joint reconstruction of $I$, $V$, $Q$, $U$. Applied to LVK O3 data, we constrain the polarized angular spectra $C^{VV}_\ell$, $C^{EE}_\ell$, $C^{BB}_\ell$ and $|C^{IV}_\ell|$. We show that an intensity-only model is biased when polarized sky components are present, since the detector-network Fisher inner product does not generally make the Stokes responses orthogonal. For transient CBC foregrounds, polarized shot noise is not parametrically suppressed relative to ordinary CBC intensity shot noise. The full Stokes framework separates the Stokes sectors while providing access to polarized anisotropies invisible to conventional intensity-only searches. less

We develop a time-domain numerical framework for global scalar wave scattering in Minkowski spacetime. The main contribution is an exact conformal matching of three compactified regions: a past hyperboloidal domain attached to $\mathscr I^-$, a Penrose domain covering a neighborhood of spatial infinity $i^0$, and a future hyperboloidal domain attached to $\mathscr I^+$. The matching surfaces are identical conformal hypersurfaces in the adjacent charts. This yields a global evolution scheme connecting $\mathscr I^-$, the neighborhood of $i^0$, and $\mathscr I^+$ without artificial timelike outer boundaries and without interpolation between scri-fixing gauges. We implement the construction for spherically symmetric scalar waves, including free propagation, localized linear scattering potentials such as the Pöschl--Teller potential, and semilinear wave equations with cubic, quintic, and septic nonlinearities. The numerical experiments demonstrate stable propagation across the matching interfaces, direct extraction of radiation at $\mathscr I^+$, and fourth-order convergence for the free and linear-potential tests. The quintic and septic nonlinear tests exhibit approximately fourth-order convergence and recover the expected late-time tail rates. The cubic case, by contrast, shows only first-order convergence, revealing a limitation of our treatment near compactified boundaries when the conformally rescaled nonlinear source remains non-vanishing. These results validate the conformal matching strategy for long-time simulations, while identifying the boundary regularity issues that must be addressed using a more robust treatment of spatial infinity. less

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