2026-03-10 | | Total: 48
Black holes are the most efficient scramblers in nature. By mapping the instantaneous mass and angular momentum of two spinless black holes in a quasi-circular binary onto those of an effective Kerr black hole, we demonstrate that the final state of the merger remnant corresponds with remarkable accuracy to the configuration that renders null geodesics unstable at the highest possible rate. This suggests a deep connection between the properties of black holes resulting from binary mergers and their unstable null orbits.
This article is devoted to the investigation of the observational properties of the Schwarzschild black hole (BH) surrounded by a dark matter (DM) halo. Our study commences with a brief review of spacetime, including the horizon structure and curvature invariants, which are the Ricci scalar, the square of the Ricci tensor, and the Kretschmann scalar. Subsequently, we explore the massive and massless particle dynamics around the Schwarzschild BH surrounded by a dark matter halo, including the innermost stable circular orbit (ISCO) and photon sphere radii. It was found that the radius of the ISCO increases under the influence of the spacetime parameters. Additionally, we investigate the weak gravitational lensing with the assumption that the BH is surrounded by a uniform and non-uniform plasma. Finally, we examine the impact of a plasma on the BH shadow and employ Event Horizon Telescope (EHT) observational data to constrain the BH's parameters.
We investigate the geodesic structure of realistic static and spherically symmetric spacetimes embedding neutron stars in metric $f(R)$ gravity, focusing on the quadratic Starobinsky model $f(R)=aR^2$ with $a<0$. Neutron-star solutions are obtained by numerically solving the modified Tolman-Oppenheimer-Volkoff system for several realistic equations of state. Such solutions are then matched consistently to the exterior vacuum geometry by enforcing the full set of junction conditions required in metric $f(R)$ theories. Using an effective potential approach, we show that stable circular orbits appear in discrete radial bands separated by forbidden regions, with a dominant principal band of stability that depends sensitively on the stellar central pressure, the equation of state, and the magnitude of the parameter $|a|$. Outside the stable bands, massive particles can have bound but unstable precessing trajectories as well as unbounded motions. On the other hand, for null geodesics, we find no evidence for photon spheres outside the neutron star within the parameter range studied.
One of the main targets for space-borne gravitational wave detectors is the detection of Extreme Mass Ratio Inspirals (EMRIs). The data analysis of EMRIs requires waveform models that are both accurate and fast. The major challenge for the fast generation of such waveforms is the generation of the Teukolsky amplitudes for generic (eccentric and inclined) Kerr orbits. The requirement for the modeling of $\sim10^5$ harmonic modes across a four-dimensional parameter space makes traditional approaches, including direct computation or dense interpolation, computationally prohibitive. To overcome this issue, we introduce a convolutional encoder-decoder architecture for a fast and end-to-end global fitting of the Teukolsky amplitudes. We also adopt a transfer learning strategy to reduce the size of the training dataset, and the model is trained gradually from the simplest Schwarzschild circular orbits to generic Kerr orbits step by step. Within this framework, we obtain a surrogate model based on a semi-analytical Post-Newtonian dataset, and the full harmonic amplitudes can be generated within milliseconds, while the median mode-distribution error for generic orbits is approximately $\sim10^{-3}$. This result indicates that the framework is viable for constructing efficient waveform models for EMRIs.
A proper-time method for constructing models of dynamic gravitational-wave fields is presented. Using the proper-time method, analytical (not numerical) models of secondary gravitational waves are constructed as perturbative solutions of linearized field equations against the background of the exact wave solution of Einstein's equations for the vacuum in the Bianchi VI universe in a privileged wave coordinate system. Relations for the proper time of test particles against the background of a strong gravitational wave are used. The analytical form of the metric components for secondary gravitational waves is found from compatibility conditions for the field equations. From the field equations, an explicit form of ordinary differential equations and their solutions is obtained for functions included in small corrections to the metric for secondary gravitational waves. It is shown that there exists a continuum of gravitational wave parameters for which the perturbative solutions are stable.
Harada's recently announced Kerr-Schild solution of Einstein's equations in general relativity [Phys. Rev. D {\bf 112}, 024020 92025)] is identified as isometric with one of the well-known NUT solutions.
This work places the invariant $ds^2$ at the center of the gravitational interaction, interpreting it not as a purely geometric object but as the differential of proper time, endowed with direct physical meaning. Starting from the extension of Fermat's principle to massive particles--namely, the requirement that freely falling bodies follow trajectories that extremize proper time, which for timelike motion corresponds to a local maximum--and invoking the universality of Galilean free fall, we derive the form of $ds^2$ in a static gravitational field. Lorentz invariance then provides the natural framework to extend this result to systems involving moving matter. The invariant derived through this procedure matches the weak-field limit of General Relativity formulated in the harmonic gauge. Within this linearized regime, we show that the structure of the theory already contains the seeds of its non-linear completion: any dynamically consistent extension to strong gravitational fields necessarily involves the Ricci tensor. From this viewpoint, Einstein's field equations appear not as a postulated geometric law, but as the unique covariant closure required to ensure energy momentum conservation and the self consistency of the gravitational interaction.
We study black hole shadows in nonminimally coupled Weyl connection gravity, a metric-affine extension of general relativity in which spacetime is described by a metric and a Weyl vector field encoding non-metricity. Despite going beyond the Riemannian framework, the presence of a non-dynamical Weyl vector ensures second-order field equations. The theory admits Schwarzschild- and Reissner--Nordström-like solutions modified by a Weyl integration constant that parametrizes deviations from General Relativity. By computing the corresponding shadow radii and confronting them with the Event Horizon Telescope constraints on Sgr A*, we place observational bounds on the Weyl parameter. Assuming an observer distance $r_O = 4.1\times 10^{10}M$ and requiring consistency at the $2σ$ level, we obtain $ω\gtrsim 10^{11.7}M$ (model I), $ω\gtrsim 10^{10.5}M$ (model II), and $ω\sim 10^{12}M$ (model III). Our results show that present horizon-scale imaging already sets meaningful limits on spacetime non-metricity. This work highlights the power of black hole shadow observations as probes of extended gravitational dynamics and establishes a direct link between Weyl-based theories and current astrophysical data.
We investigate a simple $f(R, T_φ)$ gravity model coupled to a scalar field and demonstrate that the theory admits classical degenerate metric solutions, analogous to those known in general relativity. In particular, we identify a class of solutions that exhibits a smooth transition from a Euclidean to a Lorentzian domain, thus yielding a classical dynamical realization of signature change.
We show that the hidden symmetry and separability characteristic of Kerr geometry arise as an inevitable \emph{local} consequence of the Einstein equations for rotating spacetimes. Without assuming separability, algebraic speciality, Killing-Yano symmetry, or global boundary conditions, we analyze general stationary and axisymmetric geometries in a locally non-rotating orthonormal frame. Imposing a minimal physical requirement, the local equilibrium condition, we find that the mixed Einstein equations enforce a rigid projective alignment between the radial and angular sectors. This result does not rely on vacuum: the mixed equations are insensitive to the detailed form of a comoving stress-energy tensor. Consistency then requires equality of Schwarzian derivatives, leading to a universal classification of local solutions into Möbius, exponential, and trigonometric branches. Global regularity generically excludes the trigonometric branch, leaving precisely the Kerr-type sector and the emergence of Petrov type D structure. This provides a structural precursor to Kerr uniqueness, demonstrating that the kinematical core of Kerr geometry is encoded locally within the Einstein equations.
In this paper, we extend Horndeski's theory into the Palatini approach, assuming that the metric tensor and the (symmetric) connection are a priori independent objects. We introduce an additional transformation of the connection and write down the action functional being form-invariant under both the disformal transformation of the metric and the new transformation of the connection. We show that such a theory reduces on-shell to a metric subclass of Horndeski's gravity called kinetic gravity braiding. We also introduce an invariant metric and connection, and demonstrate that quantities defined in such a way lead to a metric theory. In the second part of the paper, we consider a simple cosmological model within the theory and explore its potential links with $ k$-essence-type theories, with a non-trivial coupling between the scalar field and the matter part of the action in the Einstein frame. We show that there exists a model that reproduces late-time cosmic acceleration, approaching asymptotically the de Sitter phase, motivating further study of the theories.
We study the propagation and emission of gravitational waves in the metric-affine formulation of the bumblebee model, where spontaneous Lorentz symmetry breaking arises from a vector field acquiring a nonvanishing vacuum expectation value. Working in the geometric-optics limit of the linearized theory, we derive the modified dispersion relation governing the graviton modes and show that it depends on the orientation of the wave vector relative to the background vector. The polarization sector is examined for timelike and spacelike configurations of the Lorentz-violating vacuum. In both cases only two independent tensor modes propagate, although their propagation properties and tensor structure depend on the orientation of the background field. We then construct the retarded Green function associated with the modified wave operator and determine the radiation-zone produced by localized sources. In the timelike configuration the Lorentz-violating effects appear through a modified propagation speed and an overall amplitude renormalization, leading to a shifted retarded time while preserving the quadrupole structure of the waveform. In contrast, the spacelike sector introduces anisotropic corrections to the quadrupole amplitude together with an additional contribution proportional to the third time derivative of the quadrupole moment. As an astrophysical application, the gravitational radiation emitted by a circular binary black hole system is evaluated, allowing observational constraints on the Lorentz-violating combination $ξb^{2}$ to be estimated using multimessenger bounds from GW170817/GRB~170817A and waveform consistency requirements from gravitational wave observations.
In this work, we test the geometrical structure and thermodynamic properties of the Einstein-Maxwell-Power-Yang-Mills black hole (BH) models, which constitute a nonlinear generalization of the standard Einstein-Yang-Mills theory through the inclusion of a power-law Yang-Mills invariant. Also, we begin by analyzing the spacetime geometry via the metric function $f(r)$ and examine the modifications induced by the electromagnetic charge and nonlinear Yang-Mills parameter on the horizon structure, causal structure, and gravitational potential. Subsequently, the dynamics of photons and massive particles are explored through the study of null and timelike geodesics, allowing the determination of the effective potential, photon sphere radius, and associated BH shadow. Also, the stability of circular photon orbits is quantified using the Lyapunov exponent, which characterizes the timescale of orbital instability and provides a direct link to observable photon ring features. For massive particles, the innermost stable circular orbit (ISCO) is calculated, illustrating the influence of BH parameters on the dynamics of accretion disks. From the thermodynamic viewpoint, we compute the principal thermodynamic quantities, including the BH mass, Hawking temperature, Bekenstein-Hawking entropy, heat capacity, and Gibbs free energy, to assess both local and global stability of the system. The divergence of the heat capacity signals the occurrence of second-order phase transitions, whereas the Gibbs free energy analysis identifies possible first-order phase transitions between distinct thermodynamic configurations. In this context, our results demonstrate that the nonlinear Yang-Mills parameter strongly affects the spacetime geometry, particle dynamics, and thermodynamic phase structure, shifting the location of stability regions and critical points associated with phase transitions.
We investigate the interior structure of asymptotically flat hairy black holes (HBHs) arising in the Einstein-Klein-Gordon theory with nonpositive-definite scalar potentials, where nontrivial scalar hair exists at the event horizon. While exterior properties, including shadow imaging for HBHs supported by an inverted Higgs-like potential have been extensively investigated, their interior structure remains largely unexplored. In many gravitational theories, backreaction of classical fields can significantly eliminate the Cauchy horizon, which is known to be highly unstable due to the mass inflation effect, raising important questions regarding the validity of the Strong Cosmic Censorship conjecture. These considerations motivate us to examine the interior structure of HBHs by numerically integrating the field equations inward from the outer horizon. We find that the scalar field and the metric functions increase monotonically inside the horizon and diverge as $r \rightarrow 0$. The Ricci and Kretschmann scalars also diverge at $r=0$, confirming the presence of a genuine curvature singularity. No additional root of the metric function is observed, indicating the absence of a Cauchy horizon in the electrically neutral HBHs considered here. Furthermore, the weak energy condition is violated throughout the interior region, and the degree of violation becomes more pronounced as the scalar field at the horizon increases. These results provide new insight into the global structure of HBHs and their implications for cosmic censorship.
Motivated by the entropy-functional formulation of emergent gravity, in which spacetime is endowed with a thermodynamic entropy that, upon extremization, yields the Einstein field equations, we reformulate the onset of a neutron star as a bulk gravitational-thermodynamic threshold problem. By combining the Jacobson-Padmanabhan perspective of gravity as an equation of state with Tolman redshift, Tolman-Oppenheimer-Volkoff hydrostatics, and a bulk binding versus thermal energy balance, we derive four characteristic temperature scales for a static compact object: the Newtonian virial threshold, the redshifted threshold as observed at infinity, the degenerate-matter threshold, and the surface screen temperature.The paper supplies a mathematically expanded derivation, explicit TOV and Fermi-gas formulae, and thermodynamic consistency relations.
In this work, density perturbations are investigated within the framework of a fractional Newtonian cosmology. Focusing on the matter-dominated era and employing the fluid-flow approach, the growth equation for density perturbations is derived and solved analytically. No dynamical instability arises in the physically relevant parameter space. It is shown that both the growth equation and its solutions depend explicitly on the fractional parameter $α$, and reduce to their standard Newtonian and relativistic counterparts in the special limit $α= 1$. The existence of both growing and decaying perturbative modes is confirmed, and, in accordance with current cosmological observations, the analysis is restricted to the growing mode. Using observational relations, in particular the Sachs--Wolfe equation, an observational upper bound on the parameter $α$ is obtained, which is more restrictive than the bounds inferred from background dynamics and theoretical perturbative considerations. When combined with the independent constraints arising from the background analysis, these results confine the fractional parameter $α$ to a narrow and physically viable region of parameter space. Overall, the present study indicates that, although the background evolution of the fractional model may closely mimic that of $Λ$CDM, the associated density perturbations generically carry a distinct fractional signature that can, in principle, be tested observationally.
We present a calculation of the full retarded Green functions of the Regge-Wheeler and Teukolsky equations obeyed by gravitational field perturbations of Schwarzschild spacetime. We perform the calculations for spacetime points along: (i) a timelike circular geodesic (where null-separated points are not at caustics); and (ii) a static worldline (where null-separated points are at caustics). These Green functions show a 4-fold singularity structure away from caustics, and 2-fold at caustics (similarly to the case of scalar field perturbations, which we also reproduce). Physical oscillations near the singularities appear in the gravitational case, which were not present in the scalar case. We obtain our results by developing various numerical and analytical methods.
We investigate elastic scattering by a compact, horizonless body in curved spacetime, considering a massless scalar wave incident on a static, spherically symmetric, uniform-density star of radius $R$ and mass $M$ with a Schwarzschild exterior. We introduce an exact Debye-series decomposition of the scattering matrix, in the spirit of Debye expansions in Mie scattering. This decomposition separates direct surface reflection from contributions involving transmission into the interior and subsequent propagation, and admits a natural trajectory interpretation. We then determine the associated Regge-Debye pole spectrum in the complex angular-momentum plane. For neutron-star-like tenuities ($R>3M$), the spectrum exhibits two pole families: a surface-wave branch associated with the surface matching condition and a broad-resonance branch associated with the interior regularity condition. For ultracompact objects ($R<3M$), the surface-wave branch persists, while the interior-resonance sector splits into broad- and narrow-resonance branches. We next reconstruct the scattering amplitude from the Debye partial-wave contributions and find excellent agreement with direct partial-wave calculations. Finally, we develop complex angular-momentum representations order by order in the Debye series, making explicit how the pole families and non-pole sectors contribute to each Debye term. In the neutron-star-like regime, we find a genuine competition between Regge-Debye pole sums and branch-cut contributions, and show that, at high frequency, the rainbow-like enhancement already arises from the first interior-transmission contribution and is dominated by the interior-resonance Regge-Debye poles. By contrast, in the ultracompact regime, the Debye amplitudes are overwhelmingly pole dominated.
Previous work on the gravitational-wave background generated in a two-scalar-field cosmological model, in which warm inflation and the dark sector are unified within a single framework, has shown that the resulting spectrum could be potentially detectable by planned next-generation gravitational-wave observatories. In this work, we extend this analysis to the weak dissipation regime of warm inflation, highlighting how the features of the inflationary scenario play a crucial role in the production of gravitational waves. The full gravitational-wave energy spectrum is calculated using the formalism of continuous Bogoliubov coefficients. By comparing our results with those obtained in the strong dissipation regime and with the sensitivity curves of future detectors, we find that the weak dissipation regime improves the prospects for observational detection.
Extreme mass ratio inspirals (EMRIs), as a core target for future space-based gravitational wave detection, offer crucial observational grounds for testing strong-field gravitational theories and classifying black holes through their orbital dynamics and gravitational wave radiation characteristics. This study systematically investigates the periodic orbit characteristics and gravitational wave radiation properties of EMRIs in the spacetime of a non-minimally coupled Einstein-Yang-Mills (EYM) black hole. The results show that as the magnetic charge parameter \(Q\) and the non-minimal coupling constant \(ξ\) increase, the radii of the marginally bound orbit (\(r_{\text{MBO}}\)) and the innermost stable circular orbit (\(r_{\text{ISCO}}\)) decrease significantly, with corresponding reductions in the orbital energy \(E\) and angular momentum \(L\). Furthermore, the allowable parameter space of energy and angular momentum (\(E\)-\(L\)) for bound orbits shifts towards the left. We then plot typical periodic orbits through orbit classification, finding that the non-minimal coupling effect suppresses the orbital contraction induced by the magnetic charge, leading to degeneracy of orbits with different \(Q\) values towards the Schwarzschild case, as well as phase shifts, amplitude enhancement, and period shortening in the gravitational waveforms with increasing \(Q\) and \(ξ\). These results provide theoretical predictions for distinguishing different black hole models through future gravitational wave observations.
Black hole solutions are studied here within the symmetric teleparallel formulation of gravity, employing the $f(Q)$ model in which the gravitational dynamics are governed by the non-metricity scalar $Q$. We focus on static, circularly symmetric spacetimes in $(2+1)$-dimensions, analyzing both charged and uncharged cases. By adopting a power-law form for $f(Q)$, we derive exact black hole solutions and explore their thermodynamic and geometric properties. Curvature and non-metricity scalars reveal central singularities stronger than those in general relativity. we find that the horizon radii increase with the charge parameter while higher values of the non-metricity coefficient, $c_{4}$, or the cosmological constant $Λ$ tend to merge or eliminate horizons, reducing their total number and altering the near-origin structure of the spacetime. We perform a detailed topological analysis based on the Euler characteristic and examine the geodesic completeness of the spacetime. Our findings show that, depending on the presence of electric charge, the singularity may or may not be reachable by geodesics. The thermodynamic stability is confirmed via temperature, entropy, and heat capacity calculations. This study highlights the rich structure of $f(Q)$ gravity in lower-dimensional settings and offers new insights into the nature of singularities and black hole topologies in modified gravity theories.
We investigate Josephson phenomena in a slowly rotating stationary spacetime, emphasizing the distinct roles of gravitational redshift and rotational frame dragging motivated by [10.1007/JHEP02(2026)006]. Using a covariant formulation based on gauge-invariant phase dynamics and conserved currents within a $3+1$ decomposition, we analyze both AC and DC Josephson effects and interferometric configurations. Restricting attention to linear order in the rotation parameter $a$, and working in the Eulerian/ZAMO frame, we show that in the slow-rotation slicing adopted here the lapse function agrees with its static (Schwarzschild-type) form up to $\mathcal{O}(a^2)$, while rotational effects enter through the shift vector. Consequently, redshift effects on Josephson frequencies and DC critical currents remain unchanged relative to the non-rotating case at $\mathcal{O}(a)$. The AC Josephson relation retains its redshifted structure when expressed in terms of proper voltages and reduces to the standard flat-spacetime form when formulated in terms of asymptotic (Killing-time) observables. Likewise, the DC critical current measured at infinity scales with a single power of the lapse function and is unaffected by rotation at linear order in the absence of azimuthal condensate momentum. Rotational effects become relevant only in configurations sensitive to spatial phase transport or to synchronization with respect to the global time coordinate. In particular, RF-driven interferometric setups can acquire Sagnac-type phase offsets associated with frame dragging, whereas the DC fluxoid constraint remains unshifted at linear order in the present approximation. Our results provide a clean separation between lapse-driven redshift effects and shift-driven rotational contributions in Josephson physics and furnish a consistent framework for superconducting circuits in stationary spacetimes.
An all-orders worldline effective action for Kerr-Newman black hole is achieved in twistor particle theory. Exact hidden symmetries are identified in self-dual backgrounds.
Quadratic quasinormal modes (QQNMs) and Christodoulou memory effect are key nonlinear phenomena in gravitational wave physics. QQNMs characterize the near zone nonlinear response of a perturbed black hole, whereas the memory effect is a nonlinear remnant imprinted at null infinity by outgoing radiation. This naturally raises the question of whether and in what sense the two can be bridged. We show that they are related through bridge coefficients which depend primarily on remnant black hole parameters during ringdown. Future space-based gravitational-wave detectors can probe this relation. These results provide a new avenue for testing gravity and a fresh perspective on the nonlinear regime of general relativity.
In this work, we introduce a method for finding exact solutions to the vacuum Einstein field equations in higher dimensions from a given solution to the chiral equation. When considering a $n + 2$-dimensional spacetime with $n$ commutative Killing vectors, the metric tensor can take the form $\hat g = f ( ρ, ζ) ( d ρ^2 + d ζ^2 ) + g_{μν} ( ρ, ζ) d x^μd x^ν$. Then, the Einstein field equations in vacuum reduce to a chiral equation, $( ρg_{, z} g ^{-1} )_{, \bar z} + ( ρg_{, \bar z} g ^{-1} )_{, z} = 0$, and two differential equations, $( \ln f ρ^{1-1/n} )_{, Z} = \fracρ{2} \operatorname{tr} ( g_{, _Z} g^{-1} )^2$, where $g \in SL( n, \mathbb{R} )$ is the normalized matrix representation of $g_{μν}$, $z = ρ+ i ζ$ and $Z = z, \bar z$. We use the ansatz $g = g ( ξ^a )$, where the parameters $ξ^a$ depend on $z$ and $\bar z$ and satisfy a generalized Laplace equation, $( ρξ^a _{, z} )_{, \bar z} + ( ρξ^a _{, \bar z} )_{, z} = 0$. The chiral equation to the Killing equation, $A_{a , ξ^b} + A_{b , ξ^a} = 0$, where $A_a = g_{, ξ^a} g^{-1}$. Furthermore, we assume that the matrices $A_a$ commute with each other; in this way, they fulfill the Killing equation.